ImplicitNewmark.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/time_integrator/ImplicitTimeIntegrator.hpp>
 
namespace polyfem::time_integrator
{
    class ImplicitNewmark : public ImplicitTimeIntegrator
    {
    public:
        ImplicitNewmark() {}
 
        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 da_dx(const unsigned prev_ti = 0) const;
 
        double beta() const { return beta_; }
        double gamma() const { return gamma_; }
 
    protected:
        double beta_ = 0.25;
        double gamma_ = 0.5;
    };
} // namespace polyfem::time_integrator