68963166e5fe4e23655279b653f7cc6ed6d0d0b3,dipy/reconst/mapmri.py,MapmriModel,fit,#MapmriModel#Any#,32

Before Change

                    "CVXOPT package needed to enforce constraints")
            import cvxopt.solvers
            // rmax is linear in mu with rmax \aprox 0.3 for mu = 1/(2*pi*sqrt(700))
            rmax = 0.35 * mu.mean()* (2 * np.pi * np.sqrt(700))
            rgrad = gen_rgrid(rmax = rmax, Nstep = 10)
            K = mapmri_psi_matrix(self.radial_order,  mu, rgrad, self.tau)

After Change

                raise ValueError(
                    "CVXOPT package needed to enforce constraints")
            import cvxopt.solvers
            rmax = 2* np.sqrt(10 * evals.max()*self.tau)
            r_index, r_grad = create_rspace(11, rmax)
            K = mapmri_psi_matrix(self.radial_order,  mu, r_grad[0:len(r_grad)/2,:], self.tau)

            Q = cvxopt.matrix(np.dot(M.T,M)+ self.lambd * I)
            p = cvxopt.matrix(-1*np.dot(M.T,data))

Instances