optimtarget.hpp

Go to the documentation of this file.

#include "defs.hpp"
#include "gate.hpp"
#include "mastereq.hpp"
#pragma once
 
class OptimTarget{
    protected:
 
    PetscInt dim; 
    PetscInt dim_rho; 
    PetscInt dim_ess; 
    int noscillators; 
 
    TargetType target_type; 
    ObjectiveType objective_type; 
    Gate *targetgate; 
    double purity_rho0; 
    PetscInt purestateID; 
    Vec targetstate; 
    InitialConditionSettings initcond; 
    DecoherenceType decoherence_type; 
    int mpisize_petsc; 
    int mpirank_petsc; 
    PetscInt localsize_u; 
    PetscInt ilow; 
    PetscInt iupp; 
 
    Vec aux; 
    bool quietmode; 
 
  public:
    OptimTarget();
 
    OptimTarget(const Config& config, MasterEq* mastereq, double total_time, Vec rho_t0, bool quietmode_);
 
    ~OptimTarget();
 
    ObjectiveType getObjectiveType(){ return objective_type; };
 
    int prepareInitialState(const int iinit, const int ninit, const std::vector<size_t>& nlevels, const std::vector<size_t>& nessential,  Vec rho0);
 
    void prepareTargetState(const Vec rho);
 
    void evalJ(const Vec state, double* J_re_ptr, double* J_im_ptr);
 
    void evalJ_diff(const Vec state, Vec statebar, const double J_re_bar, const double J_im_bar);
 
    double finalizeJ(const double obj_cost_re, const double obj_cost_im); 
 
    void finalizeJ_diff(const double obj_cost_re, const double obj_cost_im, double* obj_cost_re_bar, double* obj_cost_im_bar); 
 
    double FrobeniusDistance(const Vec state);
 
    void FrobeniusDistance_diff(const Vec state, Vec statebar, const double Jbar);
 
    void HilbertSchmidtOverlap(const Vec state, const bool scalebypurity, double* HS_re_ptr, double* Hs_im_ptr );
 
    void HilbertSchmidtOverlap_diff(Vec statebar, bool scalebypurity, const double HS_re_bar, const double HS_im_bar);
};