Almost first draft for optim with W (still Compute_Omega to write)