mevals = [None, ] * len(mf)
if mevecs is None:
mevecs = [np.eye(3) for i in range(len(mf))]
for j, f in enumerate(mf):
pdf += f * single_tensor_pdf(pdf_points,
After Change
else:
sticks = [sphere2cart(1, np.deg2rad(pair[0]), np.deg2rad(pair[1]))
for pair in angles]
sticks = np.array(sticks)
pdf = np.zeros(len(pdf_points))
mevecs = []
for s in sticks:
mevecs += [all_tensor_evecs(s).T]
for j, f in enumerate(mf):
pdf += f * single_tensor_pdf(pdf_points,
evals=mevals[j], evecs=mevecs[j], tau=tau)
return pdf
