/[escript]/trunk/finley/src/Mesh_rec4.cpp
ViewVC logotype

Diff of /trunk/finley/src/Mesh_rec4.cpp

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

revision 4491 by caltinay, Fri May 31 07:09:03 2013 UTC revision 4492 by caltinay, Tue Jul 2 01:44:11 2013 UTC
# Line 41  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 41  Finley_Mesh* Finley_RectangularMesh_Rec4
41    dim_t N0,N1,NE0,NE1,i0,i1,k,Nstride0=0,Nstride1=0, local_NE0, local_NE1, local_N0=0, local_N1=0, global_i0, global_i1;    dim_t N0,N1,NE0,NE1,i0,i1,k,Nstride0=0,Nstride1=0, local_NE0, local_NE1, local_N0=0, local_N1=0, global_i0, global_i1;
42    index_t offset0=0, offset1=0, e_offset0=0, e_offset1=0;    index_t offset0=0, offset1=0, e_offset0=0, e_offset1=0;
43    dim_t totalNECount,faceNECount,NDOF0=0,NDOF1=0,NFaceElements,NN;    dim_t totalNECount,faceNECount,NDOF0=0,NDOF1=0,NFaceElements,NN;
44    Finley_ReferenceElementSet *refPoints=NULL, *refContactElements=NULL, *refFaceElements=NULL, *refElements=NULL;    ReferenceElementSet *refPoints=NULL, *refContactElements=NULL, *refFaceElements=NULL, *refElements=NULL;
45    index_t node0, myRank;    index_t node0, myRank;
46    Finley_Mesh* out;    Finley_Mesh* out;
47    Esys_MPIInfo *mpi_info = NULL;    Esys_MPIInfo *mpi_info = NULL;
# Line 71  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 71  Finley_Mesh* Finley_RectangularMesh_Rec4
71        Esys_MPIInfo_free( mpi_info );        Esys_MPIInfo_free( mpi_info );
72        return NULL;        return NULL;
73    }    }
74    refElements= Finley_ReferenceElementSet_alloc(Finley_Rec4,order,reduced_order);    refElements= ReferenceElementSet_alloc(Rec4,order,reduced_order);
75    if (useElementsOnFace) {    if (useElementsOnFace) {
76          refFaceElements=Finley_ReferenceElementSet_alloc(Finley_Rec4Face, order, reduced_order);          refFaceElements=ReferenceElementSet_alloc(Rec4Face, order, reduced_order);
77          refContactElements=Finley_ReferenceElementSet_alloc(Finley_Rec4Face_Contact, order, reduced_order);          refContactElements=ReferenceElementSet_alloc(Rec4Face_Contact, order, reduced_order);
78    } else {    } else {
79          refFaceElements=Finley_ReferenceElementSet_alloc(Finley_Line2, order, reduced_order);          refFaceElements=ReferenceElementSet_alloc(Line2, order, reduced_order);
80          refContactElements=Finley_ReferenceElementSet_alloc(Finley_Line2_Contact, order, reduced_order);          refContactElements=ReferenceElementSet_alloc(Line2_Contact, order, reduced_order);
81    }    }
82    refPoints=Finley_ReferenceElementSet_alloc(Finley_Point1, order, reduced_order);    refPoints=ReferenceElementSet_alloc(Point1, order, reduced_order);
83        
84    if ( Finley_noError()) {    if ( Finley_noError()) {
85        
# Line 289  Finley_Mesh* Finley_RectangularMesh_Rec4 Line 289  Finley_Mesh* Finley_RectangularMesh_Rec4
289    }    }
290    
291    /* free up memory */    /* free up memory */
292    Finley_ReferenceElementSet_dealloc(refPoints);    ReferenceElementSet_dealloc(refPoints);
293    Finley_ReferenceElementSet_dealloc(refContactElements);    ReferenceElementSet_dealloc(refContactElements);
294    Finley_ReferenceElementSet_dealloc(refFaceElements);    ReferenceElementSet_dealloc(refFaceElements);
295    Finley_ReferenceElementSet_dealloc(refElements);    ReferenceElementSet_dealloc(refElements);
296    Esys_MPIInfo_free( mpi_info );      Esys_MPIInfo_free( mpi_info );  
297    
298    return out;    return out;

Legend:
Removed from v.4491  
changed lines
  Added in v.4492

  ViewVC Help
Powered by ViewVC 1.1.26