Laplacian.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/assembler/Assembler.hpp>
#include <polyfem/utils/AutodiffTypes.hpp>
 
// local assembler for laplace equation
namespace polyfem
{
    namespace assembler
    {
        class Laplacian : public LinearAssembler
        {
        public:
            using LinearAssembler::assemble;
 
            std::string name() const override { return "Laplacian"; }
            std::map<std::string, ParamFunc> parameters() const override { return std::map<std::string, ParamFunc>(); }
 
            Eigen::Matrix<double, Eigen::Dynamic, 1, 0, 9, 1> assemble(const LinearAssemblerData &data) const override;
 
            VectorNd compute_rhs(const AutodiffHessianPt &pt) const override;
 
            void compute_stress_grad_multiply_mat(const OptAssemblerData &data,
                                                  const Eigen::MatrixXd &mat,
                                                  Eigen::MatrixXd &stress,
                                                  Eigen::MatrixXd &result) const override;
 
            void compute_stiffness_value(const double t,
                                         const assembler::ElementAssemblyValues &vals,
                                         const Eigen::MatrixXd &local_pts,
                                         const Eigen::MatrixXd &displacement,
                                         Eigen::MatrixXd &tensor) const override;
 
            Eigen::Matrix<AutodiffScalarGrad, Eigen::Dynamic, 1, 0, 3, 1> kernel(const int dim, const AutodiffGradPt &rvect, const AutodiffScalarGrad &r) const override;
        };
    } // namespace assembler
} // namespace polyfem