aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* remove unused (and duplicated) lables, remove unused variables, always ↵knarf2012-12-17
| | | | | | initialize 'rate' git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@122 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* remove catch of wrong function-call because 'guessed' solution passes ↵knarf2010-06-07
| | | | | | testsuites if tolerance is lowered to similar levels as other IDAxi-thorns git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@111 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* correct include file paths (the crude way(knarf2010-05-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@109 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Pass double precision instead of single precision valueschnetter2010-03-16
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@108 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Comment out call to function that has wrong arguments and thatschnetter2010-03-16
| | | | | | | (presumably) is not called anyway git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@107 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Use CCTK_REAL instead of real*8.schnetter2006-05-27
| | | | | | | | | | | | | | | Avoid coordinate singularities on the axis by adding a small epsilon. Correct a typo in a comment. Do not abort if the interpolator fails. Instead, initialise all 3D variables to Minkowski. Paper over the coordinate singularity at the origin by replacing it with Minkowski data. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@94 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Use CCTK_REAL instead of real*8schnetter2006-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@93 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Removing single quotes from comments as preprocessors don't like them.goodale2005-10-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@83 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Add an option to control whether we generate a conformal metricjthorn2005-07-08
| | | | | | | | | | (StaticConformal::conformal_state = 3) or a physical metric ( = 0), as per http://www.cactuscode.org/pipermail/developers/2005-July/001178.html and subsequent discussion git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@79 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Fortran apparently requires an explicit field width for I formatjthorn2005-05-31
| | | | | | | --> s/ write (A,I) ... / write (A,I8) ... / git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@72 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* s/1.5/1.5d0/ and similarly for a few other single-precision constantsjthorn2005-05-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@71 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Use "character(max_string_length)" instead ofschnetter2005-05-20
| | | | | | | | | "character*max_string_length", which is not legal Fortran. (The syntax "character*(max_string_length)" would also be legal, but is deprecated.) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@70 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* This is a general code cleanup, and also adds some new features tojthorn2005-05-20
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | this thorn. [Unforunately, for the reasons explained in today's E-mail discussion on the developers@cactuscode.org mailing list, it's not practical for me to break this up into multiple smaller single-purposed commits.] param.ccl * add comments * add additional interpolator parameters (described in detail below) * add "output psi on 2D grid" parameters (described in detail below) * add debugging parameters doc/documentation.tex * add "IDAxiBrillBH/" prefix to make LaTeX labels unique across thorns * correct a few typos * clarify explanation of how we solve on a 2-D grid and then interpolate to the Cactus 3-D grid * explain the new interpolator and output-psi-on-2D-grid parameters * add some comments on choosing the 2-D grid resolution parameters * mention debug parameters doc/TODO * new file noting some problems with this thorn which I found, but didn't fix src/IDAxiBrillBH.F * use new CCTK_WARN(CCTK_WARN_ABORT, ...) instead of old CCTK_WARN(0, ...) * more flexible interpolator parameters: This thorn first solves the Brill-wave equation on an internal (axisymmetric) 2D grid, then uses uses the standard Cactus CCTK_InterpLocalUniform() local-interpolator API to interpolate this to the 3D grid points. Before this patch, this thorn hard-wired the interpolation operator to "uniform cartesian", only allowed orders 1, 2, and 3, and didn't allow any other interpolator parameters to be set. This patch adds two new parameters to allow any interpolation operator and parameters to be specified. Either the old or the new parameters can be used (see doc/documentation.tex for details of how this works); the defaults leave the behavior of this thorn unchanged. * add option to output psi on 2D grid When using this thorn, I found it hard to choose the parameters defining the resolution and extent of the internal 2D grid. To help with this, I added an option to output the solution on that grid to an ASCII data file, so it can be examined and plotted. * add debugging code I've added debug options which, if enabled, print the values of various internal quantities at specified 2D and/or 3D grid points. I've left these in the final "production" code as possibly being useful for future debugging. * add many comments * correct or remove some old comments which were out of date * adjust whitespace to make the code a bit more readable: I've changed code like allocate(cc(neb,nqb),ce(neb,nqb),cw(neb,nqb),cn(neb,nqb),cs(neb,nqb), $ rhs(neb,nqb),psi2d(neb,nqb),detapsi2d(neb,nqb),dqpsi2d(neb,nqb), $ detaetapsi2d(neb,nqb),detaqpsi2d(neb,nqb),dqqpsi2d(neb,nqb), $ etagrd(neb),qgrd(nqb)) to make it clearer which arrays have the same size: allocate( cc(neb,nqb)) allocate( ce(neb,nqb)) allocate( cw(neb,nqb)) allocate( cn(neb,nqb)) allocate( cs(neb,nqb)) allocate( rhs(neb,nqb)) allocate( psi2d(neb,nqb)) allocate( detapsi2d(neb,nqb)) allocate( dqpsi2d(neb,nqb)) allocate(detaetapsi2d(neb,nqb)) allocate( detaqpsi2d(neb,nqb)) allocate( dqqpsi2d(neb,nqb)) allocate(etagrd(neb)) allocate( qgrd(nqb)) * change some Fortran write statements to CCTK_INFO [so they're properly synchronized with output from C routines, even when stdout+stderr are redirected to a log file] * include more information about what went wrong, in various error messages [eg if the interpolator returns an error code, the code now includes that error code in the error message] * rename a few variables to make their purpose clearer: ep1 --> error_at_this_grid_point ep2 --> max_error_in_grid git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@69 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* This file is unused; move it from src/ to new directory archive/ .jthorn2005-05-20
| | | | | | | | | [I chose to keep it in archive/, rather than only in the CVS attic, because the CVS attic isn't usable with normal Unix tools ("cat, awk, and grep" et al).] git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@67 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Destroy the interpolation parameter table again when all the interpolationtradke2004-11-05
| | | | | | | calls have finished. Fixes a (small) resource leak. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@66 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Use CCTK_InterpLocalUniform() instead of the deprecated CCTK_InterpLocal() API.tradke2004-05-17
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@62 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Slight formatting changes to make the red hat C preprocessor work withoutdiener2003-11-05
| | | | | | | any special options. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@56 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Started formatting info statementsallen2002-05-31
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@48 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* changing ifdefs and include filesallen2002-05-16
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@46 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* nclude files need to be from new arangement temporarilyallen2002-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@43 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Eliminate #include Einstein.hrideout2002-05-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@40 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* IDAxiBrillBH converted to CactusEinstein2rideout2002-05-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@39 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* FARGUMENTS -> ARGUMENTSallen2001-06-25
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@31 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Use the new interpolation flesh API routines and call the new interpolationtradke2001-02-15
| | | | | | | | operators from PUGHInterp (according to the setting of the IDAxiBrillBH::interpolation_order parameter). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@27 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* GetInterpHandle->InterpHandleallen2000-12-04
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@20 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Don't change parameters in place but assign their values to local variablestradke2000-11-21
| | | | | | | | instead. This fixes problems with recovery where the modified parameter values caused different results. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@18 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* GetInterpHandle -> InterpHandleallen2000-07-13
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@15 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Removing my last commitallen2000-04-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@12 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Removing warningsallen2000-04-19
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@11 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* cctk_parameters -> cctk_Parametersallen2000-01-27
| | | | | | | cctk_arguments -> cctk_Arguments git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@9 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* Fix from Thomas Radke for floating point exception on the Alphas causedallen2000-01-24
| | | | | | | | | | | | | | | | | | | | by exponention by a real number .... here's Thomas's description: The include file generated by Mathematica contains an exponentiation operation with a REAL exponent. My Fortran book says that y**x = exp(x * ln(y)) On O2K this intrinsic function checks the exponent against 0, on Alphas it doesn't and you get a FP exception if y <= 0. I just changed the type of the exponent to be INTEGER instead of REAL which was probably also the intention of the programmer. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@8 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* took out old include references to Developmentlanfer1999-11-01
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@5 0a4070d5-58f5-498f-b6c0-2693e757fa0f
* This commit was generated by cvs2svn to compensate for changes in r2, whichallen1999-11-01
included commits to RCS files with non-trunk default branches. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@3 0a4070d5-58f5-498f-b6c0-2693e757fa0f
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-21 17:19:23 +0000