[ET Trac] [Einstein Toolkit] #1599: CarpetInterp/waveinterp_2p tests constant in time initial data

Einstein Toolkit trac-noreply at einsteintoolkit.org
Tue Apr 29 00:25:37 CDT 2014


#1599: CarpetInterp/waveinterp_2p tests constant in time initial data
--------------------------+-------------------------------------------------
 Reporter:  rhaas         |       Owner:  eschnett           
     Type:  defect        |      Status:  new                
 Priority:  minor         |   Milestone:                     
Component:  Carpet        |     Version:  development version
 Keywords:  CarpetInterp  |  
--------------------------+-------------------------------------------------
 the test waveinterp_2p sets up initial data using IDWaveMoL (fom
 CactusExamples) using kx=ky=kz = 0. and slopet = 0.1. Ie the initial data
 is constant in space and has frequency 0.1. It then evolves this using
 WaveMoL and interpolates in space and time with CarpetInterp.

 Since the function is constant in space (exactly so, see eg the output in
 InterpToArray::array1d_vars[0]) the interpolation in space just tests
 roundoff errors. The second time derivative (array1d_vars[2]) should be
 zero and all we measure in the test is truncation errors on the level of
 1e-9. This makes the test very sensitive to -Ofast, -O0, gcc vs. intel
 etc. and not a very good test.

 I would rather either use a Gausian wave for initial data so that there is
 some actual data in the result and not just numerical noise, or leave out
 the 2nd time derivative from the test data.

 The test fails due to significant errors in a virtual machine using
 ubuntu.cfg (which usies -ffast-math and -O2).

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


More information about the Trac mailing list