config_defaults.hpp

Go to the documentation of this file.

#pragma once
 
#include <cstddef>
#include <string>
 
#include "defs.hpp"
 
namespace ConfigDefaults {
 
// General options
const double ROTATION_FREQUENCY = 0.0; 
const double SELFKERR = 0.0; 
const double CROSSKERR_COUPLING = 0.0; 
const double DIPOLE_COUPLING = 0.0; 
const DecoherenceType DECOHERENCE_TYPE = DecoherenceType::NONE; 
const double DECAY_TIME = 0.0; 
const double DEPHASE_TIME = 0.0; 
 
// Optimization options
const bool CONTROL_ZERO_BOUNDARY_CONDITION = true; 
const ControlType CONTROL_TYPE = ControlType::BSPLINE; 
const size_t CONTROL_SPLINE_COUNT = 10; 
const ControlInitializationType CONTROL_INIT_TYPE = ControlInitializationType::CONSTANT; 
const double CONTROL_INIT_AMPLITUDE = 0.0; 
const double CONTROL_INIT_PHASE = 0.0; 
 
const double CONTROL_AMPLITUDE_BOUND = 1e12; 
const double CARRIER_FREQ = 0.0; 
 
const TargetType OPTIM_TARGET = TargetType::NONE; 
const GateType GATE_TYPE = GateType::NONE; 
const double GATE_ROT_FREQ = 0.0; 
const ObjectiveType OPTIM_OBJECTIVE = ObjectiveType::JTRACE; 
const double OPTIM_WEIGHT = 1.0; 
 
const double OPTIM_TIKHONOV_COEFF = 1e-4; 
const bool OPTIM_TIKHONOV_USE_X0 = false; 
const double OPTIM_TOL_GRAD_ABS = 1e-4; 
const double OPTIM_TOL_GRAD_REL = 1e-4; 
const double OPTIM_TOL_FINAL_COST = 1e-8; 
const double OPTIM_TOL_INFIDELITY = 1e-5; 
const size_t OPTIM_MAXITER = 200; 
 
const double OPTIM_PENALTY_LEAKAGE = 0.0; 
const double OPTIM_PENALTY_WEIGHTEDCOST = 0.0; 
const double OPTIM_PENALTY_WEIGHTEDCOST_WIDTH = 0.5; 
const double OPTIM_PENALTY_DPDM = 0.0; 
const double OPTIM_PENALTY_ENERGY = 0.0; 
const double OPTIM_PENALTY_VARIATION = 0.01; 
 
inline const std::string OUTPUT_DIRECTORY = "./data_out"; 
 
const size_t OUTPUT_TIMESTEP_STRIDE = 1; 
const size_t OUTPUT_OPTIMIZATION_STRIDE = 10; 
 
const RunType RUNTYPE = RunType::SIMULATION; 
const bool USEMATFREE = true; 
const LinearSolverType LINEARSOLVER_TYPE = LinearSolverType::GMRES; 
const size_t LINEARSOLVER_MAXITER = 10; 
const TimeStepperType TIMESTEPPER_TYPE = TimeStepperType::IMR; 
const int RAND_SEED = 1; 
 
} // namespace ConfigDefaults
 
/* Structs to group certain configuration settings */
 
struct InitialConditionSettings {
  InitialConditionType type; 
 
  // Optional settings - populate based on type
  std::optional<std::string> filename; 
  std::optional<std::vector<size_t>> levels; 
  std::optional<std::vector<size_t>> subsystem; 
};
 
struct OptimTargetSettings {
  TargetType type = ConfigDefaults::OPTIM_TARGET; 
  std::optional<GateType> gate_type = std::nullopt; 
  std::optional<std::vector<double>> gate_rot_freq = std::nullopt; 
  std::optional<std::vector<size_t>> levels = std::nullopt; 
  std::optional<std::string> filename = std::nullopt; 
};
 
struct ControlParameterizationSettings {
  ControlType type = ConfigDefaults::CONTROL_TYPE; 
  std::optional<size_t> nspline = ConfigDefaults::CONTROL_SPLINE_COUNT; 
  std::optional<double> tstart = std::nullopt; 
  std::optional<double> tstop = std::nullopt; 
  std::optional<double> scaling = std::nullopt; 
 
};
 
struct ControlInitializationSettings {
  ControlInitializationType type = ConfigDefaults::CONTROL_INIT_TYPE; 
  std::optional<double> amplitude = ConfigDefaults::CONTROL_INIT_AMPLITUDE; 
  std::optional<double> phase = std::nullopt; 
  std::optional<std::string> filename = std::nullopt; 
};