solutionErrorMonitor.hpp

#ifndef ABLATELIBRARY_SOLUTIONERRORMONITOR_HPP
#define ABLATELIBRARY_SOLUTIONERRORMONITOR_HPP
#include <iostream>
#include <monitors/logs/log.hpp>
#include <vector>
#include "monitor.hpp"
#include "utilities/mathUtilities.hpp"
 
namespace ablate::monitors {
 
class SolutionErrorMonitor : public Monitor {
   public:
    enum class Scope { VECTOR, COMPONENT };
 
   private:
    static PetscErrorCode MonitorError(TS ts, PetscInt step, PetscReal crtime, Vec u, void* ctx);
    Scope errorScope;
    ablate::utilities::MathUtilities::Norm normType;
    const std::shared_ptr<logs::Log> log;
 
   public:
    SolutionErrorMonitor(Scope errorScope, ablate::utilities::MathUtilities::Norm normType, const std::shared_ptr<logs::Log>& log = {});
 
    PetscMonitorFunction GetPetscFunction() override { return MonitorError; }
 
    std::vector<PetscReal> ComputeError(TS ts, PetscReal time, Vec u);
};
 
std::ostream& operator<<(std::ostream& os, const SolutionErrorMonitor::Scope& v);
std::istream& operator>>(std::istream& is, SolutionErrorMonitor::Scope& v);
}  // namespace ablate::monitors
 
#endif  // ABLATELIBRARY_ERRORMONITOR_HPP