/[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 787 by bcumming, Wed Jul 26 01:46:45 2006 UTC revision 788 by bcumming, Wed Jul 26 05:12:15 2006 UTC
# Line 220  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 220  Finley_Mesh* Finley_RectangularMesh_Hex2
220        
221    /*  allocate tables: */    /*  allocate tables: */
222    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);    Finley_NodeFile_allocTable(out->Nodes,N0*N1*N2);
   Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, NDOF0*NDOF1*NDOF2, 0, 0 );  
223    Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);    Finley_ElementFile_allocTable(out->Elements,NE0*NE1*NE2);
224    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
225      Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, NDOF0*NDOF1*NDOF2, 0, 0 );
226      Finley_ElementDistribution_allocTable( out->Elements->elementDistribution, NE0*NE1*NE2, NE0*NE1*NE2);
227      Finley_ElementDistribution_allocTable( out->FaceElements->elementDistribution, NFaceElements, NFaceElements );
228  #else  #else
229    out=Finley_Mesh_alloc(name,3,order);    out=Finley_Mesh_alloc(name,3,order);
230    
# Line 821  Finley_Mesh* Finley_RectangularMesh_Hex2 Line 823  Finley_Mesh* Finley_RectangularMesh_Hex2
823    
824    /*  allocate tables: */    /*  allocate tables: */
825    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*3, 0 );  
826    Finley_ElementFile_allocTable(out->Elements,numElementsLocal*NE1*NE2);    Finley_ElementFile_allocTable(out->Elements,numElementsLocal*NE1*NE2);
827    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);    Finley_ElementFile_allocTable(out->FaceElements,NFaceElements);
828        
829      Finley_NodeDistribution_allocTable( out->Nodes->degreeOfFreedomDistribution, numDOFLocal*NDOF1*NDOF2, NDOF1*NDOF2*3, 0 );
830      Finley_ElementDistribution_allocTable( out->Elements->elementDistribution, numElementsLocal*NE1*NE2, NE1*NE2*(numElementsLocal-boundaryRight*(!periodic[1])) );
831      Finley_ElementDistribution_allocTable( out->FaceElements->elementDistribution, NFaceElements, NFaceElements-2*boundaryRight*(NE2*(!periodic[1])+NE1*(!periodic[2])) );
832    if (! Finley_noError()) {    if (! Finley_noError()) {
833        Finley_Mesh_dealloc(out);        Finley_Mesh_dealloc(out);
834        return NULL;        return NULL;

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

  ViewVC Help
Powered by ViewVC 1.1.26