Added pre()

This avoids a numerical Jacobian. It is also the same implementation as is used for the BCVTB.