| Commit message (Collapse) | Author | Age |
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1562 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
used for the metric. Setting this to three unconditionally leads to problems when running with Carpet and prolongation_order_time=1 because all GFs in Carpet are expected to have the same number of timelevels: prolongation_order_time+1.
This is not the ideal solution. That would eliminate also variables ala metric_timelevels, because they directly depend on prolongation_order_time and should not have to be set (correctly) in a parameter file. This could be done automatically.
However, this patch for now uses metric_timelevels, in order to get AHFinderDirect working quickly.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1554 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
the patch allows the possibility to schedule
AHFinderDirect also at POSTPOSTINITIAL by adding the parameter
run_at_CCTK_POSTPOSTINITIAL (default equals to false). I have used it to
be sure that an AH was found after initial data computation and before
the last call to the Con2Prim routine of Whisky before evolution was
started (in order to have hydro variables properly excised before the
beginning of the evolution).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1550 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1541 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
AHFinderDirect to look at certain grid scalars to reset the origin
used for horizon finding.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1530 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
variables are usually not yet correctly initialised in this bin.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1527 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1526 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1484 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1483 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1433 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
* decrease the default verbose_level from "algorithm highlights"
down to "physics details"
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1423 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1421 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
In a multiprocessor run, AHFinderDirect finds different horizons
in parallel on different processors, and doesn't broadcast horizon
shapes to all processors unless it thinks they'll be needed. However,
the "the information will be needed" calculation doesn't take into
account the possibility of the HorizonRadiusInDirection() aliased
function being called on another processor, and in this case this
function returns a wrong result if AHFinderDirect hasn't broadcast
the horizon shapes.
* introduce a new Boolean parameter always_broadcast_horizon_shape
(defaults to false)
* document in interface.ccl, param.ccl, and doc/documentation.tex,
that this parameter should be set to true if you expect
HorizonRadiusInDirection() to be used
* inclusive-or this parameter into the "should we broadcast the horizon shape?"
calculation in src/driver/find_horizons.cc
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1419 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1416 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
-- thanks to Frank Loeffler for the basic patches for this
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1397 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
in CCTK_ANALYSIS (which has historically been the case, and remains
the default after this commit) or CCTK_POSTSTEP (which sometimes
comes in handy).
There should be no user-visible change in behavior from this commit.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1383 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
==> With this commit, AHFinderDirect now inherits from
AEIThorns/SphericalSurface, so you must have that thorn in
your configuration to be able to compile.
* add computation of surface quadrupole moments and areal radius
* expand BH_diagnostics file format to accomodate quadrupole moments
and areal radius, and also to include not-implemented-yet columns
for 9 more diagnostics which Erik has implemented in his branch
* some other small cleanups
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1329 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
for, from 10 all the way up to 100. Hopefully that will be enough for
now...
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1234 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
|
| |
so it now tests against the angular minimum of r_inner,
not r_min (= the angular minimum of r_horizon) like it did before
... docs update will follow shortly
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1216 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1210 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
min_horizon_size_for_mask
to
min_horizon_radius_points_for_mask
to make its semantics clearer
(thanks to Erik Schnetter for suggesting the renaming!)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1209 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
| |
real min_horizon_size_for_mask \
"only set mask if min horizon radius >= this number of local grid spacings"
docs will follow shortly...
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1208 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
so it can/must be specified for each horizon
--> gives a bit more flexibility for screening out "bad" cases
(alas this isn't in the thornguide yet :()
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1205 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1203 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1194 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1193 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1188 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
because ILUCG _usually_ works ok, but sometimes blows up and is ++slow :( :(
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1182 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
| |
this will make AHFinderDirect faster (typically about 2.5 times faster),
but less accurate (on analytic data at dx=0.1m or so, typical accuracies
of 1e-4m or so, instead of 1e-5m or so for order=3). In practice, for
most of our runs the AHFinderDirect (in)accuracy is dominated by the
inaccuracy of the numerically computed g_ij and K_ij fed into it, so
the slightly lower AH-finding accuracy is unimportant.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1178 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
|
|
|
|
|
|
|
|
|
|
|
|
| |
- add new parameters
BH_diagnostics_directory
h_directory
which both default to IO::out_dir
- change code to use these for output directories
- change code to auto-create output directories if they don't already exist
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1157 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
|
|
|
|
|
|
|
| |
to bool hardwire_Schwarzschild_geometry
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1147 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
| |
CCTK_InterpGridArrays() to interpolate the Cactus geometry grid functions
to the horizon position
-- this was never used in practice, and it conflicts with the changse
I'm about to commit to make AHFinderDirect work properly under Carpet
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1146 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1141 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
parameter to (we think) fix a serious bug where if one horizon goes
off the edge of the grid, we fail to find any other horizons
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1130 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
| |
since we need some AEILocalInterp features that LocalInterp doesn't have
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1128 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
new OpenDX macros (which can read AHFinderDirect gnuplot-format output files)
> Update of /numrelcvs/AEIPhysics/Visualization/OpenDX/Macros
> In directory cvsserv:/tmp/cvs-serv29388
>
> Added Files:
> ImportAHFinderDirectGnuplot.net
> ImportAHFinderDirectGnuplotPatch.net
> Log Message:
>
> Macros which can read the AHFinderDirect gnuplot output files
> directly into OpenDX.
There is a new Boolean parameter output_OpenDX_control_files
to control whether the control files (one per horizon) should be written,
but this defaults to true, so you don't need to do anything special to get
the control files.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1125 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
| |
(the case expansion = 0 is the AH, others should be surfaces nested inside
or outside the AH)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1118 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
| |
param.ccl
schedule.ccl
drop obselete find_AHs_at{poststep,postinitial} parameters
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1117 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1095 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1081 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1079 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1072 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
| |
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1059 f88db872-0e4f-0410-b76b-b9085cfa78c5
|
|
|
|
|
|
|
|
|
|
|
| |
* 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|