| Commit message (Collapse) | Author | Age |
... | |
|
|
|
|
|
|
| |
for new style masks
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1031 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1030 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1018 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
*and* new-style (CCTK_INT) mask simultaneously -- Frank Loeffler
asked for this, since Whiskey needs it (at least until all of
Cactus is converted to the new style...)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1015 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
get rid of our own mask gridfns since SpaceMask already defines them
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@987 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
-- 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
|
|
|
|
|
|
|
| |
want to change closer to the top
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@973 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
(i.e. require the user to specify a large enough ghost size)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@955 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
| |
(18 was more than needed for most runs),
also add more comments on performance vs resolution
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@940 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@927 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
(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
|
|
|
|
|
|
|
|
|
|
|
| |
(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
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
they're out-of-date (and would quickly become out of date again if
I updated them now) -- documentation.tex now gives a summary, and
for the exact info points to either the file's header comments,
or the source code.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@895 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
add support for computing BH circumferences in xy, xz, yz planes
as requested by Ed recently
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@893 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
| |
if all AHs are entirely within a single processor's chunk of the grid,
then set AHFinderDirect::processor_number = 42
(or whatever that processor's processor number is)
and things should work ok
(this is not yet properly documented :( )
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@882 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@874 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@869 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
| |
because (for reasons I don't understand) all the Cactus arrays
show up as NULL pointers at run-time :( :(
* Add an option to also write the BH diagnostics to an output file
inside AHFinderDirect -- this works fine. Alas the output file
format is hard-wired to a simple ASCII one. :(
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@860 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
to avoid confusion with future file-format parameters for other types
of output files (eg for BH diagnostics)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@859 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
in Cactus arrays so they're accessible to other thorns and/or I/O methods
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@858 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
parameters; the top-level method parameter is now just
keyword method "what should this thorn do for each apparent horizon?"
{
# these options are mostly for testing/debugging
"horizon function" :: "evaluate the LHS function H(h)"
"test Jacobian" :: "compute/print the J[H(h)] Jacobian matrix"
# this is for normal apparent horizon finding
"find horizon" :: "find the apparent horizon"
} "find horizon"
The keyword strings have changed a bit, but since most users will
not need to change the default value, this shouldn't change their
parameter files.
* Revise driver logic to clean up handling of the different values
of the top-level method parameter.
* Slightly tweak output format for BH diagnostics, and documentation
(thorn guide) to make it explicit that what we're calculating is the
irreducible mass of the apparent horizon.
Modified Files:
param.ccl doc/documentation.tex par/misner-run.par
run/test-ahfinderdirect/misc/Kerr.par
run/test-ahfinderdirect/misc/try-5.par
run/test-ahfinderdirect/misc/try-7.5-debug.par
run/test-ahfinderdirect/misc/try-7.5.par
run/test-ahfinderdirect/misc/try-9.par src/driver/driver.hh
src/driver/find_horizons.cc src/driver/setup.cc
Added Files:
# new par file to compute/print the analytical horizon
# for this Kerr/Schild test, for comparison with numerical ones
run/test-ahfinderdirect/misc/try-7.5-horizon.par
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@855 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
from ".gnuplot" to ".gp" to avoid confusion with gnuplot command files
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@843 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
("sphere" --> "coordinate sphere", "ellipsoid" --> "coordinate ellipsoid")
* more comments in param.ccl
* add option to allow only writing "large" data files every N time steps
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@838 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
(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
|
|
|
|
|
|
|
| |
(separate parameters for checking h and the interpolated geometry fields)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@828 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
one Newton iteration to the next, Delta_h, at each Newton iteration
(for help debugging convergence problems)
also move parameters for "should we output {h,H,Delta_h} from
struct IO_info to struct solver_info
also make the instructions on setting environment variables in
make.configuration.defn even more verbose...
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@820 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
permute order of som pars
schedule.ccl
remove initial-data setting -- this is now done inside
the AHFinderDirect_find_horizons() routine
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@803 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
| |
// guess is likely to be rather inaccurate, so we may need a
// larger number of iterations. But if we've found this horizon
// before, then we have its previous position as an initial guess,
// so we shouldn't need as many iterations.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@795 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@787 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@754 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
add new "automatic choice" option for surface_integral_method
and make it the default
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@741 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
* add a bunch of comments about choosing the interpolator parameters
* change the default interpatch interpolator parameters to use
Lagrange rather than Hermite, because the current Hermite interpolator
is rather inaccurate when off-centered
* split up initial_guess__Kerr__* parameters into separate
initial_guess__Kerr_Kerr__* and initial_guess__Kerr_KerrSchild__*
parameters to match the thorn Exact conventions
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@733 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@726 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
but commented out
(it's more elegant, but appears to tickle a CST bug :( :( )
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@723 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
via surface integrals over the horizon
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@715 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@708 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
and to handle finding multiple AHs in a single slice
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@698 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@696 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
"synchronize Jacobian" to "ghost zone Jacobian" to clarify just what's
being tested here
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@682 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
--> seems to eliminate (or at least greatly reduce) convergence
problems due to non-smoothness of interpolated geometry
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@678 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
- add support for condition number estimating when solving linear system
- cleanup Schw/EF analytic geometry option
- add option to dump h and H at each Newton iteration for debuggin
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@669 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@665 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@653 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
--> 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
|