BDF.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/time_integrator/ImplicitTimeIntegrator.hpp>
 
namespace polyfem::time_integrator
{
    class BDF : public ImplicitTimeIntegrator
    {
    public:
        BDF(const int order = 1);
 
        void set_parameters(const json &params) override;
 
        void update_quantities(const Eigen::VectorXd &x) override;
 
        Eigen::VectorXd x_tilde() const override;
 
        Eigen::VectorXd compute_velocity(const Eigen::VectorXd &x) const override;
 
        Eigen::VectorXd compute_acceleration(const Eigen::VectorXd &v) const override;
 
        double acceleration_scaling() const override;
 
        double dv_dx(const unsigned prev_ti = 0) const override;
 
        double beta_dt() const;
 
        Eigen::VectorXd weighted_sum_x_prevs() const;
 
        Eigen::VectorXd weighted_sum_v_prevs() const;
 
        static const std::vector<double> &alphas(const int i);
 
        static double betas(const int i);
 
    protected:
        int max_steps() const override { return max_steps_; }
 
        int max_steps_ = 1;
    };
} // namespace polyfem::time_integrator