TimeStepCriterion.h

Go to the documentation of this file.

#pragma once
 
 
#include "common/ForwardDecl.h"
#include "objects/containers/Array.h"
#include "objects/wrappers/AutoPtr.h"
#include "quantities/QuantityIds.h"
 
NAMESPACE_SPH_BEGIN
 
class IScheduler;
 
enum class CriterionId {
    INITIAL_VALUE, 
    MAXIMAL_VALUE, 
    DERIVATIVE,    
    CFL_CONDITION, 
    ACCELERATION,  
    DIVERGENCE,    
    MAX_CHANGE,    
};
 
std::ostream& operator<<(std::ostream& stream, const CriterionId id);
 
struct TimeStep {
    Float value;
 
    CriterionId id;
};
 
class ITimeStepCriterion : public Polymorphic {
public:
    virtual TimeStep compute(IScheduler& scheduler, Storage& storage, Float maxStep, Statistics& stats) = 0;
};
 
class DerivativeCriterion : public ITimeStepCriterion {
private:
    Float factor;
 
    Float power;
 
public:
    explicit DerivativeCriterion(const RunSettings& settings);
 
    virtual TimeStep compute(IScheduler& scheduler,
        Storage& storage,
        Float maxStep,
        Statistics& stats) override;
 
private:
    template <template <typename> class Tls>
    TimeStep computeImpl(IScheduler& scheduler, Storage& storage, Float maxStep, Statistics& stats);
};
 
class AccelerationCriterion : public ITimeStepCriterion {
private:
    Float factor;
 
public:
    explicit AccelerationCriterion(const RunSettings& settings);
 
    virtual TimeStep compute(IScheduler& scheduler,
        Storage& storage,
        Float maxStep,
        Statistics& stats) override;
};
 
class DivergenceCriterion : public ITimeStepCriterion {
private:
    Float factor;
 
public:
    explicit DivergenceCriterion(const RunSettings& settings);
 
    virtual TimeStep compute(IScheduler& scheduler,
        Storage& storage,
        Float maxStep,
        Statistics& stats) override;
};
 
class CourantCriterion : public ITimeStepCriterion {
private:
    Float courant;
 
public:
    explicit CourantCriterion(const RunSettings& settings);
 
    virtual TimeStep compute(IScheduler& scheduler,
        Storage& storage,
        Float maxStep,
        Statistics& stats) override;
};
 
 
class MultiCriterion : public ITimeStepCriterion {
private:
    Array<AutoPtr<ITimeStepCriterion>> criteria;
 
    Float maxChange;
    Float lastStep;
 
public:
    explicit MultiCriterion(const RunSettings& settings);
 
    MultiCriterion(Array<AutoPtr<ITimeStepCriterion>>&& criteria, const Float maxChange, const Float initial);
 
    virtual TimeStep compute(IScheduler& scheduler,
        Storage& storage,
        const Float maxStep,
        Statistics& stats) override;
};
 
NAMESPACE_SPH_END