[triqs/dlr] adapt to new dlr mesh default symmetrize=True by explicit =False

Author: HugoStrand

c++/triqs_xca/atom_diag_utils.cpp

@@ -337,7 +337,7 @@ namespace triqs_xca::atom_diag {
 
     // Create vector of gf<dlr_imtime>
     std::vector<triqs::gfs::gf<triqs::mesh::dlr_imtime>> gf_blocks(H_block_inds.size());
-    triqs::mesh::dlr_imtime tau_mesh(beta, triqs::mesh::Fermion, Lambda / beta, eps);
+    triqs::mesh::dlr_imtime tau_mesh(beta, triqs::mesh::Fermion, Lambda / beta, eps, false);
     for (int i = 0; i < H_block_inds.size(); ++i) { gf_blocks[i] = triqs::gfs::gf<triqs::mesh::dlr_imtime>(tau_mesh, ap_blocks[i]); }
     return {gf_blocks};
   }

c++/triqs_xca/block_sparse.cpp

@@ -693,7 +693,7 @@ namespace triqs_xca::block_sparse {
     auto itops  = cppdlr::imtime_ops(Lambda, dlr_rf);
 
     // triqs gf mesh
-    auto t_mesh = triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, Lambda, eps);
+    auto t_mesh = triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, Lambda, eps, false);
     // create vector of gf
     std::vector<triqs::gfs::gf<triqs::mesh::dlr_imtime>> gf_vec(BDOF.get_num_block_cols());
 

c++/triqs_xca/block_sparse_backbone.cpp

@@ -23,7 +23,7 @@ DiagramEvaluator::DiagramEvaluator(double beta, double Lambda, double eps,
                                    nda::array_const_view<dcomplex, 3> hyb_coeffs,
                                    BlockOpSymQuartet &Fq)
    : 
-     tau_mesh(triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, Lambda / beta, eps)),
+     tau_mesh(triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, Lambda / beta, eps, false)),
      itops(tau_mesh.dlr_it()),
      dlr_it(itops.get_itnodes()),
      Fq(Fq),

c++/triqs_xca/dense_backbone.cpp

@@ -17,7 +17,7 @@ namespace triqs_xca::dense {
 
   DenseDiagramEvaluator::DenseDiagramEvaluator(double beta, double eps, imtime_ops &itops, nda::vector_const_view<double> hyb_poles,
                                                nda::array_const_view<dcomplex, 3> hyb_coeffs, DenseFSet &Fset)
-     : tau_mesh(triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, itops.lambda() / beta, eps)),
+     : tau_mesh(triqs::mesh::dlr_imtime(beta, triqs::mesh::Fermion, itops.lambda() / beta, eps, false)),
        beta(beta),
        itops(itops),
        dlr_it(itops.get_itnodes()),

python/triqs_xca/triqs_solver.py

@@ -67,8 +67,8 @@ def __init__(self, beta, gf_struct, eps, w_max, verbose=True):
         self.eps = eps
         self.w_max = w_max
 
-        self.dmesh = MeshDLR(beta=beta, statistic='Fermion', eps=eps, w_max=w_max)        
-        self.tmesh = MeshDLRImTime(beta=beta, statistic='Fermion', eps=eps, w_max=w_max)        
+        self.dmesh = MeshDLR(beta=beta, statistic='Fermion', eps=eps, w_max=w_max, symmetrize=False)
+        self.tmesh = MeshDLRImTime(beta=beta, statistic='Fermion', eps=eps, w_max=w_max, symmetrize=False)
 
         self.Delta_tau = BlockGf(mesh=self.tmesh, gf_struct=self.gf_struct)
 

test/python/block_sparse_api.py

@@ -85,7 +85,7 @@ def test_oca_diagram_cf_block_sparse_and_dense(beta=2.0, verbose=False):
 
     # -- Hybridization function and adapol fit
 
-    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps)
+    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps, symmetrize=False)
     Delta_w = Gf(mesh=mesh_w, target_shape=[2]*2)
     iwn = np.array([ complex(x) for x in mesh_w ])
     Delta_w.data[:] = make_Delta_with_cont_spec_mat(iwn, semicircular, a=a, b=b, r0=r0, eps=eps)[:, :2, :2]

test/python/block_sparse_api_one_fermion.py

@@ -139,7 +139,7 @@ def test_oca_diagram_cf_block_sparse_and_dense(
 
     print(f'conserved_operators = {conserved_operators}')    
 
-    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps)
+    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps, symmetrize=False)
     Delta_w = Gf(mesh=mesh_w, target_shape=[1]*2)
 
     Delta_w << inverse(iOmega_n - e1)

test/python/block_sparse_api_two_fermions.py

@@ -63,7 +63,7 @@ def test_diagrams_cf_block_sparse_and_dense(e1=-1.5, beta=2.0, conserved_operato
 
     from triqs.gf import MeshDLRImFreq
 
-    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps)
+    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps, symmetrize=False)
     Delta_w = Gf(mesh=mesh_w, target_shape=[2]*2)
     iwn = np.array([ complex(x) for x in mesh_w ])
 

test/python/block_sparse_bethe.py

@@ -30,7 +30,7 @@ def test_block_sparse_self_cons(verbose=False):
 
     H = -mu * N_op
 
-    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps)
+    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps, symmetrize=False)
     Delta_w = Gf(mesh=mesh_w, target_shape=[1, 1])
 
     Delta_w << 0.5 * inverse(iOmega_n - e1)

test/python/block_sparse_normalization.py

@@ -69,7 +69,7 @@ def solve_slater_condon_bethe_half_filling(
     from triqs.gf import SemiCircular
     from triqs.gf import MeshDLRImFreq, Gf, make_gf_dlr_imtime, iOmega_n, inverse
 
-    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps)
+    mesh_w = MeshDLRImFreq(beta=beta, statistic='Fermion', w_max=w_max, eps=eps, symmetrize=False)
     Delta_w = Gf(mesh=mesh_w, target_shape=[n_orb]*2)
 
     Delta_w << inverse(iOmega_n)