maubin/sbmy_control
0
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

bugfix: converting int32 to int64 to avoid overflow in products

Browse source
This commit is contained in:
Mayeul Aubin 2025-04-29 18:06:56 +02:00
parent 9352ea6f6d
commit ecf9c0416b
1 changed files with 8 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

16
analysis/power_spectrum.py
View file

@ -15,10 +15,10 @@ def get_power_spectrum(field, kmin=kmin, kmax=kmax, Nk=Nk, G=None):
field.L0, field.L0,
field.L1, field.L1,
field.L2, field.L2,
field.N0, int(field.N0),
field.N1, int(field.N1),
field.N2, int(field.N2),
k_modes=np.concat([PowerSpectrum(field.L0,field.L1,field.L2,field.N0,field.N1,field.N2,).FourierGrid.k_modes[:10],np.logspace( k_modes=np.concat([PowerSpectrum(field.L0,field.L1,field.L2,int(field.N0),int(field.N1),int(field.N2),).FourierGrid.k_modes[:10],np.logspace(
np.log10(kmin), np.log10(kmin),
np.log10(kmax), np.log10(kmax),
Nk, Nk,
@ -43,10 +43,10 @@ def get_cross_correlations(field_A, field_B, kmin=kmin, kmax=kmax, Nk=Nk, G=None
field_A.L0, field_A.L0,
field_A.L1, field_A.L1,
field_A.L2, field_A.L2,
field_A.N0, int(field_A.N0),
field_A.N1, int(field_A.N1),
field_A.N2, int(field_A.N2),
k_modes=np.concat([PowerSpectrum(field_A.L0,field_A.L1,field_A.L2,field_A.N0,field_A.N1,field_A.N2,).FourierGrid.k_modes[:10],np.logspace( k_modes=np.concat([PowerSpectrum(field_A.L0,field_A.L1,field_A.L2,int(field_A.N0),int(field_A.N1),int(field_A.N2),).FourierGrid.k_modes[:10],np.logspace(
np.log10(kmin), np.log10(kmin),
np.log10(kmax), np.log10(kmax),
Nk, Nk,