relaxation_methods.jl

abstract type AbstractRelaxationTimeIntegrationAlgorithm <: AbstractTimeIntegrationAlgorithm end

abstract type RelaxationIntegrator <: AbstractTimeIntegrator end

get_tmp_cache(integrator::RelaxationIntegrator) = (integrator.u_tmp,)

# some algorithms from DiffEq like FSAL-ones need to be informed when a callback has modified u
u_modified!(integrator::RelaxationIntegrator, ::Bool) = false

# stop the time integration
function terminate!(integrator::RelaxationIntegrator)
    integrator.finalstep = true
    empty!(integrator.opts.tstops)

    return nothing
end

@inline function update_t_relaxation!(integrator::RelaxationIntegrator)
    # Check if due to entropy relaxation the final time would not be reached
    if integrator.finalstep == true && integrator.gamma != 1
        integrator.gamma = 1
    end
    integrator.t += integrator.gamma * integrator.dt

    return nothing
end

include("entropy_relaxation.jl")
include("methods_subdiagonal.jl")
include("methods_vanderHouwen.jl")