aboutsummaryrefslogtreecommitdiff
path: root/src/gr
Commit message (Collapse)AuthorAge
* Use CCTK_VError in some placeseschnett2013-02-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1574 f88db872-0e4f-0410-b76b-b9085cfa78c5
* 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
* disable automagic template instantiation on SGI CC compilerjthorn2004-07-01
| | | | | | | | | | | | because it's totally broken: if file A contains #include "../jtutil/array.hh" the compiler tries to also #include "../jtutil/array.cc" (even though I didn't ask for that file)... then it gets fatal errors when it can't resolve relative include paths for *that* file :( git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1375 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
* Add proper CXXFLAGS to disable automatic template instantiation on DEC Alphastradke2003-12-12
| | | | | | | with native C++. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1226 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
* we don't want a (non--no-op) default target in this makefilejthorn2003-07-31
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1169 f88db872-0e4f-0410-b76b-b9085cfa78c5
* rename makefile to Makefile.standalonejthorn2003-07-31
| | | | | | | | since the latter is the Cactus standard and doesn't interfere with the Cactus build systems on Windoze systems git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1168 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
* move mask processing stuff (about the Cactus grid)jthorn2003-07-27
| | | | | | | | | | out of gr / struct cactus_grid_info into new driver / struct mask_grid_info which is private to driver/mask.cc --> also revise this info a bit to play better with Carpet git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1150 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 max_allowable_horizon_radiusjthorn2003-07-10
| | | | | | | | parameter to (we think) fix a serious bug where if one horizon goes off the edge of the grid, we fail to find any other horizons git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1130 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
* * change comments in all make.code.defn filesjthorn2003-06-04
| | | | | | | * add new sparse matrix library in src/sparse-matrix/umfpack/ git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1098 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
* rename this file so we don't have two files "misc.cc" (abeit injthorn2003-03-17
| | | | | | | | different directories) in this thorn, since that might give the linker hiccups and/or cause it to drop one file's object module git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@979 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
* add expansion_status_string()jthorn2003-03-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@969 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add gfns for skeletal patch systemjthorn2003-03-12
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@968 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
* * add support for sparse-matrix Jacobians ==> works!jthorn2003-01-22
| | | | | | | | | | * change default in param.ccl to use this * change default in src/include/config.h to default to no longer link in LAPACK routines * update documentation git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@931 f88db872-0e4f-0410-b76b-b9085cfa78c5
* start adding support for sparse matrix Jacobians (doesn't work yet)jthorn2003-01-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@926 f88db872-0e4f-0410-b76b-b9085cfa78c5
* change so that only processor #0 writes data filesjthorn2003-01-17
| | | | | | | (other processors still do all the computations, but writes are supressed) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@922 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
* oops, forgot to add this file! :)jthorn2003-01-15
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@919 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
* update to match recent change in terminology & file namesjthorn2002-11-23
| | | | | | | from H = horizon function to Theta = expansion git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@910 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
* change "apparent-horizon surface point ouside the Cactus grid"jthorn2002-10-26
| | | | | | | | | and "infinity or NaN detected in interpolated geometry" from level 1 and 0 warnings respectively, to level 2 (--> won't halt Cactus run, probably won't even print anything) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@863 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix a bug which made us use the -1 ("past") time level every 2nd timejthorn2002-10-22
| | | | | | | | | | | | | | | | | step of an evolution: Previously the code found the data pointers for the Cactus g_ij and K_ij grid arrays once in its setup phase, then used those all through the evolution. The problem is... the "current time level" data pointers change at each time step, because Cactus rotates the time-level arrays rather than copying all the (big) grid array data. This commit changes the coee to just lookup the Cactus variable indices in setup, then get a fresh set of data pointers each time we look for an apparent horizon (i.e. at each time step). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@850 f88db872-0e4f-0410-b76b-b9085cfa78c5
* reword/expand/combine error messages in geometry_is_finite()jthorn2002-10-13
| | | | | | | | --> now they give the actual interpolated geometry values (yes, all of them, 6 gij, 6 Kij, 18 partial_k gij) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@834 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * expand NaN checks in gridfns to "is finite" checksjthorn2002-10-12
| | | | | | | | | | | | | | | (also catchs +/-infinity) * move final stage of calculation of LHS function H(h) out of Maple-generated code (in src/gr.cg/) into hand-written code in src/gr/horizon_function.c (comput_H()) so we can check that HD > 0 before taking its sqrt. [HD <= 0 probably means the inerpolated g_ij isn't positive definite] * also systematize CCTK_VWarn() printing for various error conditions in horizon_function.c git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@829 f88db872-0e4f-0410-b76b-b9085cfa78c5
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-01 02:36:17 +0000