search

TRIQS / tprf

TPRF: The Two-Particle Response Function tool box for TRIQS
https://triqs.github.io/tprf
Other
14 stars 12 forks source link

Bug in screened interaction #3

Closed HugoStrand closed 5 years ago

HugoStrand commented 5 years ago

This is a bug report from M. Roesner on the screened interaction calculation https://github.com/TRIQS/tprf/blob/unstable/c%2B%2B/triqs_tprf/lattice/gw.cpp#L28

W_wk = screened_interaction_W(PI_wk, V_k)

is returning just zeros. I checked the corresponding c++ routine in gw.cpp and to me it feels like that there is something wrong with the matrix view (to_matrix).

At least I get from:

    std::cout<<V_k[k]<<"\n";
    std::cout<<V<<"\n";

something like:

[(32.5812,0) ]
[[(1.38383e-316,1.38307e-316)]]

so V_k[k] is correctly non-zero, but V = to_matrix(V_k[k]) somehow loses this value.

HugoStrand commented 5 years ago

fixed in 648f9e1c5340696b8870a216ae7d915db9dccfcb