arbitraryEulerianSource.hpp

#ifndef ABLATELIBRARY_ARBITRARYEULERIANSOURCE_HPP
#define ABLATELIBRARY_ARBITRARYEULERIANSOURCE_HPP
 
#include <string>
#include "coupledProcess.hpp"
 
namespace ablate::particles::processes {
 
class ArbitraryEulerianSource : public CoupledProcess {
   private:
    const std::string coupledFieldName;
 
    const std::shared_ptr<mathFunctions::MathFunction> sourceFunction;
 
   public:
    ArbitraryEulerianSource(std::string coupledFieldName, std::shared_ptr<mathFunctions::MathFunction> sourceFunction);
 
    void ComputeRHS(PetscReal time, accessors::SwarmAccessor& swarmAccessor, accessors::RhsAccessor& rhsAccessor, accessors::EulerianAccessor& eulerianAccessor) override {}
 
    void ComputeEulerianSource(PetscReal startTime, PetscReal endTime, accessors::SwarmAccessor& swarmAccessorPreStep, accessors::SwarmAccessor& swarmAccessorPostStep,
                               accessors::EulerianSourceAccessor& eulerianSourceAccessor) override;
};
 
}  // namespace ablate::particles::processes
#endif  // ABLATELIBRARY_INERTIAL_HPP