Method save_derivs

src/solver/solver.cxx:1269–1297

}

void Solver::save_derivs(BoutReal* dudata) {
  // Make sure vectors in correct basis
  for (const auto& v : v2d) {
    if (v.covariant) {
      v.F_var->toCovariant();
    } else {
      v.F_var->toContravariant();
    }
  }
  for (const auto& v : v3d) {
    if (v.covariant) {
      v.F_var->toCovariant();
    } else {
      v.F_var->toContravariant();
    }
  }

  // Make sure 3D fields are at the correct cell location
  for (const auto& f : f3d) {
    if (f.var->getLocation() != (f.F_var)->getLocation()) {
      throw BoutException(_("Time derivative at wrong location - Field is at {:s}, "
                            "derivative is at {:s} for field '{:s}'\n"),
                          toString(f.var->getLocation()),
                          toString(f.F_var->getLocation()), f.name);
    }
  }

  loop_vars(dudata, SOLVER_VAR_OP::SAVE_DERIVS);
}

void Solver::set_id(BoutReal* udata) { loop_vars(udata, SOLVER_VAR_OP::SET_ID); }