ELF>@@UH1EUEUH H ]UHH}HuHUHMLELMf]UHH@Eȉ}ă}EHH:HE1EYfH~EYYEfH~EYEEYfH~EYEEYEYEfH~EYEEYYEfH~EYEEYEYEYEfH~oEYEEYEEYfH~CEHx fH*HHH fH*XHEf(fHnfH~fHnUHH}HEEEimNAEE90E%EEimNAEE90e E%‹E1ЉEEimNAEE90e E%‹E1ЉEHEUE]UHH}HmHEfE HEHH HHEHH HEH HEHEHEHEHHEHEHEfHEHHHEHfHEH HHEHf]UHHH}HEH=EEEf%f ?fEHEHHEHf  %%‹E% ЉEHEHHE ЉEE \UHH`H}HuHUHEHHEHHEHEHEHEHEfE"EfEEE%EHDEЋE1ЉEHEHE}~E/HEHЋEHDEЉ1ȉ‹EHfTEЃE}~fEfEfEEPE1EHLEƋU‹EHDE1ʋEHfTEƃE}~HEH2E}~HEHUHHH}HEHu DHEHu %HEHu HEHu HEHu HEHu HEHu HEHunHHHHujlAAHǸHHHUHEEHH]UHEEHH]UHH HuHUE}u %}u&HE@Y<tHE@Y<}uHE@Y}u&HE@Y<tHE@Y<}uHE@Y<}u;HE@Y<uvHEHHMcHH7 f/L}uAHE@Y<u5HEHHMcHHfH~fHnf/UHHH(H LLH(4EH(HHEHH@HHEHH@HHEHH8 HEHH@ HEHHHEHH0 HEHH HEfEE%EHcHHHHHHEH@X<:EHcHHHHHHEH@YEHcHHHHHHEHHHi`HEHHEEHcHHHHHHEHHHHHH)HHHEHHELEH}HMHUHuHEuMIHgHfH~HEHHEHcHHHHHHEH4Hʉ@EHcHHHHHHEH@Y<u7EHcHHHHHHEH H@HHHEHcHHHHHHEH@LfZ@YE:uEYMXEEHcHHHHHHEH@YEHcHHHHHHEHHHi`HEHHEEHcHHHHHHEHHHHHH)HHHEHHxLEH}HMHUHuHEHxMIHdHfH~HEHHEHcHHHHHHEH4Hʉ8EHcHHHHHHEH@Y<u7EHcHHHHHHEH H8HHHEHcHHHHHHEH@LfZ8YEEYEMXE EE;EH HpEhHpH`HpPHPHXPXhHHHHHXH`H HPHXHPH9uUHH`H}uHUHEHEHEHHEHEH@HEHEH@ HEHEHEExEHHHEHЋHcHHHHHHHHEHHEHE@HE@HEHpLEH}HEMIHEE;E|UHHxH8H0, HHH8HXff/w\f/s YEHHH8HXff/w\f/s YEHHH8HXff/w\f/s YEMX,Ef*ME\EMX,Ef*ME\EMX,Ef*ME\EH0HE,EEEEԉEEĉEEEEEEEEf(\MY\EYxEf(Y\pMXYUXYhf(\MY\EY`Ef(Y\XMXYUXYPf(\MY\EYHEf(Y\@MXYUXY8EYxY`YHEUM40,E(4(Й((0(Й(Ѝ(ȋ,(Й(HHHEHH H HH HHXHHHHH HHHH9uEYxY`Y@EU EE ЙЙЍȋЙHHHEHHHHHHHXHHHHH HHHH9uEYxY`Y8EUMEЙЙЍȋЙHHHEHHHHHHHXHHHHH HHHH9uEYxYXYHEPEEEЙЙЍȋЙHHHEHHHHHxHxHxXpHpHHHH HxHHxH9uEYxYXY@EEEEЙЙЍȋЙHHHEHHxHHpH`H`Hh`XxXHXHHhHpH H`HhH`H9uEYxYXY8EPElEhdE`l`Й``h`Й`Ѝ`ȋd`Й`HHHEHHXPHXHHHXHHHHPHXP@H@HHPHHH HHHPHHH9uEYxYPYHEUMD@{>@B>{>=,?{>@B>{>@B>{>=,?{>=,?    X?zŎ?`?E??䙂?O?f?Ѐ?1?u?O?p?}?E$?w ?A=?_?>?j?-2?5??` ?o?#?"? ?!??? ????:???\ ??? ?? ? ?? ????( ? ??K ? ? ?????s?5??? ??@?H??%???? ??5???????h??????S?_???????`??-??z?K???G?t???'?8????????v?T??F????s?:?x?6?b?mattercdmgasstarBHpressureneutrinoneutrino0neutrino1power_spectrum.c[%04i] %s %s:%s():%i: Do not recognize the power spectrum type '%s'. # %s-%s cross-spectrum # %s power spectrum ###################################################### # # (0) Redshift # Unit = dimensionless # (1) Fourier scale k # Unit = %e cm**-1 # Unit = %e pc**-1 # Unit = %e Mpc**-1 # (2) The shot-noise subtracted power # Unit = %e cm**3 # Unit = %e pc**3 # Unit = %e Mpc**3 # Unit = %e Mpc**3 * eV**2 * cm**-6 # Unit = %e Mpc**3 * eV * cm**-3 # (3) The shot-noise power (scale-independant) (3) (2) (1) (0) # %13s %15s %15s %15s P_noiseP(k)kz[%04i] %s %s: Preparing to calculate the %s-%s power spectrum. [%04i] %s %s:%s():%i: Cell infrastructure is not in place for power spectra. [%04i] %s %s: Calculating the shot noise. power_%s-%s[%04i] %s %s: Calculating the power for folding num. %d. power_spectra/foldings%s/%s_%04d_%d.txtwMpc^3Mpc^3 (eV cm^(-3))^2Mpc^3 eV cm^(-3)# Folding %d, all lengths/volumes are comoving. k-bin centres are not corrected for the weights of the modes. # Shotnoise [%s] %g # Redshift [dimensionless] # k [Mpc^(-1)] p [%s] %g %g power_spectra%s/%s_%04d.txt%15.8f %15.8e %15.8e %15.8e [%04i] %s %s: Calculated the power! Cleaning up and leaving. PowerSpectrum:grid_side_lengthPowerSpectrum:num_foldsPowerSpectrum:fold_factorPowerSpectrum:window_order[%04i] %s %s:%s():%i: Power spectrum calculation is not implemented for %dth order! [%04i] %s %s: WARNING: window order is recommended to be at least 2 (CIC). [%04i] %s %s: WARNING: fold factor is recommended not to exceed 4 for a mass assignment order of 2 (CIC) or below. [%04i] %s %s: WARNING: fold factor is recommended not to exceed 6 for a mass assignment order of 3 (TSC) or below. PowerSpectrum:shift_centre_small_k_bins[%04i] %s %s:%s():%i: Combination of power grid size and fold factor do not allow for enough overlap between foldings! PowerSpectrum:requested_spectra-[%04i] %s %s:%s():%i: Requested power spectra are not in the format type1-type2! [%04i] %s %s: took %.3f %s. power spectrum datapower types 1power types 2power_spectrumMatterCDMstars/BHselectron pressureneutrino (even)neutrino (odd)power_initpower_spectrum_get_typecalc_all_power_spectra???.A@-DT! @p@-DT!@GCC: (GNU) 14.1.0zRx AC V <"AC ] \tAC o |AC  AC  AC  <AC 7 qAC l AC T <AC T \LAC G |AC  AC  $'AC ' n AC i AC  AC  <9AC 4 \AC  |iAC d  AC E   QAC JB  AC H  AC G ,TAC O LAC  lAC  0)"G@ [ 8h lt l 4@ =tI< q#K4@DdY@@m\$'5n BPx    4\}#Ow$+3GQ}Let EI9JKi.VQ,*:Veup  2JhovwxTx4yRpower_spectrum.cgettickskernel_coeffscooling_get_electron_pressurelightcone_map_typesrunner_shiftrunner_flipsortlistIDsid_scalesub_sid_flaginteger_powinl_rand_rinl_drand48_iterateinl_erand48random_unit_intervalcorrection_shift_k_valuespower_spectrum_get_type__FUNCTION__.1get_powtype_namepowtype_names.4get_powtype_filenamepowtype_filenames.3gpart_to_grid_TSCgpart_to_grid_CICgpart_to_grid_NGPpower_init_output_file__FUNCTION__.5__FUNCTION__.2__FUNCTION__.0lightcone_map_total_mass_get_valuelightcone_map_total_mass_type_contributeslightcone_map_total_mass_baseline_valuelightcone_map_gas_mass_get_valuelightcone_map_gas_mass_type_contributeslightcone_map_dark_matter_mass_get_valuelightcone_map_dark_matter_mass_type_contributeslightcone_map_stellar_mass_get_valuelightcone_map_stellar_mass_type_contributeslightcone_map_black_hole_mass_get_valuelightcone_map_black_hole_mass_type_contributeslightcone_map_sfr_get_valuelightcone_map_sfr_type_contributesstrcasecmpstdoutfflushclocks_get_timesincestartengine_rankstderrfprintfompi_mpi_comm_worldMPI_Abortshould_collect_massshotnoisetermsgpart_neutrino_weight_mesh_onlyshotnoise_mappercell_to_powgridcell_to_powgrid_mappermass_to_contrast_mapperpow_from_grid_mappersinfwritepower_spectrumgather_neutrino_constsfftw_alloc_realunits_cgs_conversion_factorthreadpool_mapompi_mpi_op_sumompi_mpi_doubleMPI_AllreducemallocsprintfstrlenmemsetMPI_Reducefftw_execute_dft_r2cfopenfclosefftw_freepower_initparser_get_opt_param_intparser_get_param_intfftw_init_threadsfftw_plan_with_nthreadsparser_get_param_string_arraystrtokstrcpyfftw_plan_dft_r2c_3dsafe_checkdircalc_all_power_spectraclocks_getunitclocks_from_tickspower_cleanfftw_destroy_planfftw_cleanup_threadspower_spectrum_struct_dumprestart_write_blockspower_spectrum_struct_restorerestart_read_blocks^  g2 3  3) 13H P3g o3 3 3 345678    9)815; :@;^ w @k < s > < v >H=wPam$$$=IZr$77#8DDDG<G H>J@LDLMD1MMDMcP pP9P .P9P HPEP PEP PEP QEQ QE.QEQ RQ9cQsQQ Q9QQ,QQ Q9Q RE'R45R2HR @UR9fRwR4R2R XR9RR,R4R2R pS9#S/S,@S4NS2aS nS9SS,S4S2S S9S SE T4T2,T @9T9JT[T4iT2|T XT9TT,T4T2T pT9UU,$U42U2EU RU9eUqU,U4U2U U9U UEU U !U )U 1U 9U9V PV X V ]V _V 9#V9@W<WGW6W7W xW hWHW4X5 X6X7X8+X x1X 6X CX9NX8VX5`X :eX;XIXI!ZJ-\66\7;\ xB\ L\H.] ?6]KK] :Q] LV] Mh]N^O^O^O_D_D``O{`Oc #cP=cQkc ,xcPc6c7c xc 0cHgdLdRdRe AeK]e AeeKe :e Le MeSe :e Le MfSPf :\f Laf MnfSf :f Lf MfSgBRg BZgKgBg BgK-hThhThRhRiC?i CGiKi ji iPi iUj $jPDj QjP_j ljP{j j9j 7j9j Ij9j MjEk Ik92k i?k9k k9kV