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;