Basically all working (with basic QL). Checked against matlab code without noise, with noise giving different results. Need to implement MF save and AM/dissipation save.