fieldMonitor.hpp

#ifndef ABLATELIBRARY_FIELDMONITOR_HPP
#define ABLATELIBRARY_FIELDMONITOR_HPP
 
#include <petsc.h>
#include <vector>
#include "io/serializable.hpp"
#include "monitor.hpp"
 
namespace ablate::monitors {
 
class FieldMonitor : public Monitor, public io::Serializable {
   private:
    // Reuse the domain object to set up a domain to hold the monitor vector/dm
    std::shared_ptr<ablate::domain::Domain> monitorDomain = nullptr;
 
   protected:
    // Hold onto the subdomain for this monitor.  It should be over the entire monitorDomain
    std::shared_ptr<ablate::domain::SubDomain> monitorSubDomain = nullptr;
 
   public:
    [[nodiscard]] const std::string& GetId() const override { return monitorDomain->GetName(); }
 
    void Register(std::shared_ptr<solver::Solver> solverIn) override { Monitor::Register(solverIn); }
 
    void Register(std::string id, std::shared_ptr<solver::Solver> solverIn, std::vector<std::shared_ptr<domain::FieldDescriptor>> fields);
 
    PetscErrorCode Save(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
 
    PetscErrorCode Restore(PetscViewer viewer, PetscInt sequenceNumber, PetscReal time) override;
};
 
}  // namespace ablate::monitors
 
#endif  // ABLATELIBRARY_FIELDMONITOR_HPP