blob: 9f0fb6380e513f1c580f3fbfa5dbc1811559476c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
small things
set centroid variables for drift correction
I think we work properly with FishEye, but I should verify that
read IO parameters to set (default) directory in which we write output files
medium things
compute Gaussian curvature, cf AHFinder/src/AHFinder_gau.F
lib/make/extras/LAPACK/ shell scripts to replace environment variables
LAPACK_DIR, LAPACK_EXTRA_LIBS, and LAPACK_EXTRA_LIBDIRS
there should be a Cactus test suite
HDF5 data files (simple format and/or Werner Benger's fancy one)
move origin point to track moving BHs
read AHFinder Fourier coeffs for initial guess
write data files giving intersection of AH with a Cactus output plane
(or more generally, with a hyperslab)
export a "compute h(theta,phi)" function
so other thorns could test point membership in/outside AH etc
prevent h from going negative if the iteration isn't converging
(maybe use log(h) in Newton iteration?)
handle excision properly
(detect when interpolator tries to use data from excised region
and treat this as failure to evalute H(h) --> failure to find horizon
checkpoint/recover doesn't work
(needs storage of h gridfn in Cactus arrays,
or else some equivalent mechanism)
if two apparent horizons appear together, then move apart, right now
even with N_horizons > 1 we will probably have each of our AHs track
the same one; there should be some way to say "use AH 1 + something
as the initial guess for AH2" or suchlike, so we could find both AHs
large things
multi-processor (interpolator, overall control)
sparse matrix elliptic solver
handle quadrant/octant grids with rotating BCs
(needs Jacobian of one ghost zone to cover two different ghost zones)
"isolated horizons" computation of BH mass and spin
detect Cartoon and rotate geometry interpolation into Cartoon plane
|