Move call to optimizers in mode-computing to separate file to allow other function to access all integrated optimizers

- also moves several options to subfields of options structure
- allows setting options of newrat
- always outputs hessian returned by optimizer (empty matrix if not computed) and deletes subsequent overwriting if cova_compute=0