dmipy icon indicating copy to clipboard operation
dmipy copied to clipboard

Refactor spherical mean estimation of data

Open rutgerfick opened this issue 6 years ago • 0 comments

Current implementation in dmipy/dmipy/utils/spherical_mean.py is very slow and ugly with shell-wise for-loops. In spherical mean model fitting it is often half the estimation time.

  • precalculate and pre-invert multishell spherical mean observation matrix A_sm_inv in acquisition scheme
  • multi-shell spherical mean estimation should be nothing more than sm = np.dot(A_sm_inv, data)

rutgerfick avatar Jun 28 '18 11:06 rutgerfick