Create a spherical harmonic ODF that works the same way as the one in MTEX

[argerlt]

Here is a link to NFFT:
https://www-user.tu-chemnitz.de/~potts/nfft/download.php

NFFT has a function called NFSOFT, which is the thing we need to make ODFs work. it, in turn, uses FFTW (another C code)

In MTEX, Ralph wrapped these files as a series of .mex objects, and that was basically the core of MTEX

we need to do the same thing, but wrapping it with python/cython/whatever

Here is an example of one of the files that does the actual call to NFSOFT from MTEX:
https://github.com/mtex-toolbox/mtex/blob/develop/ODFAnalysis/%40unimodalComponent/calcFourier.m
I'd recommend just downloading a version of MTEX to see how it does it, pausing it at somewhere in this point. the command you want is either "eval()" or "calcDensity" (depends on the MTEX version)

Ask steve for more details