aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* Use "epsilon" also for Brill-LindquistHEADmastereschnett2011-10-27
| | | | | | | | | | 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
* The Kerr initial data (ineschnett2011-03-08
| | | | | | | | | | | | | | | | | | 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
* Segfault with kerr initial data using IDAnalyticBHbmundim2011-02-28
| | | | | | | | | | | | | | | | | | | | 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
* correct computation of shift vector in Kerr metric, along with the ↵knarf2010-10-07
| | | | | | 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
* Decrease the fudge factor which avoids nans again. Rewriting theschnetter2005-04-14
| | | | | | | 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
* Increase the fudge factor which avoids nans.schnetter2005-04-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@158 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Produce good initial data even on the axis. Patch by Bernard Kelly.schnetter2004-12-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@157 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* assert fixyye002004-08-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@154 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Use a fudge factor to remove singularities when grid points fall ontoschnetter2004-07-06
| | | | | | | an axis. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@152 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Do not produce nans on the origin.schnetter2004-06-23
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@151 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Add lapse "schwarzschild" as an alias for "schwarz".schnetter2004-06-23
| | | | | | | Do not produce nans on the origin. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@150 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Untabified.goodale2004-05-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@148 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Untabified. Please try to adhere to the Cactus coding guidelines whengoodale2004-05-17
| | | | | | | modifying or adding files. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@147 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Bug-fix from Frank Herrman.goodale2004-05-13
| | | | | | | Fixes PR 1692. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@146 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* * fix mismatched CCTK_VWarn() formats and arguments forjthorn2003-06-04
| | | | | | | | | "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
* since I haven't tracked down the infinite memory leak, at least putjthorn2003-05-06
| | | | | | | | | | | | 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
* add CCTK_VInfo() calls to tell the user that we're in this thornjthorn2003-05-06
| | | | | | | setting up initial data git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@131 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* fix grdoc comment for iso(),jthorn2003-05-06
| | | | | | | add const qualifier for eval_bh_psi() git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@130 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* add prototypes for a bunch of un-prototyped functionsjthorn2003-05-06
| | | | | | | | | | | 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
* fix a bug where settingjthorn2003-05-06
| | | | | | | | | | | | | | | 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
* This patch fixes a bug injthorn2003-05-02
| | | | | | | | | | | | | | | | | 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
* Change from using memset() to using a new functionjthorn2003-03-24
| | | | | | | | | 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
* Fix to a bug in the Schwarzschild metric for the physical case (Note: this ↵allen2002-11-01
| | | | | | | | | | | 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
* Adding missing @history tags.goodale2002-08-14
| | | | | | | Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@121 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* tiding and changing routine name to avoid clashesallen2002-04-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@104 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Set conformal_state to 1 if only the conformal factor is being initialised.goodale2002-04-26
| | | | | | | Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@101 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Removing some excess code I accidently left in.goodale2002-04-26
| | | | | | | Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@100 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Just spotted that Schwarzschild does calculate the second derivatives, sogoodale2002-04-26
| | | | | | | | | now set the conformal_state appropriately. Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@99 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Converted to use new metric_type and conformal_storage variables.goodale2002-04-26
| | | | | | | | | | | | The thorns calculate the appropriate number of derivatives and set the conformal_state variable appropriately - note that currently there is no code to do second derivatives of Schwarzschild, so this sets the conformal_state to a maximum of 2 (i.e. factor and first derivatives). Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@98 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Remove compiler warning on originallen2002-01-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@94 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Formatting changesallen2001-10-21
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@92 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Cannot write to parameters since these are read-only now.tradke2001-05-22
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@90 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Fixed CCTK_FILEVERSIONallen2001-05-21
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@89 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* CCTK_FILEVERSIONallen2001-05-11
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@87 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Making rcsid strings constant - PR 669.goodale2001-05-10
| | | | | | | Tom git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@86 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* CCTK_FILEVERSIONallen2001-05-09
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@85 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Fixed wrong argument types to pow(3) which caused the testsuites on DEC Alphastradke2001-05-07
| | | | | | | | (with native C compiler) to fail if the second argument to this function was given as an integer. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@84 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Removing compiler warningsallen2001-02-25
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@82 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Fixing warningsallen2001-02-24
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@81 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Changed C++ comments back to C comments.miguel2000-10-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@75 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* John Bakers initial data for Kerr black holesallen2000-09-15
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@74 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Some optimizations: use single loop with linear indices, replace divs by mults.tradke2000-06-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@70 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* 20% speedup by taking the sinh and tanh calculations out of the loop.tradke2000-06-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@69 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Loop optimizations and divs replaced by mults make calculation of ID fortradke2000-06-22
| | | | | | | two Misner BHs 25% faster (on origin). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@67 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Change CARGUMENTS to ARGUMENTSallen2000-03-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@65 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Accidently added a return code to a void function yesturdayallen2000-03-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@64 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Added return code so that automatic warning removal worksallen2000-03-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@63 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Removing a print Hello statement ...allen2000-03-23
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@62 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* Removing USE_CCTK_ARGUMENTS and USE_CCTK_PARAMETERS since they areallen2000-02-05
| | | | | | | automatically included now git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@60 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
* cctk_parameters -> cctk_Parametersallen2000-01-27
| | | | | | | cctk_arguments -> cctk_Arguments git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@59 6a3ddf76-46e1-4315-99d9-bc56cac1ef84
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-18 02:28:50 +0000