solver.h

Go to the documentation of this file.

 
#pragma once
 
#include "eddy.h"
 
class domain;
 
using namespace std;
 
 
class solver {
 
 
 
    public:
 
        domain         *domn;          
 
        eddy           *ed3;           
        domain         *eddl3;         
 
        double         time;           
        double         t0;             
        double         dtSmean;        
        double         dtCUmax;        
 
        bool           LeddyAccepted;  
        int            iEtrials;       
 
        double         PaSum;          
        int            nPaSum;         
        int            neddies;        
        double         PaSumC;         
        int            nPaSumC;        
 
 
    public:
 
        virtual void calculateSolution();
 
    private:
 
        bool   sampleEddyAndImplementIfAccepted();
        void   computeDtSmean();
        void   computeDtCUmax();
        double sampleDt();
        void   diffusionCatchUpIfNeeded(bool Ldoit=false);
        void   raiseDtSmean();
        void   lowerDtSmean();
        bool   testLES_elapsedTime(const double time, const double tauEddy);
        bool   testLES_fracDomain( const double eSize);
        bool   testLES_integralLength(const double time, const double eSize);
        bool   testLES_thirds();
 
 
 
    public:
 
        solver(){ ed3=0; eddl3=0; }
        virtual void init(domain *p_domn);
        virtual ~solver();
 
};