RootfindingIVP.h

Go to the documentation of this file.

#ifndef ROOTFINDINGIVP_H_
#define ROOTFINDINGIVP_H_
 
#include "MUQ/Modeling/ODEBase.h"
 
namespace muq {
  namespace Modeling {
    class RootfindingIVP : public ODEBase {
    public:
 
      RootfindingIVP(std::shared_ptr<WorkPiece> rhs, std::shared_ptr<WorkPiece> root, boost::property_tree::ptree const& pt, std::shared_ptr<muq::Utilities::AnyAlgebra> algebra = std::make_shared<muq::Utilities::AnyAlgebra>());
 
      virtual ~RootfindingIVP();
      
    private:
 
 
      virtual void EvaluateImpl(ref_vector<boost::any> const& inputs) override;
 
 
      virtual void JacobianImpl(unsigned int const wrtIn, unsigned int const wrtOut, ref_vector<boost::any> const& inputs) override;
 
 
      Eigen::VectorXi FindRoot(ref_vector<boost::any> const& inputs, int const wrtIn = -1, int const wrtOut = -1, DerivativeMode const& mode = DerivativeMode::Jac);
 
 
      void UpdateInputOutputTypes();
 
 
      static int EvaluateRoot(realtype t, N_Vector state, realtype *root, void *user_data);
 
      std::shared_ptr<WorkPiece> root;
 
      const unsigned int maxSteps;
 
      const double maxTime;
 
      const unsigned int maxErrorTests;
 
    };
  } // namespace Modeling
} // namespace muq
 
#endif