Parametrization.hpp

Go to the documentation of this file.

#pragma once
 
#include <memory>
#include <vector>
 
#include <Eigen/Core>
 
namespace polyfem::solver
{
    class Parametrization
    {
    public:
        Parametrization() {}
        virtual ~Parametrization() {}
 
        virtual Eigen::VectorXd inverse_eval(const Eigen::VectorXd &y);
 
        virtual int size(const int x_size) const = 0; // just for verification
        virtual Eigen::VectorXd eval(const Eigen::VectorXd &x) const = 0;
        virtual Eigen::VectorXd apply_jacobian(const Eigen::VectorXd &grad_full, const Eigen::VectorXd &x) const = 0;
    };
 
    class CompositeParametrization : public Parametrization
    {
    public:
        CompositeParametrization() : parametrizations_(std::vector<std::shared_ptr<Parametrization>>()) {}
        CompositeParametrization(std::vector<std::shared_ptr<Parametrization>>&& parametrizations) : parametrizations_(parametrizations) {}
        virtual ~CompositeParametrization() {}
 
        Eigen::VectorXd inverse_eval(const Eigen::VectorXd &y) override;
 
        int size(const int x_size) const override;
        Eigen::VectorXd eval(const Eigen::VectorXd &x) const override;
        Eigen::VectorXd apply_jacobian(const Eigen::VectorXd &grad_full, const Eigen::VectorXd &x) const override;
 
    private:
        const std::vector<std::shared_ptr<Parametrization>> parametrizations_;
    };
} // namespace polyfem::solver