gravity.hpp

#ifndef ABLATELIBRARY_GRAVITY_HPP
#define ABLATELIBRARY_GRAVITY_HPP
 
#include "flowProcess.hpp"
 
namespace ablate::finiteVolume::processes {
 
class Gravity : public FlowProcess {
   private:
    const std::vector<PetscReal> gravityVector;
 
    static PetscErrorCode ComputeGravitySource(PetscInt dim, PetscReal time, const PetscFVCellGeom* cg, const PetscInt uOff[], const PetscScalar u[], const PetscInt aOff[], const PetscScalar a[],
                                               PetscScalar f[], void* ctx);
 
   public:
    explicit Gravity(std::vector<double> gravityVector);
    void Setup(ablate::finiteVolume::FiniteVolumeSolver& fv) override;
};
 
}  // namespace ablate::finiteVolume::processes
#endif  // ABLATELIBRARY_GRAVITY_HPP