base_parameterization.hpp

#ifndef __DESIGN_PARAMETERIZATION_BASE_H__
#define __DESIGN_PARAMETERIZATION_BASE_H__

#include "mesh/high_order_grid.h"

namespace PHiLiP {
    
template<int dim>
class BaseParameterization {

    using VectorType = dealii::LinearAlgebra::distributed::Vector<double>; 
    using MatrixType = dealii::TrilinosWrappers::SparseMatrix; 

public:
    explicit BaseParameterization(std::shared_ptr<HighOrderGrid<dim,double>> _high_order_grid);
    
    virtual ~BaseParameterization() = default;
    
    virtual void initialize_design_variables(VectorType &design_var) = 0; 
    
    virtual bool update_mesh_from_design_variables( 
        const MatrixType &dXv_dXp,
        const VectorType &design_var) = 0;

    virtual void compute_dXv_dXp(MatrixType &dXv_dXp) const = 0;

    virtual void output_design_variables(const unsigned int /*iteration_no*/) const;
    
    virtual unsigned int get_number_of_design_variables() const = 0;

    bool has_design_variable_been_updated(const VectorType &current_design_var, const VectorType &updated_design_var) const;

    std::shared_ptr<HighOrderGrid<dim,double>> high_order_grid;
    
    MPI_Comm mpi_communicator;

    dealii::ConditionalOStream pcout;

    int mpi_rank;

    int n_mpi;
};

} // PHiLiP namespace

#endif