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')