params.hpp

// Copyright (c) 2022--present, The Simons Foundation
// Copyright (c) 2022--present, Max Planck Institute for Polymer Research, Mainz, Germany
// This file is part of TRIQS/ctseg and is licensed under the terms of GPLv3 or later.
// SPDX-License-Identifier: GPL-3.0-or-later
// See LICENSE in the root of this distribution for details.
 
#pragma once
 
#include <triqs/gfs.hpp>
#include <triqs/operators/many_body_operator.hpp>
using namespace triqs::gfs;
 
namespace triqs_ctseg {

  struct constr_params_t {

    double beta;

    gf_struct_t gf_struct;

    int n_tau = 10001;

    int n_tau_bosonic = 10001;
  };
 
  //---------------------------------------------

  struct solve_params_t {

    triqs::operators::many_body_operator h_int;

    triqs::operators::many_body_operator h_loc0;

    int n_tau_G = 0;

    int n_tau_chi2 = 0;

    double dlr_omega_max = 100;

    double dlr_epsilon = 1e-8;

    int n_w_b_vertex = 10;

    int n_w_f_vertex = 10;

    int n_cycles;

    int length_cycle = 50;

    int length_cycle_g3 = 0;

    int n_warmup_cycles = 5000;

    int random_seed = 34788 + 928374 * mpi::communicator().rank();

    std::string random_name = "";

    int max_time = -1;

    int verbosity = mpi::communicator().rank() == 0 ? 3 : 0;
 
    // -------- Move control --------------

    bool move_insert_segment = true;

    bool move_remove_segment = true;

    bool move_double_insert_segment = true;

    bool move_double_remove_segment = true;

    bool move_move_segment = true;

    bool move_split_segment = true;

    bool move_regroup_segment = true;

    bool move_insert_spin_segment = true;

    bool move_remove_spin_segment = true;

    bool move_split_spin_segment = true;

    bool move_regroup_spin_segment = true;

    bool move_swap_spin_lines = true;
 
    // -------- Measure control --------------

    bool measure_pert_order = true;

    bool measure_G_tau = true;

    bool measure_F_tau = false;

    bool measure_densities = true;

    bool measure_average_sign = true;

    bool measure_nn_static = false;

    bool measure_nn_tau = false;

    bool measure_nn_nu_dlr = false;

    bool measure_Sperp_tau = false;

    bool measure_state_hist = false;

    bool measure_g2w = false;

    bool measure_g3w = false;
 
    // -------- Misc parameters --------------

    double imag_threshold = 1.e-13;

    int det_init_size = 100;

    int det_n_operations_before_check = 100;

    double det_precision_warning = 1.e-8;

    double det_precision_error = 1.e-5;

    double det_singular_threshold = -1;

    int histogram_max_order = 1000;

    bool visualize_config = false;
  };

  struct params_t : constr_params_t, solve_params_t {
    params_t(constr_params_t const &constr_params_, solve_params_t const &solve_params_)
       : constr_params_t{constr_params_}, solve_params_t{solve_params_} {}
  };

  void h5_write(h5::group h5group, std::string subgroup_name, constr_params_t const &c);

  void h5_read(h5::group h5group, std::string subgroup_name, constr_params_t &c);

  void h5_write(h5::group h5group, std::string subgroup_name, solve_params_t const &c);

  void h5_read(h5::group h5group, std::string subgroup_name, solve_params_t &c);
 
} // namespace triqs_ctseg