I still have negative definite hessians every time I have \sigma above 1.0  This is not good.  I set up natural coordinate descent, but it is still buggy.