read_arrays()

read_arrays

(

)

Actual file reading. Reads arrays from different files of file_info struct

                 {
  FILE* fchief = NULL;
  FILE* fmics;
  float ver;
  long double x,y,z;
  int i,j, imic, ivbase, icount;
  int ntmp,netmp;
  char line[100];
  int i1,i2,i3;
  long double d1,d2,d3;
  int j1,j2,j3,j4,ftype,datasize,eltype;
  int ielg,ntag,nphy,ngeo,npart;

  // CHIEF points
  if( 0 < modgeominfo->nchief){
    logger(LOG_DEBUG,"reading file %s\n",fileinfo->chief_file);
    fchief = fopen(fileinfo->chief_file,"r");
    if(NULL == fchief){
      printf(MSG_GEOM_FILE_NOTFOUND,fileinfo->chief_file);
      return ACOUSTO_ERROR_GEOMETRY;
    }
  }

  // Microphones
  logger(LOG_DEBUG,"reading file %s\n",fileinfo->mics_file);
  fmics = fopen(fileinfo->mics_file,"r");
  if(NULL == fmics){
    printf(MSG_GEOM_FILE_NOTFOUND,fileinfo->mics_file);
    return ACOUSTO_ERROR_GEOMETRY;
  }


  /* Reading chief */
  if(0 < modgeominfo->nchief){
   for(i=0;i<modgeominfo->nchief;i++){
     fgets(line,100,fchief);
     sscanf(line,VECTOR_SCANF,&x,&y,&z);

     geometry->chief[i].x = x;
     geometry->chief[i].y = y;
     geometry->chief[i].z = z;

     /* Adding CHIEF points into geometry */
     vec_copy(&geometry->cntrc[modgeominfo->ncntr+i], geometry->chief[i]);

     _DEBUG_IDX_AND_VECTOR("CHIEF",i,geometry->chief[i]);
   }
  }
  /* Reading mics */
  fgets(line,100,fmics);
  if(0 == strcmp("$MeshFormat\n", line)){
      /* Reading header */
      fgets(line,100,fmics);
      sscanf(line,GMSH_VER_SCANF,&ver,&ftype,&datasize);
      if(ver < 2.2){
        logger(LOG_ERROR,"Gmsh file version %f not supported\n",ver);
        return ACOUSTO_ERROR_GEOMETRY;
      }
      if(0 != ftype){
        logger(LOG_ERROR,"Gmsh file type %d not supported\n",ftype);
        return ACOUSTO_ERROR_GEOMETRY;
      }

      icount=0;
      while (fgets(line,100,fmics)!=NULL){
        if(0 == strcmp("$Nodes\n",line)){
          fgets(line,100,fmics);
          sscanf(line,"%d",&ntmp);
          for(imic=0;imic<ntmp;imic++){
              fgets(line,100,fmics);
              sscanf(line,GMSH_NOD_SCANF,&i,&x,&y,&z);
              geometry->mics[imic].x = x;
              geometry->mics[imic].y = y;
              geometry->mics[imic].z = z;
          }
        }
        if(0 == strcmp("$Elements\n",line)){
          fgets(line,100,fmics);
          sscanf(line,"%d",&netmp);
          for(j=0;j<netmp;j++){
              fgets(line,100,fmics);
              sscanf(line,"%d %d %d",&ielg,&eltype,&ntag);
              if(2==eltype){         // Triangle
                ivbase=icount*4;
                if(2==ntag){
                  sscanf(line,"%d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&j1,&j2,&j3);
                  geometry->jnodm[ivbase   ] = j1-1;
                  geometry->jnodm[ivbase+1 ] = j2-1;
                  geometry->jnodm[ivbase+2 ] = j3-1;
                  geometry->jnodm[ivbase+3 ] = j3-1;
                }
                if(3==ntag){
                  sscanf(line,"%d %d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&npart,&j1,&j2,&j3);
                  geometry->jnodm[ivbase   ] = j1-1;
                  geometry->jnodm[ivbase+1 ] = j2-1;
                  geometry->jnodm[ivbase+2 ] = j3-1;
                  geometry->jnodm[ivbase+3 ] = j3-1;
                }
                geometry->kelmm[icount] = 2;
                icount++;
              }
              if(3==eltype){         // Quadrangle
                ivbase=icount*4;
                if(2==ntag){
                 sscanf(line,"%d %d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&j1,&j2,&j3,&j4);
                 geometry->jnodm[ivbase   ] = j1-1;
                 geometry->jnodm[ivbase+1 ] = j2-1;
                 geometry->jnodm[ivbase+2 ] = j3-1;
                 geometry->jnodm[ivbase+3 ] = j4-1;
                }
                if(3==ntag){
                  sscanf(line,"%d %d %d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&npart,&j1,&j2,&j3,&j4);
                  geometry->jnodm[ivbase   ] = j1-1;
                  geometry->jnodm[ivbase+1 ] = j2-1;
                  geometry->jnodm[ivbase+2 ] = j3-1;
                  geometry->jnodm[ivbase+3 ] = j4-1;
                }
                geometry->kelmm[icount] = 3;
                icount++;
              }
              if(4==eltype){         // Tetrahedra
                ivbase=icount*4;
                if(2==ntag){
                  sscanf(line,"%d %d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&j1,&j2,&j3,&j4);
                  geometry->jnodm[ivbase   ] = j1-1;
                  geometry->jnodm[ivbase+1 ] = j2-1;
                  geometry->jnodm[ivbase+2 ] = j3-1;
                  geometry->jnodm[ivbase+3 ] = j4-1;
                }
                if(3==ntag){
                  sscanf(line,"%d %d %d %d %d %d %d %d %d %d\n",&ielg,&eltype,&ntag,&nphy,&ngeo,&npart,&j1,&j2,&j3,&j4);
                  geometry->jnodm[ivbase   ] = j1-1;
                  geometry->jnodm[ivbase+1 ] = j2-1;
                  geometry->jnodm[ivbase+2 ] = j3-1;
                  geometry->jnodm[ivbase+3 ] = j4-1;
                }
                geometry->kelmm[icount] = 4;
                icount++;
              }
          }
          if(icount != modgeominfo->nemics){
            logger(LOG_ERROR,"Number of read elements %d not equal to mics elements %d \n",icount,modgeominfo->nemics);
            return ACOUSTO_ERROR_GEOMETRY;
          }
        }
      }
  }else{
      sscanf(line,f3I3FSCAN,&i1,&i2,&i3,&d1,&d2,&d3);
      for(i=0;i<modgeominfo->nmics;i++){
        fgets(line,100,fmics);
        if(NULL != line){ // if EOF it means we declared less microphones  than node numbers
          sscanf(line,VECTOR_SCANF,&x,&y,&z);
        }
        geometry->mics[i].x = x;
        geometry->mics[i].y = y;
        geometry->mics[i].z = z;
        _DEBUG_IDX_AND_VECTOR("MICS",i,geometry->mics[i]);
      }
  }

  if(modgeominfo->nchief > 0) fclose(fchief);
  fclose(fmics);

  if(modgeominfo->nchief > 0 )logger(LOG_DEBUG,"file %s read\n",fileinfo->chief_file);
  if(fileinfo->mics_file)     logger(LOG_DEBUG,"file %s read\n",fileinfo->mics_file);
  return 1;

}