def _get_ps_func_for_field(input_field, box_dims):
'''
Return ps(k) for the specified field. For internal use.
'''
ps_input, k_input, n_modes = power_spectrum_1d(input_field, \
box_dims=box_dims, kbins=10, return_n_modes=True)
ps_k = interp1d(k_input[n_modes>0], ps_input[n_modes>0], kind='linear', \
bounds_error=False, fill_value=0.)
return ps_k
def _get_ps_func_for_field(input_field, box_dims):
'''
Return ps(k) for the specified field. For internal use.
'''
ps_input, k_input, n_modes = power_spectrum_1d(input_field, \
box_dims=box_dims, kbins=10, return_n_modes=True)
ps_k = interp1d(k_input[n_modes>0], ps_input[n_modes>0], kind='linear', \
bounds_error=False, fill_value=0.)
return ps_k
def do_powerspectrum(filename,doubleflag):
xfrac = read_xfrac(filename,doubleflag).data
ps,bins = power_spectrum.power_spectrum_1d(xfrac, [47.,47.,47.], kbins=153)
return ps
def power_spectrum(id=''):
for i in range(len(redshifts)):
print "Doing redshift: " + str(redshifts[i])
data=np.load(setup_dirs.resultsdir()+"map_dbt_"+id+str('%.3f' % redshifts[i])+".bin")
powerspec=ps.power_spectrum_1d(data,100)
IO.write2data(powerspec[0],powerspec[1],setup_dirs.resultsdir()+'/powerSpectra_100b_'+id+str('%.3f' % redshifts[i])+'.dat',setup_dirs.resultsdir()+'/powerSpectraFrequencies_dbt_100b_'+str('%.3f' % redshifts[i])+'.dat')
def temp_power_spectrum():
for i in range(len(redshifts)):
print "Doing redshift: " + str(redshifts[i])
data=c2t.TemperFile(setup_dirs.path()+"Temper3D_"+str('%.3f' % redshifts[i])+".bin").temper
powerspec=ps.power_spectrum_1d(data,100)
IO.write2data(powerspec[0],powerspec[1],setup_dirs.resultsdir()+'/powerSpectra_temp_'+str('%.3f' % redshifts[i])+'.dat',setup_dirs.resultsdir()+'/powerSpectraFrequencies_dbt_100b_'+str('%.3f' % redshifts[i])+'.dat')
def xfrac_power_spectrum(id):
for i in range(len(redshifts)):
print "Doing redshift: " + str(redshifts[i])
data=c2t.XfracFile(setup_dirs.path()+"xfrac3d"+id+'_'+str('%.3f' % redshifts[i])+".bin").xi
powerspec=ps.power_spectrum_1d(data,100)
IO.write2data(powerspec[0],powerspec[1],setup_dirs.resultsdir()+'/powerSpectra_xfrac_'+id+'_'+str('%.3f' % redshifts[i])+'.dat',setup_dirs.resultsdir()+'/powerSpectraFrequencies_dbt_100b_'+str('%.3f' % redshifts[i])+'.dat')
