aboutsummaryrefslogtreecommitdiff
path: root/src/elliptic
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
* 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
* 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
* many changes tojthorn2003-06-04
| | | | | | | | | | * use UMFPACK sparse matrix library * use new ILUCG sparse matrix library API * allow passing Cactus parameters down into the linear solvers * fix assorted small bugs git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1094 f88db872-0e4f-0410-b76b-b9085cfa78c5
* tiny wording tweak in an error message;jthorn2003-06-02
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1089 f88db872-0e4f-0410-b76b-b9085cfa78c5
* 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
* 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
* 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
* 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 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
* fix typo in commentjthorn2003-05-26
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1064 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
* 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
* change Fortran commentsjthorn2003-03-24
| | | | | | | | | | | | | c blah blah's blah to c blah blah is blah to avoid confusing the C preprocessor on these files git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1011 f88db872-0e4f-0410-b76b-b9085cfa78c5
* the Compaq f90 compiler complains about empty (or all-comment)jthorn2003-03-24
| | | | | | | | files, so we still have to define empty subroutines even if they're never called. :( Oh well, memory's cheap... git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1010 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 copies of header files that have to be included verbatim here :( :(jthorn2003-01-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@943 f88db872-0e4f-0410-b76b-b9085cfa78c5
* oops, forgot to add this file to cvsjthorn2003-01-30
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@942 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
* must define HAVE_DENSE_JACOBIAN or HAVE_ROW_SPARSE_JACOBIANjthorn2003-01-23
| | | | | | | | | in include/config.h rather than in elliptic/Jacobian.hh, since we want these to be usable in #ifdefs in elliptic/*_wrapper.F77 (which only include config.h) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@939 f88db872-0e4f-0410-b76b-b9085cfa78c5
* the cactus preprocessor requires that all preprocessor commandsjthorn2003-01-23
| | | | | | | | | start with # in column 1 :( :( :( ==> this change makes the code uglier and harder to read, but otherwise we can't compile on some platforms (eg my laptop) :( :( git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@938 f88db872-0e4f-0410-b76b-b9085cfa78c5
* this file renamed to 'lapack_wrapper.F77'jthorn2003-01-22
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@937 f88db872-0e4f-0410-b76b-b9085cfa78c5
* only call the single precision lapack/ilucg routines if #ifdef FP_IS_FLOAT,jthorn2003-01-22
| | | | | | | | | | only call the double precision lapack/ilucg routines if #ifdef FP_IS_DOUBLE this way we can link ok even on systems that have (eg) only double precision LAPACK, so long as that's allwe're using git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@936 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
* add more Jacobian-matrix stuffjthorn2003-01-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@930 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 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
* remove unneeded ; injthorn2002-10-07
| | | | | | | | | | | | namespace jtutil { ... }; (icc warns about this) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@810 f88db872-0e4f-0410-b76b-b9085cfa78c5
* This change fixes problems with some compilers recognizing only <vector>,jthorn2002-10-07
| | | | | | | | | | | | | | | | | and others recognizing only <vector.h>; change from #include <vector.h> to #include "stl_vector.hh" (this now after includeing "cctk.h"), where "stl_vector.hh" is a new header file that looks at the symbols HAVE_VECTOR and HAVE_VECTOR_H (defined by Cactus at configuration time) to see what to #include and whether or not to use an explicit using declaration to bring vector into the global namespace git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@809 f88db872-0e4f-0410-b76b-b9085cfa78c5
* - rename "config.hh" to "config.h" since it's now C++/C/Fortranjthorn2002-09-29
| | | | | | | | | | | | minor code tweaks for better portability (std:: namespace stuff) - switch all (new-style) #include <cstdio> style includes back to (old-style) #include <stdio.h> style - drop all using std::printf; declarations - change all uses of std::min and std::max to jtuti:: instead - drop all explicit std:: qualifiers git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@777 f88db872-0e4f-0410-b76b-b9085cfa78c5
* comment and whitespace changes onlyjthorn2002-09-20
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@758 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * add #include of "config.hh" so we know what Jacobian types are configuredjthorn2002-09-20
| | | | | | | * change C++ to C-style comments on #endif lines git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@757 f88db872-0e4f-0410-b76b-b9085cfa78c5
* * change from $Id: to $Header:jthorn2002-09-16
| | | | | | | * put LAPACK wrapper functions inside #ifdef HAVE_DENSE_JACOBIAN git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@747 f88db872-0e4f-0410-b76b-b9085cfa78c5
* change $Id:$ --> $Header:$jthorn2002-09-13
| | | | | | | | | rename src/util/ --> src/patch/ add Jacobian code for computing surface integrals add support for choosing ASCII or HDF5 data files (HDF5 not implemented yet) git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@727 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add #define HAVE_DENSE_JACOBIAN for dense Jacobiansjthorn2002-09-11
| | | | | | | | | -- in the future there will be such a #define for each Jacobian type, so we can easily disable ones we don't want to avoid having to link against the corresponding linear-solver libraries git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@712 f88db872-0e4f-0410-b76b-b9085cfa78c5
* move decoding of the Jacobian type into this directory (used to be in ../driver)jthorn2002-09-03
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@702 f88db872-0e4f-0410-b76b-b9085cfa78c5
* fix a bug where we took A(i,j) as the infinity-norm instead of |A(i,j)|jthorn2002-07-31
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@676 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add support for estimating condition number when we solve a linear systemjthorn2002-07-29
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@668 f88db872-0e4f-0410-b76b-b9085cfa78c5
* re-sync changes from laptopjthorn2002-07-26
| | | | | | | | - parameter to control how Jacobian is computed - can hardwire geometry to Schwarzschild/EF git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@661 f88db872-0e4f-0410-b76b-b9085cfa78c5
* various changes including d/dr terms in Jacobian by numerical perturbation,jthorn2002-07-22
| | | | | | | | | tweak I/O parameters, move printing Jacobian out of Jacobian class into test driver, drop unused array BLAS routines in jtutil:: git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@654 f88db872-0e4f-0410-b76b-b9085cfa78c5
* prototypes for LAPACK routinesjthorn2002-07-22
| | | | git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@649 f88db872-0e4f-0410-b76b-b9085cfa78c5
* add a whole bunch of changes from working at homejthorn2002-07-22
| | | | | | | --> AHFinderDirect now finds AH correctly for Kerr/offset!!! git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@648 f88db872-0e4f-0410-b76b-b9085cfa78c5