eschnett/cactusamrex - Change 4LCLQYKMQSLNHY47BPLR5BULPHFPN3KR3RPK2YIORUKFWIVKJOXQC

CarpetX: More run-time options for linear solver

Created by Erik Schnetter on June 9, 2020

4LCLQYKMQSLNHY47BPLR5BULPHFPN3KR3RPK2YIORUKFWIVKJOXQC

Dependencies

In channels

master

main

Change contents

Replacement in CarpetX/interface.ccl at line 32 [3.17976]

B:BD[2.72] → [2.72:94]

  CCTK_INT IN gi_res)

[2.72]

[4.29]

  CCTK_INT IN gi_res,
  CCTK_REAL IN reltol,
  CCTK_REAL IN abstol,
  CCTK_REAL OUT res_initial,
  CCTK_REAL OUT res_final)

Insertion in CarpetX/src/linsolve.cxx at line 11 [4.115]
[4.294]
[4.294]
```
#include <cmath>
```

Replacement in CarpetX/src/linsolve.cxx at line 19 [4.115]

B:BD[2.215] → [2.215:277]

                                     const CCTK_INT gi_res) {

[2.215]

[2.277]

                                     const CCTK_INT gi_res,
                                     const CCTK_REAL reltol,
                                     const CCTK_REAL abstol,
                                     CCTK_REAL *restrict const res_initial,
                                     CCTK_REAL *restrict const res_final) {

Replacement in CarpetX/src/linsolve.cxx at line 87 [4.115]
B:BD[4.2682] → [2.602:618]
```
  if (have_res)
```
[4.2682]
[2.618]
```
  if (have_res) {
```

Insertion in CarpetX/src/linsolve.cxx at line 89 [4.115]

[2.659]

[5.590]

    *res_initial = 0;
    for (int level = 0; level < int(ghext->leveldata.size()); ++level)
      *res_initial =
          fmax(*res_initial, ress.at(level)->norminf(vi, 0, false, true));
  } else {
    *res_initial = NAN;
  }

Replacement in CarpetX/src/linsolve.cxx at line 112 [4.115]

B:BD[4.3271] → [4.3271:3335]

B:BD[4.3335] → [6.687:750]

  const CCTK_REAL rtol = 0.0;
  const CCTK_REAL atol = 1.0e-12;
  const CCTK_REAL maxerr = mlmg.solve(sols, rhss, rtol, atol);

[4.3271]

[6.750]

  // const CCTK_REAL rtol = 0.0;
  // const CCTK_REAL atol = 1.0e-12;
  const CCTK_REAL maxerr = mlmg.solve(sols, rhss, reltol, abstol);

Replacement in CarpetX/src/linsolve.cxx at line 118 [4.115]
B:BD[6.835] → [2.875:891]
```
  if (have_res)
```
[6.835]
[2.891]
```
  if (have_res) {
```

Insertion in CarpetX/src/linsolve.cxx at line 120 [4.115]

[2.932]

[6.874]

    *res_final = 0;
    for (int level = 0; level < int(ghext->leveldata.size()); ++level)
      *res_final =
          fmax(*res_final, ress.at(level)->norminf(vi, 0, false, true));
  } else {
    *res_final = NAN;
  }