parsedSeries.hpp

#ifndef ABLATELIBRARY_PARSEDSERIES_HPP
#define ABLATELIBRARY_PARSEDSERIES_HPP
#include <muParser.h>
#include "formulaBase.hpp"
#include "parameters/parameters.hpp"
 
namespace ablate::mathFunctions {
class ParsedSeries : public FormulaBase {
   private:
    // The local count for the series
    mutable double i = 0;
 
    const int lowerBound;
 
    const int upperBound;
 
   private:
    static PetscErrorCode ParsedPetscSeries(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nf, PetscScalar* u, void* ctx);
 
   public:
    ParsedSeries(const ParsedSeries&) = delete;
    void operator=(const ParsedSeries&) = delete;
 
    explicit ParsedSeries(std::string functionString, int lowerBound = 1, int upperBound = 1000, const std::shared_ptr<ablate::parameters::Parameters>& constants = {});
 
    double Eval(const double& x, const double& y, const double& z, const double& t) const override;
 
    double Eval(const double* xyz, const int& ndims, const double& t) const override;
 
    void Eval(const double& x, const double& y, const double& z, const double& t, std::vector<double>& result) const override;
 
    void Eval(const double* xyz, const int& ndims, const double& t, std::vector<double>& result) const override;
 
    void* GetContext() override { return this; }
 
    PetscFunction GetPetscFunction() override { return ParsedPetscSeries; }
};
}  // namespace ablate::mathFunctions
 
#endif  // ABLATELIBRARY_PARSEDFUNCTION_HPP