ENH: Use analytic scores and Hessians
ENH: Change default solver back to bfgs
TST: Use BFGS in test