From ecf9c0416b366ea8ae00cac99ae73c20fc884e5e Mon Sep 17 00:00:00 2001 From: Mayeul Aubin Date: Tue, 29 Apr 2025 18:06:56 +0200 Subject: [PATCH] bugfix: converting int32 to int64 to avoid overflow in products --- analysis/power_spectrum.py | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/analysis/power_spectrum.py b/analysis/power_spectrum.py index 64dacb5..93bb8ee 100644 --- a/analysis/power_spectrum.py +++ b/analysis/power_spectrum.py @@ -15,10 +15,10 @@ def get_power_spectrum(field, kmin=kmin, kmax=kmax, Nk=Nk, G=None): field.L0, field.L1, field.L2, - field.N0, - field.N1, - 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( + int(field.N0), + int(field.N1), + int(field.N2), + 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(kmax), 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.L1, field_A.L2, - field_A.N0, - field_A.N1, - 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( + int(field_A.N0), + int(field_A.N1), + int(field_A.N2), + 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(kmax), Nk,