inner_vol_parameterization.hpp

#ifndef __DESIGN_PARAMETERIZATION_INNERVOL_H__
#define __DESIGN_PARAMETERIZATION_INNERVOL_H__

#include "base_parameterization.hpp"

namespace PHiLiP {

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

public:
    explicit InnerVolParameterization(
        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;
    

    void compute_inner_vol_index_to_vol_index();

private:
    VectorType current_design_var;
    unsigned int n_inner_nodes;
    dealii::IndexSet inner_vol_range;
    dealii::LinearAlgebra::distributed::Vector<int> inner_vol_index_to_vol_index;
};

} // namespace PHiLiP
#endif