aboutsummaryrefslogtreecommitdiff
path: root/m/McLachlan_BSSN.m
Commit message (Collapse)AuthorAge
* McLachlan_BSSN.m: Correct harmonic shift for conformalMethod = phiBernard Kelly2013-04-15
| | | | | | | The harmonic shift equations implemented in 7f81c9901489f8a94493320a27437bd125b69010 are correct only for the W conformal method, and lead to wrong results if the phi conformal method is used. This commit corrects this.
* McLachlan_BSSN.m: Add comment to explain confusing situation with advection ↵Ian Hinder2013-04-02
| | | | term for \tilde \Gamma^i in \partial_t B^i equation
* Fix the sign for 2nd and 6th order dissipation operators.Peter Diener2013-01-29
| | | | | | Fix the sign for 2nd and 6th order dissipation operators. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
* McLachlan_BSSN: fix the Hamiltonian constraint for CCZ4.Anton Khirnov2012-12-18
| | | | | | A typo in a8794cd2227c87c7457ba51d040e863be2f981e7 has caused the additional CCZ4 terms to overwrite the BSSN value instead of adding to it.
* Change the dafault value of apply_dissipation to "never".Peter Diener2012-10-23
| | | | | | | | Change the default value of apply_dissipation from "always" to "never". Since the default dissipation strength is zero this should not change results for users using default values but should save time. Signed-off-by: Peter Diener <diener@cct.lsu.edu>
* Correct some errors in CCZ4, in particular the GammaShift parameterErik Schnetter2012-05-29
|
* Do not generate *_O* and _MP thornsErik Schnetter2012-05-02
|
* Rename CMBSSNW and CMBSSNphi to CMW and CMphi.Barry Wardell2012-05-01
| | | | They can be used with both BSSN and CCZ4 formulations.
* Generate a separate thorn for CCZ4.Barry Wardell2012-05-01
| | | | This still uses the same Kranc script as the standard BSSN code, but having a separate thorn means that the BSSN is not negatively affected performance-wise.
* Add formulation parameter to select between BSSN and CCZ4.Barry Wardell2012-04-30
| | | | The default selects BSSN for backwards compatibility.
* Some superficial changes:Erik Schnetter2012-04-27
| | | | | | * Rename Tet -> Theta and Zet -> Z * Introduce an enum to select between CMBSSNphi, CMBSSNW and CMCCZ4. * Rename gamashift parameter to GammaShift
* Merge branch 'master' into CCZ4Barry Wardell2012-04-27
|\ | | | | | | | | Conflicts: m/McLachlan_BSSN.m
| * Correct etaExpr to avoid (inconsequential) division by zeroErik Schnetter2012-02-05
| |
| * Apply advection terms to A or B^i only if they are evolvedErik Schnetter2012-01-04
| |
| * Correct definition of dot[A] advection termsErik Schnetter2012-01-02
| | | | | | | | | | There terms were implemented wrong, probably when the advection terms were split out of the main RHS routine. Without this correction, Kerr-Schild is not manifestly stationary any more. Some of these terms are in code paths that are (to my knowledge) unused, and so it is not clear what to do. I've tried to make these code paths as similar to the others as possible.
| * Remove unused function definition SetGroupTimelevelsErik Schnetter2012-01-02
| |
| * Don't replace Pi by its numeric equivalent -- Kranc does this alreadyErik Schnetter2012-01-02
| |
* | Merge branch 'master' into CCZ4Barry Wardell2012-01-01
|\| | | | | | | | | Conflicts: m/McLachlan_BSSN.m
| * Correct description of LapseAdvectionCoeffErik Schnetter2011-12-23
| |
| * Implement harmonic shift conditionErik Schnetter2011-12-15
| |
| * Fix inconsistent index positionsIan Hinder2011-11-20
| | | | | | No change to generated code
| * Initialise boundaries of A and B^i together with boundaries of Gamma^iErik Schnetter2011-10-02
| |
| * Use kranc script instead of calling Mathematica directlyIan Hinder2011-09-24
| | | | | | This allows use of better error handling provided by Kranc
| * Remove debugging codeIan Hinder2011-09-23
| |
| * McLachlan_BSSN.m: Make finite difference order a run-time parameterIan Hinder2011-09-19
| |
| * Enable Vectorisation.Barry Wardell2011-08-18
| |
| * Correct schedule, in particular for checkpoint/recoveryErik Schnetter2011-08-02
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Do not mark ADMBase variables for non-checkpointing if they have multiple timelevels. (Variables with multiple timelevels must always be checkpointed, because the past timelevels cannot be regenerated after recovery.) Finally remove all perl post-processing of the auto-generated code; instead, use proper Kranc mechanisms. Schedule the ADM constraints and ADM quantities after MoL_PostStep, since this is where the ADMBase variables are set. Schedule enforcing the BSSN constraints in the new schedule group MoL_PostStepModify, since they should not be enforced after recovery. (This would lead to inconsistencies at floating-point round-off level.) Regenerate all thorns.
* | Implement CCZ4 formulation.Daniela Alic2012-01-01
| | | | | | | | | | The conformal and covariant formulation of the Z4 system with constraint-violation damping is described in detail in arXiv:1106.2254.
* | Enable Vectorisation.Barry Wardell2012-01-01
|/
* Correct scheduling orderErik Schnetter2011-06-23
|
* Use new Jacobian mechanism for all thornsErik Schnetter2011-06-23
|
* Set RHS variables to zero in ANALYSIS before computing themIan Hinder2011-06-13
| | | | The RHS variables can only be computed on the interior in ANALYSIS, so initialise the entire grid to zero before computing the variable. An alternative would be to create a new looping method "Exterior" which looped over all the non-interior points. That would avoid setting the points twice, though setting to zero is inexpensive so this is not really necessary.
* Undo "Revert recent commits"Ian Hinder2011-06-13
| | | | | | It was only the tests which were wrong. This reverts commit 05347a08d0c9bd2a87846ab4ad8990fe26274a4a.
* Revert recent commitsIan Hinder2011-06-12
| | | | | | | | These are causing NaNs when run with poison. I don't know if this is due to the tests or the code, so I am reverting the commits from 863a3e5b25e7150148f9d2b60b4b362628c675f7 to 2725eb1eb32525486df76a3686f8e550155c8e0c while the problem is being diagnosed.
* McLachlan_BSSN.m: Add InitGamma calculationIan Hinder2011-06-11
| | | | This is necessary to temporarily work around a problem with multipatch (points which are both interpatch and outer boundary are not initialised in the correct order)
* McLachlan_BSSN.m: Change from old TensorTools to new Kranc Jacobians for MPIan Hinder2011-06-11
| | | | | | | The Jacobian gridfunction groups are no longer hard-coded into McLachlan, and the Coordinates thorn is no longer inherited. Instead, the user must use GenericFD parameters to specify the Jacobian to use. The new Kranc Jacobian implementation allows a single map to be specified on which the Jacobian is not applied, which greatly improves performance in typical multipatch simulations. To use the new Jacobian method with Llama, set GenericFD::jacobian_group = "Coordinates::jacobian" GenericFD::jacobian_derivative_group = "Coordinates::jacobian2" GenericFD::jacobian_identity_map = 0
* McLachlan_BSSN.m: Split advection terms into a separate calculationIan Hinder2011-06-09
| | | This helps to fit the RHS loops in the Intel instruction cache and benefits performance.
* McLachlan_BSSN.m: Split dissipation into a separate calculationIan Hinder2011-06-09
| | | Now controlled by parameter "apply_dissipation = always/never"
* McLachlan_BSSN.m: Apply dissipation consistentlyIan Hinder2011-06-09
| | | Previously, the dissipation term for Xt^i was being added both in the Xt^i RHS and in the B^i RHS, due to the way the temporary variable for dot[X^i} was used in the Gamma driver shift.
* Add comment about constraint projectionErik Schnetter2011-01-03
|
* Split the BSSN constraint calculations differently into two calculations to ↵Erik Schnetter2011-01-03
| | | | reduce code size
* 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 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.
* Change AssertSymmetricDecreasing to AssertSymmetricIncreasingErik Schnetter2010-10-11
|
* Update path to KrancErik Schnetter2010-08-04
|
* 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>
* Add missing definition of temporary variableErik Schnetter2010-05-03
|
* 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.
* 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-31 23:56:07 +0000