finiteElementSolver.hpp

#ifndef ABLATELIBRARY_FINITEELEMENTSOLVER_HPP
#define ABLATELIBRARY_FINITEELEMENTSOLVER_HPP
 
#include <solver/timeStepper.hpp>
#include <string>
#include <vector>
#include "boundaryConditions/boundaryCondition.hpp"
#include "eos/eos.hpp"
#include "mathFunctions/fieldFunction.hpp"
#include "solver/solver.hpp"
 
namespace ablate::finiteElement {
 
class FiniteElementSolver : public solver::Solver, public solver::IFunction, public solver::BoundaryFunction {
   private:
    const std::vector<std::shared_ptr<boundaryConditions::BoundaryCondition>> boundaryConditions;
    const std::vector<std::shared_ptr<mathFunctions::FieldFunction>> auxiliaryFieldsUpdaters;
 
   public:
    FiniteElementSolver(std::string solverId, std::shared_ptr<domain::Region> region, std::shared_ptr<parameters::Parameters> options,
                        std::vector<std::shared_ptr<boundaryConditions::BoundaryCondition>> boundaryConditions, std::vector<std::shared_ptr<mathFunctions::FieldFunction>> auxiliaryFields);
 
    void Register(std::shared_ptr<ablate::domain::SubDomain> subDomain) override;
    void Setup() override;
    void Initialize() override;
 
    virtual void CompleteFlowInitialization(DM, Vec) = 0;
 
    static void UpdateAuxFields(TS ts, FiniteElementSolver& fe);
 
    PetscErrorCode ComputeIFunction(PetscReal time, Vec locX, Vec locX_t, Vec locF) override;
    PetscErrorCode ComputeIJacobian(PetscReal time, Vec locX, Vec locX_t, PetscReal X_tShift, Mat Jac, Mat JacP) override;
    PetscErrorCode ComputeBoundary(PetscReal time, Vec locX, Vec locX_t) override;
};
}  // namespace ablate::finiteElement
 
#endif  // ABLATELIBRARY_FINITEVOLUME_HPP