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@HHEHHp HEHHx HEHHHEHHh 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?mattercdmgasstarBHpressureneutrinoneutrino0neutrino1../../src/power_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 BP    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^  g$2$ 3  3) 13H P3g o3 3 3 345678    9)815; :@;^ w @k,, < $s > < $v >H=w,,,P,a,m,4,,,,,,$4=,I,Z,r,,,4,,,7$7$#8$D,D,D,G<G$ H>J@LDL$MD1M$MDM$cP pP9P 6P9P PPEP PEP PEP QEQ QE.Q$EQ RQ9cQ$sQ$Q Q9Q$Q<Q$Q Q9Q  RE'RD5R2HR HUR9fR$wRDR2R `R9R$R<RDR2R xS9#S$/S<@SDNS2aS nS9S$S<SDS2S S9S SE TDT2,T H9T9JT$[TDiT2|T `T9T$T<TDT2T xT9U$U<$UD2U2EU RU9eU$qU<UDU2U U9U UEU !U )U 1U 9U AU9V XV ` V eV gV A#V9@WLWGW6W7W W pWHW4X5 X6X7X8+X 1X 6X CX9NX8VX5`X :eX;XIXI!ZJ-\66\7;\ B\ L\H.] ?6]KK] :Q] LV] Mh]N^O^O^O_T_T``O{`Oc +cP=cQkc 4xcPc6c7c c 8cHgd\dRdRe AeK]e AeeKe :e Le MeSe :e Le MfSPf :\f Laf MnfSf :f Lf MfSgBRg BZgKgBg BgK-hThhThRhRiC?i CGiKi ri iPi iUj $jPDj QjP_j ljP{j j9j ?j9j Qj9j UjEk Qk92k q?k9k k9kV