<div dir="ltr"><div>Dear Taiki,<br><br>This is an experimental response prepared using a set of LLM agents and as a<br>test of the newly created GRHayLET knowledge base:<br><a href="https://github.com/GRHayL/GRHayLET/blob/main/AGENTS.md">https://github.com/GRHayL/GRHayLET/blob/main/AGENTS.md</a><br><br>Zach will be traveling for the next month and hopes that this response will<br>be helpful in the meantime.<br><br>The single Psi4 value is consistent with the simulation stopping at iteration<br>0. In kuibit, `SimDir.gws[100]` selects the extraction radius r=100, and<br>`[(2,2)]` selects one spherical-harmonic mode. The resulting `.x` and `.y`<br>arrays are the time and complex data samples for that particular radius and<br>mode; they are not containers for all radii or modes. You can inspect what was<br>read with<br><br> print(SimDir.gws.radii)<br> print(SimDir.gws[100].available_lm)<br> print(SimDir.gws[100][(2,2)].x)<br><br>Since the last command returns only t=0 and the Psi4 ASCII files contain only<br>one row, kuibit appears to be reading the available output correctly. With<br>`Multipole::out_every = 128`, the next sample would be written at iteration<br>128; the attached parameter file instead uses 64, so its next sample would be<br>at iteration 64. Neither setting can produce another sample if the evolution<br>never advances past iteration 0. Files containing spatial profiles can have<br>many rows while still representing only iteration 0.<br><br>The IllinoisGRMHD message does report C2P recovery failures, but it is an INFO<br>message and does not by itself explain why Cactus stopped:<br><br>* `Averaged pts = 38250 Font1D 38250` says that 38,250 points first required<br> neighbor-averaging attempts and then reached the Font1D fallback.<br>* `Failures: 2250` says that recovery still failed at 2,250 points and those<br> points were reset to atmosphere.<br>* These are counters for the local grid data handled by that printed C2P<br> invocation, not necessarily global counts over the whole simulation.<br>* `tau nan` and the NaN momentum diagnostics mean that an aggregate<br> post-recovery error calculation encountered non-finite values. They do not<br> imply that every value of tau or every momentum component is NaN.<br><br>Thus, there are genuine C2P failures, but C2P may be reacting to invalid<br>metric, primitive, EOS, or conservative input rather than creating the first<br>NaN. The first priority is to distinguish the reason for termination from the<br>first source of non-finite data.<br><br>First, please inspect the batch-system exit status and the final 100--200 lines<br>of both stdout and stderr. Search backward for the termination request,<br>especially `WARNING level 0`, `CCTK_ERROR`, NaNChecker, TerminationTrigger,<br>MPI abort, or scheduler kill messages; an unrelated earlier warning need not<br>explain the stop. The attached parameter file also enables<br>`TerminationTrigger::termination_from_file` with the file name `TERMINATE`.<br>Please check for a stale `TERMINATE` file in the run/output directory. Such a<br>file can request a clean stop at iteration 0 even if the C2P INFO message did<br>not terminate the run.<br><br>There also appears to be a mismatch between the quoted diagnostic and the<br>requested recovery configuration. If the quoted text is the complete<br>diagnostic block, a line containing `Font1D` but neither entropy nor Y_e<br>resembles the Hybrid non-entropy path in current IllinoisGRMHD. In contrast,<br>the attached file requests a tabulated EOS with entropy evolution and three<br>backup routines. The TabulatedEntropy diagnostic includes entropy and Y_e and<br>does not print that `Font1D` field:<br><br><a href="https://github.com/GRHayL/GRHayLET/blob/301648c62bb87544adc7399d3239fef97ae0e5d4/IllinoisGRMHD/src/Hybrid/conservs_to_prims.c#L334-L354">https://github.com/GRHayL/GRHayLET/blob/301648c62bb87544adc7399d3239fef97ae0e5d4/IllinoisGRMHD/src/Hybrid/conservs_to_prims.c#L334-L354</a><br><br><a href="https://github.com/GRHayL/GRHayLET/blob/301648c62bb87544adc7399d3239fef97ae0e5d4/IllinoisGRMHD/src/TabulatedEntropy/conservs_to_prims.c#L337-L361">https://github.com/GRHayL/GRHayLET/blob/301648c62bb87544adc7399d3239fef97ae0e5d4/IllinoisGRMHD/src/TabulatedEntropy/conservs_to_prims.c#L337-L361</a><br><br>Please check the complete diagnostic block, the startup schedule, and the<br>parameter file copied into the simulation output directory. In particular:<br><br>1. Run once in a fresh output directory with `IOUtil::recover = "no"`, or<br> otherwise disable recovery explicitly, and ensure that no stale termination<br> file is present. A fresh output directory alone is insufficient here<br> because the attached `IOUtil::recover_dir` names a separate<br> `output_directory_01/checkpoints` directory.<br>2. Confirm that the runtime parameter file contains<br><br> GRHayLib::EOS_type = "Tabulated"<br> GRHayLib::evolve_entropy = yes<br> GRHayLib::con2prim_routine = "Palenzuela1D"<br><br>3. Check whether `ID_converter_ILGRMHD` is active and whether startup prints<br> `Running GRHayLib with backward-compatible settings`. In current<br> IllinoisGRMHD, activating that thorn adds legacy Hybrid scheduling and<br> initializes Noble2D/Font1D compatibility settings independently of the<br> modern EOS-selected schedule, potentially creating conflicting paths. For<br> a modern GRHayLib/GRHayLID setup, remove it and confirm that only the<br> TabulatedEntropy C2P routine is scheduled. The attached file does not list<br> it, making this check especially useful for detecting submission of a<br> different generated parameter file or executable.<br>4. Keep `Carpet::ghost_size = 3`, as already specified in the attached file.<br> IllinoisGRMHD requires at least three ghost zones, and increasing this to<br> four will not repair invalid C2P input.<br>5. Leave the Multipole interpolation settings aside until the evolution<br> advances. Multipole analysis cannot repair an iteration-0 Prim2Con/C2P<br> problem.<br><br>For a more direct NaN check, add `NaNChecker` to `ActiveThorns` and use<br>something like<br><br> NaNChecker::check_after = 0<br> NaNChecker::check_every = 1<br> NaNChecker::check_for = "both"<br> NaNChecker::action_if_found = "just warn"<br> NaNChecker::verbose = "all"<br> NaNChecker::check_vars = "<br> ADMBase::metric ADMBase::lapse ADMBase::shift<br> HydroBase::rho HydroBase::press HydroBase::eps<br> HydroBase::entropy HydroBase::Y_e HydroBase::temperature<br> IllinoisGRMHD::grmhd_velocities<br> IllinoisGRMHD::u0<br> IllinoisGRMHD::Ax IllinoisGRMHD::Ay IllinoisGRMHD::Az<br> IllinoisGRMHD::phitilde<br> IllinoisGRMHD::grmhd_B_stagger<br> IllinoisGRMHD::grmhd_B_center<br> IllinoisGRMHD::grmhd_conservatives<br> IllinoisGRMHD::ent_star IllinoisGRMHD::Ye_star<br> "<br><br>NaNChecker reports fields at its scheduled check point; that timing depends on<br>the NaNChecker version and does not prove which routine first produced a<br>non-finite value. To distinguish invalid initial primitives from Prim2Con<br>output and C2P effects, check immediately before Prim2Con, immediately after<br>Prim2Con but before C2P, and immediately after C2P. The NaNChecker<br>documentation is:<br><br><a href="https://einsteintoolkit.org/thornguide/CactusUtils/NaNChecker/documentation.html">https://einsteintoolkit.org/thornguide/CactusUtils/NaNChecker/documentation.html</a><br><br>A finite rho diagnostic together with non-finite tau and momentum only weakly<br>narrows the cause. GRHayL computes the conserved density from fewer inputs,<br>whereas conserved energy and momentum use additional thermodynamic, velocity,<br>metric, and magnetic quantities:<br><br><a href="https://github.com/GRHayL/GRHayL/blob/8ef7e60884360ed4e6bef52c1be8fcff4558f3f5/GRHayL/Con2Prim/compute_conservs.c#L92-L98">https://github.com/GRHayL/GRHayL/blob/8ef7e60884360ed4e6bef52c1be8fcff4558f3f5/GRHayL/Con2Prim/compute_conservs.c#L92-L98</a><br><br>I would therefore check lapse and metric components, velocity/Lorentz factor,<br>pressure, specific internal energy, temperature, Y_e, and vector-potential or<br>magnetic-field components at iteration 0. Although the attached file sets<br>`Seed_Magnetic_Fields::A_b = 0`, that setting alone does not prove that every<br>magnetic array is initialized and finite, so I would not rank the EOS and<br>magnetic possibilities until those checks are available.<br><br>Please also report the Einstein Toolkit release and the exact GRHayL and<br>GRHayLET commit IDs used to build the executable. Older combinations may<br>predate relevant fixes. For example, GRHayLET pull request 14 fixed missing<br>error-control flow in the HybridEntropy, Tabulated, and TabulatedEntropy<br>recovery branches (not the Hybrid branch):<br><br><a href="https://github.com/GRHayL/GRHayLET/pull/14">https://github.com/GRHayL/GRHayLET/pull/14</a><br><br>Finally, the current Gallery `bns.par` uses GRHydro with an ideal-fluid EOS,<br>not IllinoisGRMHD with an SFHo table. Changing the hydrodynamics thorn and EOS<br>is therefore a substantial new configuration rather than a small change to a<br>validated Gallery setup:<br><br><a href="https://einsteintoolkit.org/gallery/bns/index.html">https://einsteintoolkit.org/gallery/bns/index.html</a><br><br><a href="https://einsteintoolkit.org/gallery/bns/bns.par">https://einsteintoolkit.org/gallery/bns/bns.par</a><br><br>The most useful items to send next would be the complete iteration-0 C2P<br>diagnostic block, the first warning/error and final output lines, the batch<br>exit status, the copied runtime parameter file, and the build commit IDs.<br><br>Best regards,<br></div><div><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:12.8px">-Zach</div><div style="font-size:12.8px"><br></div><span style="font-size:12.8px">* * *</span><br style="font-size:12.8px"><span style="font-size:12.8px">Zachariah Etienne</span></div><div><span style="font-size:12.8px">Prof. of Physics, U. of Idaho</span></div><div><span style="font-size:12.8px">Adjunct Prof. of Physics & Astronomy, West Virginia U.</span></div><div dir="ltr"><div><a href="https://etienneresearch.com" target="_blank">https://etienneresearch.com</a></div><div><a href="https://blackholesathome.net/" target="_blank">https://blackholesathome.net</a><br></div></div></div></div></div></div></div></div></div></div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 3, 2026 at 3:10 AM MAKINO Taiki via Users <<a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Dear Jordan</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you for your quick response.</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I understand what I should do. I will wait for feedback from the others, and I consider attending the weekly meeting.</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you again for your advice.</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Best regards,</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Taiki Makino</div>
<div style="text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Graduate School of Science, Nagoya University</div>
<div style="text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,120,212)">
<a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" title="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" style="color:rgb(0,120,212);margin-top:0px;margin-bottom:0px" target="_blank"><u>makino</u>.taiki.g6@s.mail.nagoya-u.ac.jp</a></div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="m_8705149526771489587m_3375230692233155978ms-outlook-mobile-signature">
<div class="MsoNormal" style="direction:ltr">取得 <a href="https://aka.ms/GetOutlookForMac" target="_blank">
Outlook for Mac</a></div>
</div>
<div id="m_8705149526771489587m_3375230692233155978mail-editor-reference-message-container">
<div style="direction:ltr">
</div>
<div style="padding:3pt 0in 0in;border-width:1pt medium medium;border-style:solid none none;border-color:rgb(181,196,223) currentcolor currentcolor">
<div style="text-align:left;font-family:Aptos;font-size:12pt;color:black">
<b>差出人: </b>Jordan Nicoules <<a href="mailto:jnicoules@ua.pt" target="_blank">jnicoules@ua.pt</a>><br>
<b>日付: </b>金曜日, 2026年7月3日 18:38<br>
<b>宛先: </b>MAKINO Taiki <<a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" target="_blank">makino.taiki.g6@s.mail.nagoya-u.ac.jp</a>><br>
<b>CC: </b><a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a> <<a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a>><br>
<b>件名: </b>Re: Question about IllinoisGRMHD C2P failure and kuibit Psi4 extraction<br>
<br>
</div>
</div>
<div id="m_8705149526771489587m_3375230692233155978divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<p style="margin-top:0px;margin-bottom:0px">Hi Taiki,</p>
<p style="margin-top:0px;margin-bottom:0px">(apologies for my previous typo)</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">The next step would be to understand why the simulation crashes and doesn't go past t=0. This is beyond what I can quickly hotfix I think, so
<b>I'm calling to other members to answer that.</b> Potentially, if you have the chance, the weekly calls can also be an opportunity to interact with people directly and put forth your questions.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Regarding files that have more than one line of data, that would make sense for the 1D profiles, but they still probably only correspond to iteration 0 (t=0).</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Best,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Jordan</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_8705149526771489587m_3375230692233155978divRplyFwdMsg" dir="ltr"><span style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(0,0,0)"><b>From:</b> MAKINO Taiki <<a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" target="_blank">makino.taiki.g6@s.mail.nagoya-u.ac.jp</a>><br>
<b>Sent:</b> Friday, July 3, 2026 10:14:57 AM<br>
<b>To:</b> Jordan Nicoules<br>
<b>Cc:</b> <a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a><br>
<b>Subject:</b> Re: Question about IllinoisGRMHD C2P failure and kuibit Psi4 extraction</span>
<div> </div>
</div>
<table style="width:100%">
<tbody>
<tr>
<td style="background-color:rgb(255,185,0);width:5px">
</td>
<td style="background-color:rgb(255,255,224);padding:5px">
<div style="font-family:Arial;font-size:10pt">
<b>CUIDADO:</b> Email de um sistema externo. Cuidado com links, anexos e pedidos de dados/senhas.<br>
<b>CAUTION:</b> Email from an external system. Be careful with links, attachments, and requests for data/passwords.</div>
</td>
</tr>
</tbody>
</table>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Dear Jordan</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you for your response.</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif">When</span><span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:16px"> I checked the result of my simulation, SimDir.gws[100][(2,2)].x was indeed 0 as you said and the simulation
stopped at iteration=0. </span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:16px"><br>
</span></p>
<p style="direction:ltr;line-height:normal;margin:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:16px">Also, I checked the psi4 ASCII files and all of them contained only one line. (However, the ASCII files that do not contain information of psi4 have several lines of data.)</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif">I</span><span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:16px"> am not sure what is causing this. What should I do next ?</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif;font-size:16px"><br>
</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif">Best regards,</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Arial,Helvetica,sans-serif"><br>
</span></p>
<p style="direction:ltr;text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px">
<span style="font-family:"Yu Gothic Regular";color:rgb(0,0,0)">Taiki Makino</span></p>
<p style="text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px">
<span style="font-family:"Yu Gothic Regular";color:rgb(0,0,0)">Graduate School of Science, Nagoya University</span></p>
<p style="text-align:left;text-indent:0px;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;color:rgb(0,120,212)"><a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" title="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" style="color:rgb(0,120,212);margin-top:0px;margin-bottom:0px" target="_blank"><u>makino</u>.taiki.g6@s.mail.nagoya-u.ac.jp</a></span></p>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="m_8705149526771489587m_3375230692233155978ms-outlook-mobile-signature">
<div class="MsoNormal" style="direction:ltr">取得 <a href="https://aka.ms/GetOutlookForMac" target="_blank">
Outlook for Mac</a></div>
</div>
<div id="m_8705149526771489587m_3375230692233155978mail-editor-reference-message-container">
<div style="padding:3pt 0in 0in;border-width:1pt medium medium;border-style:solid none none;border-color:rgb(181,196,223) currentcolor currentcolor">
<div style="text-align:left;font-family:Aptos;font-size:12pt;color:black">
<b>差出人: </b>Jordan Nicoules <<a href="mailto:jnicoules@ua.pt" target="_blank">jnicoules@ua.pt</a>><br>
<b>日付: </b>木曜日, 2026年7月2日 19:31<br>
<b>宛先: </b>MAKINO Taiki <<a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" target="_blank">makino.taiki.g6@s.mail.nagoya-u.ac.jp</a>>; <a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a> <<a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a>><br>
<b>件名: </b>Re: Question about IllinoisGRMHD C2P failure and kuibit Psi4 extraction<br>
<br>
</div>
</div>
<div id="m_8705149526771489587m_3375230692233155978divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif">
<div id="m_8705149526771489587m_3375230692233155978divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<p style="margin-top:0px;margin-bottom:0px">Hi Taiko,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">I will (partly) skip the first part of your message and let others reply on that, as I don't have experience with NS simulations.</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Regarding the gravitational wave extraction, I think you are accessing the data correctly (your example gives the time series for one mode at one radius), and there should be more than one element in general.
From your parameter file, it seems the output should go to an ascii file (it's the default), so you can double check in that file too.</p>
<p style="margin-top:0px;margin-bottom:0px">However, since you have `Multipole::out_every = 128`, does your simulation reach iteration 128 and later? In particular, because you're getting nans at iteration 0, does the simulation actually go past t=0? Given
what you're saying, I guess your single element of `SimDir.gws[100][(2,2)].x` is 0 (i.e t=0).</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Best,</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px">Jordan</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
<p style="margin-top:0px;margin-bottom:0px"><br>
</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_8705149526771489587m_3375230692233155978divRplyFwdMsg" dir="ltr"><span style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(0,0,0)"><b>From:</b> Users <<a href="mailto:users-bounces@einsteintoolkit.org" target="_blank">users-bounces@einsteintoolkit.org</a>> on behalf of MAKINO Taiki via Users <<a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a>><br>
<b>Sent:</b> Thursday, July 2, 2026 12:39:00 AM<br>
<b>To:</b> <a href="mailto:users@einsteintoolkit.org" target="_blank">users@einsteintoolkit.org</a><br>
<b>Subject:</b> [Users] Question about IllinoisGRMHD C2P failure and kuibit Psi4 extraction</span>
<div> </div>
</div>
<table style="width:100%">
<tbody>
<tr>
<td style="background-color:rgb(255,185,0);width:5px"></td>
<td style="background-color:rgb(255,255,224);padding:5px"><span style="font-family:Arial;font-size:10pt"><b>CUIDADO:</b> Email de um sistema externo. Cuidado com links, anexos e pedidos de dados/senhas.<br>
<b>CAUTION:</b> Email from an external system. Be careful with links, attachments, and requests for data/passwords.</span></td>
</tr>
</tbody>
</table>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span style="text-transform:none">Dear Einstein Toolkit Users,</span></div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I am a new Einstein Toolkit user, and I have some questions about a BNS simulation.</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt">
<span style="color:rgb(0,0,0)">I am trying to run a BNS simulation using bns.par, which I downloaded from Einstein Toolkit Gallery (</span><span style="color:rgb(0,120,212)"><a href="https://einsteintoolkit.org/gallery/bns/index.html" title="https://einsteintoolkit.org/gallery/bns/index.html" style="color:rgb(0,120,212)" target="_blank">https://einsteintoolkit.org/gallery/bns/index.html</a></span><span style="color:rgb(0,0,0)">)
and slightly modified.</span></div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I can successfully create and submit the simulation but the following message appears during the simulation:</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:40px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
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</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:40px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
To investigate this problem, I added the following parameters:</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:40px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Multipole::interpolator_name = "Lagrange polynomial interpolation"<br>
Multipole::interpolator_pars = "order=4"</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I also changed the parameter Driver::ghost_size from 3 to 4, but the same error occurred.</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;margin-left:0px;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I would like to ask:</div>
<ul style="direction:ltr;text-align:left;margin-top:0px;margin-bottom:0px;list-style-type:disc">
<li style="font-family:"Yu Gothic Regular";font-size:12pt;color:rgb(0,0,0);text-indent:0px">
<div role="presentation"><span style="text-transform:none">What could cause tau, Sx, Sy, and Sz to become NaN in this situation?</span></div>
</li><li style="font-family:"Yu Gothic Regular";font-size:12pt;color:rgb(0,0,0);text-indent:0px">
<div role="presentation"><span style="text-transform:none">Does this indicate a failure of C2P conversion?</span></div>
</li><li style="font-family:"Yu Gothic Regular";font-size:12pt;color:rgb(0,0,0);text-indent:0px">
<div role="presentation"><span style="text-transform:none">Are there any parameter settings or initial data that I should check first?</span></div>
</li></ul>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
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, </div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
SimDir.gws[100][(2,2)].y=[-9.79586237e-15-5.27476073e-14j].</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Other attributes such as</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
SimDir.gws[100][(2,2)].x </div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
also contain only one element. Is this a correct behavior or could it indicate a problem with my .par file? </div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
For reference, I have attached my parameter file.</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Thank you very much for your time and any suggestions.</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Best regards,</div>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<p style="text-align:left;text-indent:0px;text-transform:none;margin-top:0px;margin-bottom:0px">
<span style="font-family:"Yu Gothic Regular";font-size:12pt;color:rgb(0,0,0);text-transform:none">Taiki Makino</span></p>
<p style="text-align:left;text-indent:0px;text-transform:none;margin-top:0px;margin-bottom:0px">
<span style="font-family:"Yu Gothic Regular";font-size:12pt;color:rgb(0,0,0)">Graduate School of Science, Nagoya University</span></p>
<p style="text-align:left;text-indent:0px;text-transform:none;margin-top:0px;margin-bottom:0px">
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,120,212)"><a href="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" title="mailto:makino.taiki.g6@s.mail.nagoya-u.ac.jp" style="color:rgb(0,120,212);margin-top:0px;margin-bottom:0px" target="_blank"><u>makino</u>.taiki.g6@s.mail.nagoya-u.ac.jp</a></span></p>
<div style="direction:ltr;text-align:left;text-indent:0px;text-transform:none;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,120,212)">
<br>
</div>
<div style="font-family:"Yu Gothic Regular";color:rgb(33,33,33)">取得 <span style="color:rgb(0,120,212)">
<a href="https://aka.ms/GetOutlookForMac" title="https://aka.ms/GetOutlookForMac" style="color:rgb(0,120,212)" target="_blank">Outlook for Mac</a></span></div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="direction:ltr;font-family:Aptos,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="m_8705149526771489587m_3375230692233155978ms-outlook-mobile-signature">
<div class="MsoNormal" style="direction:ltr">取得 <a href="https://aka.ms/GetOutlookForMac" target="_blank">
Outlook for Mac</a></div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div>