Method loadMatrix

Source/radiation/HypreExtMultiABec.cpp:148–837

}

void HypreExtMultiABec::loadMatrix()
{
  BL_PROFILE("HypreExtMultiABec::loadMatrix");

  HypreMultiABec::loadMatrix();

  if (0 && verbose >= 1 && ParallelDescriptor::IOProcessor()) {
    std::cout << "In HypreExtMultiABec::loadMatrix(), the multipliers are:" << std::endl;
    std::cout << "  HypreExtMultiABec::alpha2 = " << alpha2 << std::endl;
    std::cout << "  HypreExtMultiABec::gamma  = " << gamma  << std::endl;
    std::cout << "  HypreExtMultiABec::delta1 = " << delta1 << std::endl;
    std::cout << "  HypreExtMultiABec::delta2 = " << delta2 << std::endl;
  }

  // These really ought to be members and already defined:
#if (AMREX_SPACEDIM == 1)
  // if we were really 1D:
/*
  int offsets[3][1] = {{ 0}
                       {-1},
                       { 1}};
*/
  // fake 1D as a 2D problem:
  int offsets[3][2] = {{ 0,  0},
                       {-1,  0},
                       { 1,  0}};
#elif (AMREX_SPACEDIM == 2)
  int offsets[5][2] = {{ 0,  0},
                       {-1,  0},
                       { 1,  0},
                       { 0, -1},
                       { 0,  1}};
#elif (AMREX_SPACEDIM == 3)
  int offsets[7][3] = {{ 0,  0,  0},
                       {-1,  0,  0},
                       { 1,  0,  0},
                       { 0, -1,  0},
                       { 0,  1,  0},
                       { 0,  0, -1},
                       { 0,  0,  1}};
#endif

  const int size = 2 * AMREX_SPACEDIM + 1;

  int stencil_indices[size];

  for (int i = 0; i < size; i++) {
    stencil_indices[i] = i;
  }

  FArrayBox matfab;
  FArrayBox mat_tmpfab;
  for (int level = crse_level; level <= fine_level; level++) {
    int part = level - crse_level;

    for (MFIter mfi(*acoefs[level]); mfi.isValid(); ++mfi) {
      int i = mfi.index();
      const Box &reg = grids[level][i];