create_power_operator() returns operator with incorrect data type
If I understood the purpose of
create_power_operator() correctly, it should produce a
DiagonalOperator which imprints the supplied power spectrum onto a given field in harmonic space.
I would expect that this returned operator contains values of real type, but its data type is complex, even if I explicitly ask for real numbers:
from nifty import * import numpy as np s=RGSpace(10,harmonic=True) def ps(x): return x x=create_power_operator(s,ps,dtype=np.float64) print x.diagonal().dtype
What am I missing?