Still trying to fix the bugs.  I ran the solver using a simulated B Table and low sigmas, but still no luck.  I might need a new way to choose step sizes, stochastic gradient descent, step-wise descent, or some other change.  I'm also considering refactoring entirely into something that could run for both the isotropic model and (multiple parameterizations) of the DTI model.