eschnett/cactusamrex - Change 5A7OHW6V23YT7O3PJKS7J6IKPLXIMCRTW6VIWPZKYMEFRYQX5VZQC

CarpetX: More work on linear solver

Created by Erik Schnetter on June 8, 2020

5A7OHW6V23YT7O3PJKS7J6IKPLXIMCRTW6VIWPZKYMEFRYQX5VZQC

Dependencies

In channels

master

main

Change contents

Insertion in CarpetX/src/linsolve.cxx at line 30 [3.115]
[3.906]
[3.906]
```
  // TODO: Look at Geometry for this
```

Replacement in CarpetX/src/linsolve.cxx at line 32 [3.115]

B:BD[3.937] → [2.13227:13438]

      {amrex::LinOpBCType::Periodic, amrex::LinOpBCType::Periodic,
       amrex::LinOpBCType::Periodic},
      {amrex::LinOpBCType::Periodic, amrex::LinOpBCType::Periodic,
       amrex::LinOpBCType::Periodic});

[3.937]

[3.1092]

      {amrex::LinOpBCType::Dirichlet, amrex::LinOpBCType::Dirichlet,
       amrex::LinOpBCType::Dirichlet},
      {amrex::LinOpBCType::Dirichlet, amrex::LinOpBCType::Dirichlet,
       amrex::LinOpBCType::Dirichlet});

Deletion in CarpetX/src/linsolve.cxx at line 76 [3.115]

B:BD[3.2536] → [3.2536:2669]


  // Initial guess
  for (int level = 0; level < int(ghext->leveldata.size()); ++level) {
    sols.at(level)->setVal(0.0, vi);
  }

Insertion in CarpetX/src/linsolve.cxx at line 82 [3.115]

[4.615]

    CCTK_VINFO("Before solving:");
    for (int i = -1; i < 10; ++i) {
      const int j = 4, k = 4;
      CCTK_VINFO("phi[%d,%d,%d]=%g", i, j, k,
                 double(sols.at(0)->array(0)(i, j, k, vi)));
    }

Replacement in CarpetX/src/linsolve.cxx at line 90 [3.115]

B:BD[4.734] → [4.734:804]

                 double(rhss.at(level)->norminf(0, 0, false, true)));

[4.734]

[4.804]

                 double(rhss.at(level)->norminf(vi, 0, false, true)));

Replacement in CarpetX/src/linsolve.cxx at line 93 [3.115]

B:BD[4.923] → [4.923:993]

                 double(sols.at(level)->norminf(0, 0, false, true)));

[4.923]

[4.993]

                 double(sols.at(level)->norminf(vi, 0, false, true)));

Replacement in CarpetX/src/linsolve.cxx at line 96 [3.115]

B:BD[4.1112] → [4.1112:1182]

                 double(ress.at(level)->norminf(0, 0, false, true)));

[4.1112]

[4.1182]

                 double(ress.at(level)->norminf(vi, 0, false, true)));

Replacement in CarpetX/src/linsolve.cxx at line 101 [3.115]

B:BD[3.3335] → [3.3335:3373]

  mlmg.solve(sols, rhss, rtol, atol);

[3.3335]

[3.3373]

  const CCTK_REAL maxerr = mlmg.solve(sols, rhss, rtol, atol);
#pragma omp critical
  CCTK_VINFO("Solution error (norm_inf): %g", double(maxerr));
  mlmg.compResidual(ress, sols, rhss);
#pragma omp critical
  {
    CCTK_VINFO("After solving:");
    for (int i = -1; i < 10; ++i) {
      const int j = 4, k = 4;
      CCTK_VINFO("phi[%d,%d,%d]=%g", i, j, k,
                 double(sols.at(0)->array(0)(i, j, k, vi)));
    }
    for (int level = 0; level < int(ghext->leveldata.size()); ++level)
      CCTK_VINFO("norm_inf rhs[%d]: %g", level,
                 double(rhss.at(level)->norminf(vi, 0, false, true)));
    for (int level = 0; level < int(ghext->leveldata.size()); ++level)
      CCTK_VINFO("norm_inf sol[%d]: %g", level,
                 double(sols.at(level)->norminf(vi, 0, false, true)));
    for (int level = 0; level < int(ghext->leveldata.size()); ++level)
      CCTK_VINFO("norm_inf res[%d]: %g", level,
                 double(ress.at(level)->norminf(vi, 0, false, true)));
  }