read_solution()

void read_solution

(

char *

fname

)

Reads an externally evaluated solution on Master node. Physical meaning of read quantity depends on Dirichlet or Neumann BC. If dirneu = 1 External file "fname" contains the normal derivatives of primary, secondary and total field on the boundary. If dirneu = 2 External file "fname" contains the primary, secondary and total field on the boundary.

Parameters

[in]

fname

File name.

                               {
  FILE *file;
  int mpi_ncntr;
  //Vec locSol;
  int i ;
  int     idummy;
  long double  x;
  long double  y;
  long double  z;
  long double  rphiinc;
  long double  iphiinc;
  long double  phinctot;
  long double  rphisc;
  long double  iphisc;
  long double  phisctot;
  long double  rphitot;
  long double  iphitot;
  long double  phitot;
  char s[255];

#define Rf1I12F "%d %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf %Lf\n"
#define Pf1I12F "%d %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE %+10.8LE\n"
  mpi_ncntr  = modgeominfo->ncntr + modgeominfo->nchief;


    file = fopen(fname,"r");
    // reading header
    for(i=0;i<7;i++){
      fgets(s,255,file);
    }
    for(i=0;i<mpi_ncntr;i++){
      fscanf(file,
           Rf1I12F,
           &idummy,   // CONTROL POINT INFO
           &x,
           &y,
           &z,
           &rphiinc, // INCIDENT
           &iphiinc,
           &phinctot,
           &rphisc, // SCATTERING
           &iphisc,
           &phisctot,
           &rphitot, // TOTAL
           &iphitot,
           &phitot
           );

      if(1 == modsolinfo->knw){
        solution->sol[i] = (DOUBLE)rphitot + I*(DOUBLE)iphitot;
      }else{
        solution->sol[i] = (DOUBLE)rphisc + I*(DOUBLE)iphisc;
      }
    }
    
    fclose(file);
}