<div dir="ltr"><div><div>Hello Vassili,<br>No, I don&#39;t think I&#39;m using it. I&#39;ll try it to see if it fixes the problem. <br></div>Thank you. <br></div>Petr.<br><div><div><div><div id="dc_vk_code" style></div></div>
</div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Nov 26, 2013 at 12:18 PM, Vassilios Mewes <span dir="ltr">&lt;<a href="mailto:vassilios.mewes@uv.es" target="_blank">vassilios.mewes@uv.es</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Petr,<div><br></div><div>are you using the </div><div><br></div><div>GRHydro::reconstruct_Wv = yes<br>
</div><div><br></div><div>parameter? </div><div>I need to use it in my accretion disk runs when using the Ideal_Fluid EOS...Using the polytrope this parameter wasn&#39;t needed if I recall correctly..</div>
<div><br></div><div>hope this helps..</div><div><br></div><div>best wishes,</div><div><br></div><div>Vassili</div></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On Tue, Nov 26, 2013 at 6:00 PM, Petr Tsatsin <span dir="ltr">&lt;<a href="mailto:ptsatsin@my.fau.edu" target="_blank">ptsatsin@my.fau.edu</a>&gt;</span> wrote:<br>

</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr"><div>Hello,<br></div>I&#39;ve recently tried to perform simulation of the binary neutron stars started from Meudon initial data. I&#39;ve encountered a following problem. I can run the simulations without any problems when I use polytropic EOS. However, when I switch to Gamma-Law simulations crushed  soon after start somewhere in the conversion of primitive quantities to conserved quantities. The typical log is look like this one: <br>


 <br clear="all"><div><div>WARNING[L1,P28] (GRHydro): Con2Prim failed; stopping the code.<br>WARNING[L1,P28] (GRHydro): Even with mesh refinement, this point is not restricted from a finer level, so this is really an error<br>


WARNING[L1,P28] (GRHydro):             on carpet reflevel:  2<br>WARNING[L1,P28] (GRHydro):                   xyz location:     12.00000        0.000000       -2.000000<br>WARNING[L1,P28] (GRHydro):                         radius:     12.16553<br>


WARNING[L1,P28] (GRHydro): hydro_excision_mask, C2Pfailed:   0   1.0000000<br>WARNING[L1,P28] (GRHydro):        rho eps press w_lorentz:    0.1047685E-08   0.5280317E-01   0.5532109E-10    1.013257<br>WARNING[L1,P28] (GRHydro):                     velocities:   -0.3371566E-02   0.1176264      -0.3492353E-01<br>


WARNING[L1,P28] (GRHydro):                       dens tau:    0.2037191E-08   0.1388378E-09<br>WARNING[L1,P28] (GRHydro):                           scon:   -0.2835004E-10   0.4723716E-09  -0.1051106E-09<br>WARNING[L1,P28] (GRHydro):                       3-metric:     1.420906      -0.7104738E-01  -0.2077218E-01    1.770598       0.4372303E-01    1.468061<br>


WARNING[L1,P28] (GRHydro):                   lapse, shift:    0.8118283       0.1691328E-01  -0.3824454E-01  -0.5210642E-03<br>WARNING level 0 in thorn GRHydro processor 28 host <a href="http://c559-402.stampede.tacc.utexas.edu" target="_blank">c559-402.stampede.tacc.utexas.edu</a><br>


  (line 2404 of GRHydro_Con2Prim.F90):<br>  -&gt; Aborting.<br>cactus_SBNS: /work/02566/tsatsin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:314: int Carpet::Abort(const _cGH *, int): Assertion `0&#39; failed.<br><br>


</div><div>The parameter files for &quot;polytropic&quot; runs and &quot;Gamma-Law&quot; are identical except, obviously, the parameters controlling the EOS.<br></div><div>I was wondering if somebody experienced similar problems before? I would be very grateful for any help and suggestions.<br>


</div><div>I also include the some relevant parameters I used for &quot;Gamma-law&quot; EOS below. May be I&#39;m missing some important parameters or values are not properly set?<br><br> # Hydro paramters<br>ActiveThorns = &quot;EOS_Omni&quot;<br>


EOS_Omni::poly_k = 123.6<br>EOS_Omni::poly_gamma = 2.0<br>EOS_Omni::gl_k = 123.6<br>EOS_Omni::gl_gamma = 2.0<br><br>ActiveThorns = &quot;GRHydro&quot;<br><br>GRHydro::initial_Gamma = 2.0<br><br>HydroBase::evolution_method      = &quot;GRHydro&quot;<br>


<br>#GRHydro::riemann_solver            = &quot;Marquina&quot;<br>GRHydro::riemann_solver     = &quot;HLLE&quot;<br>GRHydro::GRHydro_eos_type   = &quot;General&quot;<br>GRHydro::GRHydro_eos_table  = &quot;Ideal_Fluid&quot;<br>


GRHydro::recon_method       = &quot;ppm&quot;<br>GRHydro::GRHydro_stencil    = 3<br>GRHydro::bound              = &quot;flat&quot;<br>GRHydro::rho_abs_min        = 1.e-10<br>#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau<br>


#GRHydro::GRHydro    = 10    # gij(6), alpha, beta(3)<br><br>GRHydro::c2p_reset_pressure        = &quot;yes&quot;<br>grhydro::grhydro_c2p_warnlevel         =  0<br>grhydro::grhydro_c2p_warn_from_reflevel      =  2<br><br>


<br>ActiveThorns = &quot;GenericFD NewRad&quot;<br>ActiveThorns = &quot;ML_BSSN ML_BSSN_Helper&quot;<br>ActiveThorns = &quot;ML_ADMConstraints&quot;<br><br><br></div><div>Thank you and Happy Thanksgiving.<span><font color="#888888"><br>

</font></span></div><span><font color="#888888"><div>Petr.<br>
</div><div>-- <br><div dir="ltr"><div>Petr Tsatsin</div><div>Graduate Student - Department of Physics</div><div>Charles E. Schmidt College of Science</div><div>Florida Atlantic University</div><div><a href="tel:%28561%29-297-3386" value="+15612973386" target="_blank">(561)-297-3386</a></div>

<div>
cell: <a href="tel:1-561-299-8048" value="+15612998048" target="_blank">1-561-299-8048</a></div><div><a href="mailto:ptsatin@fau.edu" target="_blank">ptsatin@fau.edu</a>, <a href="mailto:ptsatsin@gmail.com" target="_blank">ptsatsin@gmail.com</a></div>

<div></div><div></div><div></div><div></div>
<div></div></div>
</div></font></span></div></div>
<br></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" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div>Petr Tsatsin</div><div>Graduate Student - Department of Physics</div><div>Charles E. Schmidt College of Science</div><div>Florida Atlantic University</div>
<div>(561)-297-3386</div><div>cell: 1-561-299-8048</div><div><a href="mailto:ptsatin@fau.edu" target="_blank">ptsatin@fau.edu</a>, <a href="mailto:ptsatsin@gmail.com" target="_blank">ptsatsin@gmail.com</a></div><div></div>
<div></div><div></div><div></div><div></div></div>
</div>