georeader()

georeader

(

)

Main Geometry function, reads data, builds geometry

               {

 int nvrtp,nreal,nintt;
#ifdef _ACOUSTO_DEBUG 
 FILE *f;
 char fname[20];
#endif

 geometry = malloc(sizeof(struct geometry_struct)); // allocating structure

 runinfo->ndim = 3; /* built-in */
 nvrtp = 4;         // this version works for quadrangular panels 

 modgeominfo->nvelb = modgeominfo->nelmb*nvrtp; // number of vertices

 allocgeom(); // allocation
 geometry->ncntrs   = modgeominfo->ncntr; // number of control points

 // calculating allocation size
 nreal=runinfo->ndim*(modgeominfo->nnodb+
                      modgeominfo->ncntr+
                      modgeominfo->nelmb+
                      modgeominfo->nchief)+
                      modgeominfo->nelmb +
       runinfo->ndim* modgeominfo->nmics;

 nintt= modgeominfo->nvelb+5*modgeominfo->nelmb+4*modgeominfo->nemics;


 // Geometry is read only on master node and sent to every other node
 if(0 == rank ){
   if(1 != build_geoms()){
     logger(LOG_ERROR,MSG_GEOM_READ_ERROR);
     return ACOUSTO_ERROR_GEOMETRY;
   }
   if(1 != read_arrays()){
     logger(LOG_ERROR,MSG_GEOM_READ_ERROR);
     return ACOUSTO_ERROR_GEOMETRY;
   }
 }
 logger(LOG_DEBUG," Total n.of nodes          = %d\n",modgeominfo->nnodb);
 logger(LOG_DEBUG," Total n.of elements       = %d\n",modgeominfo->nelmb);
 logger(LOG_DEBUG," Total n.of control points = %d\n",modgeominfo->ncntr);

 // Broadcasting geometry
 logger(LOG_DEBUG," Broadcasting geometry...\n");
 MPI_Bcast(geometry->nodes    ,modgeominfo->nnodb ,mpi_vector_type     ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->jnodb    ,modgeominfo->nvelb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->jphib    ,modgeominfo->nelmb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->jpsib    ,modgeominfo->nelmb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->mvrtb    ,modgeominfo->nelmb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->kcreb    ,modgeominfo->nelmb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->kelmb    ,modgeominfo->nelmb ,MPI_INT             ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->mics     ,modgeominfo->nmics ,mpi_vector_type     ,0,MPI_COMM_WORLD);
 if(modgeominfo->nchief > 0) MPI_Bcast(geometry->chief    ,modgeominfo->nchief,mpi_vector_type     ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->cntrc    ,modgeominfo->ncntr+
                               modgeominfo->nchief,mpi_vector_type     ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->cntr     ,modgeominfo->ncntr ,mpi_vector_type     ,0,MPI_COMM_WORLD);
 MPI_Bcast(geometry->elements ,modgeominfo->nelmb ,mpi_element_type    ,0,MPI_COMM_WORLD);
 if(modgeominfo->nemics > 0) MPI_Bcast(geometry->jnodm    ,4*modgeominfo->nemics ,MPI_INT          ,0,MPI_COMM_WORLD);
 if(modgeominfo->nemics > 0) MPI_Bcast(geometry->kelmm    ,modgeominfo->nemics ,MPI_INT          ,0,MPI_COMM_WORLD);
 logger(LOG_DEBUG," Broadcasting geometry done\n");
 
 return 1;

}