#2119: The binary neutron star gallery example gives different results with the release candiate. ---------------------------+--------------------------------- Reporter: Peter Diener | Owner: Peter Diener Type: defect | Status: assigned Priority: blocker | Milestone: ET_2018_02 Component: Other | Version: development version Resolution: | Keywords: ---------------------------+---------------------------------
Comment (by Peter Diener):
A part of the difference comes from a change to Meudon_Bin_NS to use EOS_Omni in commit: 1b9d6a3eff34224a357a61e3b74d5468f8efa82b. After this change some variables have slightly different profiles on the initial data slice. This affects: hydrobase::press, hydrobase::eps, grhydro::tau and grhydro::scon. The following plot shows the maximum density as function of time for 3 different runs. The red curve is for the original code at the time of the creation of the gallery example. The green curve is for the current development version of the toolkit. The blue curve is for the current development version with Meudon_Bin_NS reverted to the commit before the EOS_Omni commit.
[[Image(dens_compare.png)]]
As can be seen the blue and red curves agrees much better up to the merger, but still some differences appear after that.
The affect on the Psi_4 waveform can be seen here where I have zoomed in to the merger and post-merger phase:
[[Image(psi4_compare.png)]]
Also in this case, the agreement is better between the red and blue curves, but significant differences still appear late in the waveform.
I don't know if Meudon_Bin_NS is more correct before or after the commit mentioned above. And I don't know what then causes the additional differences.