Method HypreMultiABec

Source/radiation/HypreMultiABec.cpp:613–687

}

HypreMultiABec::HypreMultiABec(int _crse_level, int _fine_level,
                               int _solver_flag)
  : crse_level(_crse_level), fine_level(_fine_level),
    solver_flag(_solver_flag),
    geom(fine_level+1),
    grids(fine_level+1),
    dmap(fine_level+1),
    fine_ratio(fine_level+1),
    bd(fine_level+1),
    subgrids(fine_level+1),
    acoefs(fine_level+1),
    bcoefs(fine_level+1),
    SPa(fine_level+1),
    cintrp(fine_level+1),
    ederiv(fine_level+1),
    c_cintrp(fine_level+1),
    c_ederiv(fine_level+1),
    c_entry(fine_level+1),
    hgrid(NULL), stencil(NULL), graph(NULL),
    A(NULL), A0(NULL), b(NULL), x(NULL),
    sstruct_solver(NULL), solver(NULL), precond(NULL)
{
  ParmParse pp("hmabec");

  verbose = 1;       pp.query("v", verbose); pp.query("verbose", verbose);
  verbose_threshold = 0; pp.query("verbose_threshold", verbose_threshold);
  bho = 0;           pp.query("bho", bho);
  use_subgrids = 0;  pp.query("use_subgrids", use_subgrids);

  static int first = 1;
  if (verbose >= 1 && first && ParallelDescriptor::IOProcessor()) {
    first = 0;
    std::cout << "hmabec.bho               = " << bho << std::endl;
    std::cout << "hmabec.use_subgrids      = " << use_subgrids << std::endl;
    std::cout << "hmabec.verbose           = " << verbose << std::endl;
    std::cout << "hmabec.verbose_threshold = " << verbose_threshold << std::endl;
  }

  if (solver_flag == 100 || solver_flag == 102 || solver_flag == 104 ||
      solver_flag == 105 ||
      solver_flag == 150 || solver_flag == 151 || solver_flag == 153 ||
      solver_flag == 1002) {
    ObjectType = HYPRE_PARCSR;
  }
  else if (solver_flag == 101 || solver_flag == 103 ||
           solver_flag == 106 || solver_flag == 107 ||
           solver_flag == 152 ||
           solver_flag == 1003) {
    ObjectType = HYPRE_SSTRUCT;
  }
  else if (solver_flag == 108 || solver_flag == 109) {
    ObjectType = HYPRE_STRUCT;
  }
  else {
    std::cout << "HypreMultiABec: no such solver" << std::endl;
    exit(1);
  }