eschnett/cactusamrex - Change QGS6SF4SPCD2A7EINWQPICO4YXIM5YL5KIENFPF6KG77AB7KRJSQC

ODESolvers: Don't special-case first RHS evaluation

Created by Erik Schnetter on July 4, 2020

QGS6SF4SPCD2A7EINWQPICO4YXIM5YL5KIENFPF6KG77AB7KRJSQC

Dependencies

In channels

main

Change contents

Deletion in ODESolvers/src/solve.cxx at line 263 [6.541]

B:BD[5.191] → [5.191:250]

∅:D[5.250] → [7.239:264]

B:BD[6.4334] → [7.239:264]

B:BD[7.264] → [2.260:348]

∅:D[2.348] → [7.335:382]

B:BD[7.335] → [7.335:382]

  *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
  // Calculate first RHS
  if (verbose)
    CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
  CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Insertion in ODESolvers/src/solve.cxx at line 274 [6.541]

[8.172]

[6.4491]


    // Calculate first RHS
    *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
    if (verbose)
      CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
    CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Insertion in ODESolvers/src/solve.cxx at line 291 [6.541]
[6.4645]
[6.4645]
```
    *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
```

Insertion in ODESolvers/src/solve.cxx at line 293 [6.541]

[6.4646]

[6.4696]

    // Step 1
    if (verbose)
      CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
    CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Insertion in ODESolvers/src/solve.cxx at line 303 [6.541]
[6.4836]
[2.349]
```
    // Step 2
```

Insertion in ODESolvers/src/solve.cxx at line 319 [6.541]

[3.243]

    *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
    // Step 1
    if (verbose)
      CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
    CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Insertion in ODESolvers/src/solve.cxx at line 367 [6.541]

[8.447]

    *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
    // Step 1
    if (verbose)
      CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
    CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Insertion in ODESolvers/src/solve.cxx at line 416 [6.541]

[5.377]

    *const_cast<CCTK_REAL *>(&cctkGH->cctk_time) = old_time;
    // Step 1
    if (verbose)
      CCTK_VINFO("Calculating RHS #1 at t=%g", double(cctkGH->cctk_time));
    CallScheduleGroup(cctkGH, "ODESolvers_RHS");

Replacement in ODESolvers/src/solve.cxx at line 534 [6.541]

B:BD[4.3847] → [4.3847:3962]

    // Special-case step 1 (why?)
    ks.push_back(rhs.copy());
    for (size_t step = 1; step < nsteps; ++step) {

[4.3847]

[4.3962]

    for (size_t step = 0; step < nsteps; ++step) {

Replacement in ODESolvers/src/solve.cxx at line 626 [6.541]

B:BD[4.8273] → [4.8273:8388]

    // Special-case step 1 (why?)
    ks.push_back(rhs.copy());
    for (size_t step = 1; step < nsteps; ++step) {

[4.8273]

[4.8388]

    for (size_t step = 0; step < nsteps; ++step) {

Replacement in ODESolvers/src/solve.cxx at line 665 [6.541]
B:BD[6.5163] → [4.9339:9364]
```
} // namespace ODESolver
```
[6.5163]
[6.5165]
```
}
```