Tags tg_n1 and tg_n3 added to config header for passing messages.

MPI_Send calls added to Config::sendResults for n1 and n3

Arrays added to main to hold n1 and n3 for each replica. Recieve calls added.

prob modified to use new probability.
Changed code so that swaps happen during equlibration. Now equilibration is just the ‘zeroth’ bin, in which everything is done as usual, except the cutoff is updated after each MC sweep and data is collected, but not written to file.
This version produces segfaults. I’m not sure why.
Still segfaulting, fixed error where swaps didn’t occur for large systems
Segfault was coming from the fact that nreps was unset in slave nodes (only set in master). Different implementations disagreed about what to do in this situation, so the result was undefined behavior. Now everything should work as expected.
Cleaned up code, got rid of unnecessary repMags array, commented out mc step progress reporting.
Added a section that does isteps for equilibration and mcsteps for all the other bins.
Turned off reporting to stdout of the number of successful swaps.
Writes out the number of type 1 and type 3 operators in each replica as a testing procedure.
Added debugging print statements to check the math on the probability calculation and print the number of operators of each type for each replica.
Print the number of successful swaps for each hnum at the end of each bin.
It turns out my acceptance probability is fine, it just simulations can’t easily pass from the high magnetization regime to the low-magnetization regime.

I also have to make sure that I can’t divide by zero by accident.
Fixed probability to allow for epsilon != 0 and to avoid dividing by zero when h could be zero.
Fixed error discovered in swap acceptance for h_A = 0.