aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* oops, mistakenly committed this with lots of run-time checking ofjthorn2003-06-02
| | | | | | | | internal data structures ==> made this thorn run O(10) times slower! :( (this is ok for debugging, but not for production!) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1088 f88db872-0e4f-0410-b76b-b9085cfa78c5
* import of these files into AHFinderDirect cvsjthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1087 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add 1st draft of code for UMFPACK sparse-matrix solverjthorn2003-06-02
| | | | | | | | | | -- this is not tested yet, by default this code is *not* compiled at present (#undef HAVE_ROW_SPARSE_JACOBIAN__UMFPACK in ../incloude/config.h) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1086 f88db872-0e4f-0410-b76b-b9085cfa78c5
* define fp datatype roundoff threshold,jthorn2003-06-02
| | | | | | | add definition for UMFPACK sparse solver (currently *not* compiled by default) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1085 f88db872-0e4f-0410-b76b-b9085cfa78c5
* Jacobian::solve_linear_system() now returns *reciprocal* condition numberjthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1084 f88db872-0e4f-0410-b76b-b9085cfa78c5
* solve_linear_system() now returns *reciprocal* condition numberjthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1083 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * solve_linear_system() now returns *reciprocal* condition numberjthorn2003-06-02
| | | | | | | * add hooks for UMFPACK linear solver git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1082 f88db872-0e4f-0410-b76b-b9085cfa78c5
* comment tweakjthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1081 f88db872-0e4f-0410-b76b-b9085cfa78c5
* no more ILUCG wrapper routines!jthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1080 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add parameters for convergence criteria for ILUCG sparse linear solverjthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1079 f88db872-0e4f-0410-b76b-b9085cfa78c5
* set up and pass parameter structure down to Jacobian::solve_linear_system()jthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1078 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * add explicit parameter structure for solve_linear_system()jthorn2003-06-02
| | | | | | | | | * change row_sparse_Jacobian__ILUCG::solve_linear_system() to use ILUCG parameters in parameter structure for convergence tolerance and iteration limit git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1077 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * remove ilucg.h (moved to ../sparse-matrix/ilucg/)jthorn2003-06-02
| | | | | | | * remove ilucg_wrapper.F77 (don't need this any more) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1076 f88db872-0e4f-0410-b76b-b9085cfa78c5
* don't need this any more -- ILUCG subroutines modified to have directlyjthorn2003-06-02
| | | | | | | C-callable API git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1075 f88db872-0e4f-0410-b76b-b9085cfa78c5
* moved to ../sparse-matrix/ilucg/jthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1074 f88db872-0e4f-0410-b76b-b9085cfa78c5
* ilucg.h move here from ../../elliptic/jthorn2003-06-02
| | | | | | | | | | and update prototypes for remaining code [sd]ilucg.f77 convert from being Fortran logical functions to subroutines, add explicit error return argument git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1073 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add option for UMFPACK sparse matrix library (code isn't there yet)jthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1072 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add an example par file showing how to set the interpolator parametersjthorn2003-05-27
| | | | | | | to allow running with driver::ghost_size = 2 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1071 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix a bug where the old-style mask was set to the buffer valuejthorn2003-05-27
| | | | | | | | in the outside region :( :( -- thanks to Michael Koppitz for spotting this! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1070 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix output file/directory namesjthorn2003-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1069 f88db872-0e4f-0410-b76b-b9085cfa78c5
* new par files with shorter & more mnemonic namesjthorn2003-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1068 f88db872-0e4f-0410-b76b-b9085cfa78c5
* get ride of these out-of-date par filesjthorn2003-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1067 f88db872-0e4f-0410-b76b-b9085cfa78c5
* get rid of these out-of-date par filesjthorn2003-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1066 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix IOASCII parameters for gnuplot outputjthorn2003-05-27
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1065 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix typo in commentjthorn2003-05-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1064 f88db872-0e4f-0410-b76b-b9085cfa78c5
* clarify pseudocode algorithm to explain how noshrink option worksjthorn2003-05-26
| | | | | | | for excision mask git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1063 f88db872-0e4f-0410-b76b-b9085cfa78c5
* explicitly give what the horizon shape is when you specifyjthorn2003-05-23
| | | | | | | | AHFinderDirect::initial_guess_method[n] = "Kerr/Kerr" or "Kerr/Kerr-Schild" git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1062 f88db872-0e4f-0410-b76b-b9085cfa78c5
* - update to use new return codes from SpaceMask (-1=error instead of 0=error)jthorn2003-05-23
| | | | | | | | | | - split up a single huge function into several smaller/more readable subfunctions - general cleanup of the logic for the noshrink option - more documentation on same git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1061 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add to "small things" sectionjthorn2003-05-16
| | | | | | | | | | | there should be an option to specify using the parameter-file initial guess at selected times (or cactus iteration numbers) instead of always using the AH posn from the previous time step -- this would be useful if there are an inner & outer horizon which start out coincident and you want to track them separately git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1060 f88db872-0e4f-0410-b76b-b9085cfa78c5
* allow searching for up to 10 horizons in a slicejthorn2003-05-14
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1059 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
* expand/clarify some commentsjthorn2003-05-11
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1055 f88db872-0e4f-0410-b76b-b9085cfa78c5
* tiny code tweakjthorn2003-05-09
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1054 f88db872-0e4f-0410-b76b-b9085cfa78c5
* When unpacking the interprocessor synchronization results from ourjthorn2003-05-09
| | | | | | | | | | | (CCTK_REAL) buffer back to user-typed variables, the code used to do static_cast<enum blah>(a CCTK_REAL value). It seems the Hitachi C++ compiler doesn't like this, so change the code to static_cast<enum blah>( static_cast<int>(a CCTK_REAL value) ) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1053 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
* add #ifndef include guards to a bunch of files, even though they shouldn'tjthorn2003-05-06
| | | | | | | | | | | | be necessary (I never include a file twice!)... since certain lame/broken/dumb compilers (like DEC/Compaq/HP/whatever-they-call-themselves-this-week C++ version 6.something on Alpha Linux) still give multiple inclusions even with automagic template instantiation turned off :( :( :( -- thanks to Frank Loeffler for helping track this problem down! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1050 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add /*NOTREACHED*/ comments after some error_exit() callsjthorn2003-05-06
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1049 f88db872-0e4f-0410-b76b-b9085cfa78c5
* Harwell sparse matrix library, as obtained from netlib,jthorn2003-05-05
| | | | | | | plus my README file explaining this, plus my www file giving url where I got it git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1048 f88db872-0e4f-0410-b76b-b9085cfa78c5
* detailed description moved down to subdirectoriesjthorn2003-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1047 f88db872-0e4f-0410-b76b-b9085cfa78c5
* explain the origin of these routinesjthorn2003-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1046 f88db872-0e4f-0410-b76b-b9085cfa78c5
* move sparse matrix routines from elliptic/ to new directory sparse-matrix/jthorn2003-05-05
| | | | | | | and subdirectories git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1045 f88db872-0e4f-0410-b76b-b9085cfa78c5
* sparse-matrix routines are now moved to ../sparse-matrix/ subdirectoriesjthorn2003-05-05
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1044 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add ILUCG routines in a separate directoryjthorn2003-05-05
| | | | | | | (split up single/double precision into separate files) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1043 f88db872-0e4f-0410-b76b-b9085cfa78c5
* split former initial-data optionjthorn2003-05-02
| | | | | | | | | | | | | | | | | | | keyword initial_guess_method[5] \ "method used to set up initial guess for apparent horizon shape" { "read from file" :: "read from input file" into two options keyword initial_guess_method[5] \ "method used to set up initial guess for apparent horizon shape" { "read from named file" :: "read from explicitly-named input file" "read from h file" :: \ "read from input file named the same as the (later) h output file" git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1042 f88db872-0e4f-0410-b76b-b9085cfa78c5
* Change from using a private Strlcat() function to the Util_Strlcat()jthorn2003-05-02
| | | | | | | | | | | | | one that's now in the Cactus flesh. ***** ***** After this change, you will need to cvs-update your copy of the ***** Cactus flesh (files src/util/String.c and src/util/util_String.h) ***** in order to recompile this thorn. ***** git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1041 f88db872-0e4f-0410-b76b-b9085cfa78c5
* small tweaksjthorn2003-05-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1040 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix typo in comment -- no executable code changesjthorn2003-04-10
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1039 f88db872-0e4f-0410-b76b-b9085cfa78c5