<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd">
<html><head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="generator" content="Osso Notes">
<title></title></head>
<body>
<p>Hi Yosef,
<br>
<br>I believe the code was reorganised in carpet and this needed the algorithm to change in a way that meant numerical results would change.  I can't check right now, but this should have been mentioned in the announcement of the new version of carpet as well as the ET release notes.
<br>
<br><font color="#999999">-- </font>
<br><font color="#999999">Ian Hinder</font>
<br><font color="#999999"><a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a></font>
<br><font color="#999999"></font>
<br>
<br>----- Original message -----
<br>> While trying to understand why the results from a BBH run
<br>> differed when using the 2011_05 and 2011_10 version of ET
<br>> I noticed that the two versions gave different results on every
<br>> other point of the buffer region (in a 2 level run) on the 2nd
<br>> time level.
<br>>
<br>> I reran the test using McLachlan and I found a similar
<br>> behavior, but for fewer points (here, only 3, rather
<br>> than 8) and only the innermost points.
<br>>
<br>>
<br>> diff -w -u ~/ET/Old/Cactus/test/mc/gxx.x.asc ~/ET/Develop/Cactus/test/mc
<br>> | grep "^[+-]2"
<br>> -2  0 0 0 0 26 8 8  0.1  1.8 0 0 5.86163790879802
<br>> +2  0  0 0 0  26 8 8  0.1  1.8 0 0 5.86163790879801
<br>> -2  0 1 0 0 26 8 8  0.1  1.8 0 0 5.86163790879802
<br>> +2  0  1 0 0  26 8 8  0.1  1.8 0 0 5.86163790879801
<br>> -2  0 1 0 0 59 8 8  0.1  5.1 0 0 2.06954055131151
<br>> +2  0  1 0 0  59 8 8  0.1  5.1 0 0 2.06940807802512
<br>> -2  0 1 0 0 61 8 8  0.1  5.3 0 0 2.01672019300839
<br>> +2  0  1 0 0  61 8 8  0.1  5.3 0 0 2.01674238678588
<br>> -2  0 1 0 0 63 8 8  0.1  5.5 0 0 1.96875965038995
<br>> +2  0  1 0 0  63 8 8  0.1  5.5 0 0 1.96875664800624
<br>>
<br>>
<br>>
<br>> One of the major differences between these two versions of ET is the
<br>> switch from the git to hg version of carpet.
<br>>
<br>> Here is the full parfile that I used.
<br>>
<br>> #==============================================================================
<br>>
<br>> ActiveThorns = "admanalysis admbase admcoupling admmacros coordgauge
<br>> spacemask StaticConformal  boundary time cartgrid3d ioutil iobasic
<br>> nanchecker CoordBase aeilocalinterp Slab SphericalSurface LocalReduce
<br>> MoL Carpet CarpetInterp CarpetIOASCII CarpetLib CarpetReduce CarpetSlab
<br>> CarpetRegrid2 SphericalSurface CarpetIOHDF5 CarpetIOScalar InitBase
<br>> SymBase LoopControl GSL ReflectionSymmetry"
<br>>
<br>>
<br>> ActiveThorns = "TimerReport"
<br>>
<br>> TimerReport::out_every  = 512
<br>> TimerReport::out_filename = "TimerReport"
<br>> Carpet::output_timers_every = 512
<br>> #------------------------------------------------------------------------------
<br>>
<br>> # Run parameters
<br>> #------------------------------------------------------------------------------
<br>>
<br>> # Grid
<br>>
<br>> mol::ode_method  = "RK4"
<br>> mol::mol_intermediate_steps  = 4
<br>> mol::mol_num_scratch_levels  = 1
<br>>
<br>>
<br>>
<br>> CoordBase::domainsize  = "minmax"
<br>> CoordBase::xmin  = 0
<br>> CoordBase::ymin  = 0
<br>> CoordBase::zmin  = 0
<br>> CoordBase::xmax  = 10
<br>> CoordBase::ymax  = 10
<br>> CoordBase::zmax  = 10
<br>> CoordBase::spacing  = "numcells"
<br>> CoordBase::ncells_x  = 50
<br>> CoordBase::ncells_y  = 50
<br>> CoordBase::ncells_z  = 50
<br>> CartGrid3D::type  = "coordbase"
<br>> CartGrid3D::domain  = "full"
<br>> CartGrid3D::bitant_plane  = "xy"
<br>> CartGrid3D::avoid_originx  = no
<br>> CartGrid3D::avoid_originy  = no
<br>> CartGrid3D::avoid_originz  = no
<br>>
<br>> CoordBase::boundary_size_x_lower  = 4
<br>> CoordBase::boundary_size_y_lower  = 4
<br>> CoordBase::boundary_size_z_lower  = 4
<br>> CoordBase::boundary_shiftout_x_lower  = 1
<br>> CoordBase::boundary_shiftout_y_lower  = 1
<br>> CoordBase::boundary_shiftout_z_lower  = 1
<br>>
<br>> CoordBase::boundary_size_x_upper  = 4
<br>> CoordBase::boundary_size_y_upper  = 4
<br>> CoordBase::boundary_size_z_upper  = 4
<br>>
<br>> ReflectionSymmetry::reflection_x = yes
<br>> ReflectionSymmetry::reflection_y = yes
<br>> ReflectionSymmetry::reflection_z = yes
<br>> ReflectionSymmetry::avoid_origin_x = no
<br>> ReflectionSymmetry::avoid_origin_y = no
<br>> ReflectionSymmetry::avoid_origin_z = no
<br>>
<br>>
<br>>
<br>> driver::ghost_size  = 4
<br>>
<br>> #---------------------------------------------------------------------------
<br>>
<br>> cactus::cctk_itlast  = 16
<br>>
<br>>
<br>>
<br>> time::dtfac  = .5
<br>> Carpet::domain_from_coordbase  = yes
<br>> Carpet::prolongation_order_space  = 5
<br>> Carpet::prolongation_order_time  = 2
<br>> Carpet::max_refinement_levels  = 2
<br>> Carpet::use_buffer_zones = yes
<br>> Carpet::verbose  = no
<br>> Carpet::schedule_barriers = no
<br>> Carpet::veryverbose  = no
<br>>
<br>> CarpetLib::output_bboxes  = no
<br>> CarpetLib::check_bboxes  = no
<br>> CarpetLib::interleave_communications = yes
<br>> CarpetLib::combine_sends  = yes
<br>>
<br>> #Carpet::time_refinement_factors = "[1,1,1,1,2,2,4,8,16,32,64]"
<br>> Carpet::init_3_timelevels = "no"
<br>> Carpet::init_each_timelevel = "no"
<br>> Carpet::init_fill_timelevels = "yes"
<br>> Carpet::enable_all_storage  = no
<br>> Carpet::regrid_during_recovery = no
<br>>
<br>> ActiveThorns="CarpetRegrid2"
<br>> CarpetRegrid2::regrid_every = 16
<br>> CarpetRegrid2::num_centres = 1
<br>> CarpetRegrid2::num_levels_1 = 2
<br>> CarpetRegrid2::radius_1[ 1] = 5
<br>>
<br>>
<br>> # MODEL
<br>> ActiveThorns="IDAnalyticBH"
<br>> ADMBase::initial_data="bl_bh"
<br>> IDAnalyticBH::bl_nbh = 2
<br>> IDAnalyticBH::bl_x0_1 = 1.01001
<br>> IDAnalyticBH::bl_y0_1 = 1.01001
<br>> IDAnalyticBH::bl_x0_2 = -1.01001
<br>> IDAnalyticBH::bl_y0_2 = -1.01001
<br>> initbase::initial_data_setup_method = init_some_levels
<br>>
<br>>
<br>>
<br>>
<br>> admbase::metric_type  = "Physical"
<br>> #staticconformal::conformal_storage  = "factor+derivs+2nd derivs"
<br>>
<br>> #--- numerics
<br>>
<br>>
<br>> ActiveThorns = "ML_BSSN ML_BSSN_Helper TmunuBase GenericFD"
<br>>
<br>> ADMBase::evolution_method  = "ML_BSSN"
<br>> ADMBase::lapse_evolution_method = "ML_BSSN"
<br>> ADMBase::shift_evolution_method = "ML_BSSN"
<br>>
<br>> ML_BSSN::timelevels = 3
<br>>
<br>> ML_BSSN::harmonicN  = 1  # 1+log
<br>> ML_BSSN::harmonicF  = 2.0  # 1+log
<br>> ML_BSSN::ShiftGammaCoeff  = 0.75
<br>> ML_BSSN::BetaDriver  = 1.0
<br>> ML_BSSN::LapseAdvectionCoeff = 1
<br>> ML_BSSN::ShiftAdvectionCoeff = 1
<br>>
<br>> ML_BSSN::my_boundary_condition = "Minkowski"
<br>>
<br>> ML_BSSN::ML_log_confac_bound = "none"
<br>> ML_BSSN::ML_metric_bound  = "none"
<br>> ML_BSSN::ML_Gamma_bound  = "none"
<br>> ML_BSSN::ML_trace_curv_bound = "none"
<br>> ML_BSSN::ML_curv_bound  = "none"
<br>> ML_BSSN::ML_lapse_bound  = "none"
<br>> ML_BSSN::ML_dtlapse_bound  = "none"
<br>> ML_BSSN::ML_shift_bound  = "none"
<br>> ML_BSSN::ML_dtshift_bound  = "none"
<br>>
<br>> admbase::initial_lapse  = "one"
<br>>
<br>>
<br>> ADMBase::initial_shift  = "zero"
<br>> ADMBase::initial_dtlapse = "zero"
<br>> ADMBase::initial_dtshift = "zero"
<br>>
<br>> #------------------------------------------------------------------------------
<br>> # Output
<br>> #------------------------------------------------------------------------------
<br>> IO::out_dir  = $parfile
<br>> IO::out_fileinfo   = "none"
<br>>
<br>>
<br>>
<br>>
<br>> CarpetIOASCII::out1D_every  = 1
<br>> CarpetIOASCII::out1D_vars  = "admbase::alp
<br>> admbase::metric  admbase::curv admbase::shift "
<br>> CarpetIOASCII::out1D_d=no
<br>> _______________________________________________
<br>> Users mailing list
<br>> <a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a>
<br>> <a href="http://lists.einsteintoolkit.org/mailman/listinfo/users">http://lists.einsteintoolkit.org/mailman/listinfo/users</a>
<br><br></p>
</body>
</html>