| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
|
|
| |
IDAnalyticBH has an "epsilon" parameter to smooth out singularites.
This parameter should also be used in Brill-Lindquist initial data.
Otherwise, moving puncture evolution of these data is not possible
without an explicit smoothing step.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@183 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
EinsteinInitialData?/IDAnalyticBH/src/Kerr.c) has a defect along the
z-axis (that is, the spin axis of the hole). The three-metric (for
instance) is not smooth across this axis, leading to discontinuities
in second spatial derivatives near there.
This appears to be due to overly crude safety measures protecting from
division by zero when calculating this data. Tweaking these measures
for more self-consistency leads to much smoother behaviour. I'm
attaching a simple patch that fixes the problem at least for
metric_type = physical.
[Patch from Bernard Kelly]
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@179 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Yosef's patch:
If initial_data=kerr and metric_type=physical, IDAnalyticBH
writes to psi[x-z] and psi[x-z][x-z] (it tries to set these
GFs to zero) even if no storage is allocated. The problem
is with the code:
IDAnalyticBH_zero_CCTK_REAL_array(npoints, psix);
etc.., which is only executed when the metric is
physical.
I was able to use the Kerr data with this patch, but just
deleting the offending lines should be sufficient.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@178 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
| |
documentation and a proof that this is right
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@175 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@169 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@168 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@167 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
data.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@166 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@165 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@164 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
| |
the testsuite for testcomplex because of a new reduction scheme for complex numbers. WaveToy1DF77 also fails with inconsistency, failing at 3,7 processors with varying files and differences between them. Currently investigating these failurse in detail.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@163 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@162 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
equation to make it manifestly non-negative is enough.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@159 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@158 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@157 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
consistent.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@155 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@154 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@153 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
an axis.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@152 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@151 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
Do not produce nans on the origin.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@150 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@149 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@148 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
modifying or adding files.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@147 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
Fixes PR 1692.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@146 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@144 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@143 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@141 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@140 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@139 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@138 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@137 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@136 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
| |
"impossible value for conformal_storage" error msgs
* also slight tweaks to surrounding code so the compiler knows that
make_conformal_derivs is always set before it's used
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@135 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
for Misner 3-BH data
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@134 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
|
| |
in a level 1 warning to tell users that they're in trouble:
WARNING[L1,P0] (IDAnalyticBH):
Misner_init(): warning: about to call fill_iso() ; at present
this routine seems to leak quasi-infinite amounts of
memory at the rate of O(100 megabytes/second) :( :( :(
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@133 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@132 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
setting up initial data
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@131 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
add const qualifier for eval_bh_psi()
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@130 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
| |
some code cleanups in Misner_points.c
*** I think there is still a quasi-infinite memory leak
*** (c. 100 megabytes/second on a xeon) somewhere around the function
*** fill_iso() in Misner_points.c
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@129 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
ADMBase::metric_type = "static conformal"
StaticConformal::conformal_storage = "factor"
led to the local flag make_conformal_derivs never being initialized.
Alas, later code then tested this flag, and if the random garbage that's
in uninitialized memory happened to be nonzero, the later code would then
proceed to try to initialize the 1st derivatives of the conformal factor.
Since these don't have storage allocated, this core-dumps. :( :(
The bug is fixed by changing the logic to always set make_conformal_derivs .
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@128 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
void Misner_standard(CCTK_ARGUMENTS)
where setting
ADMBase::metric_type = "static conformal"
StaticConformal::conformal_storage = "factor"
led to the local flag make_conformal_derivs never being initialized.
Alas, later code then tested this flag, and if the random garbage that's
in uninitialized memory happened to be nonzero, the later code would then
proceed to try to initialize the 1st derivatives of the conformal factor.
Since these don't have storage allocated, this core-dumps. :( :(
The bug is fixed by changing the logic to always set make_conformal_derivs .
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@127 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
| |
IDAnalyticBH_zero_CCTK_REAL_array() when setting a grid array to all zeros.
This is necessary because memset() gives all bits 0, but we want
floating point (CCTK_REAL) 0.0, which might be a different bit pattern.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@126 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@125 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
Fixed broken references.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@124 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
|
|
|
|
| |
is just
in the physical metric, the conformal (usual) metric is fine!).
Spotted and fixed by Frank Herrmann.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@123 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
Fixes Cactus/1031
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@122 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
|
|
|
| |
Tom
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@121 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@120 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@118 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
|