10#include "cantera/oneD/refine.h"
33 if (!sol || !(sol->thermo())) {
35 "Missing or incomplete Solution object.");
48 if (nv != m_nv || !m_refiner) {
50 m_refiner = make_unique<Refiner>(*
this);
53 m_name.resize(m_nv,
"");
54 m_max.resize(m_nv, 0.0);
55 m_min.resize(m_nv, 0.0);
57 m_rtol_ss.resize(m_nv, 1.0e-4);
58 m_atol_ss.resize(m_nv, 1.0e-9);
59 m_rtol_ts.resize(m_nv, 1.0e-4);
60 m_atol_ts.resize(m_nv, 1.0e-11);
63 m_slast.resize(m_nv *
m_points, 0.0);
69 if (m_name[n] !=
"") {
72 return fmt::format(
"component {}", n);
84 "no component named "+name);
90 for (n = 0; n < m_nv; n++) {
103 for (n = 0; n < m_nv; n++) {
123 auto wrap_tols = [
this](
const vector<double>& tols) {
126 set<double> unique_tols(tols.begin(), tols.end());
127 if (unique_tols.size() == 1) {
131 for (
size_t i = 0; i < tols.size(); i++) {
138 state[
"type"] =
type();
139 state[
"points"] =
static_cast<long int>(
nPoints());
141 state[
"tolerances"][
"transient-abstol"] = wrap_tols(m_atol_ts);
142 state[
"tolerances"][
"steady-abstol"] = wrap_tols(m_atol_ss);
143 state[
"tolerances"][
"transient-reltol"] = wrap_tols(m_rtol_ts);
144 state[
"tolerances"][
"steady-reltol"] = wrap_tols(m_rtol_ss);
152 "Domain needs to be installed in a container before calling asArray.");
165 "Domain needs to be installed in a container before calling fromArray.");
175 auto set_tols = [&](
const AnyValue& tols,
const string& which, vector<double>& out)
177 if (!tols.
hasKey(which)) {
180 const auto& tol = tols[which];
181 if (tol.isScalar()) {
186 if (tol.hasKey(name)) {
187 out[i] = tol[name].asDouble();
189 warn_user(
"Domain1D::setMeta",
"No {} found for component '{}'",
196 if (meta.
hasKey(
"tolerances")) {
197 const auto& tols = meta[
"tolerances"];
198 set_tols(tols,
"transient-abstol", m_atol_ts);
199 set_tols(tols,
"transient-reltol", m_rtol_ts);
200 set_tols(tols,
"steady-abstol", m_atol_ss);
201 set_tols(tols,
"steady-reltol", m_rtol_ss);
229 for (
size_t j = 0; j <
m_points; j++) {
238 for (
size_t i = 0; i < nn; i++) {
239 writeline(
'-', 79,
false,
true);
241 for (
size_t n = 0; n < 5; n++) {
244 writeline(
'-', 79,
false,
true);
245 for (
size_t j = 0; j <
m_points; j++) {
247 for (
size_t n = 0; n < 5; n++) {
248 double v = value(x, i*5+n, j);
254 size_t nrem = m_nv - 5*nn;
255 writeline(
'-', 79,
false,
true);
257 for (
size_t n = 0; n < nrem; n++) {
260 writeline(
'-', 79,
false,
true);
261 for (
size_t j = 0; j <
m_points; j++) {
263 for (
size_t n = 0; n < nrem; n++) {
264 double v = value(x, nn*5+n, j);
271void Domain1D::setProfile(
const string& name,
double* values,
double* soln)
273 for (
size_t n = 0; n < m_nv; n++) {
275 for (
size_t j = 0; j <
m_points; j++) {
276 soln[index(n, j) +
m_iloc] = values[j];
281 throw CanteraError(
"Domain1D::setProfile",
"unknown component: "+name);
286 for (
size_t j = 0; j <
m_points; j++) {
287 for (
size_t n = 0; n < m_nv; n++) {
Header file for class ThermoPhase, the base class for phases with thermodynamic properties,...
A map of string keys to values whose type can vary at runtime.
bool hasKey(const string &key) const
Returns true if the map contains an item named key.
A wrapper for a variable whose type is determined at runtime.
bool hasKey(const string &key) const
Returns true if this AnyValue is an AnyMap and that map contains a key with the given name.
Base class for exceptions thrown by Cantera classes.
void setTransientTolerances(double rtol, double atol, size_t n=npos)
Set tolerances for time-stepping mode.
size_t lastPoint() const
The index of the last (that is, right-most) grid point belonging to this domain.
size_t m_iloc
Starting location within the solution vector for unknowns that correspond to this domain.
shared_ptr< Solution > m_solution
Composite thermo/kinetics/transport handler.
size_t nComponents() const
Number of components at each grid point.
double rtol(size_t n)
Relative tolerance of the nth component.
virtual void setMeta(const AnyMap &meta)
Retrieve meta data.
virtual void _finalize(const double *x)
In some cases, a domain may need to set parameters that depend on the initial solution estimate.
size_t nPoints() const
Number of grid points in this domain.
shared_ptr< vector< double > > m_state
data pointer shared from OneDim
virtual void resize(size_t nv, size_t np)
Resize the domain to have nv components and np grid points.
void setSteadyTolerances(double rtol, double atol, size_t n=npos)
Set tolerances for steady-state mode.
virtual shared_ptr< SolutionArray > asArray(const double *soln) const
Save the state of this domain as a SolutionArray.
virtual string componentName(size_t n) const
Name of the nth component. May be overloaded.
void setSolution(shared_ptr< Solution > sol)
Set the solution manager.
double atol(size_t n)
Absolute tolerance of the nth component.
shared_ptr< SolutionArray > toArray(bool normalize=false) const
Save the state of this domain to a SolutionArray.
size_t m_points
Number of grid points.
string type() const
String indicating the domain implemented.
virtual double initialValue(size_t n, size_t j)
Initial value of solution component n at grid point j.
virtual size_t componentIndex(const string &name) const
index of component with name name.
virtual void fromArray(SolutionArray &arr, double *soln)
Restore the solution for this domain from a SolutionArray.
void needJacUpdate()
Set this if something has changed in the governing equations (for example, the value of a constant ha...
virtual void _getInitialSoln(double *x)
Writes some or all initial solution values into the global solution array, beginning at the location ...
Domain1D(size_t nv=1, size_t points=1, double time=0.0)
Constructor.
virtual size_t loc(size_t j=0) const
Location of the start of the local solution vector in the global solution vector,.
void locate()
Find the index of the first grid point in this domain, and the start of its variables in the global s...
virtual AnyMap getMeta() const
Retrieve meta data.
virtual void show(std::ostream &s, const double *x)
Print the solution.
virtual void setupGrid(size_t n, const double *z)
called to set up initial grid, and after grid refinement
void setAge(int age)
Set the Jacobian age.
An error indicating that an unimplemented function has been called.
virtual void resize()
Call after one or more grids has changed size, for example after being refined.
void saveStats()
Save statistics on function and Jacobian evaluation, and reset the counters.
MultiJac & jacobian()
Return a reference to the Jacobian evaluator of an OneDim object.
void writelog(const string &fmt, const Args &... args)
Write a formatted message to the screen.
void warn_user(const string &method, const string &msg, const Args &... args)
Print a user warning raised from method as CanteraWarning.
Namespace for the Cantera kernel.
const size_t npos
index returned by functions to indicate "no position"