AdaptiveDMD.h

 
// Description: Computes the AdaptiveDMD algorithm on the given snapshot matrix.
//              The AdaptiveDMD algorithm should be used if the dt changes
//              between the samples. This algorithm uniformly interpolates the samples
//              that may have been taken with variable steps, using the constant step dt
//              (a prequisite of the DMD algorithm). The smaller dt is, the finer
//              the fidelity of the interpolation.
 
#ifndef included_AdaptiveDMD_h
#define included_AdaptiveDMD_h
 
#include "DMD.h"
#include <vector>
 
namespace CAROM {
 
class Vector;
 
class AdaptiveDMD : public DMD
{
public:
 
    AdaptiveDMD(int dim, double desired_dt = -1.0, std::string rbf = "G",
                std::string interp_method = "LS", double closest_rbf_val = 0.9,
                bool alt_output_basis = false,
                std::shared_ptr<Vector> state_offset = nullptr);
 
    ~AdaptiveDMD() { }
 
    void train(double energy_fraction, const Matrix* W0 = NULL,
               double linearity_tol = 0.0);
 
    void train(int k, const Matrix* W0 = NULL, double linearity_tol = 0.0);
 
    double getTrueDt() const;
 
    std::unique_ptr<const Matrix> getInterpolatedSnapshots();
 
private:
 
    AdaptiveDMD();
 
    AdaptiveDMD(
        const AdaptiveDMD& other);
 
    AdaptiveDMD&
    operator = (
        const AdaptiveDMD& rhs);
 
    std::string d_rbf;
 
    std::string d_interp_method;
 
    std::vector<std::shared_ptr<Vector>> d_interp_snapshots;
 
    void interpolateSnapshots();
 
    double d_closest_rbf_val;
};
 
}
 
#endif