KernelProblem.hpp

Go to the documentation of this file.

#pragma once
 
#include "ProblemWithSolution.hpp"
 
#include <Eigen/Dense>
 
#include <vector>
#include <string>
 
namespace polyfem
{
    namespace problem
    {
 
        class KernelProblem : public ProblemWithSolution
        {
        public:
            KernelProblem(const std::string &name, const assembler::Assembler &assembler);
 
            VectorNd eval_fun(const VectorNd &pt, const double t) const override;
            AutodiffGradPt eval_fun(const AutodiffGradPt &pt, const double t) const override;
            AutodiffHessianPt eval_fun(const AutodiffHessianPt &pt, const double t) const override
            {
                assert(false);
                return AutodiffHessianPt(1);
            }
 
            void rhs(const assembler::Assembler &assembler, const Eigen::MatrixXd &pts, const double t, Eigen::MatrixXd &val) const override;
 
            void set_parameters(const json &params) override;
            bool is_scalar() const override;
 
        private:
            const assembler::Assembler &assembler_;
            int n_kernels_ = 5;
            double kernel_distance_ = 0.05;
            Eigen::VectorXd kernel_weights_;
        };
    } // namespace problem
} // namespace polyfem