arbitrarySource.hpp

#ifndef ABLATELIBRARY_ARBITRARYSOURCE_HPP
#define ABLATELIBRARY_ARBITRARYSOURCE_HPP
 
#include "process.hpp"
namespace ablate::finiteVolume::processes {
class ArbitrarySource : public Process {
    const std::map<std::string, std::shared_ptr<ablate::mathFunctions::MathFunction>> functions;
 
    static PetscErrorCode ComputeArbitrarySource(PetscInt dim, PetscReal time, const PetscFVCellGeom* cg, const PetscInt uOff[], const PetscScalar u[], const PetscInt aOff[], const PetscScalar a[],
                                                 PetscScalar f[], void* ctx);
 
    struct PetscFunctionStruct {
        mathFunctions::PetscFunction petscFunction;
        void* petscContext;
        PetscInt fieldSize;
    };
 
    std::vector<PetscFunctionStruct> petscFunctions;
 
   public:
    explicit ArbitrarySource(std::map<std::string, std::shared_ptr<ablate::mathFunctions::MathFunction>> functions);
 
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& fvmSolver) override;
};
 
}  // namespace ablate::finiteVolume::processes
 
#endif  // ABLATELIBRARY_ARBITRARYSOURCE_HPP