aboutsummaryrefslogtreecommitdiff
path: root/src/gr/expansion.cc
Commit message (Collapse)AuthorAge
* Make "Erik" branch the default trunk for AHFinderDirect.svnadmin2008-06-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1526 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix a C++ namespaces bug found by Erik:jthorn2005-02-01
| | | | | | | | | | | | | | | error_exit() was declared extern "C" --> it wasn't namespace-mangled --> it could conflict with a similar error_exit() in another thorn although this fix touches a *lot* of files, it's conceptually simple: * error_exit() is not extern "C" any more * everything in src/jtutil/ (which used to be in the jtutil:: namespace) is now in the AHFinderDirect::jtutil:: namespace --> all C++ code in this thorn is now inside the AHFinderDirect:: namespace git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1395 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * add interface to Erik's SphericalSurface thornjthorn2004-05-12
| | | | | | | | | | | | | | ==> With this commit, AHFinderDirect now inherits from AEIThorns/SphericalSurface, so you must have that thorn in your configuration to be able to compile. * add computation of surface quadrupole moments and areal radius * expand BH_diagnostics file format to accomodate quadrupole moments and areal radius, and also to include not-implemented-yet columns for 9 more diagnostics which Erik has implemented in his branch * some other small cleanups git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1329 f88db872-0e4f-0410-b76b-b9085cfa78c5
* remove unneededjthorn2004-05-06
| | | | | | | #include "cctk_Arguments.h" git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1322 f88db872-0e4f-0410-b76b-b9085cfa78c5
* move lots and lots of code into C++ namespace AHFinderDirectjthorn2003-08-18
| | | | | | | so it can't conflict with stuff in other thorns git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1184 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix a typo in checking the geometry interpolator statusjthorn2003-07-28
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1152 f88db872-0e4f-0410-b76b-b9085cfa78c5
* change enum geometry_methodjthorn2003-07-27
| | | | | | | to bool hardwire_Schwarzschild_geometry git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1147 f88db872-0e4f-0410-b76b-b9085cfa78c5
* drop support for (optionally) using CCTK_InterpLocalUniform() instead ofjthorn2003-07-27
| | | | | | | | | | CCTK_InterpGridArrays() to interpolate the Cactus geometry grid functions to the horizon position -- this was never used in practice, and it conflicts with the changse I'm about to commit to make AHFinderDirect work properly under Carpet git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1146 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * add support for CCTK_InterpGridArrays() returning per-pointjthorn2003-07-24
| | | | | | | | | | | | | | interpolator status (use it if it's available, otherwise fall back to CCTK_InterpGridArrays() function result) * drop support for interpolators returning status info on which side of the grid an outside-the-grid point is outside -- this turned out to be awkward to implement in combinatino with per-point status, I only used it to produce slightly more informative error messages in a few cases (ditto in ../patch/) * convert some commented-out debugging code to #ifdef on various DEBUG symbols git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1138 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add new parameter and code to allow computing surfaces of constant expansionjthorn2003-07-04
| | | | | | | | (the case expansion = 0 is the AH, others should be surfaces nested inside or outside the AH) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1118 f88db872-0e4f-0410-b76b-b9085cfa78c5
* set parameter "supress_warnings" in parameter table before callingjthorn2003-06-04
| | | | | | | | CCTK_InterpGridArrays() so it doesn't print warning messages about points outside of the grid git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1097 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix error-handling code for when one processor gets ↵jthorn2003-05-14
| | | | | | CCTK_INTERP_ERROR_POINT_OUTSIDE git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1058 f88db872-0e4f-0410-b76b-b9085cfa78c5
* if we get a CCTK_ERROR_INTERP_POINT_OUTSIDE error code (point outside grid),jthorn2003-05-13
| | | | | | | | then make sure the local interpolator also returned this error code before we use its error status git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1057 f88db872-0e4f-0410-b76b-b9085cfa78c5
* redo excision mask:jthorn2003-05-13
| | | | | | | | | | | * there are now 3 states (inside/buffer/outside) * there is now a "noshrink" option * many other changes too -- see thorn guide for details the new excision support isn't very well tested yet :( :( git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1056 f88db872-0e4f-0410-b76b-b9085cfa78c5
* change error handling if geometry interpolatorjthorn2003-05-07
| | | | | | | | | | | | | | returns CCTK_ERROR_INTERP_POINT_OUTSIDE: if the interpolator sets the error_{pt,axis,direction} variables in the parameter table, fine --> use them in our message but if not, don't abort, rather just produce a slightly less informative msg also, only bother doing all this if we're going to print a msg -- if not, skip testing the error_* variables entirely git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1052 f88db872-0e4f-0410-b76b-b9085cfa78c5
* if we get the error code CCTK_ERROR_INTERP_POINT_OUTSIDE back fromjthorn2003-05-06
| | | | | | | | | the geometry interpolator, add text to the error message to say that in a uniprocessor run, this can also be caused by too small a driver::ghost_size git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1051 f88db872-0e4f-0410-b76b-b9085cfa78c5
* clarify error messages now that "point outside grid"jthorn2003-03-23
| | | | | | | and "ghost size too small" are distinct interpolator error codes git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1005 f88db872-0e4f-0410-b76b-b9085cfa78c5
* catch "ghost size too small" error code from CCTK_InterpGridArrays()jthorn2003-03-22
| | | | | | | and print a suitable error message git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1002 f88db872-0e4f-0410-b76b-b9085cfa78c5
* clarify local vs global Cactus coordinate origin,jthorn2003-03-17
| | | | | | | | add member functions to struct cactus_grid_info to do coordinate transformations git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@980 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add 1st draft of support for setting an excision maskjthorn2003-03-17
| | | | | | | -- alas this doesn't work properly (yet) for multiple processors git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@975 f88db872-0e4f-0410-b76b-b9085cfa78c5
* print a bit more info on point-outside-grid errorsjthorn2003-03-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@970 f88db872-0e4f-0410-b76b-b9085cfa78c5
* expansion() and expansion_Jacobian() now return an enum expansion_statusjthorn2003-03-02
| | | | | | | | | instead of just a Boolean success/fail flag --> this will let higher-level code know why a failing evaluation failed, so it can potentially print a message to tell the user that (eg) the ghost size is too small git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@959 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * general tidying of multiprocessor codejthorn2003-02-16
| | | | | | | | | | | | | | * add current iteration number to what's sent to processor #0 for printing at each iteration * add a new verbosity level "algorithm debug" for what used to be #ifdef DEBUG debugging code * add new parameters & logic to set CCTK_Warn() levels for various warnings (eg AH is outside the gid, metric isn't +ve definite, etc) * restore the ability to output the horizon shape to a data file, which I mistakenly deleted a few days ago :( git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@947 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * major changes to driver routines to find multiple horizons in paralleljthorn2003-02-15
| | | | | | | | | | | | across multiple processors -- see src/driver/README.parallel for details * drop convergence checks on ||Delta_h|| in param.ccl because they don't fit well with parallelization changes ==> With this changes, AHFinderDirect is now (I think) multiprocessor-ready!! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@946 f88db872-0e4f-0410-b76b-b9085cfa78c5
* update interpolator error return codesjthorn2003-02-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@944 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * print a bit more info about setting up data structuresjthorn2003-01-30
| | | | | | | | * fix a few bugs where error messages had mismatched printf (CCTK_VWarn) formats and arguments git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@941 f88db872-0e4f-0410-b76b-b9085cfa78c5
* Changed the order of function parameters in CCTK_InterpGridArrays().tradke2003-01-22
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@932 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * change to not use STL vector class any morejthorn2003-01-16
| | | | | | | | | | | (it caused portability problems on platinum) ==> move stl_vector.hh from src/include/ to archive/ directory * convert patch_system_type and initial_guess_method parameters to be array parameters, i.e. they're now set individually for each AH we want to find git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@920 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add beginnings of support for using new CCTK_InterpGridArrays()jthorn2003-01-13
| | | | | | | | | | | global interpolator & multiprocessor operation ... at present the CCTK_InterpGridArrays() support works fine, but only on 1 processor, I seem to get an infinite loop on multiple processors git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@917 f88db872-0e4f-0410-b76b-b9085cfa78c5
* change in terminology/notation:jthorn2002-11-21
LHS of apparent horizon equation was formerly called $H$, changed to $\Theta$ because it is in fact precisely the expansion of the surface r = h(angle). This implies renaming a lot of variables & functions & a few parameters (which aren't specified in most par files outside my own for testing this thorn) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@901 f88db872-0e4f-0410-b76b-b9085cfa78c5