42
42
Finley_ReferenceElementSet *refPoints=NULL, *refContactElements=NULL, *refFaceElements=NULL, *refElements=NULL;
43
43
index_t node0, myRank;
45
Paso_MPIInfo *mpi_info = NULL;
45
Esys_MPIInfo *mpi_info = NULL;
47
47
#ifdef Finley_TRACE
48
48
double time0=Finley_timer();
51
51
/* get MPI information */
52
mpi_info = Paso_MPIInfo_alloc( MPI_COMM_WORLD );
52
mpi_info = Esys_MPIInfo_alloc( MPI_COMM_WORLD );
53
53
if (! Finley_noError()) {
66
66
sprintf(name,"Rectangular %d x %d mesh",N0,N1);
67
67
out=Finley_Mesh_alloc(name,DIM, mpi_info);
68
68
if (! Finley_noError()) {
69
Paso_MPIInfo_free( mpi_info );
69
Esys_MPIInfo_free( mpi_info );
72
refElements= Finley_ReferenceElementSet_alloc(Rec4,order,reduced_order);
72
refElements= Finley_ReferenceElementSet_alloc(Finley_Rec4,order,reduced_order);
73
73
if (useElementsOnFace) {
74
refFaceElements=Finley_ReferenceElementSet_alloc(Rec4Face, order, reduced_order);
75
refContactElements=Finley_ReferenceElementSet_alloc(Rec4Face_Contact, order, reduced_order);
74
refFaceElements=Finley_ReferenceElementSet_alloc(Finley_Rec4Face, order, reduced_order);
75
refContactElements=Finley_ReferenceElementSet_alloc(Finley_Rec4Face_Contact, order, reduced_order);
77
refFaceElements=Finley_ReferenceElementSet_alloc(Line2, order, reduced_order);
78
refContactElements=Finley_ReferenceElementSet_alloc(Line2_Contact, order, reduced_order);
77
refFaceElements=Finley_ReferenceElementSet_alloc(Finley_Line2, order, reduced_order);
78
refContactElements=Finley_ReferenceElementSet_alloc(Finley_Line2_Contact, order, reduced_order);
80
refPoints=Finley_ReferenceElementSet_alloc(Point1, order, reduced_order);
80
refPoints=Finley_ReferenceElementSet_alloc(Finley_Point1, order, reduced_order);
82
82
if ( Finley_noError()) {
95
Paso_MPIInfo_Split(mpi_info,NE1,&local_NE1,&e_offset1);
95
Esys_MPIInfo_Split(mpi_info,NE1,&local_NE1,&e_offset1);
99
Paso_MPIInfo_Split(mpi_info,NE0,&local_NE0,&e_offset0);
99
Esys_MPIInfo_Split(mpi_info,NE0,&local_NE0,&e_offset0);
291
291
Finley_ReferenceElementSet_dealloc(refContactElements);
292
292
Finley_ReferenceElementSet_dealloc(refFaceElements);
293
293
Finley_ReferenceElementSet_dealloc(refElements);
294
Paso_MPIInfo_free( mpi_info );
294
Esys_MPIInfo_free( mpi_info );