Multimodal simulations of fibre optical parametric amplifiers and oscillators