Modify

Opened 6 months ago

Last modified 8 weeks ago

#2119 assigned defect

The binary neutron star gallery example gives different results with the release candiate.

Reported by: Peter Diener Owned by: Peter Diener
Priority: blocker Milestone: ET_2018_02
Component: Other Version: development version
Keywords: Cc: frank.loeffler@…

Description (last modified by Peter Diener)

The late part of the waveform plot looks very different in the release
candidate than on the gallery page. For l=2, m=2 mode of psi_4 extracted
at R=300 (that seems to be the correct data file as otherwise the scale on the y-axis don't match) the merger seem to happen slightly later and the waveform after merger is completely different.

The gallery page plot:

Gallery plot of psi_4

The new plot: New plot of psi_4

Attachments (3)

mp_Psi4_l2_m2_r300.00.png (11.2 KB) - added by Peter Diener 6 months ago.
Gallery plot of psi_4
mp_Psi4_l2_m2_r300.00_new.png (1.8 KB) - added by Peter Diener 6 months ago.
New plot of psi_4
mp_Psi4_l2_m2_r300.00_newer.png (2.1 KB) - added by Peter Diener 5 months ago.
With GRHydro::sources_spatial_order = 4

Download all attachments as: .zip

Change History (12)

Changed 6 months ago by Peter Diener

Attachment: mp_Psi4_l2_m2_r300.00.png added

Gallery plot of psi_4

Changed 6 months ago by Peter Diener

New plot of psi_4

comment:1 Changed 6 months ago by Peter Diener

Description: modified (diff)

comment:2 Changed 6 months ago by Peter Diener

I will try to add "GRHydro::sources_spatial_order = 4" and see if this fixes this problem as well.

comment:3 Changed 6 months ago by Roland Haas

This may well be the same issue as with the single NS #2118, namely having to set GRHydro::sources_spatial_order = 4 since http://einsteintoolkit.org/gallery/bns/nsnstohmns.par does not set it and the default is 2. Admittedly normally a too low order should be more dissipative and should lead to an earlier merger but still, having and order of 2 makes little sense. One would have to go back and check whether at the time the example was created (by me) GRHydro used ADMBase::spatial_order or GRHydro::sources_spatial_order.

Changed 5 months ago by Peter Diener

With GRHydro::sources_spatial_order = 4

comment:4 Changed 5 months ago by Peter Diener

With GRHydro::sources_spatial_order = 4 the plot of psi_4 look more similar to the gallery page than before. However,there are still differences after merger as can be seen from this plot:
With GRHydro::sources_spatial_order = 4

comment:5 Changed 5 months ago by roberto.depetri@…

Quite likely the new run use the c++ code instead of the fortran one (that was the default at that time) and small difference can be present. A better comparison (for the post-merger phase) should be to consider (also) the module of psi_4 instead of just the real part. In the post merger phase there are more than one modes with diffrent constructive-distructive interference) and a small phase differnce on the two-three excited mode may show up. Running the Gallery example on Marconi.

comment:6 Changed 5 months ago by Steven R. Brandt

Owner: set to Peter Diener
Status: newassigned

comment:7 Changed 4 months ago by Steven R. Brandt

Peter, what is the status of this ticket?

comment:8 Changed 2 months ago by Frank Löffler

Cc: frank.loeffler@… added

Putting in my 2 cents, but no solution. I am currently also investigating a difference in our NSNS production runs when we move from the 2015_05 release to the current one. Details are of course different, but the most notable one is also a difference in the post BH-formation GW form, so CCing me here.

Last edited 2 months ago by Frank Löffler (previous) (diff)

comment:9 Changed 8 weeks ago by Steven R. Brandt

Peter, any progress?

Modify Ticket

Change Properties
Set your email in Preferences
Action
as assigned The owner will remain Peter Diener.
Next status will be 'review'.
as The resolution will be set.
to The owner will be changed from Peter Diener to the specified user.
The owner will be changed from Peter Diener to anonymous.

Add Comment


E-mail address and name can be saved in the Preferences.

 
Note: See TracTickets for help on using tickets.