<div dir="ltr">Dear Lorenzo,<div><br></div><div>It is to be expected that the region with the highest density and gravity is the one with the highest constraint violation.</div><div><br></div><div>I think that the problem you are seeing is an output one. </div><div><br></div><div>(Disclaimer, this is my understanding; I might very well be wrong)</div><div><br></div><div>You have several refinement levels. Different refinement levels evolve at a different pace (determined by your</div><div>time_refinement_factors, which by default is essentially proportionally to your local dx), and iteration #4 (that you</div><div>are plotting) is not one where all the refinement levels are synced. </div><div><br></div><div>My guess is that the stars will appear if you plot iterations 128.</div><div><br></div><div>You might also be able use kuibit to inspect your HDF5 files and see what refinement levels are in the data.</div><div><br></div><div>That would be something like:</div><div>```</div><div>import kuibit</div><div>print(SimDir("PATH_OF_DATA_FOLDER").gf.xy["rho"][4].available_refinement_levels)</div><div>```</div><div><br></div><div>Best,</div><div>Gabriele</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jan 18, 2024 at 3:30 PM Lorenzo Cipriani <<a href="mailto:lorenzo.cipriani@graduate.univaq.it">lorenzo.cipriani@graduate.univaq.it</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 class="msg8950940236650797097">
<div dir="ltr">
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">I noticed that the Hamiltonian constraint violation along the x-axis (y=0, z=0) (in
the attached image, taken from the H.xy.h5 file), has a very peculiar shape, aligned with the stars profile. This is very wrong, am I right?</span></div>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">So, the problem is either in the EoS, more likely, or Lorene.</span></div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div id="m_2829804400576971955Signature">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Lorenzo</div>
</div>
<div id="m_2829804400576971955appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="m_2829804400576971955divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>Da:</b> Bruno Giacomazzo <<a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a>><br>
<b>Inviato:</b> giovedì 18 gennaio 2024 15:13<br>
<b>A:</b> Lorenzo Cipriani <<a href="mailto:lorenzo.cipriani@graduate.univaq.it" target="_blank">lorenzo.cipriani@graduate.univaq.it</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>Oggetto:</b> Re: [Users] Star disappear from grid</font>
<div> </div>
</div>
<div>
<div dir="ltr">Lorenzo,
<div> I opened the hdf5 file you sent me and I indeed see the same problem (I attached a couple of snapshots). So it is not a problem of resolution. Did you check that hydro quantities and spacetime variables are initialized correctly (just in case there
is a problem with the initial data file)?</div>
<div><br>
</div>
<div>Cheers,</div>
<div>Bruno</div>
<div><br>
</div>
</div>
<br>
<div>
<div dir="ltr">Il giorno gio 18 gen 2024 alle ore 14:55 Lorenzo Cipriani <<a href="mailto:lorenzo.cipriani@graduate.univaq.it" target="_blank">lorenzo.cipriani@graduate.univaq.it</a>> ha scritto:<br>
</div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div dir="ltr">
<div><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)">Dear Bruno,</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)"><br>
</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)">I increased
the resolution so that the finest level now has a resolution of 7/2^5 = 0.21875 using a smaller coarse grid, but the output remains unchanged.</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)"><br>
</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)">Please
let me know your thoughts on this!</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)"><br>
</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)">Best,</span></div>
<div style="text-align:left;margin:0px;background-color:rgb(255,255,255)"><span style="letter-spacing:normal;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;font-weight:400;color:rgb(0,0,0)">Lorenzo</span></div>
<div id="m_2829804400576971955x_m_7766451353610490557Signature"></div>
<div id="m_2829804400576971955x_m_7766451353610490557appendonsend"></div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<hr style="display:inline-block;width:98%">
<div id="m_2829804400576971955x_m_7766451353610490557divRplyFwdMsg" dir="ltr"><span style="font-family:Calibri,sans-serif;font-size:11pt;color:rgb(0,0,0)"><b>Da:</b> Bruno Giacomazzo <<a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a>><br>
<b>Inviato:</b> mercoledì 17 gennaio 2024 16:44<br>
<b>A:</b> Lorenzo Cipriani <<a href="mailto:lorenzo.cipriani@graduate.univaq.it" target="_blank">lorenzo.cipriani@graduate.univaq.it</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>Oggetto:</b> Re: [Users] Star disappear from grid</span>
<div> </div>
</div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Lorenzo,</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"> the perfectly cut away may be due to the fact that you may have had failures in con2prim on
the finest levels (and maybe points got set to atmosphere), but the coarser levels were not computed yet. Your standard output reports indeed errors in c2p. You use 6 refinement levels and therefore all levels are at the same time every 2^(6-1)=32 iterations.
At iteration 4 only the 3 most inner levels would be at the same.</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"> I would increase the resolution since your finest grid is using a resolution of 12/2^5 = 0.375
so that you have a resolution of at least 0.24 on the finest level (level 5) and check if this solves the problem.</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Cheers,</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Bruno</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Il giorno mer 17 gen 2024 alle ore 15:41 Lorenzo Cipriani <<a href="mailto:lorenzo.cipriani@graduate.univaq.it" id="m_2829804400576971955x_m_7766451353610490557OWA1bf6e79b-cf0f-ea66-91b3-84899bea220e" target="_blank">lorenzo.cipriani@graduate.univaq.it</a>>
ha scritto:</span></div>
<blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left:1px solid rgb(204,204,204)">
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Dear all,</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">my name is Lorenzo Cipriani, PhD student at the University of L'Aquila.</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">My goal is to simulate a BNS merger using initial data computed from Lorene and evolved according
to the parameter file I have attached.</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Cactus reads the initial data correctly, as you can see in the image "rhoxy_initial.png" but, immediately
after, the stars practically disappear from the grid, perfectly cut away (see image "rhoxy_it4.png").</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Does anybody have a hint on what might be happening? I have never encountered a behavior like this.</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Thank you for your assistance,</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">best regards,</span></div>
<div style="direction:ltr"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">Lorenzo</span></div>
<div id="m_2829804400576971955x_m_7766451353610490557x_m_-5433044787563394965Signature"></div>
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" id="m_2829804400576971955x_m_7766451353610490557OWAdcb2ba5f-8108-93f0-948a-9c16014ff442" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" id="m_2829804400576971955x_m_7766451353610490557OWA7c208f9c-8eb4-f34b-a956-1c8e1ae3582d" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</span></blockquote>
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span>
<div><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span></div>
<span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">--<br>
</span>
<p style="direction:ltr;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,0,0)">Pr</span><span style="font-family:arial,sans-serif;font-size:12pt;color:rgb(0,0,0)">of.
Bruno Giacomazzo<br>
Department of Physics<br>
University of Milano-Bicocca<br>
Piazza della Scienza 3<br>
20126 Milano<br>
</span><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12.8px;color:rgb(0,0,0)">Italy</span></p>
<p style="direction:ltr;margin-top:0px;margin-bottom:0px"><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12.8px;color:rgb(0,0,0)">email:
<a href="mailto:bruno.giacomazzo@unimib.it" id="m_2829804400576971955x_m_7766451353610490557OWA54944286-b4bc-c311-8dd5-59da98cf5643" style="margin-top:0px;margin-bottom:0px" target="_blank">
bruno.giacomazzo@unimib.it</a><br>
phone: (+39) 02 6448 2321</span><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)"><br>
</span><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12.8px;color:rgb(0,0,0)">web:
<a href="http://www.brunogiacomazzo.org/" id="m_2829804400576971955x_m_7766451353610490557OWAa5dd11a2-0918-4efa-230b-1890352c693d" style="margin-top:0px;margin-bottom:0px" target="_blank">
http://www.brunogiacomazzo.org</a></span></p>
<p style="direction:ltr;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,0,0)">---------------------------------------------------------------------<br>
</span><span style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12.8px;color:rgb(0,0,0)">There are only 10 types of people in the world:<br>
Those who understand binary, and those who don't<br>
----------------------------------------------------------------------</span></p>
</div>
</div>
</blockquote>
</div>
<br clear="all">
<div><br>
</div>
<span>-- </span><br>
<div dir="ltr">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<p><font color="#000000">Pr<font face="arial, sans-serif">of. Bruno Giacomazzo<br>
Department of Physics<br>
University of Milano-Bicocca<br>
</font></font><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">Piazza della Scienza 3<br>
</font></span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">20126 Milano<br>
</font></span><span style="font-size:12.8px;color:rgb(0,0,0)">Italy</span></p>
<p><font color="#000000"><span style="font-size:12.8px">email: </span><span style="font-size:12.8px"><a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a><br>
</span><span style="font-size:12.8px">phone: (+39) 02 6448 2321</span></font><br>
<font color="#000000"><span style="font-size:12.8px">web: </span></font><span style="color:rgb(0,0,0);font-size:12.8px"><a href="http://www.brunogiacomazzo.org/" style="font-size:12.8px" target="_blank">http://www.brunogiacomazzo.org</a></span></p>
<p><font color="#000000">---------------------------------------------------------------------<br>
<span style="font-size:12.8px">There are only 10 types of people in the world:<br>
</span><span style="font-size:12.8px">Those who understand binary, and those who don't<br>
</span><span style="font-size:12.8px">----------------------------------------------------------------------<br>
</span></font></p>
</div>
</div>
</div>
</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>
</div></blockquote></div>