compute two-point correlation function between quark loops
loop_2pt(loop_snk, loop_src, random_vectors_outer_product = FALSE,
nstoch_to_avg = "all")'raw_cf' container with spin-projected quark loop at sink (see loop_spin_project)
'raw_cf' container with spin-projected quark loop at source (see loop_spin_project)
Boolean. For testing purposes, the average over all random sample combinations can be carried out explicitly as \( \sum_{i \neq j} Tr[ \Gamma_snk M_i ] Tr[ \Gamma_src M_j ] \) instead of the (much faster) equivalent \( (\sum_i Tr[ \Gamma_snk M_i ])*(\sum_j Tr[ \Gamma_src M_j ]) - \sum_i ( Tr[ \Gamma_snk M_i ] Tr[ \Gamma_src M_i ] ) \).
String or integer, how many of the available stochastic samples should be averaged over. See loop_stochav for details.
'raw_cf' container with two-point function of these two quark loops. In the calculation, both averaging over all source locations and the average over all stochastic sample combinations are performed.