gaussianNormalPair does not work with _use_same_rnd_numbers (Defect #472)
In e7b7aafe I added test for volume sources with use_same_rnd_numbers.
However, I observed nans on CPU and GPU with this setting.
Different to #470, this was not caused by the looping over the sites, all sites are reached. It has to do with the gaussian normal pair.
In fact, the kernel works with a z4 sourcecontent.
Also if one excludes the for-loop by inserting
hmc_complex tmp2 = gaussianNormalPair(&rnd); printf("%f %f\n", tmp2.re, tmp2.im); return;
before the loop one can see that the output of gaussianNormalPair is undefined.
I checked the compile parameters of the two kernels "spinorfield_volume_source" and "spinorfield_gaussian", they differ by "KAPPA=0.125 -D MKAPPA=-0.125 -D NUM_SOURCES=12 -D SOURCE_CONTENT=2". The first 3 are not needed in the kernel, but commenting them out does not change the result.