LCPQ/dft_tools
3
0
Fork 0
mirror of https://github.com/triqs/dft_tools synced 2025-04-16 13:39:25 +02:00
Code Issues Releases Wiki Activity

test of re_im_freq_time.cpp ok

Browse Source
This commit is contained in:
Laura Messio 2013-07-19 11:33:39 +02:00 committed by Olivier Parcollet
parent 91432bfb80
commit 5f133154e8
1 changed files with 25 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

39
test/triqs/gf/gf_re_im_freq_time.cpp
View File

@ -8,18 +8,11 @@
#include <triqs/arrays.hpp>
namespace tql= triqs::clef;
using triqs::gf::Fermion;
using triqs::gf::make_gf;
using triqs::gf::retime;
using triqs::gf::refreq;
using triqs::gf::refreq_imtime;
using triqs::gf::re_im_time;
using triqs::gf::re_im_freq;
using triqs::arrays::make_shape;
using triqs::gf::scalar_valued;
using namespace triqs::gf;
int main() {
double precision=10e-9;
double beta =1.;
double tmin=0.;
@ -35,13 +28,28 @@ int main() {
auto G_t_tau= make_gf<re_im_time, scalar_valued>( tmin, tmax, n_re_time, beta, Fermion, n_im_time);
auto G_w_wn = make_gf<re_im_freq, scalar_valued>( wmin, wmax, n_re_freq, beta, Fermion, n_im_freq);
auto G_w_tau= make_gf<refreq_imtime, scalar_valued>(wmin, wmax, n_re_freq, beta, Fermion, n_im_time);
auto G_w= make_gf<refreq, scalar_valued>(wmin, wmax, n_re_freq);
triqs::clef::placeholder<0> w_;
triqs::clef::placeholder<1> wn_;
triqs::clef::placeholder<2> tau_;
G_w_wn(w_,wn_)<<1/(wn_-1)/( pow(w_,3) );
G_w_tau(w_,tau_)<< exp( -2*tau_ ) / (w_*w_ +1 );
G_w_tau(w_,tau_)<< exp( -2*tau_ ) / (w_*w_ + 1 );
int index = n_re_freq/3;
double tau = std::get<1>(G_w_tau.mesh().components())[index];
//identical functions
G_w(w_) << exp( -2*tau ) / (w_*w_ + 1 );
//the singularity must be removed as it is inexistent in re_im_time, to give the same TF.
G_w.singularity()(0)=triqs::arrays::matrix<double>{{0}};
G_w.singularity()(1)=triqs::arrays::matrix<double>{{0}};
G_w.singularity()(2)=triqs::arrays::matrix<double>{{0}};
auto G_w2 = slice_mesh_imtime(G_w_tau, index);
for(auto& w:G_w.mesh())
if ( std::abs(G_w(w)-G_w2(w)) > precision) TRIQS_RUNTIME_ERROR<<" fourier_slice error : w="<< w <<" ,G_w="<< G_w(w)<<" ,G_w2="<< G_w2(w) <<"\n";
//test of the interpolation
std::cout << G_t_tau(0.789,0.123) << std::endl;
std::cout << G_w_wn( 0.789,0.123) << std::endl;
std::cout << G_w_tau(0.789,0.123) << std::endl;
@ -54,11 +62,14 @@ int main() {
// try to slice it
auto gt = slice_mesh_imtime(G_t_tau, 1);
std::cout << gt.data()<< std::endl ;
h5_write(file, "gt0", gt);
auto gw = slice_mesh_imtime(G_w_tau, 1);
std::cout << gw.data()<< std::endl ;
h5_write(file, "gw0", gw);
//comparison of the TF of the one time and sliced two times GF's
auto G_t = inverse_fourier(G_w);
auto G_t2 = inverse_fourier(slice_mesh_imtime(G_w_tau, index) );
for(auto& t:G_t.mesh())
if ( std::abs(G_t(t)-G_t2(t)) > precision) TRIQS_RUNTIME_ERROR<<" fourier_slice_re_time error : t="<< t <<" ,G_t="<< G_t(t) <<" ,G_t2="<< G_t2(t) <<"\n";
}