def call_phys_space_run(cube1_file, cube2_file,
inifile=None):
"""Directly call the power spectral estimation on some physical vol"""
params_init = {
"unitless": True,
"return_3d": False,
"truncate": False,
"window": "blackman",
"bins": [0.00765314, 2.49977141, 35]
}
prefix = 'xs_'
params = parse_ini.parse(inifile, params_init, prefix=prefix)
if inifile is None:
print "WARNING: no ini file for pwrspec estimation"
bparam = params['bins']
bins = np.logspace(math.log10(bparam[0]),
math.log10(bparam[1]),
num=bparam[2], endpoint=True)
retval = pe.calculate_xspec_file(cube1_file, cube2_file, bins,
weight1_file=None, weight2_file=None,
truncate=params['truncate'], window=params['window'],
return_3d=params['return_3d'], unitless=params['unitless'])
return retval
def call_phys_space_run(cube1_file, cube2_file):
"""Directly call the power spectral estimation on some physical vol"""
# define the 1D spectral bins
#nbins=40
#bins = np.logspace(math.log10(0.00702349679605685),
# math.log10(2.81187396154818),
# num=(nbins + 1), endpoint=True)
bins = np.logspace(math.log10(0.00765314),
math.log10(2.49977141),
num=35, endpoint=True)
retval = pe.calculate_xspec_file(cube1_file, cube2_file, bins,
weight1_file=None, weight2_file=None)
