ffd_parameterization.cpp

#include "ffd_parameterization.hpp"

namespace PHiLiP {

template<int dim>
FreeFormDeformationParameterization<dim> :: 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)
    : BaseParameterization<dim>(_high_order_grid)
    , ffd(_ffd)
    , ffd_design_variables_indices_dim(_ffd_design_variables_indices_dim)
    {}

template<int dim>
void FreeFormDeformationParameterization<dim> :: initialize_design_variables(
    VectorType &ffd_des_var) 
{
    const unsigned int n_design_variables = ffd_design_variables_indices_dim.size();
    const dealii::IndexSet row_part = dealii::Utilities::MPI::create_evenly_distributed_partitioning(MPI_COMM_WORLD,n_design_variables);
    dealii::IndexSet ghost_row_part(n_design_variables);
    ghost_row_part.add_range(0,n_design_variables);
    ffd_des_var.reinit(row_part, ghost_row_part, MPI_COMM_WORLD);

    ffd.get_design_variables(ffd_design_variables_indices_dim, ffd_des_var);
    
    initial_ffd_des_var = ffd_des_var;
    initial_ffd_des_var.update_ghost_values();
}

template<int dim>
void FreeFormDeformationParameterization<dim> :: compute_dXv_dXp(MatrixType &dXv_dXp) const
{
    ffd.get_dXvdXp(*(this->high_order_grid), ffd_design_variables_indices_dim, dXv_dXp);
}

template<int dim>
bool FreeFormDeformationParameterization<dim> :: update_mesh_from_design_variables(
    const MatrixType &dXv_dXp,
    const VectorType &ffd_des_var)
{
    AssertDimension(ffd_des_var.size(), initial_ffd_des_var.size());
    VectorType current_ffd_des_var = ffd_des_var;
    ffd.get_design_variables( ffd_design_variables_indices_dim, current_ffd_des_var);
    
    bool design_variable_has_changed = this->has_design_variable_been_updated(current_ffd_des_var, ffd_des_var);
    
    bool mesh_updated;
    if(!design_variable_has_changed)
    {
        mesh_updated = false;
        return mesh_updated;
    }
    
    ffd.set_design_variables( ffd_design_variables_indices_dim, ffd_des_var);

    VectorType dXp = ffd_des_var;
    dXp -= initial_ffd_des_var;
    dXp.update_ghost_values();
    VectorType dXv = this->high_order_grid->volume_nodes;
    dXv_dXp.vmult(dXv, dXp);
    dXv.update_ghost_values();
    this->high_order_grid->volume_nodes = this->high_order_grid->initial_volume_nodes;
    this->high_order_grid->volume_nodes += dXv;
    this->high_order_grid->volume_nodes.update_ghost_values();
    mesh_updated = true;
    return mesh_updated;
}

template<int dim>
void FreeFormDeformationParameterization<dim> :: output_design_variables(const unsigned int iteration_no) const
{
    ffd.output_ffd_vtu(iteration_no);
}

template<int dim>
unsigned int FreeFormDeformationParameterization<dim> :: get_number_of_design_variables() const
{
    return ffd_design_variables_indices_dim.size();
}

template class FreeFormDeformationParameterization<PHILIP_DIM>;
} // PHiLiP namespace