debugBoundarySolver.hpp

#ifndef ABLATELIBRARY_DEBUGBOUNDARYSOLVER_HPP
#define ABLATELIBRARY_DEBUGBOUNDARYSOLVER_HPP
 
#include "boundarySolver.hpp"
 
namespace ablate::boundarySolver {
 
class DebugBoundarySolver : public BoundarySolver {
   public:
    DebugBoundarySolver(std::string solverId, std::shared_ptr<domain::Region> region, std::shared_ptr<domain::Region> fieldBoundary, std::vector<std::shared_ptr<BoundaryProcess>> boundaryProcesses,
                        std::shared_ptr<parameters::Parameters> options, bool mergeFaces = false);
 
    void Setup() override;
 
    PetscErrorCode ComputeRHSFunction(PetscReal time, Vec locXVec, Vec locFVec) override;
 
   private:
    void OutputStencilCellLocation(std::ostream& stream, PetscInt cell);
};
 
}  // namespace ablate::boundarySolver
#endif  // ABLATELIBRARY_DEBUGBOUNDARYSOLVER_HPP