Use Powell to fit parameters rather than L-BFGS-B
Remove anneal from optimization and fix tests
Return m(Z)^2 when K_l is singular