This function evaluates Kennedy's \(\prod\) product, but with the
additional feature that \(\prod_i^j=0\) if \(i>j+1\).
This seems to work in practice.
References
M. C. Kennedy and A. O'Hagan 2000. “Predicting the output from
a complex computer code when fast approximations are available”
Biometrika, 87(1): pp1-13