G4-NGENIC: gsl: qag.c:257: ERROR: maximum number of subdivisions reached

From: Tiago Castro <tiagobscastro_at_gmail.com>
Date: Fri, 30 Apr 2021 19:31:46 +0200

Dear list,

  I have tried to run an Einstein de Sitter/power-law power-spectrum nbody
simulation with G4. However, during IC creation, gadget returns the error
bellow if the spectrum is steeper than k**-2. For the specific log pasted,
the power spectrum was k**-2.5. Have someone faced a similar issue?


Kind regards,
Tiago.

----
MALLOC: Allocation of shared memory took 0.00993109 sec
found 63 times in output-list.
BEGRUN: Hubble (internal units)  = 0.1
BEGRUN: h                        = 0.6766
BEGRUN: G (internal units)       = 43018.7
BEGRUN: UnitMass_in_g            = 1.989e+43
BEGRUN: UnitLenth_in_cm          = 3.08568e+21
BEGRUN: UnitTime_in_s            = 3.08568e+16
BEGRUN: UnitVelocity_in_cm_per_s = 100000
BEGRUN: UnitDensity_in_cgs       = 6.76991e-22
BEGRUN: UnitEnergy_in_cgs        = 1.989e+53
NGENIC: generated grid of size 1024
NGENIC: computing displacement fields...
NGENIC: vel_prefac1= 5.57763  hubble_a=55.7764   fom1=0.999999
NGENIC: vel_prefac2= 11.1553  hubble_a=55.7764   fom2=2
found 800 rows in input spectrum table
Normalization of spectrum in file:  Sigma8 = 0.8277
Normalization adjusted to  Sigma8=0.83   (Normfac=1.00556)
NGENIC: Dplus=78.5555
NGENIC_2LPT: Computing secondary source term, derivatices 0 0
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Computing secondary source term, derivatices 1 1
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Computing secondary source term, derivatices 2 2
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Computing secondary source term, derivatices 0 1
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Computing secondary source term, derivatices 0 2
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Computing secondary source term, derivatices 1 2
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Secondary source term computed in real space
NGENIC_2LPT: Done transforming it to k-space
NGENIC_2LPT: Obtaining second order displacements for axes=0
NGENIC_2LPT: Obtaining second order displacements for axes=1
NGENIC_2LPT: Obtaining second order displacements for axes=2
NGENIC_2LPT: Obtaining Zeldovich displacements for axes=0
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Obtaining Zeldovich displacements for axes=1
NGENIC: setting up modes in kspace...
NGENIC_2LPT: Obtaining Zeldovich displacements for axes=2
NGENIC: setting up modes in kspace...
NGENIC: Maximum displacement: 1061.45, in units of the part-spacing=
1.08692
NGENIC: Maximum velocity component: 5919.66
gsl: qag.c:257: ERROR: maximum number of subdivisions reached
Default GSL error handler invoked.

Received on 2021-04-30 19:32:29