GCPContactForceDerivative.hpp

Go to the documentation of this file.

#pragma once
 
#include <Eigen/Core>
#include <polyfem/solver/forms/GCPContactForm.hpp>
#include <ipc/gcp/gcp_collisions.hpp>
 
namespace polyfem::solver
{
    class GCPContactForceDerivative
    {
    public:
        static void force_shape_derivative(
            const GCPContactForm &form,
            const ipc::GCPCollisions &collision_set,
            const Eigen::MatrixXd &solution,
            const Eigen::VectorXd &adjoint_sol,
            Eigen::VectorXd &term);
    };
} // namespace polyfem::solver