| Commit message (Collapse) | Author | Age |
|
|
|
| |
They can be used with both BSSN and CCZ4 formulations.
|
| |
|
|
|
|
| |
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.
|
| |
|
|
|
|
| |
The default selects BSSN for backwards compatibility.
|
|
|
|
|
|
| |
* Rename Tet -> Theta and Zet -> Z
* Introduce an enum to select between CMBSSNphi, CMBSSNW and CMCCZ4.
* Rename gamashift parameter to GammaShift
|
|\
| |
| |
| |
| | |
Conflicts:
m/McLachlan_BSSN.m
|
| | |
|
| | |
|
| |
| |
| |
| |
| | |
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.
|
| | |
|
| | |
|
|\|
| |
| |
| |
| | |
Conflicts:
m/McLachlan_BSSN.m
|
| | |
|
| | |
|
| |
| |
| | |
No change to generated code
|
| | |
|
| | |
|
| | |
|
| |
| |
| | |
This allows use of better error handling provided by Kranc
|
| | |
|
| | |
|
| |
| |
| |
| |
| |
| |
| | |
Don't use "ADM*" as glob, since this refers to ADM, ADMConstraints,
and ADMQuantities.
Add ADMQuantities to cleanup rule.
|
| | |
|
| | |
|
| | |
|
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| |
| | |
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.
|
| |
| |
| |
| |
| | |
The conformal and covariant formulation of the Z4 system with constraint-violation
damping is described in detail in arXiv:1106.2254.
|
|/ |
|
| |
|
| |
|
| |
|
|
|
|
|
|
| |
Tensor indices with numbers are not supported any more in Kranc;
rename them to use letters instead.
(originally by Erik Schnetter)
|
|
|
|
| |
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.
|
|
|
|
|
|
| |
It was only the tests which were wrong.
This reverts commit 05347a08d0c9bd2a87846ab4ad8990fe26274a4a.
|
|
|
|
| |
This reverts commit a8b69dc06e565427bc96d16c443bb4b3392c4365.
|
|
|
|
|
| |
Tensor indices with numbers are not supported any more in Kranc;
rename them to use letters instead.
|
|
|
|
|
|
|
|
| |
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.
|
|
|
|
| |
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)
|
|
|
|
|
|
|
| |
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
|
|
|
| |
This helps to fit the RHS loops in the Intel instruction cache and benefits performance.
|
|
|
| |
Now controlled by parameter "apply_dissipation = always/never"
|
|
|
| |
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.
|
| |
|
| |
|
|
|
|
| |
reduce code size
|
|
|
|
| |
constraints
|
|
|
|
| |
the time evolution; this simplifies the generated code significantly.
|
|
|
|
| |
are scheduled in MoL_PostStep, since MoL_PostStep is already scheduled in post_recover_variables.
|
| |
|