Helper functions implementing bits of Mk

These are just useful for comparison and not as part of forest
Basic Mk2 model.

Forward and backward calculations done, but no tree integration yet