[ET Trac] [Einstein Toolkit] #1199: [PATCH]McLachlan_BSSN: fix the Hamiltonian constraint for CCZ4

[Einstein Toolkit]

#1199: [PATCH]McLachlan_BSSN: fix the Hamiltonian constraint for CCZ4
---------------------------------+------------------------------------------
 Reporter:  anton@…              |        Type:  defect               
   Status:  new                  |    Priority:  major                
Milestone:                       |   Component:  EinsteinToolkit thorn
  Version:  development version  |    Keywords:                       
---------------------------------+------------------------------------------
 A typo in a8794cd2227c87c7457ba51d040e863be2f981e7 has caused the
 additional CCZ4 terms to overwrite the BSSN value instead of adding to it.

 Tested to work on TwoPunctures initial data:

 h  = h5py.File('/tmp/testbssn/H.h5')  # BSSN reference
 h1 = h5py.File('/tmp/testccz4/H.h5')  # with the patch
 h2 = h5py.File('/tmp/testccz42/H.h5') # without the patch

 In [89]: numpy.sum(numpy.abs(h.values()[0][:] - h1.values()[0][:]))
 Out[89]: 8.0082674998918897e-15

 In [90]: numpy.sum(numpy.abs(h.values()[0][:] - h2.values()[0][:]))
 Out[90]: 7.9536720382719803

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1199>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit