identity_parameterization.hpp

#ifndef __DESIGN_PARAMETERIZATION_IDENTITY_H__
#define __DESIGN_PARAMETERIZATION_IDENTITY_H__

#include "base_parameterization.hpp"

namespace PHiLiP {

template<int dim>
class IdentityParameterization : public BaseParameterization<dim> {
    
    using VectorType = dealii::LinearAlgebra::distributed::Vector<double>; 
    using MatrixType = dealii::TrilinosWrappers::SparseMatrix; 

public:
    explicit IdentityParameterization(
        std::shared_ptr<HighOrderGrid<dim,double>> _high_order_grid); 

    void initialize_design_variables(VectorType &design_var) override;
    

    void compute_dXv_dXp(MatrixType &dXv_dXp) const override;
    
    bool update_mesh_from_design_variables(
        const MatrixType &dXv_dXp,
        const VectorType &design_var) override;

    unsigned int get_number_of_design_variables() const override;

private:
    VectorType current_volume_nodes;
};

} // namespace PHiLiP
#endif