/[escript]/trunk/finley/src/CPPAdapter/MeshAdapterFactory.cpp
ViewVC logotype

Diff of /trunk/finley/src/CPPAdapter/MeshAdapterFactory.cpp

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

revision 3258 by jfenwick, Tue Sep 14 05:01:00 2010 UTC revision 3259 by jfenwick, Mon Oct 11 01:48:14 2010 UTC
# Line 12  Line 12 
12  *******************************************************/  *******************************************************/
13    
14    
15  #ifdef PASO_MPI  #ifdef ESYS_MPI
16  #include <mpi.h>  #include <mpi.h>
17  #endif  #endif
18  #ifdef USE_NETCDF  #ifdef USE_NETCDF
# Line 52  namespace finley { Line 52  namespace finley {
52    Domain_ptr loadMesh(const std::string& fileName)    Domain_ptr loadMesh(const std::string& fileName)
53    {    {
54  #ifdef USE_NETCDF  #ifdef USE_NETCDF
55      Paso_MPIInfo *mpi_info = Paso_MPIInfo_alloc( MPI_COMM_WORLD );      Esys_MPIInfo *mpi_info = Esys_MPIInfo_alloc( MPI_COMM_WORLD );
56      AbstractContinuousDomain* temp;      AbstractContinuousDomain* temp;
57      Finley_Mesh *mesh_p=NULL;      Finley_Mesh *mesh_p=NULL;
58      char error_msg[LenErrorMsg_MAX];      char error_msg[LenErrorMsg_MAX];
59    
60      char *fName = Paso_MPI_appendRankToFileName(fileName.c_str(),      char *fName = Esys_MPI_appendRankToFileName(fileName.c_str(),
61                                                  mpi_info->size,                                                  mpi_info->size,
62                                                  mpi_info->rank);                                                  mpi_info->rank);
63    
# Line 75  namespace finley { Line 75  namespace finley {
75      if (!dataFile.is_valid()) {      if (!dataFile.is_valid()) {
76        sprintf(error_msg,"loadMesh: Opening file NetCDF %s for reading failed.", fName);        sprintf(error_msg,"loadMesh: Opening file NetCDF %s for reading failed.", fName);
77        Finley_setError(IO_ERROR,error_msg);        Finley_setError(IO_ERROR,error_msg);
78        Paso_MPIInfo_free( mpi_info );        Esys_MPIInfo_free( mpi_info );
79        throw DataException(error_msg);        throw DataException(error_msg);
80      }      }
81    
# Line 194  namespace finley { Line 194  namespace finley {
194    
195          /* read elements */          /* read elements */
196          if (Finley_noError()) {          if (Finley_noError()) {
197            Finley_ReferenceElementSet  *refElements= Finley_ReferenceElementSet_alloc((ElementTypeId)Elements_TypeId,order, reduced_order);            Finley_ReferenceElementSet  *refElements= Finley_ReferenceElementSet_alloc((Finley_ElementTypeId)Elements_TypeId,order, reduced_order);
198            if (Finley_noError())  {            if (Finley_noError())  {
199                mesh_p->Elements=Finley_ElementFile_alloc(refElements, mpi_info);                mesh_p->Elements=Finley_ElementFile_alloc(refElements, mpi_info);
200            }            }
# Line 267  namespace finley { Line 267  namespace finley {
267    
268          /* get the face elements */          /* get the face elements */
269          if (Finley_noError()) {          if (Finley_noError()) {
270            Finley_ReferenceElementSet *refFaceElements=  Finley_ReferenceElementSet_alloc((ElementTypeId)FaceElements_TypeId ,order, reduced_order);            Finley_ReferenceElementSet *refFaceElements=  Finley_ReferenceElementSet_alloc((Finley_ElementTypeId)FaceElements_TypeId  ,order, reduced_order);
271            if (Finley_noError())  {            if (Finley_noError())  {
272                mesh_p->FaceElements=Finley_ElementFile_alloc(refFaceElements, mpi_info);                mesh_p->FaceElements=Finley_ElementFile_alloc(refFaceElements, mpi_info);
273            }            }
# Line 338  namespace finley { Line 338  namespace finley {
338    
339          /* get the Contact elements */          /* get the Contact elements */
340          if (Finley_noError()) {          if (Finley_noError()) {
341            Finley_ReferenceElementSet *refContactElements=   Finley_ReferenceElementSet_alloc((ElementTypeId)ContactElements_TypeId,order, reduced_order);            Finley_ReferenceElementSet *refContactElements=   Finley_ReferenceElementSet_alloc((Finley_ElementTypeId)ContactElements_TypeId,order, reduced_order);
342            if (Finley_noError())  {            if (Finley_noError())  {
343                mesh_p->ContactElements=Finley_ElementFile_alloc(refContactElements, mpi_info);                mesh_p->ContactElements=Finley_ElementFile_alloc(refContactElements, mpi_info);
344            }            }
# Line 410  namespace finley { Line 410  namespace finley {
410    
411          /* get the Points (nodal elements) */          /* get the Points (nodal elements) */
412          if (Finley_noError()) {          if (Finley_noError()) {
413            Finley_ReferenceElementSet *refPoints=    Finley_ReferenceElementSet_alloc((ElementTypeId)Points_TypeId,order, reduced_order);            Finley_ReferenceElementSet *refPoints=    Finley_ReferenceElementSet_alloc((Finley_ElementTypeId)Points_TypeId,order, reduced_order);
414            if (Finley_noError())  {            if (Finley_noError())  {
415                mesh_p->Points=Finley_ElementFile_alloc(refPoints, mpi_info);                mesh_p->Points=Finley_ElementFile_alloc(refPoints, mpi_info);
416            }            }

Legend:
Removed from v.3258  
changed lines
  Added in v.3259

  ViewVC Help
Powered by ViewVC 1.1.26