ffd_parameterization.hpp

#ifndef __DESIGN_PARAMETERIZATION_FFD_H__
#define __DESIGN_PARAMETERIZATION_FFD_H__

#include "base_parameterization.hpp"
#include "mesh/free_form_deformation.h"

namespace PHiLiP {

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

public:
    FreeFormDeformationParameterization(
        std::shared_ptr<HighOrderGrid<dim,double>> _high_order_grid,
        const FreeFormDeformation<dim> &_ffd,
        std::vector< std::pair< unsigned int, unsigned int > > &_ffd_design_variables_indices_dim);
    
    void initialize_design_variables(VectorType &ffd_des_var) override;
    
    void compute_dXv_dXp(MatrixType &dXv_dXp) const override;
    
    bool update_mesh_from_design_variables(
        const MatrixType &dXv_dXp,
        const VectorType &ffd_des_var) override;

    void output_design_variables(const unsigned int iteration_no) const override;
    
    unsigned int get_number_of_design_variables() const override;

private:
    FreeFormDeformation<dim> ffd;

    const std::vector< std::pair< unsigned int, unsigned int > > ffd_design_variables_indices_dim;

    VectorType initial_ffd_des_var;
};

} // namespace PHiLiP
#endif