Computes the quark-line disconnected contribution to a three-point function of the form $$ C_3(t, \Delta t = t_{snk} - t_{src}) = C_2(t_{snk}, t_{src}) * L(t) $$ \(\forall t\) considering only the case t_snk > t_src.
disc_3pt(cf_2pt, loop, src_ts, dt, reim_loop = "both", reim_2pt = "both",
vev_subtract = FALSE)'raw_cf' container holding two-point part of three-point function in lattice-absolute coordinates (not relative to source!)
'raw_cf' container holding loop contribution, suitably spin-projected and averaged over stochastic samples.
Integer vector, the source time slices that were used for the computation
of the two-point function in lattice-absolute coordinates.
Must be of the same length as the number of measurements in cf_2pt.
Integer, the source-sink separation that should be computed.
String, one of 'real', 'imag' or 'both'. Specifies whether just the real or imaginary part should be considered when constructing the correlation with the two-point function.
String, same as reim_loop but for the two-point
contribution to the three-point function.
Boolean, whether the loop contains a vev which should be subtracted.
raw_cf container with the product of loop and 2pt function, shifted
in time to be relative to source using the info from src_ts