formulaBase.hpp

#ifndef ABLATELIBRARY_FORMULABASE_HPP
#define ABLATELIBRARY_FORMULABASE_HPP
 
#include <muParser.h>
#include <random>
#include "mathFunction.hpp"
#include "parameters/parameters.hpp"
 
namespace ablate::mathFunctions {
 
class FormulaBase : public MathFunction {
   private:
    std::minstd_rand0 pseudoRandomEngine{0};
 
    std::default_random_engine randomEngine{0};
 
   protected:
    mutable double coordinate[3] = {0, 0, 0};
 
    mutable double time = 0.0;
 
    mu::Parser parser;
 
    const std::string formula;
 
    explicit FormulaBase(std::string functionString, const std::shared_ptr<ablate::parameters::Parameters>& constants);
 
    static std::invalid_argument ConvertToException(mu::Parser::exception_type& exception);
 
   public:
    FormulaBase(const FormulaBase&) = delete;
    void operator=(const FormulaBase&) = delete;
 
   private:
    static mu::value_type PowerFunction(mu::value_type a, mu::value_type b);
 
    static mu::value_type RandomFunction(void* data, mu::value_type lowerBound, mu::value_type upperBound);
 
    static mu::value_type PseudoRandomFunction(void* data, mu::value_type lowerBound, mu::value_type upperBound);
 
    static mu::value_type ModulusOperator(mu::value_type left, mu::value_type right);
};
 
}  // namespace ablate::mathFunctions
#endif  // ABLATELIBRARY_FORMULABASE_HPP