Parametrization.hpp

Go to the documentation of this file.

#pragma once
 
#include <memory>
#include <utility>
#include <vector>
 
#include <Eigen/Core>
 
namespace polyfem::solver
{
    class Parametrization
    {
    public:
        virtual ~Parametrization() = default;
 
        virtual int inverse_size(int y_size) const = 0;
 
        virtual Eigen::VectorXd inverse_eval(const Eigen::VectorXd &y) const = 0;
 
        virtual int size(const int x_size) const = 0;
 
        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() = default;
        CompositeParametrization(std::vector<std::shared_ptr<Parametrization>> parametrizations) : parametrizations_(std::move(parametrizations)) {}
 
        int inverse_size(int y_size) const override;
        Eigen::VectorXd inverse_eval(const Eigen::VectorXd &y) const 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