added "AddJacobianDiagonals", successfully computed f using J*v = f. Next need to remove all instances of MI and MR.