self.model.cache_set("interp_coords", sphere, interp_coords)
Pr = self.pdf()
//calculate the orientation distribution function
return pdf_odf(Pr, sphere, self.model.qradius, interp_coords)
def pdf_interp_coords(sphere, rradius, origin):
Precompute coordinates for ODF calculation from the PDF
After Change
Pr = self.pdf()
//calculate the orientation distribution function
odf = pdf_odf(Pr, sphere, self.model.qradius, interp_coords)
// We compute the gfa here, since we have the ODF available and
// the storage requirements are minimal. Otherwise, we"d need to
// recompute the ODF at significant computational expense.
self._gfa = gfa(odf)
return odf
@property
def gfa(self):
if self._gfa is None: