<div dir="ltr">Roland,<div>  thank you very much for the information. Sorry that I was not able to join the telecon, but Jay informed me about the discussion.</div><div><br></div><div>Thanks,</div><div>Bruno</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 12 nov 2020 alle ore 16:38 Roland Haas &lt;<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>&gt; ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Bruno, Jay, (to mark this email as &quot;responded&quot;)<br>
<br>
Jay brought this up in today&#39;s ET call:<br>
<br>
<a href="http://lists.einsteintoolkit.org/pipermail/users/2020-November/007670.html" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/pipermail/users/2020-November/007670.html</a><br>
<br>
and Peter, Zach, I responded:<br>
<br>
--8&lt;--<br>
Jay Kalinani: TOV simulation using new Spritz code. When adding a dipolar<br>
magnetic field to the star, there is some confusion about the behavior of<br>
cctk_time. There&#39;s a desire to not store RK intermediate values of<br>
cctk_time to cctk_time itself.<br>
Zach suggested that a local variable could be defined that simply outputs<br>
cctk_iteration*dt if desired.<br>
Roland had some comments about scheduling such a routine.<br>
--8&lt;--<br>
<br>
basically: &quot;this is as designed in case the RHS in has an explicit<br>
dependency on time&quot;.<br>
<br>
Yours,<br>
Roland<br>
<br>
&gt; Hi Erik,<br>
&gt;   those should be internal steps of the Runge-Kutta integrator when going<br>
&gt; from t to t+dt. Why is the cctk_time changed in those intermediate steps?<br>
&gt; Shouldn&#39;t cctk_time be set only at each iteration in time?<br>
&gt; <br>
&gt; Thanks,<br>
&gt; Bruno<br>
&gt; <br>
&gt; <br>
&gt; Il giorno mer 11 nov 2020 alle ore 15:11 Erik Schnetter &lt;<br>
&gt; <a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>&gt; ha scritto:  <br>
&gt; <br>
&gt; &gt; Jay<br>
&gt; &gt;<br>
&gt; &gt; This intermediate step comes from the time integrator. The Einstein<br>
&gt; &gt; Toolkit usually uses Runge-Kutta type time integrators that perform<br>
&gt; &gt; multiple sub-steps for each time step. These time integrators are<br>
&gt; &gt; provided by the &quot;Methods of Lines&quot; thorn (&quot;MoL&quot;). The respective<br>
&gt; &gt; documentation should have further pointers, and the respective<br>
&gt; &gt; parameter settings should also tell you more. (See also e.g.<br>
&gt; &gt; &lt;<a href="https://en.wikipedia.org/wiki/Runge%E2%80%93Kutta_methods" rel="noreferrer" target="_blank">https://en.wikipedia.org/wiki/Runge–Kutta_methods</a>&gt;).<br>
&gt; &gt;<br>
&gt; &gt; -erik<br>
&gt; &gt;<br>
&gt; &gt; On Tue, Nov 10, 2020 at 12:20 PM Jay Vijay Kalinani<br>
&gt; &gt; &lt;<a href="mailto:jayvijay.kalinani@phd.unipd.it" target="_blank">jayvijay.kalinani@phd.unipd.it</a>&gt; wrote:  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Dear all,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I am trying to perform a magnetised TOV simulation using Einstein  <br>
&gt; &gt; Toolkit and Spritz GRMHD code. In particular, currently I have a subroutine<br>
&gt; &gt; that adds a dipolar magnetic field to the TOV star after its initial data<br>
&gt; &gt; has been generated via TOVSolver thorn.  <br>
&gt; &gt; &gt; Within this subroutine, I have added the following lines that print the  <br>
&gt; &gt; &#39;cctk_time&#39; value as:  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; character(len=200) warnline<br>
&gt; &gt; &gt; WRITE( warnline, &#39;( &quot;cctk_time: &quot;, g16.7)&#39; ) cctk_time<br>
&gt; &gt; &gt; call CCTK_INFO(warnline)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; This subroutine is scheduled in HydroBase_Con2Prim after Con2Prim step.  <br>
&gt; &gt; And I get the following lines in the .out file as:  <br>
&gt; &gt; &gt; .........<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (TOVSolver): Integrated TOV equation<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (TOVSolver): Information about the TOVs used:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (): TOV    radius    mass  bary_mass mass(g) cent.rho rho(cgi)  <br>
&gt; &gt;   K   K(cgi)    Gamma  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO ():   1   8.12502  1.40016  1.50618 2.78e+33  0.00128 7.92e+14  <br>
&gt; &gt; 100 1.45e+05        2  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (TOVSolver): Not using old matter initial data<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (TOVSolver): Done interpolation.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:     0.000000<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): Adding a poloidal magnetic field to the TOV  <br>
&gt; &gt; star on reflevel 0  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): Adding atmo to the TOV star on reflevel =  <br>
&gt; &gt; &#39;0&#39;  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:     0.000000<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  <br>
&gt; &gt; -----------------------------------------------------------------------------------------------------------------------------------  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Iteration      Time |              ADMBASE::alp |  <br>
&gt; &gt; HYDROBASE::rho |          HYDROBASE::press |            SPRITZ::B_norm  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                     |      minimum      maximum |      minimum  <br>
&gt; &gt; maximum |      minimum      maximum |      minimum      maximum  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  <br>
&gt; &gt; -----------------------------------------------------------------------------------------------------------------------------------  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;         0     0.000 |    0.6698612    0.9603815 | 5.827798e-11  <br>
&gt; &gt; 0.0012800 |    0.0000000    0.0001638 | 8.399814e-09    0.0000399  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.3906250E-01<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.3906250E-01<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.7812500E-01<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.7812500E-01<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;         1     0.078 |    0.6698612    0.9603833 | 5.827814e-11  <br>
&gt; &gt; 0.0012799 |    0.0000000    0.0001638 | 8.294987e-09    0.0000395  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.1171875<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.1171875<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.1562500<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; INFO (Spritz_DipoleAndAtmo): cctk_time:    0.1562500<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;         2     0.156 |    0.6698612    0.9603850 | 5.827829e-11  <br>
&gt; &gt; 0.0012797 |    0.0000000    0.0001638 | 8.294994e-09    0.0000390  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ..........<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; The simulation is performed in uniform grid with dx=dy=dz=0.3125,  <br>
&gt; &gt; time::dtfac = 0.25 and HydroBase::timelevels = 2.  <br>
&gt; &gt; &gt; As you may see from the output lines, in between iteration 0 which is at  <br>
&gt; &gt; cctk_time=0.0 and iteration 1 which is at cctk_time=0.078125, there is an<br>
&gt; &gt; intermediate time step with &#39;cctk_time=0.0390625&#39; when the subroutine is<br>
&gt; &gt; called.  <br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I was hoping to understand better why this intermediate time step exists  <br>
&gt; &gt; between two consecutive iterations, and if there is any documentation I<br>
&gt; &gt; could look at.  <br>
&gt; &gt; &gt; Any help would be much appreciated.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Thank you very much.<br>
&gt; &gt; &gt; Sincerely,<br>
&gt; &gt; &gt; Jay Kalinani<br>
&gt; &gt; &gt; _______________________________________________<br>
&gt; &gt; &gt; Users mailing list<br>
&gt; &gt; &gt; <a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
&gt; &gt; &gt; <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a>  <br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; Erik Schnetter &lt;<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>&gt;<br>
&gt; &gt; <a href="http://www.perimeterinstitute.ca/personal/eschnetter/" rel="noreferrer" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; Users mailing list<br>
&gt; &gt; <a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
&gt; &gt; <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
&gt; &gt;  <br>
&gt; <br>
&gt; <br>
<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><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&#39;t<br></span><span style="font-size:12.8px">----------------------------------------------------------------------<br></span></font></p></div></div></div></div></div></div></div>