estimate_memory2()

MEMORY estimate_memory2

(

)

                          {
  MEMORY ESTIMATED_MEMORY = 0L;
  int lrow, lcol, lrowM, lcolE;
  int izero = 0;

  lrow = numroc_(&modgeominfo->ncntr, &runinfo->row_block_size,
      &runinfo->myrow, &izero, &runinfo->nprows);
  lrowM = numroc_(&modgeominfo->nmics, &runinfo->row_block_size,
      &runinfo->myrow, &izero, &runinfo->nprows);
  lcol = numroc_(&modgeominfo->ncntr, &runinfo->col_block_size,
      &runinfo->mycol, &izero, &runinfo->npcols);
  lcolE = numroc_(&modgeominfo->nelmb, &runinfo->col_block_size,
      &runinfo->mycol, &izero, &runinfo->npcols);

  if (modsolinfo->pre_calculate_coefs) {
    ESTIMATED_MEMORY += (MEMsizeof(struct bodycoefs_struct));
    ESTIMATED_MEMORY += 2L* TOMEMORY(lrow) * TOMEMORY(lcol) * TOMEMORY(runinfo->nsymm) * MEMsizeof(DOUBLE); // B Body Coefs
    ESTIMATED_MEMORY +=
        TOMEMORY(4L*( TOMEMORY(lrow) ) * TOMEMORY(lcol)* TOMEMORY(runinfo->nsymm) * MEMsizeof(DOUBLE)); // C Body Coefs

    ESTIMATED_MEMORY += TOMEMORY(lrow)
        * MEMsizeof(DOUBLE); // temporary vector
    ESTIMATED_MEMORY += MEMsizeof(struct micscoefs_struct);
    ESTIMATED_MEMORY += 2L * TOMEMORY(lrowM)
        * TOMEMORY(lcolE)* TOMEMORY(runinfo->nsymm) * MEMsizeof(DOUBLE); // 2 mics B coefficients
    ESTIMATED_MEMORY += 3L * TOMEMORY(lrowM)
        * TOMEMORY(lcolE)* TOMEMORY(runinfo->nsymm) * MEMsizeof(DOUBLE); // 3 mics C coefficients
  }
  ESTIMATED_MEMORY +=
      (TOMEMORY(modgeominfo->nnodb) * MEMsizeof(struct vector)); // Nodes
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nvelb) * MEMsizeof(int)); // Vertex Indices
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(int)); // Index Vector
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(int)); // Index Vector
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(int)); // mvtrb
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(int)); // kcreb
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(int)); // kelmb
  ESTIMATED_MEMORY +=
      (TOMEMORY(modgeominfo->nmics) * MEMsizeof(struct vector)); // microphones
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->nchief)
      * MEMsizeof(struct vector)); // chief
  ESTIMATED_MEMORY += (TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)
      * MEMsizeof(struct vector)); // ctrl + chief
  ESTIMATED_MEMORY +=
      (TOMEMORY(modgeominfo->ncntr) * MEMsizeof(struct vector)); // cntr
  ESTIMATED_MEMORY +=
      (TOMEMORY(modgeominfo->nelmb) * MEMsizeof(struct panel4)); // elements
    if(PSEUDOINVERSE == modsolinfo->solver){
      ESTIMATED_MEMORY+= (2L*MEMsizeof(COMPLEX)); // work in linsys
      ESTIMATED_MEMORY+= (TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)*MEMsizeof(COMPLEX));// work in linsys
    }else if (GMRES == modsolinfo->solver){
      ESTIMATED_MEMORY+= (4*TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)*MEMsizeof(COMPLEX));// work in linsys
      ESTIMATED_MEMORY+= (3*TOMEMORY(lrow)*MEMsizeof(COMPLEX));// work in linsys
      ESTIMATED_MEMORY+= (TOMEMORY(lrow)*TOMEMORY(modsolinfo->restart)*MEMsizeof(COMPLEX));// work in linsys
      ESTIMATED_MEMORY+= (TOMEMORY(modsolinfo->restart+1)*TOMEMORY(lrow)*MEMsizeof(COMPLEX));// work in linsys
    }
  ESTIMATED_MEMORY += MEMsizeof(struct run_info);
  ESTIMATED_MEMORY += MEMsizeof(struct run_details);
  ESTIMATED_MEMORY += MEMsizeof(struct modgeom_info);
  ESTIMATED_MEMORY += MEMsizeof(struct modcoefac_info);
  ESTIMATED_MEMORY += MEMsizeof(struct modcoemic_info);
  ESTIMATED_MEMORY += MEMsizeof(struct modsol_info);
  ESTIMATED_MEMORY += MEMsizeof(struct solution_struct);
  ESTIMATED_MEMORY += ((strlen(datfile) + 1) * MEMsizeof(char)); // datfile
  ESTIMATED_MEMORY += (TOMEMORY(modsolinfo->n_freqresp_cnt) * MEMsizeof(int)); // ctrl freq response
  ESTIMATED_MEMORY += (TOMEMORY(modsolinfo->n_freqresp_mic) * MEMsizeof(int)); // mics freq response

  if (modsolinfo->pre_calculate_coefs) {
    ESTIMATED_MEMORY+= (TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)*TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX)); // matB
    ESTIMATED_MEMORY+= (TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)*TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX));// matBH
    ESTIMATED_MEMORY+= (TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)*TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX));// matBK
  }
  ESTIMATED_MEMORY+= (TOMEMORY(lrow)*TOMEMORY(lcol)*MEMsizeof(COMPLEX)); // matY
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);//
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->nsig ) * MEMsizeof(DOUBLE);
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->nome ) * MEMsizeof(DOUBLE);
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->nsourc )
      * MEMsizeof(struct acoustic_point);
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->ndipoles )
      * MEMsizeof(struct dipole);

  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->nome )*TOMEMORY(3*modgeominfo->ncntr) * MEMsizeof(COMPLEX); // custombc
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->n_freqresp_mic) * MEMsizeof(int);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(5*modgeominfo->nmics)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY+= TOMEMORY(modgeominfo->ncntr)*MEMsizeof(COMPLEX);
  ESTIMATED_MEMORY += (MEMsizeof(struct geometry_struct)); // allocating structure
  ESTIMATED_MEMORY += (MEMsizeof(struct file_info));
  ESTIMATED_MEMORY += (MEMsizeof(struct sphere_geom));
  ESTIMATED_MEMORY += (MEMsizeof(struct cylinder_geom));
  ESTIMATED_MEMORY += (MEMsizeof(struct nodes_geom));
  ESTIMATED_MEMORY += (MEMsizeof(struct plate_geom));

  return ESTIMATED_MEMORY;
}