3 Methods to Compute Meson 2-Pt Function

[chrisculver]

Hello! Recently I've been trying to better understand Hadrons and was looking at the Meson module. The code & documentation say the meson 2-pt function can be computed with three options for the sink parameter. I have attempted to do exactly this and am currently confused by the results. A summary of the problem is:

Using an MSink object (I believe a ScalarSink) reproduces the pion correlation function in Test_hadrons_spectrum.cpp. Using the same MSource object as used to compute the propagator, produces a correlator with only one non-zero entry at t=0. Using a SlicedPropagator produces a symmetric correlation function, but not one that decays exponentially.

Firstly here are a few time slices of each correlator to showcase the problem

meson correlator - ScalarSink
t=0 (0.47442148, 0.)
t=1 (0.14026629, 0.)
t=2 (0.01591819, 0.)

meson correlator - sources
t=0 (0.00053426, 0.)
t=1 (0., 0.)
t=2 (0., 0.)

meson correlator - Slice Propagator

t=0 (0.16742971, 0.)
t=1 (5.06398922, 0.)
t=2 (0.4132842, 0.)
t=3 (1.03891438, 0.)
t=4 (0.40092583, 0.)
t=5 (1.03891438, 0.)

I've also attached the main.cpp.txt main file which computes the correlator each way.

I appreciate I may just be misusing the other sink methods, and would appreciate any help.