| Commit message (Collapse) | Author | Age |
... | |
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1229 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
| |
switch to UMFPACK 4.1 and ATLAS BLAS
(http://www.cise.ufl.edu/research/sparse/umfpack/ and links therein)
--> should make AHFinderDirect considerablyfaster
(optionally) remember the "sparse structure decompositon" for all
Newton solves for a given AH within a Cactus run
--> uses more memory, but should make AHFinderDirect much faster
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1228 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
#
# To detect some divergence cases faster, we also deem any horizon where
# the infinity-norm of Theta exceeds this parameter, to be "not found".
# The default is effectively "no limit".
#
real max_allowable_Theta \
"max ||Theta||_infinity allowed for any trial surface \
before we give up and say we can't find this horizon"
{
(0.0:* :: "any positive real number"
} 1.0e10
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1225 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
export a "compute h(theta,phi)" function
so other thorns could test point membership in/outside AH etc
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1224 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
on a per-horizon basis:
The old parameter set_mask is now replaced by two new parameters:
set_mask_for_all_horizons
set_mask_for_individual_horizon[n]
For most purposes, you can just change set_mask to
set_mask_for_all_horizons in your parameter files. If you
really want to control the setting of the mask on a per-horizon basis,
read section 4.6 in the thorn guide to find out how the new parameters
work.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1222 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
give a cogent error message (instead of core-dumping like we do right now)
if we try to set a mask and nobody's registered the bitfield yet
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1220 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1217 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
option
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1214 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1213 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1212 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
min_horizon_radius_points_for_maks
parameter
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1211 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
max_allowable_horizon_radius
option
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1206 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1201 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1192 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1191 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1189 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
maybe try the TATElliptic interface for solving the AH equation
as an alternative to the home-grown solvers I use now???
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1186 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
now that order=2 is the default geometry interpolation
(--> can use driver::ghost_size=2)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1183 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1181 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
desribe this and the tradeoffs for this vs order=3
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1179 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
which we want to find a given horizon(s)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1177 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1176 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1175 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
specified AH centroid -- this uses the new function-aliasing DriftCorrect
interface developed by Peter Diener and myself (with much assistance
from Ian Hawke)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1174 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1172 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
this directory; since this directory is never seen by the Cactus build
system there's no problem with either name...
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1170 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1163 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1161 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1159 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
< read IO parameters to set (default) directory in which we write output files
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1158 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
to CCTK_ANALYSIS (which does, yea!)
...
Alas, due to the @*(#&%@*&(%#$ "trigger" behavior of Cactus,
the only way to force a routine to run in CCTK_ANALYSIS is
to use a dummy I/O method to fake doing output at every time
step. Fortunately, Erik Schnetter has already written the code
for this, so AHFinderDirect now just used that. Thus AHFinderDirect
now depends on thorn TAT/always. All thornlists and par files
will have to be changed to compile in always and activate it
before activating AHFinderDirect.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1156 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
| |
how_often_to_output_h --> output_h_every
how_often_to_output_Theta --> output_Theta_every
to match the usual Cactus convention for "*_every" parameters
(including our own find_every parameter)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1154 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
with new integer parameter find_every
so user can choose how often to find AHs
(with Carpet, at present you can only find AHs on a time level
where the coarsest grid is defined ==> you must use find_every )
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1153 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
HIGH PRIORITY: once Thomas Radke fixes bug CactusPUGH/1563 and I make
some corresponding changes in AEIThorns/AEILocalInterp, I need to
change AHFinderDirect to properly handle the case where you're finding
AHs on multiple processors, some surfaces are outside the grid, and
other surfaces are all inside the grid: right now AHFinderDirect
gives up finding *all* the AHs, even the ones which are entirely
within the grid.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1145 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1140 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
< Peter Diener would like a "+xz quadrant (mirrored)" patch system type
< (eg for BBH0 with BHs on the +/- y axis)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1139 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1137 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
Peter Diener would like a "+xz quadrant (mirrored)" patch system type
(eg for BBH0 with BHs on the +/- y axis)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1132 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
add an option to generate OpenDX headers for the ASCII output files,
cf Peter's perl script
numrelcvs:AEIPhysics/Utilities/ahd_dx
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1131 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
HIGH PRIORITY: once Thomas Radke fixes bug CactusPUGH/1563 and I make
some corresponding changes in AEIThorns/AEILocalInterp, I need to
change AHFinderDirect to properly handle the case where you're finding
AHs on multiple processors, some surfaces are outside the grid, and
other surfaces are all inside the grid: right now AHFinderDirect
gives up finding *all* the AHs, even the ones which are entirely
within the grid.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1129 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
OpenDX macros
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1126 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1124 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1122 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
document new AHFinderDirect option to find surfaces of constant expansion
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1121 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1120 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1119 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
also document move to AEIThorns and s/LocalInterp/AEILocalInterp/g
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1116 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
(it was formerly in AEIDevelopment)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1113 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1112 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
change arrangement from AEIDevelopment to AEIThorns
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1111 f88db872-0e4f-0410-b76b-b9085cfa78c5
|