File resolutionpreserving.h

File List > code_source > cosmolattice > include > CosmoInterface > definitions > defectsmodule > resolutionpreserving.h

Go to the documentation of this file

#ifndef COSMOINTERFACE_DEFINITIONS_DEFECTSMODULE_RESOLUTIONPRESERVING_H
#define COSMOINTERFACE_DEFINITIONS_DEFECTSMODULE_RESOLUTIONPRESERVING_H

/* This file is part of CosmoLattice, available at www.cosmolattice.net .
 *   Copyright Daniel G. Figueroa, Adrien Florio, Francisco Torrenti and Wessel Valkenburg.
 *   Released under the MIT license, see LICENSE.md. */

// File info: Main contributor(s): Jorge Baeza-Ballesteros,  Year: 2026


#include "TempLat/util/exception.h"
#include "CosmoInterface/evolvers/evolver.h"
#include "TempLat/util/almostequal.h"

namespace TempLat
{

    class ResolutionPreserving
    {
    public:
        // Put public methods here. These should change very little over time.

        ResolutionPreserving() = delete;

        template<typename Model, typename T>
        static void updateResolutionPreservingFactor(Model &model, T tminust0, T t0, T tRP0, T tRPMax, T sRP)
        {
            if(tminust0 + t0 < tRPMax && tminust0 + t0 >= tRP0 - model.dt) {
                model.resolutionPreservingFactor *= pow(model.aI / model.aIM, 2. * (sRP - 1.));
            }
        }

        template<typename Model, typename T>
        static bool updateTimeStep(Model &model, RunParameters<T> &rPar, T t)
        {
            if (t > (rPar.tRPMax - rPar.dt) && t < rPar.tRPMax && !AlmostEqual(rPar.tRPMax, t) && !AlmostEqual(rPar.tRPMax, t + rPar.dt)) {
                model.dt = rPar.tRPMax - t;
                return true;
            }
            else if ( AlmostEqual(rPar.tRPMax, t) && !AlmostEqual(model.dt, rPar.dt) ) model.dt = rPar.dt - model.dt;
            else if ( !AlmostEqual(model.dt, rPar.dt) ) model.dt = rPar.dt;
            return false;
        }

    };


} // namespace TempLat

#endif