eval_nw()

void eval_nw	(	COMPLEX	freq,
		int	icounter
	)

This routine evaluates the boundary conditions, according to the formulation to be used.

Parameters

[in]	freq	COMPLEX frequency
[in]	icounter	frequency index

                                       {

  int i;
  int ipsiu,npsiu;
  COMPLEX psiuval;


  int locRows;
  int gRows;
  int izero = 0;
#ifdef _ACOUSTO_DEBUG
  FILE *f;
  char fname[MAX_PATH];
#endif


  npsiu = modgeominfo->ncntr;
  gRows = npsiu;

  locRows  = numroc_(&gRows,&runinfo->row_block_size,&runinfo->myrow,&izero,&runinfo->nprows);


  if(icounter == 0){
   CREATE_VEC(solution->vecPsiu    ,locRows);
  }
  VecZeroEntries(solution->vecPsiu,locRows);

  if(runinfo->mycol != 0){ // vectors are distributed only on first col of grid
    return;
  }

  for(i=0;i<locRows;i++){
    sc_l2g(i,runinfo->myrow,gRows,runinfo->nprows,runinfo->row_block_size,&ipsiu);
    if(2 == modsolinfo->knw){
      psiuval = solution->vecPsinc[i];
      solution->vecPsiu[i] -= psiuval;
    }
  }


#ifdef _ACOUSTO_DEBUG
  sprintf(fname,"psiu.%d%d",runinfo->myrow,runinfo->mycol);
  f = fopen(fname,"w");
  for(i=0;i<locRows;i++){
    sc_l2g(i,runinfo->myrow,gRows,runinfo->nprows,runinfo->row_block_size,&ipsiu);
    fprintf(f," %d %f %f\n",ipsiu,
                           CREAL(solution->vecPsiu[i]),CIMAG(solution->vecPsiu[i])
                           );
  }
  fclose(f);
#endif
  

}