Dear Einstein Toolkit Users,

I am a new Einstein Toolkit user, and I have some questions about a BNS simulation.

I am trying to run a BNS simulation using bns.par, which I downloaded from Einstein Toolkit Gallery (https://einsteintoolkit.org/gallery/bns/index.html) and slightly modified.

I can successfully create and submit the simulation but the following message appears during the simulation:

INFO (IllinoisGRMHD): C2P: Iter. # 0, Lev: 5 NumPts= 132300 | Backups: 0 0 0 | Fixes: VL= 0 rho*= 0 Averaged pts = 38250 Font1D 38250 | Failures: 2250 InHoriz= 0 / 0 | 1.00 iters/gridpt Error, Sum: rho 6.378e-08, 2.433e-05 | tau nan, -nan | Sx nan, nan | Sy nan, nan | Sz nan, nan

To investigate this problem, I added the following parameters:

Multipole::interpolator_name = "Lagrange polynomial interpolation"
Multipole::interpolator_pars = "order=4"

I also changed the parameter Driver::ghost_size  from 3 to 4, but the same error occurred.

I would like to ask:

Additionally, I have another question. When I use kuibit to extract Psi4, I expected to obtain the full time series (or all available modes) but the returned object contains only one element, for example, 

      SimDir.gws[100][(2,2)].y=[-9.79586237e-15-5.27476073e-14j].

Other attributes such as

       SimDir.gws[100][(2,2)].x 

also contain only one element. Is this a correct behavior or could it indicate a problem with my .par file? 

For reference, I have attached my parameter file.

Thank you very much for your time and any suggestions.

Best regards,

取得 Outlook for Mac