/[escript]/branches/domexper/dudley/src/Mesh_tet4.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Mesh_tet4.c

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

revision 787 by bcumming, Wed Jul 26 01:46:45 2006 UTC revision 788 by bcumming, Wed Jul 26 05:12:15 2006 UTC
# Line 225  Finley_Mesh* Finley_RectangularMesh_Hex8 Line 225  Finley_Mesh* Finley_RectangularMesh_Hex8
225        
226    /*  allocate tables: */    /*  allocate tables: */
227    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);
228      Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);
229      Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
230  #ifdef PASO_MPI  #ifdef PASO_MPI
231    Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, NDOF0*NDOF1*NDOF2, 0, 0 );    Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, NDOF0*NDOF1*NDOF2, 0, 0 );
232      Finley_ElementDistribution_allocTable( out->Elements->elementDistribution, NE0*NE1*NE2, NE0*NE1*NE2);
233      Finley_ElementDistribution_allocTable( out->FaceElements->elementDistribution, NFaceElements, NFaceElements );
234  #endif  #endif
   Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);  
   Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);  
235    if (! Finley_noError()) {    if (! Finley_noError()) {
236        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
237        return NULL;        return NULL;
# Line 666  Finley_Mesh* Finley_RectangularMesh_Hex8 Line 668  Finley_Mesh* Finley_RectangularMesh_Hex8
668        
669    /*  allocate tables: */    /*  allocate tables: */
670    Finley_NodeFile_allocTable(out->Nodes,N0t*N1*N2);    Finley_NodeFile_allocTable(out->Nodes,N0t*N1*N2);
   Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, numDOFLocal*NDOF1*NDOF2, NDOF1*NDOF2*2, 0 );  
671    Finley_ElementFile_allocTable(out->Elements,(numElementsLocal)*NE1*NE2);    Finley_ElementFile_allocTable(out->Elements,(numElementsLocal)*NE1*NE2);
672    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
673    
674      Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, numDOFLocal*NDOF1*NDOF2, NDOF1*NDOF2*2, 0 );
675      Finley_ElementDistribution_allocTable( out->Elements->elementDistribution, numElementsLocal*NE1*NE2, NE1*NE2*(numElementsLocal-boundaryRight*(!periodic[1])) );
676      Finley_ElementDistribution_allocTable( out->FaceElements->elementDistribution, NFaceElements, NFaceElements-2*boundaryRight*(NE2*(!periodic[1])+NE1*(!periodic[2])) );
677    if (! Finley_noError()) {    if (! Finley_noError()) {
678        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
679        return NULL;        return NULL;

Legend:
Removed from v.787  
changed lines
  Added in v.788

  ViewVC Help
Powered by ViewVC 1.1.26