Method solve

src/oneD/Sim1D.cpp:374–420

}

void Sim1D::solve(int loglevel, bool refine_grid)
{
    int new_points = 1;
    m_attempt_counter = 0;
    finalize();
    if (loglevel > 6) {
        clearDebugFile();
    }

    while (new_points > 0) {
        SteadyStateSystem::solve(loglevel);
        if (loglevel > 0) {
            writelog("\nNewton steady-state solve succeeded.\n\n");
            writelog("Problem solved on [");
            for (size_t mm = 1; mm < nDomains(); mm+=2) {
                writelog("{}", domain(mm).nPoints());
                if (mm + 2 < nDomains()) {
                    writelog(", ");
                }
            }
            writelog("] point grid(s).\n");
            if (loglevel > 3) {
                show();
            }
        }
        if (m_steady_callback) {
            m_steady_callback->eval(0);
        }
        if (refine_grid) {
            new_points = refine(loglevel);
            writeDebugInfo("Regridding", "After regridding", loglevel,
                           m_attempt_counter);
        } else {
            debuglog("grid refinement disabled.\n", loglevel);
            new_points = 0;
        }
    }
    if (new_points < 0) {
        // If the solver finished after removing grid points, do one final evaluation
        // of the governing equations to update internal arrays in each domain that may
        // be used for data saved in output files.
        for (auto dom : m_dom) {
            dom->eval(npos, m_state->data() + dom->loc(), m_xnew.data() + dom->loc(),
                      m_mask.data());
        }
    }
}

int Sim1D::refine(int loglevel)
{