MultiJac::MultiJac(OneDim& r) : BandMatrix(r.size(),r.bandwidth(),r.bandwidth()) { m_size = r.size(); m_points = r.points(); m_resid = &r; m_r1.resize(m_size); m_ssdiag.resize(m_size); m_mask.resize(m_size); m_elapsed = 0.0; m_nevals = 0; m_age = 100000; doublereal ff = 1.0; while (1.0 + ff != 1.0) { ff *= 0.5; } m_atol = sqrt(ff); m_rtol = 1.0e-5; }