Bugs in Reverse

As I noted, there were problems with backward coupling. Oh, it ran, but it wouldn’t couple. Going back over the .backstep code revealed a small bug. The stop condition I had originally coded was this

checkCondition <- (cauchy/(proposal*truncEff)) - (comp/inf$infimumK);
STOP <- if(checkCondition >= 0) {TRUE} else {FALSE};

but the cauchy variable is on a log scale, so it needed to be written this way:

checkCondition <- (exp(cauchy)/(proposal*truncEff)) - (comp/inf$infimumK);
STOP <- if(checkCondition >= 0) {TRUE} else {FALSE};

Maybe this will fix the issue with backward coupling. Of course, it could and I still never see it, so I’m hoping for a quick, successful coupling so I can at least know it works in theory.