rol_objective.hpp

#ifndef __ROLOBJECTIVESIMOPT_H__
#define __ROLOBJECTIVESIMOPT_H__

#include "ROL_Objective_SimOpt.hpp"

#include "mesh/free_form_deformation.h"

#include "functional/functional.h"

#include "design_parameterization/base_parameterization.hpp"

namespace PHiLiP {

using ROL_Vector = ROL::Vector<double>;


template <int dim, int nstate>
class ROLObjectiveSimOpt : public ROL::Objective_SimOpt<double> {
private:
    Functional<dim,nstate,double> &functional;

    std::shared_ptr<BaseParameterization<dim>> design_parameterization;

    dealii::LinearAlgebra::distributed::Vector<double> design_var;

public:

    dealii::TrilinosWrappers::SparseMatrix dXvdXp;

    ROLObjectiveSimOpt(
        Functional<dim,nstate,double> &_functional,
        std::shared_ptr<BaseParameterization<dim>> _design_parameterization,
        std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> precomputed_dXvdXp = nullptr);
  
    using ROL::Objective_SimOpt<double>::value;
    using ROL::Objective_SimOpt<double>::update;
  
    void update(
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        bool flag = true,
        int iter = -1) override;
  
    double value(
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  
    void gradient_1(
        ROL::Vector<double> &g,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  
    void gradient_2(
        ROL::Vector<double> &g,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  

    void hessVec_11(
        ROL::Vector<double> &output_vector,
        const ROL::Vector<double> &input_vector,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  

    void hessVec_12(
        ROL::Vector<double> &output_vector,
        const ROL::Vector<double> &input_vector,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  

    void hessVec_21(
        ROL::Vector<double> &output_vector,
        const ROL::Vector<double> &input_vector,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;
  

    void hessVec_22(
        ROL::Vector<double> &output_vector,
        const ROL::Vector<double> &input_vector,
        const ROL::Vector<double> &des_var_sim,
        const ROL::Vector<double> &des_var_ctl,
        double &/*tol*/ ) override;

}; // ROLObjectiveSimOpt

} // PHiLiP namespace

#endif