estimate_memory()

MEMORY estimate_memory

(

)

                         {
  MEMORY ESTIMATED_MEMORY = 0L;

  if (modsolinfo->pre_calculate_coefs) {
    ESTIMATED_MEMORY += (MEMsizeof(struct bodycoefs_struct));
    ESTIMATED_MEMORY += 2L
        * (TOMEMORY(modgeominfo->ncntr) + TOMEMORY(modgeominfo->nchief))
        * TOMEMORY(modgeominfo->ncntr) * MEMsizeof(DOUBLE); // B Body Coefs
    ESTIMATED_MEMORY +=
        TOMEMORY(4L*( TOMEMORY(modgeominfo->ncntr) + TOMEMORY(modgeominfo->nchief) ) * TOMEMORY(modgeominfo->ncntr) * MEMsizeof(DOUBLE)); // C Body Coefs

    ESTIMATED_MEMORY += TOMEMORY(modgeominfo->ncntr+modgeominfo->nchief)
        * MEMsizeof(DOUBLE); // temporary vector
    ESTIMATED_MEMORY += MEMsizeof(struct micscoefs_struct);
    ESTIMATED_MEMORY += 2L * TOMEMORY(modgeominfo->nmics)
        * TOMEMORY(modgeominfo->nelmb) * MEMsizeof(DOUBLE); // 2 mics B coefficients
    ESTIMATED_MEMORY += 3L * TOMEMORY(modgeominfo->nmics)
        * TOMEMORY(modgeominfo->nelmb) * 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
  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(modgeominfo->ncntr+modgeominfo->nchief)*TOMEMORY(modgeominfo->ncntr)*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+= (modsolinfo->nimped > 0 ) ? TOMEMORY(modsolinfo->nimped) * MEMsizeof(COMPLEX) : 0L;
  ESTIMATED_MEMORY +=
      (modsolinfo->nimped > 0) ?
          TOMEMORY(modsolinfo->nimped) * MEMsizeof(int) : 0L;
  ESTIMATED_MEMORY+= (modsolinfo->nradian > 0 ) ? TOMEMORY(modsolinfo->nradian) * MEMsizeof(COMPLEX) : 0L;
  ESTIMATED_MEMORY +=
      (modsolinfo->nradian > 0) ?
          TOMEMORY(modsolinfo->nradian) * MEMsizeof(int) : 0L;
  ESTIMATED_MEMORY += TOMEMORY(modsolinfo->nplanw)
      * MEMsizeof(struct acoustic_point);
  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(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;
}