14 MultiJac::MultiJac(OneDim& r)
15 : BandMatrix(r.size(),r.bandwidth(),r.bandwidth())
18 m_points = r.points();
21 m_ssdiag.resize(m_size);
22 m_mask.resize(m_size);
26 m_atol = sqrt(std::numeric_limits<double>::epsilon());
30 void MultiJac::updateTransient(doublereal rdt, integer* mask)
32 for (
size_t n = 0; n < m_size; n++) {
33 value(n,n) = m_ssdiag[n] - mask[n]*rdt;
37 void MultiJac::incrementDiagonal(
int j, doublereal d)
40 value(j,j) = m_ssdiag[j];
50 for (
size_t j = 0; j < m_points; j++) {
52 for (
size_t n = 0; n < nv; n++) {
54 double xsave = x0[ipt];
57 dx = xsave*m_rtol + m_atol;
59 dx = xsave*m_rtol - m_atol;
69 for (
size_t i = j - 1; i != j+2; i++) {
70 if (i !=
npos && i < m_points) {
73 for (
size_t m = 0; m < mv; m++) {
74 value(m+iloc,ipt) = (m_r1[m+iloc] - resid0[m+iloc])*rdx;
83 for (
size_t n = 0; n < m_size; n++) {
84 m_ssdiag[n] =
value(n,n);
87 m_elapsed += double(clock() - t0)/CLOCKS_PER_SEC;
void eval(size_t j, double *x, double *r, doublereal rdt=-1.0, int count=1)
Evaluate the multi-domain residual function.
const size_t npos
index returned by functions to indicate "no position"
size_t loc(size_t jg)
Location in the solution vector of the first component of global point jg.
OneDim * m_resid
Residual evaluator for this Jacobian.
doublereal & value(size_t i, size_t j)
Return a changeable reference to element (i,j).
void eval(doublereal *x0, doublereal *resid0, double rdt)
Evaluate the Jacobian at x0.
void bfill(doublereal v=0.0)
Fill or zero the matrix.
size_t nVars(size_t jg)
Number of solution components at global point jg.
Namespace for the Cantera kernel.