/[escript]/trunk/finley/src/Mesh_hex20.c
ViewVC logotype

Diff of /trunk/finley/src/Mesh_hex20.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 616 by elspeth, Wed Mar 22 02:46:56 2006 UTC revision 730 by bcumming, Mon May 15 04:03:49 2006 UTC
# Line 99  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 99  Finley_Mesh* Finley_RectangularMesh_Hex2
99    /*  allocate mesh: */    /*  allocate mesh: */
100        
101    sprintf(name,"Rectangular %d x %d x %d mesh",N0,N1,N2);    sprintf(name,"Rectangular %d x %d x %d mesh",N0,N1,N2);
102      /* TEMPFIX */
103    #ifndef PASO_MPI
104    out=Finley_Mesh_alloc(name,3,order);    out=Finley_Mesh_alloc(name,3,order);
105    
106    if (! Finley_noError()) return NULL;    if (! Finley_noError()) return NULL;
107    
108    out->Elements=Finley_ElementFile_alloc(Hex20,out->order);    out->Elements=Finley_ElementFile_alloc(Hex20,out->order);
# Line 111  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 114  Finley_Mesh* Finley_RectangularMesh_Hex2
114       out->ContactElements=Finley_ElementFile_alloc(Rec8_Contact,out->order);       out->ContactElements=Finley_ElementFile_alloc(Rec8_Contact,out->order);
115    }    }
116    out->Points=Finley_ElementFile_alloc(Point1,out->order);    out->Points=Finley_ElementFile_alloc(Point1,out->order);
117    #else
118      /* TODO */
119      PASO_MPI_TODO;
120      out = NULL;
121    #endif
122    if (! Finley_noError()) {    if (! Finley_noError()) {
123        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
124        return NULL;        return NULL;
# Line 118  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 126  Finley_Mesh* Finley_RectangularMesh_Hex2
126    
127        
128    /*  allocate tables: */    /*  allocate tables: */
129      #ifndef PASO_MPI
130    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);
131    #else
132      /* TODO */
133    #endif
134    Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);    Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);
135    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
136    if (! Finley_noError()) {    if (! Finley_noError()) {

Legend:
Removed from v.616  
changed lines
  Added in v.730

  ViewVC Help
Powered by ViewVC 1.1.26