aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* Don't define the (unused) upwinding differencing operators in the ADM ↵Erik Schnetter2011-01-03
| | | | constraints
* Calculate the Ricci scalar for the BSSN constrains in the same way as for ↵Erik Schnetter2011-01-03
| | | | the time evolution; this simplifies the generated code significantly.
* Do not schedule the conversion to ADMBase at post_recover_variables if they ↵Erik Schnetter2011-01-03
| | | | are scheduled in MoL_PostStep, since MoL_PostStep is already scheduled in post_recover_variables.
* Regenerate codeErik Schnetter2010-12-08
|
* Correct boundary conditions for ADMConstraints and ADMQuantitiesErik Schnetter2010-12-07
|
* Do not schedule the RHS calculation at postinitial any moreErik Schnetter2010-12-07
|
* Do not schedule a group for the constraint boundary conditions any moreErik Schnetter2010-12-07
|
* Do not fix up the boundaries of analysis quantities any more; Kranc now ↵Erik Schnetter2010-12-07
| | | | handles this correctly.
* ML_BSSN_Test: Update testsuites to use 4th order differencing in ExactIan Hinder2010-12-07
| | | | | | | | | | | | | | | | | | | | | | | | | The Exact thorn uses finite differencing to compute the extrinsic curvature rather than using the analytic result. The step size used is independent of the Cactus grid. Using a very small step size for accuracy can lead to the subtraction of floating point numbers which are very close to each other, leading to a large loss of precision. This loss of precision can amplify differences between compilers and optimisation settings. A recent Exact commit implements 4th order differencing and a configurable step size. Hence for a given finite differencing error, a larger step size can be used, reducing the loss of precision. This commit to the McLachlan testsuites (parameter files and data) implements this change by enabling 4th order differencing and a step size of 2e-4. A convergence test in the step size has shown that this leads to a finite differencing error of ~1e-14. It also leads to more consistent results across compilers and optimisation settings. The data was generated with gcc using -O2 and the testsuites pass with the Intel compiler also using -O2. This testsuite data also takes into account the modification to the constraint boundary condition in the previous commit.
* Do not apply constraint boundaries manually; Kranc does this now automatically.Erik Schnetter2010-12-06
|
* ML_BSSN: Simplify handling of upwinding derivativesErik Schnetter2010-12-06
| | | | | | | | | Use a new function Upwind to calculate upwind derivatives. Introduce a new (Mathematica) parameter that decides whether to implement upwind derivatives the old way (using IfThen depending on the direction) or the new way (using Abs and the symmetric/antisymmetric parts of the stencils). Improve how to create partial calculations automatically. Evaluate the constraints in two partial calculations.
* par: Correct output variables in example parameter filesErik Schnetter2010-12-06
|
* Regenerate code after setting the default tensor weight in Kranc to 0Erik Schnetter2010-11-23
|
* ML_BSSN_Test: Remove testsuite output for constraintsIan Hinder2010-11-01
| | | The current Cactus testsuite mechanism does not allow you to set per-variable absolute and relative tolerances. The constraints are typically very small in comparison to the evolved variables used to compute them. This means that roundoff errors in the evolved variables can lead to large relative errors in the constraints. The correct solution would be to specify only an absolute tolerance for the constraints, but this is not possible at the moment. In future, a separate test could be written to test that the constraints are being correctly computed on constraint-violating initial data. For now, we omit the tests of the constraints from the testsuites.
* Regenerate codeErik Schnetter2010-10-11
|
* Schedule ADMBase variable calculations in MoL_PseudoEvolutionErik Schnetter2010-10-11
|
* Use boundary_width instead of 1 when selecting boundary conditionsErik Schnetter2010-10-11
|
* Change AssertSymmetricDecreasing to AssertSymmetricIncreasingErik Schnetter2010-10-11
|
* Update path to KrancErik Schnetter2010-08-04
|
* Merge branches 'master', 'master', 'master', 'master', 'master', 'master', ↵Erik Schnetter2010-06-18
|\ | | | | | | 'master', 'master', 'master', 'master', 'master', 'master', 'master', 'master' and 'master' of carpetcode.dyndns.org:McLachlan
| * Add test cases for BSSN - 3D diagonal gauge waveIan Hinder2010-06-16
| | | | | | | | | | | | | | For some reason, this test only passes when run with OMP_NUM_THREADS=1. If you use OMP_NUM_THREADS=2, or leave it unset on a multi-core machine, there are differences above tolerance in the Hamiltonian constraint, though the evolved variables are OK.
* | Merge branches 'master' and 'master' of carpetcode.dyndns.org:McLachlanErik Schnetter2010-06-01
|/
* Reintroduce an optimization.Peter Diener2010-05-05
| | | | | | Reintroduce an optimization in the evolution calculation, that was previously only introduced in the hand split optimised calculation. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
* Correct error in Makefile.Erik Schnetter2010-05-03
| | | | Regenerate code.
* Regenerate codeErik Schnetter2010-05-03
|
* Regenerate codeErik Schnetter2010-05-03
|
* Add missing definition of temporary variableErik Schnetter2010-05-03
|
* Add benchmark parameter fileErik Schnetter2010-05-02
|
* Update example parameter filesErik Schnetter2010-05-02
|
* Regenerate codeErik Schnetter2010-05-02
|
* Add dissipation terms to the RHS.Erik Schnetter2010-05-02
| | | | | | Create a new function PartialCalculation, which splits a calculation, retaining only those parts evaluating certain variables. This is useful to more easily split calculations. Modify the way upwind derivatives are calculated. Instead of using integer dir^i selectors, split them into an antisymmetric and a symmetric part, which are multiplied by the shift and its absolute value, respectively.
* Don't use explicit path name when including Slicing.hErik Schnetter2010-04-21
|
* Merge branch 'master' of /Users/eschnett/Cbeta/arrangements/McLachlanErik Schnetter2010-04-07
|\ | | | | | | | | | | | | | | | | Conflicts: ML_BSSN/src/ML_BSSN_RHS1.c ML_BSSN_MP/src/ML_BSSN_MP_RHS1.c ML_BSSN_O2/src/ML_BSSN_O2_RHS1.c ML_BSSN_O8/src/ML_BSSN_O8_RHS1.c m/McLachlan_BSSN.m
| * Move RHS evaluations around to optimize for 32Kb instruction cache machines.Peter Diener2010-03-30
| | | | | | | | | | | | | | | | | | | | | | | | Move RHS evaluations for trK and the lapse to RHS1 to better even out the instruction count between RHS1 and RHS2 leading to better performance on 32Kb instruction cache machines. Also replace covariant derivatives with respect to the physical metric with covariant derivatives with respect to the conformal metric. This avoids having to compute the physical Christoffel symbols. The performance on 64Kb instruction cache machines is unchanged. Also get rid of some unused shorthands and generate code for 3 timelevels. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
| * Fix and clarify parameter conversions.Peter Diener2010-03-29
| | | | | | | | | | | | Fix and (hopefully) clarify parameter conversions. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
* | Add SpatialShiftGammaCoeffRadius to reduce ShiftGammaCoeff outside of a ↵Erik Schnetter2010-04-07
|/ | | | certain radius
* Don't build or schedule old RHS routine any moreErik Schnetter2010-03-21
|
* Use Ian's version of Kranc, which is found atErik Schnetter2010-03-18
| | | | <git@github.com:ianhinder/Kranc.git>.
* Correct boundary conditionsErik Schnetter2010-03-14
|
* Schedule boundary conditions in ML_PseudoEvolutionErik Schnetter2010-03-14
|
* Generate and add 8th order version as wellErik Schnetter2010-03-12
|
* Make eta (the Gamma driver damping parameter) a local variable only, not a ↵Erik Schnetter2010-03-09
| | | | grid function
* Add second order accurate versions of all thornsErik Schnetter2010-02-19
|
* Use new schedule group MoL_PseudoEvolution for the ADM constraints and the ↵Erik Schnetter2010-02-19
| | | | | | ADM quantities. Re-generate all code.
* Regenerate code after correcting error in Kranc.Erik Schnetter2010-02-13
|
* Remove useMatter parameters which are now unused.Erik Schnetter2010-02-13
|
* Treate Tmunu variables as tensors, instead of accessing them via local ↵Erik Schnetter2010-02-13
| | | | helper scalars.
* MergeErik Schnetter2010-02-13
|
* Merge /Users/eschnett/Cbeta/arrangements/McLachlanErik Schnetter2010-02-13
|\ | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Conflicts: ML_BSSN/src/ML_BSSN_Minkowski.c ML_BSSN/src/ML_BSSN_convertFromADMBase.c ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.c ML_BSSN/src/ML_BSSN_setBetaDriverSpatial.c ML_BSSN/src/make.code.defn ML_BSSN_Helper/src/SetGroupTags.c ML_BSSN_M/param.ccl ML_BSSN_M/schedule.ccl ML_BSSN_M/src/ML_BSSN_M_Minkowski.c ML_BSSN_M/src/ML_BSSN_M_convertFromADMBase.c ML_BSSN_M/src/ML_BSSN_M_convertFromADMBaseGamma.c ML_BSSN_M/src/ML_BSSN_M_setBetaDriverSpatial.c ML_BSSN_M/src/make.code.defn ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c ML_BSSN_MP/src/ML_BSSN_MP_setBetaDriverSpatial.c ML_BSSN_MP/src/make.code.defn ML_BSSN_MP_Helper/src/SetGroupTags.c ML_BSSN_M_Helper/schedule.ccl ML_BSSN_M_Helper/src/SetGroupTags.c m/McLachlan_BSSN.m m/prototype/ML_BSSN_Helper/src/SetGroupTags.c
| * Schedule ML_BSSN_convertToADMBaseGroupWrapper also at ↵Frank Loeffler2010-02-08
| | | | | | | | | | | | | | | | CCTK_POST_RECOVER_VARIABLES -> recompute ADM variables after recovery before first timestep Signed-off-by: Frank Loeffler <knarf@cct.lsu.edu>
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-15 01:44:28 +0000