diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-10-17 22:27:44 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-10-17 22:27:44 +0000 |
commit | f6d399914423ae47efcad3beb3c9ca6e2e96488f (patch) | |
tree | 4d1cad667a9c5e580c2de2b543ac2f9f17bee546 | |
parent | 992b4ad347b6a9e861dc81bf651b64bef60325b3 (diff) |
expand documentation to describe the multipatch system,
and add a figure illustrating this
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@847 f88db872-0e4f-0410-b76b-b9085cfa78c5
-rw-r--r-- | doc/3patch.gnuplot | 52 | ||||
-rw-r--r-- | doc/3patch.par | 50 | ||||
-rw-r--r-- | doc/AEIDevelopment_AHFinderDirect_3patch.eps | 1723 | ||||
-rw-r--r-- | doc/documentation.tex | 79 | ||||
-rw-r--r-- | doc/makefile | 28 |
5 files changed, 1920 insertions, 12 deletions
diff --git a/doc/3patch.gnuplot b/doc/3patch.gnuplot new file mode 100644 index 0000000..dc7ec05 --- /dev/null +++ b/doc/3patch.gnuplot @@ -0,0 +1,52 @@ +# gnuplot command file to produce postscript file for 3patch figure + +set term postscript eps color solid +set size 1.5, 2.25 +set output 'AEIDevelopment_AHFinderDirect_3patch.eps' + +unset key +set view 60, 111 +set border 0 + +unset xtics +unset ytics +unset ztics + +splot '-' w lines -1, \ + '-' w lines -1, \ + '-' w lines -1, \ + '<../src/misc/select.patch +z <3patch.h.t0.ah1.gp' \ + using 4:5:6 w lines 1, \ + '<../src/misc/select.patch +x <3patch.h.t0.ah1.gp' \ + using 4:5:6 w lines 2, \ + '<../src/misc/select.patch +y <3patch.h.t0.ah1.gp' \ + using 4:5:6 w lines 3, \ + "<../src/misc/select.patch +z <3patch.h.t0.ah1.gp \ + | gawk '($4 == 0) && ($5 >= 0.0)' -" \ + using 4:5:6 w lines -1, \ + "<../src/misc/select.patch +z <3patch.h.t0.ah1.gp \ + | gawk '($4 >= 0) && ($5 == 0.0)' -" \ + using 4:5:6 w lines -1, \ + "<../src/misc/select.patch +x <3patch.h.t0.ah1.gp \ + | gawk '($5 == 0) && ($6 >= 0.0)' -" \ + using 4:5:6 w lines -1, \ + "<../src/misc/select.patch +x <3patch.h.t0.ah1.gp \ + | gawk '($5 >= 0) && ($6 == 0.0)' -" \ + using 4:5:6 w lines -1, \ + "<../src/misc/select.patch +y <3patch.h.t0.ah1.gp \ + | gawk '($4 == 0) && ($6 >= 0.0)' -" \ + using 4:5:6 w lines -1, \ + "<../src/misc/select.patch +y <3patch.h.t0.ah1.gp \ + | gawk '($4 >= 0) && ($6 == 0.0)' -" \ + using 4:5:6 w lines -1 +0.0 0.0 0.0 +1.3 0.0 0.0 +eof +0.0 0.0 0.0 +0.0 1.3 0.0 +eof +0.0 0.0 0.0 +0.0 0.0 1.3 +eof + +set output diff --git a/doc/3patch.par b/doc/3patch.par new file mode 100644 index 0000000..7905a3e --- /dev/null +++ b/doc/3patch.par @@ -0,0 +1,50 @@ +# par file to write out data used to produce 3patch figure + +ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling StaticConformal CoordGauge Exact AHFinderDirect" + +# PUGH +Driver::ghost_size = 0 +Driver::global_nx = 31 +Driver::global_ny = 31 +Driver::global_nz = 31 + +# CartGrid3D +Grid::domain = "full" +Grid::type = "byspacing" +Grid::dxyz = 0.1 +Grid::xyzmin = -3.0 +Grid::xyzmax = 3.0 + +# initial data only +cactus::cctk_itlast = 0 + +# ADMBase +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_data = "exact" +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" +ADMBase::metric_type = "physical" + +# Exact +Exact::exact_model = "Schwarzschild/EF" +Exact::Schwarzschild_EF__mass = 1.0 + +AHFinderDirect::find_AHs_at_poststep = "false" # initial data only +AHFinderDirect::method = "horizon function" # just evaluate H(h) function +AHFinderDirect::output_ghost_zones_for_h = "true" +AHFinderDirect::h_base_file_name = "3patch.h" + +AHFinderDirect::patch_system_type = "+xyz octant" +AHFinderDirect::delta_drho_dsigma = 4.5 + +AHFinderDirect::N_horizons = 1 +AHFinderDirect::origin_x[1] = 0.0 +AHFinderDirect::origin_y[1] = 0.0 +AHFinderDirect::origin_z[1] = 0.0 + +AHFinderDirect::initial_guess_method = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 1.0 diff --git a/doc/AEIDevelopment_AHFinderDirect_3patch.eps b/doc/AEIDevelopment_AHFinderDirect_3patch.eps new file mode 100644 index 0000000..be02964 --- /dev/null +++ b/doc/AEIDevelopment_AHFinderDirect_3patch.eps @@ -0,0 +1,1723 @@ +%!PS-Adobe-2.0 EPSF-2.0 +%%Title: AEIDevelopment_AHFinderDirect_3patch.eps +%%Creator: gnuplot 3.8h patchlevel 0 +%%CreationDate: Thu Oct 17 17:38:51 2002 +%%DocumentFonts: (atend) +%%BoundingBox: 50 50 590 617 +%%Orientation: Portrait +%%EndComments +/gnudict 256 dict def +gnudict begin +/Color true def +/Solid true def +/gnulinewidth 5.000 def +/userlinewidth gnulinewidth def +/vshift -46 def +/dl {10.0 mul} def +/hpt_ 31.5 def +/vpt_ 31.5 def +/hpt hpt_ def +/vpt vpt_ def +/M {moveto} bind def +/L {lineto} bind def +/R {rmoveto} bind def +/V {rlineto} bind def +/N {newpath moveto} bind def +/f {rlineto fill} bind def +/vpt2 vpt 2 mul def +/hpt2 hpt 2 mul def +/Lshow { currentpoint stroke M + 0 vshift R show } def +/Rshow { currentpoint stroke M + dup stringwidth pop neg vshift R show } def +/Cshow { currentpoint stroke M + dup stringwidth pop -2 div vshift R show } def +/UP { dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def + /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def } def +/DL { Color {setrgbcolor Solid {pop []} if 0 setdash } + {pop pop pop Solid {pop []} if 0 setdash} ifelse } def +/BL { stroke userlinewidth 2 mul setlinewidth } def +/AL { stroke userlinewidth 2 div setlinewidth } def +/UL { dup gnulinewidth mul /userlinewidth exch def + dup 1 lt {pop 1} if 10 mul /udl exch def } def +/PL { stroke userlinewidth setlinewidth } def +/LTb { BL [] 0 0 0 DL } def +/LTa { AL [1 udl mul 2 udl mul] 0 setdash 0 0 0 setrgbcolor } def +/LT0 { PL [] 1 0 0 DL } def +/LT1 { PL [4 dl 2 dl] 0 1 0 DL } def +/LT2 { PL [2 dl 3 dl] 0 0 1 DL } def +/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def +/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def +/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def +/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def +/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def +/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def +/Pnt { stroke [] 0 setdash + gsave 1 setlinecap M 0 0 V stroke grestore } def +/Dia { stroke [] 0 setdash 2 copy vpt add M + hpt neg vpt neg V hpt vpt neg V + hpt vpt V hpt neg vpt V closepath stroke + Pnt } def +/Pls { stroke [] 0 setdash vpt sub M 0 vpt2 V + currentpoint stroke M + hpt neg vpt neg R hpt2 0 V stroke + } def +/Box { stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M + 0 vpt2 neg V hpt2 0 V 0 vpt2 V + hpt2 neg 0 V closepath stroke + Pnt } def +/Crs { stroke [] 0 setdash exch hpt sub exch vpt add M + hpt2 vpt2 neg V currentpoint stroke M + hpt2 neg 0 R hpt2 vpt2 V stroke } def +/TriU { stroke [] 0 setdash 2 copy vpt 1.12 mul add M + hpt neg vpt -1.62 mul V + hpt 2 mul 0 V + hpt neg vpt 1.62 mul V closepath stroke + Pnt } def +/Star { 2 copy Pls Crs } def +/BoxF { stroke [] 0 setdash exch hpt sub exch vpt add M + 0 vpt2 neg V hpt2 0 V 0 vpt2 V + hpt2 neg 0 V closepath fill } def +/TriUF { stroke [] 0 setdash vpt 1.12 mul add M + hpt neg vpt -1.62 mul V + hpt 2 mul 0 V + hpt neg vpt 1.62 mul V closepath fill } def +/TriD { stroke [] 0 setdash 2 copy vpt 1.12 mul sub M + hpt neg vpt 1.62 mul V + hpt 2 mul 0 V + hpt neg vpt -1.62 mul V closepath stroke + Pnt } def +/TriDF { stroke [] 0 setdash vpt 1.12 mul sub M + hpt neg vpt 1.62 mul V + hpt 2 mul 0 V + hpt neg vpt -1.62 mul V closepath fill} def +/DiaF { stroke [] 0 setdash vpt add M + hpt neg vpt neg V hpt vpt neg V + hpt vpt V hpt neg vpt V closepath fill } def +/Pent { stroke [] 0 setdash 2 copy gsave + translate 0 hpt M 4 {72 rotate 0 hpt L} repeat + closepath stroke grestore Pnt } def +/PentF { stroke [] 0 setdash gsave + translate 0 hpt M 4 {72 rotate 0 hpt L} repeat + closepath fill grestore } def +/Circle { stroke [] 0 setdash 2 copy + hpt 0 360 arc stroke Pnt } def +/CircleF { stroke [] 0 setdash hpt 0 360 arc fill } def +/C0 { BL [] 0 setdash 2 copy moveto vpt 90 450 arc } bind def +/C1 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 90 arc closepath fill + vpt 0 360 arc closepath } bind def +/C2 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 90 180 arc closepath fill + vpt 0 360 arc closepath } bind def +/C3 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 180 arc closepath fill + vpt 0 360 arc closepath } bind def +/C4 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 180 270 arc closepath fill + vpt 0 360 arc closepath } bind def +/C5 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 90 arc + 2 copy moveto + 2 copy vpt 180 270 arc closepath fill + vpt 0 360 arc } bind def +/C6 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 90 270 arc closepath fill + vpt 0 360 arc closepath } bind def +/C7 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 270 arc closepath fill + vpt 0 360 arc closepath } bind def +/C8 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 270 360 arc closepath fill + vpt 0 360 arc closepath } bind def +/C9 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 270 450 arc closepath fill + vpt 0 360 arc closepath } bind def +/C10 { BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill + 2 copy moveto + 2 copy vpt 90 180 arc closepath fill + vpt 0 360 arc closepath } bind def +/C11 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 180 arc closepath fill + 2 copy moveto + 2 copy vpt 270 360 arc closepath fill + vpt 0 360 arc closepath } bind def +/C12 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 180 360 arc closepath fill + vpt 0 360 arc closepath } bind def +/C13 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 0 90 arc closepath fill + 2 copy moveto + 2 copy vpt 180 360 arc closepath fill + vpt 0 360 arc closepath } bind def +/C14 { BL [] 0 setdash 2 copy moveto + 2 copy vpt 90 360 arc closepath fill + vpt 0 360 arc } bind def +/C15 { BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill + vpt 0 360 arc closepath } bind def +/Rec { newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto + neg 0 rlineto closepath } bind def +/Square { dup Rec } bind def +/Bsquare { vpt sub exch vpt sub exch vpt2 Square } bind def +/S0 { BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare } bind def +/S1 { BL [] 0 setdash 2 copy vpt Square fill Bsquare } bind def +/S2 { BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare } bind def +/S3 { BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare } bind def +/S4 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare } bind def +/S5 { BL [] 0 setdash 2 copy 2 copy vpt Square fill + exch vpt sub exch vpt sub vpt Square fill Bsquare } bind def +/S6 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare } bind def +/S7 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill + 2 copy vpt Square fill + Bsquare } bind def +/S8 { BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare } bind def +/S9 { BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare } bind def +/S10 { BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill + Bsquare } bind def +/S11 { BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill + Bsquare } bind def +/S12 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare } bind def +/S13 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill + 2 copy vpt Square fill Bsquare } bind def +/S14 { BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill + 2 copy exch vpt sub exch vpt Square fill Bsquare } bind def +/S15 { BL [] 0 setdash 2 copy Bsquare fill Bsquare } bind def +/D0 { gsave translate 45 rotate 0 0 S0 stroke grestore } bind def +/D1 { gsave translate 45 rotate 0 0 S1 stroke grestore } bind def +/D2 { gsave translate 45 rotate 0 0 S2 stroke grestore } bind def +/D3 { gsave translate 45 rotate 0 0 S3 stroke grestore } bind def +/D4 { gsave translate 45 rotate 0 0 S4 stroke grestore } bind def +/D5 { gsave translate 45 rotate 0 0 S5 stroke grestore } bind def +/D6 { gsave translate 45 rotate 0 0 S6 stroke grestore } bind def +/D7 { gsave translate 45 rotate 0 0 S7 stroke grestore } bind def +/D8 { gsave translate 45 rotate 0 0 S8 stroke grestore } bind def +/D9 { gsave translate 45 rotate 0 0 S9 stroke grestore } bind def +/D10 { gsave translate 45 rotate 0 0 S10 stroke grestore } bind def +/D11 { gsave translate 45 rotate 0 0 S11 stroke grestore } bind def +/D12 { gsave translate 45 rotate 0 0 S12 stroke grestore } bind def +/D13 { gsave translate 45 rotate 0 0 S13 stroke grestore } bind def +/D14 { gsave translate 45 rotate 0 0 S14 stroke grestore } bind def +/D15 { gsave translate 45 rotate 0 0 S15 stroke grestore } bind def +/DiaE { stroke [] 0 setdash vpt add M + hpt neg vpt neg V hpt vpt neg V + hpt vpt V hpt neg vpt V closepath stroke } def +/BoxE { stroke [] 0 setdash exch hpt sub exch vpt add M + 0 vpt2 neg V hpt2 0 V 0 vpt2 V + hpt2 neg 0 V closepath stroke } def +/TriUE { stroke [] 0 setdash vpt 1.12 mul add M + hpt neg vpt -1.62 mul V + hpt 2 mul 0 V + hpt neg vpt 1.62 mul V closepath stroke } def +/TriDE { stroke [] 0 setdash vpt 1.12 mul sub M + hpt neg vpt 1.62 mul V + hpt 2 mul 0 V + hpt neg vpt -1.62 mul V closepath stroke } def +/PentE { stroke [] 0 setdash gsave + translate 0 hpt M 4 {72 rotate 0 hpt L} repeat + closepath stroke grestore } def +/CircE { stroke [] 0 setdash + hpt 0 360 arc stroke } def +/Opaque { gsave closepath 1 setgray fill grestore 0 setgray closepath } def +/DiaW { stroke [] 0 setdash vpt add M + hpt neg vpt neg V hpt vpt neg V + hpt vpt V hpt neg vpt V Opaque stroke } def +/BoxW { stroke [] 0 setdash exch hpt sub exch vpt add M + 0 vpt2 neg V hpt2 0 V 0 vpt2 V + hpt2 neg 0 V Opaque stroke } def +/TriUW { stroke [] 0 setdash vpt 1.12 mul add M + hpt neg vpt -1.62 mul V + hpt 2 mul 0 V + hpt neg vpt 1.62 mul V Opaque stroke } def +/TriDW { stroke [] 0 setdash vpt 1.12 mul sub M + hpt neg vpt 1.62 mul V + hpt 2 mul 0 V + hpt neg vpt -1.62 mul V Opaque stroke } def +/PentW { stroke [] 0 setdash gsave + translate 0 hpt M 4 {72 rotate 0 hpt L} repeat + Opaque stroke grestore } def +/CircW { stroke [] 0 setdash + hpt 0 360 arc Opaque stroke } def +/BoxFill { gsave Rec 1 setgray fill grestore } def +/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont +dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall +currentdict end definefont pop +end +%%EndProlog +gnudict begin +gsave +50 50 translate +0.050 0.050 scale +0 setgray +newpath +(Helvetica) findfont 140 scalefont setfont +1.000 UL +LTb +1.000 UP +1.000 UL +LTb +4199 5261 M +2498 2952 L +1.000 UL +LTb +4199 5261 M +8629 4374 L +1.000 UL +LTb +4199 5261 M +0 4284 V +1.000 UL +LT0 +7065 6805 M +-52 -78 V +-57 -81 V +-63 -83 V +-70 -87 V +-77 -92 V +-84 -96 V +-94 -101 V +6465 6078 L +6349 5962 L +6221 5837 L +6076 5702 L +5913 5555 L +5728 5395 L +5520 5221 L +6910 7056 M +-56 -85 V +-64 -89 V +-70 -92 V +-77 -98 V +-85 -103 V +-94 -108 V +6361 6365 L +6248 6242 L +6123 6111 L +5984 5969 L +5830 5817 L +5660 5653 L +5470 5476 L +5259 5286 L +6739 7295 M +-61 -90 V +-69 -96 V +-77 -102 V +-84 -107 V +-92 -115 V +6255 6664 L +6144 6535 L +6023 6397 L +5890 6250 L +5745 6093 L +5585 5924 L +5411 5744 L +5221 5551 L +5013 5346 L +6553 7521 M +-67 -95 V +-74 -102 V +-82 -110 V +-90 -117 V +-99 -126 V +6034 6838 L +5917 6695 L +5790 6543 L +5653 6380 L +5504 6208 L +5343 6023 L +5169 5828 L +4982 5620 L +4782 5402 L +6352 7733 M +-71 -99 V +-79 -108 V +-87 -117 V +-96 -127 V +5916 7146 L +5804 7000 L +5683 6844 L +5552 6678 L +5413 6501 L +5263 6314 L +5103 6115 L +4933 5905 L +4754 5684 L +4564 5453 L +6139 7930 M +-76 -102 V +-83 -114 V +-92 -124 V +5788 7455 L +5681 7308 L +5566 7151 L +5442 6983 L +5311 6803 L +5171 6613 L +5024 6411 L +4868 6198 L +4705 5975 L +4535 5742 L +4358 5500 L +5914 8110 M +-80 -105 V +-88 -118 V +-95 -130 V +5548 7613 L +5438 7458 L +5321 7290 L +5197 7110 L +5067 6918 L +4930 6715 L +4787 6500 L +4638 6275 L +4484 6040 L +4325 5796 L +4163 5544 L +5679 8272 M +-85 -106 V +-91 -122 V +-98 -136 V +5300 7758 L +5189 7594 L +5072 7416 L +4949 7226 L +4822 7023 L +4689 6808 L +4553 6582 L +4412 6346 L +4269 6100 L +4124 5846 L +3977 5585 L +5435 8416 M +-88 -107 V +-95 -125 V +5151 8043 L +5045 7887 L +4934 7715 L +4819 7529 L +4699 7330 L +4576 7117 L +4450 6892 L +4322 6656 L +4192 6410 L +4061 6155 L +3929 5892 L +3799 5624 L +5184 8541 M +-92 -108 V +-98 -127 V +4892 8162 L +4786 8000 L +4676 7822 L +4563 7629 L +4448 7422 L +4331 7201 L +4213 6968 L +4094 6723 L +3975 6469 L +3857 6206 L +3741 5936 L +3627 5661 L +4928 8646 M +-96 -109 V +4732 8409 L +4629 8262 L +4523 8097 L +4416 7914 L +4307 7716 L +4197 7502 L +4087 7274 L +3978 7034 L +3869 6783 L +3763 6522 L +3659 6252 L +3558 5976 L +-97 -281 V +4667 8730 M +-99 -108 V +4466 8493 L +4363 8344 L +4259 8176 L +4154 7990 L +4050 7788 L +3947 7570 L +3845 7338 L +3745 7093 L +-97 -257 V +-94 -266 V +-89 -275 V +-85 -281 V +-81 -285 V +4404 8794 M +4302 8685 L +4199 8556 L +4096 8407 L +3994 8238 L +3894 8050 L +-99 -204 V +-97 -220 V +-93 -235 V +-90 -249 V +-85 -259 V +-81 -270 V +-75 -278 V +-70 -285 V +-63 -289 V +998 3075 R +4035 8728 L +3932 8600 L +3830 8450 L +-99 -168 V +-97 -188 V +-93 -204 V +-88 -221 V +-85 -236 V +-79 -249 V +-74 -261 V +-68 -272 V +-61 -280 V +-54 -287 V +-48 -292 V +890 3064 R +3769 8750 L +3666 8622 L +3566 8475 L +-96 -168 V +-92 -185 V +-86 -203 V +-82 -220 V +-75 -234 V +-69 -248 V +-63 -261 V +-55 -271 V +-48 -281 V +-39 -288 V +-32 -294 V +5520 5221 M +-261 65 V +-246 60 V +-231 56 V +-218 51 V +-206 47 V +-195 44 V +-186 41 V +-178 39 V +-172 37 V +-166 34 V +-162 34 V +-158 32 V +-157 31 V +-155 30 V +5728 5395 M +-258 81 V +-249 75 V +-239 69 V +-228 64 V +-219 58 V +-210 54 V +-201 50 V +-195 46 V +-188 44 V +-183 40 V +-178 38 V +-176 36 V +-172 34 V +-171 32 V +5913 5555 M +-253 98 V +-249 91 V +-242 84 V +-236 77 V +-228 70 V +-221 65 V +-215 60 V +-208 55 V +-204 51 V +-198 46 V +-194 43 V +-191 40 V +-188 36 V +-186 33 V +6076 5702 M +-246 115 V +-245 107 V +-242 99 V +-240 92 V +-235 83 V +-230 77 V +-226 71 V +-220 64 V +-217 59 V +-212 53 V +-209 48 V +-205 43 V +-202 38 V +-199 34 V +6221 5837 M +-237 132 V +-239 124 V +-241 115 V +-241 106 V +-239 97 V +-237 89 V +-234 82 V +-231 74 V +-228 67 V +-225 60 V +-221 53 V +-218 47 V +-215 40 V +-212 33 V +6349 5962 M +-226 149 V +-233 139 V +-237 130 V +-240 121 V +-242 112 V +-241 102 V +-241 93 V +-239 84 V +-237 76 V +-235 66 V +-233 59 V +-230 49 V +-226 42 V +-223 33 V +6465 6078 M +-217 164 V +-225 155 V +-233 146 V +-238 135 V +-241 125 V +-244 115 V +-245 105 V +-246 94 V +-245 84 V +-244 73 V +-242 64 V +-240 53 V +-237 42 V +-233 32 V +6568 6187 M +-207 178 V +-217 170 V +-227 160 V +-234 149 V +-241 139 V +-245 127 V +-248 116 V +-250 104 V +-251 92 V +-251 80 V +-250 68 V +-249 56 V +-245 43 V +-243 30 V +6662 6288 M +-198 193 V +-209 183 V +-221 174 V +-230 162 V +-238 151 V +-245 139 V +-249 126 V +-253 113 V +-256 100 V +-256 87 V +-257 72 V +-255 58 V +-254 44 V +-249 29 V +6746 6384 M +-188 205 V +-202 196 V +-215 186 V +-225 175 V +-235 162 V +-243 150 V +-249 136 V +-255 121 V +-258 107 V +-260 92 V +-262 76 V +-260 60 V +-260 44 V +-256 28 V +6823 6476 M +-180 216 V +-195 208 V +-208 197 V +-221 185 V +-231 173 V +-240 158 V +-248 145 V +-255 129 V +-259 113 V +-263 97 V +-264 79 V +-265 62 V +-263 44 V +-261 25 V +6893 6563 M +-173 227 V +-188 217 V +-202 207 V +-215 195 V +-227 181 V +-237 167 V +-246 151 V +-254 135 V +-259 119 V +currentpoint stroke M +-263 100 V +-266 82 V +-267 63 V +-266 43 V +-264 25 V +6956 6646 M +-166 236 V +-181 227 V +-197 215 V +-210 202 V +-222 188 V +-234 173 V +-243 157 V +-251 140 V +-258 122 V +-262 103 V +-266 84 V +-267 63 V +-267 44 V +-266 22 V +7013 6727 M +-159 244 V +-176 234 V +-192 221 V +-205 208 V +-218 194 V +-229 177 V +-240 161 V +-247 143 V +-255 124 V +-260 104 V +-264 85 V +-266 63 V +-267 43 V +-266 22 V +7065 6805 M +-155 251 V +-171 239 V +-186 226 V +-201 212 V +-213 197 V +-225 180 V +-235 162 V +-244 144 V +-251 125 V +-256 105 V +-261 84 V +-263 64 V +-265 42 V +-265 20 V +1.000 UL +LT1 +2396 3123 M +251 -73 V +260 -59 V +268 -46 V +274 -32 V +279 -18 V +282 -4 V +284 11 V +283 25 V +281 39 V +278 53 V +273 67 V +265 81 V +258 93 V +246 106 V +-3799 -3 R +253 -72 V +263 -60 V +270 -47 V +277 -34 V +282 -22 V +284 -8 V +286 5 V +286 18 V +283 30 V +279 44 V +273 57 V +266 68 V +257 81 V +247 91 V +-3811 99 R +253 -70 V +264 -59 V +271 -48 V +278 -36 V +283 -26 V +285 -13 V +287 -2 V +286 9 V +283 22 V +280 33 V +273 44 V +267 56 V +257 66 V +246 76 V +2379 3873 M +253 -67 V +263 -58 V +270 -48 V +277 -39 V +282 -30 V +284 -19 V +286 -9 V +286 1 V +283 11 V +279 21 V +273 31 V +266 41 V +257 51 V +247 60 V +2396 4142 M +251 -64 V +260 -56 V +268 -49 V +274 -42 V +279 -33 V +282 -25 V +284 -17 V +283 -8 V +281 0 V +278 9 V +273 17 V +265 26 V +258 35 V +246 42 V +2423 4418 M +248 -60 V +256 -55 V +264 -49 V +270 -44 V +274 -37 V +278 -31 V +279 -25 V +280 -17 V +279 -11 V +275 -4 V +272 3 V +265 11 V +257 17 V +248 25 V +2461 4699 M +243 -57 V +251 -52 V +258 -49 V +263 -45 V +269 -41 V +271 -37 V +274 -31 V +275 -27 V +275 -22 V +273 -17 V +269 -11 V +264 -5 V +257 0 V +249 7 V +2510 4983 M +236 -52 V +244 -50 V +251 -49 V +256 -46 V +260 -44 V +264 -41 V +267 -39 V +269 -35 V +269 -33 V +269 -29 V +266 -25 V +263 -21 V +257 -17 V +250 -12 V +2569 5270 M +229 -48 V +235 -48 V +242 -47 V +246 -47 V +251 -47 V +255 -45 V +258 -45 V +261 -43 V +263 -42 V +263 -40 V +263 -39 V +261 -37 V +257 -33 V +250 -31 V +2638 5557 M +220 -44 V +226 -45 V +230 -46 V +236 -47 V +240 -48 V +244 -48 V +248 -50 V +252 -50 V +255 -50 V +256 -51 V +258 -51 V +258 -51 V +256 -50 V +252 -50 V +2716 5843 M +210 -40 V +215 -42 V +218 -45 V +223 -46 V +228 -48 V +232 -51 V +236 -53 V +241 -55 V +245 -57 V +249 -60 V +252 -62 V +255 -63 V +254 -66 V +253 -68 V +2804 6126 M +198 -36 V +202 -40 V +206 -42 V +209 -45 V +213 -48 V +219 -51 V +223 -55 V +229 -58 V +234 -63 V +239 -66 V +245 -71 V +249 -75 V +252 -79 V +254 -84 V +2900 6404 M +186 -33 V +188 -36 V +191 -40 V +194 -43 V +198 -46 V +204 -51 V +208 -55 V +215 -60 V +221 -65 V +228 -70 V +236 -77 V +242 -84 V +249 -91 V +253 -98 V +3005 6676 M +172 -30 V +172 -33 V +175 -37 V +179 -40 V +182 -44 V +186 -49 V +193 -54 V +199 -59 V +207 -65 V +215 -72 V +224 -80 V +234 -90 V +242 -99 V +252 -109 V +3117 6939 M +156 -26 V +157 -30 V +158 -33 V +161 -37 V +165 -41 V +169 -46 V +175 -51 V +182 -56 V +190 -63 V +200 -71 V +211 -81 V +222 -90 V +235 -104 V +247 -117 V +6178 3366 M +7 148 V +2 151 V +-2 154 V +-7 158 V +-10 164 V +-16 171 V +-21 178 V +-28 188 V +-34 198 V +-42 211 V +-51 226 V +-63 242 V +-76 260 V +-92 278 V +5932 3260 M +6 163 V +3 166 V +-3 170 V +-6 176 V +-12 181 V +-17 189 V +-22 197 V +-28 207 V +-36 217 V +-43 229 V +-52 242 V +-62 256 V +-75 271 V +-87 286 V +5674 3167 M +7 175 V +3 181 V +-3 185 V +-7 192 V +-11 199 V +-17 206 V +-22 214 V +-28 223 V +-35 234 V +-41 245 V +-50 255 V +-59 268 V +-68 279 V +-80 291 V +5409 3086 M +6 188 V +2 193 V +-2 200 V +-6 207 V +-11 214 V +-16 222 V +-21 230 V +-26 239 V +-32 248 V +-38 257 V +-44 267 V +-52 277 V +-60 285 V +-68 291 V +95 -3385 R +6 198 V +2 206 V +-2 213 V +-6 221 V +-10 228 V +-13 236 V +-18 244 V +-23 253 V +-27 260 V +-32 268 V +-37 276 V +-43 283 V +-48 288 V +-55 292 V +28 -3519 R +5 207 V +1 217 V +-1 225 V +-5 233 V +-7 241 V +-11 249 V +-14 256 V +-17 264 V +-20 271 V +-25 277 V +-27 282 V +-32 287 V +-35 290 V +-40 291 V +4577 2927 M +3 216 V +1 225 V +-1 236 V +-3 244 V +-5 252 V +-7 260 V +-8 267 V +-11 273 V +-12 279 V +-15 284 V +-16 288 V +-19 289 V +-21 290 V +-23 289 V +4294 2902 M +0 223 V +1 234 V +-1 244 V +0 253 V +-2 261 V +-2 270 V +-2 275 V +-3 281 V +-3 286 V +-4 289 V +-4 291 V +-5 291 V +-5 289 V +-6 286 V +4010 2891 M +-2 229 V +0 241 V +0 251 V +2 261 V +3 269 V +3 276 V +5 283 V +6 287 V +7 291 V +8 292 V +9 293 V +10 291 V +10 288 V +12 283 V +3728 2895 M +-4 233 V +-1 246 V +1 257 V +4 267 V +7 275 V +10 282 V +12 287 V +15 291 V +18 294 V +20 295 V +22 293 V +25 291 V +28 286 V +29 280 V +3449 2913 M +-7 237 V +-2 250 V +2 261 V +7 270 V +12 279 V +15 286 V +21 290 V +24 294 V +29 295 V +32 295 V +37 293 V +40 289 V +44 284 V +46 277 V +3175 2945 M +-10 239 V +-3 252 V +3 264 V +10 273 V +16 281 V +22 287 V +28 291 V +34 295 V +39 295 V +currentpoint stroke M +45 294 V +51 292 V +55 287 V +59 281 V +64 274 V +2907 2991 M +-12 240 V +-4 253 V +4 264 V +12 274 V +20 281 V +28 287 V +35 291 V +43 293 V +51 294 V +57 293 V +63 289 V +70 285 V +75 278 V +81 270 V +2647 3050 M +-15 241 V +-5 252 V +5 263 V +15 272 V +24 280 V +33 284 V +42 289 V +52 291 V +60 291 V +68 290 V +76 287 V +84 281 V +91 275 V +96 267 V +2396 3123 M +-17 240 V +-5 250 V +5 260 V +17 269 V +27 276 V +38 281 V +49 284 V +59 287 V +69 287 V +78 286 V +88 283 V +96 278 V +105 272 V +112 263 V +1.000 UL +LT2 +5140 3125 M +271 14 V +263 28 V +254 40 V +243 52 V +230 65 V +216 76 V +202 88 V +185 98 V +168 108 V +150 117 V +131 127 V +112 134 V +91 141 V +71 148 V +5143 3273 M +273 28 V +265 41 V +257 53 V +245 64 V +233 75 V +219 86 V +205 95 V +188 106 V +170 114 V +153 122 V +132 129 V +113 136 V +92 142 V +71 147 V +5144 3423 M +273 44 V +267 56 V +257 66 V +246 76 V +234 86 V +220 96 V +206 104 V +189 112 V +171 120 V +153 126 V +134 133 V +113 137 V +92 141 V +70 145 V +5143 3577 M +273 61 V +265 70 V +257 80 V +245 90 V +233 97 V +219 106 V +205 112 V +188 120 V +170 125 V +153 130 V +132 134 V +113 138 V +92 140 V +71 142 V +5140 3736 M +271 78 V +263 86 V +254 95 V +243 102 V +230 109 V +216 115 V +202 121 V +185 126 V +168 129 V +150 133 V +131 136 V +112 137 V +91 138 V +71 138 V +5134 3901 M +269 95 V +260 103 V +249 109 V +239 115 V +225 120 V +211 124 V +197 129 V +181 131 V +163 133 V +147 135 V +128 136 V +110 136 V +90 135 V +70 134 V +5127 4074 M +264 112 V +255 119 V +244 123 V +232 127 V +218 131 V +205 133 V +190 135 V +174 136 V +158 136 V +142 136 V +124 135 V +107 134 V +89 131 V +70 129 V +5117 4255 M +259 130 V +248 134 V +236 137 V +224 139 V +210 140 V +196 140 V +182 140 V +166 140 V +151 137 V +136 136 V +119 133 V +103 131 V +87 127 V +70 123 V +5104 4446 M +252 147 V +240 149 V +227 150 V +214 149 V +200 148 V +186 147 V +172 143 V +157 141 V +142 138 V +129 134 V +113 131 V +99 126 V +84 122 V +69 118 V +5088 4648 M +243 165 V +230 163 V +217 161 V +202 158 V +188 154 V +175 150 V +160 146 V +147 141 V +133 136 V +120 131 V +107 127 V +95 121 V +81 117 V +67 111 V +5068 4865 M +233 180 V +219 176 V +203 170 V +190 164 V +175 158 V +161 152 V +148 145 V +136 139 V +123 133 V +111 127 V +100 121 V +89 116 V +77 110 V +67 105 V +5043 5097 M +221 194 V +206 185 V +190 177 V +175 168 V +160 159 V +148 151 V +135 143 V +123 135 V +112 128 V +102 121 V +92 115 V +83 109 V +74 104 V +64 97 V +5013 5346 M +208 205 V +190 193 V +174 180 V +160 169 V +145 157 V +133 147 V +121 138 V +111 129 V +101 121 V +92 115 V +84 107 V +77 102 V +69 96 V +61 90 V +4977 5615 M +192 213 V +174 195 V +157 180 V +142 166 V +130 152 V +117 140 V +107 131 V +98 121 V +90 113 V +83 105 V +75 99 V +70 94 V +64 88 V +58 83 V +4933 5905 M +174 215 V +156 194 V +139 174 V +125 158 V +113 144 V +102 131 V +94 121 V +85 111 V +79 103 V +72 96 V +67 90 V +63 84 V +58 80 V +55 75 V +7727 4361 M +32 250 V +10 254 V +-10 257 V +-32 257 V +-54 257 V +-74 255 V +-95 251 V +-116 247 V +-135 240 V +-153 232 V +-172 222 V +-189 212 V +-205 200 V +-219 186 V +7656 4213 M +32 251 V +11 256 V +-11 260 V +-32 261 V +-53 261 V +-74 260 V +-95 257 V +-115 252 V +-133 247 V +-153 238 V +-169 230 V +-186 219 V +-202 207 V +-216 194 V +7565 4072 M +31 250 V +11 257 V +-11 261 V +-31 263 V +-52 264 V +-73 264 V +-93 261 V +-112 257 V +-130 252 V +-149 245 V +-166 236 V +-181 227 V +-197 215 V +-210 202 V +7453 3938 M +30 248 V +11 256 V +-11 260 V +-30 264 V +-50 265 V +-70 266 V +-89 264 V +-108 262 V +-126 257 V +-143 250 V +-160 243 V +-175 234 V +-190 223 V +-203 212 V +7322 3811 M +29 246 V +9 252 V +-9 259 V +-29 262 V +-47 265 V +-66 267 V +-84 266 V +-102 264 V +-120 261 V +-136 256 V +-152 249 V +-167 242 V +-181 231 V +-195 221 V +7172 3694 M +26 241 V +9 248 V +-9 255 V +-26 259 V +-44 263 V +-61 266 V +-78 266 V +-95 266 V +-111 264 V +-127 260 V +-143 255 V +-157 248 V +-172 241 V +-184 230 V +7004 3586 M +24 235 V +8 242 V +-8 250 V +-24 255 V +-39 259 V +-56 263 V +-71 265 V +-86 265 V +-102 266 V +-117 263 V +-132 260 V +-146 255 V +-161 249 V +-173 240 V +6819 3488 M +21 227 V +7 236 V +-7 242 V +-21 249 V +-35 254 V +-49 258 V +-63 261 V +-77 264 V +-92 266 V +-106 265 V +-119 264 V +-134 261 V +-148 257 V +-160 250 V +6617 3400 M +18 220 V +6 227 V +-6 234 V +-18 240 V +-30 246 V +-42 252 V +-55 256 V +-67 261 V +-80 263 V +-94 266 V +-106 266 V +-120 266 V +-134 264 V +-147 260 V +6401 3324 M +15 210 V +5 217 V +-5 224 V +-15 231 V +-25 237 V +-36 243 V +-46 249 V +-57 254 V +-69 260 V +-80 264 V +-93 267 V +-105 270 V +-118 271 V +-132 269 V +6171 3259 M +12 200 V +4 206 V +-4 213 V +-12 219 V +-20 226 V +-29 232 V +-38 240 V +-47 246 V +-57 254 V +-67 260 V +-78 266 V +-90 272 V +-103 276 V +-115 277 V +5928 3207 M +10 188 V +3 194 V +-3 199 V +-10 207 V +-16 213 V +-22 220 V +-30 228 V +-37 236 V +-45 245 V +currentpoint stroke M +-55 254 V +-63 262 V +-75 271 V +-85 279 V +-98 285 V +5674 3167 M +7 175 V +3 181 V +-3 185 V +-7 192 V +-11 199 V +-17 206 V +-22 214 V +-28 223 V +-35 234 V +-41 245 V +-50 255 V +-59 268 V +-68 279 V +-80 291 V +5411 3139 M +5 162 V +1 166 V +-1 171 V +-5 176 V +-8 182 V +-12 190 V +-15 199 V +-20 208 V +-25 220 V +-30 232 V +-37 246 V +-43 260 V +-52 277 V +-62 292 V +33 -2995 R +3 148 V +1 150 V +-1 154 V +-3 159 V +-6 165 V +-7 173 V +-10 181 V +-13 191 V +-16 202 V +-20 217 V +-25 232 V +-30 249 V +-36 269 V +-44 290 V +1.000 UL +LTb +4199 8556 M +267 -63 V +266 -84 V +262 -103 V +258 -122 V +251 -140 V +243 -157 V +234 -173 V +222 -188 V +210 -202 V +197 -215 V +181 -227 V +166 -236 V +1.000 UL +LTb +4199 8556 M +4096 8407 L +3994 8238 L +3894 8050 L +-99 -204 V +-97 -220 V +-93 -235 V +-90 -249 V +-85 -259 V +-81 -270 V +-75 -278 V +-70 -285 V +-63 -289 V +1.000 UL +LTb +3430 6883 M +-81 -270 V +-75 -278 V +-70 -285 V +-63 -289 V +-57 -293 V +-51 -294 V +-43 -293 V +-35 -291 V +-28 -287 V +-20 -281 V +-12 -274 V +-4 -264 V +1.000 UL +LTb +2891 3484 M +271 -48 V +278 -36 V +283 -26 V +285 -13 V +287 -2 V +286 9 V +283 22 V +280 33 V +273 44 V +267 56 V +257 66 V +246 76 V +1.000 UL +LTb +6202 7526 M +210 -202 V +197 -215 V +181 -227 V +166 -236 V +149 -245 V +130 -252 V +112 -257 V +93 -261 V +73 -264 V +52 -264 V +31 -263 V +11 -261 V +1.000 UL +LTb +5144 3423 M +273 44 V +267 56 V +257 66 V +246 76 V +234 86 V +220 96 V +206 104 V +189 112 V +171 120 V +153 126 V +134 133 V +113 137 V +1.000 UP +stroke +grestore +end +showpage +%%Trailer +%%DocumentFonts: Helvetica diff --git a/doc/documentation.tex b/doc/documentation.tex index c2a888f..68fcbe2 100644 --- a/doc/documentation.tex +++ b/doc/documentation.tex @@ -201,11 +201,6 @@ equations for $h$ at the angular grid points, and solve this system of equations by a global Newton's method (or a variant with improved convergence). -I use a multiple grid patch system to cover $S^2$ without coordinate -singularities. In general there are 6~patches, one each covering a -neighborhood of the $\pm z$, $\pm x$, and $\pm y$ axes, but if the slice -has suitable symmetries, this may be reduced to 5, 4, or 3~patches. - Computationally, this algorithm has 3 main parts: \begin{itemize} \item Computation of the ``horizon function'' $H(h)$ given a trial @@ -233,6 +228,70 @@ Computationally, this algorithm has 3 main parts: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +\section{The Multipatch System} + +Perhaps the most unusual feature of \thorn{AHFinderDirect} is the +``multipatch'' system used to cover $S^2$ without coordinate singularities. +In general there are 6~patches, one each covering a neighborhood of +the $\pm z$, $\pm x$, and $\pm y$ axes, but this may be reduced in +the presence of suitable symmetries. For example, +figure~\ref{AEIDevelopment_AHFinderDirect_fig-3patch} +on page~\pageref{AEIDevelopment_AHFinderDirect_fig-3patch} shows a system +of 3~patches covering the $+xyz$~octant of $S^2$. This would be +suitable for finding an apparent horizon with mirror symmetry about +the (local) $z=0$~plane, and either 90~degree periodic rotation symmetry +about the (local) $z$~axis, or mirror symmetry about each of the (local) +$x$~and $y$~axes. + +\begin{figure}[htbp] +\begin{center} +\includegraphics{AEIDevelopment_AHFinderDirect_3patch.eps} +\end{center} +\caption[Illustration of the Multipatch System] + { + This figure shows a multipatch system covering the + $(+,+,+)$~octant of the unit sphere~$S^2$ with 3~patches. + The angular resolution is 5~degrees. Notice that the + patches overlap by several ``ghost zone'' grid points. + } +\label{AEIDevelopment_AHFinderDirect_fig-3patch} +\end{figure} + +To allow easy angular finite differencing within the patch system, +each patch is extended beyond its nominal extent by a ``ghost zone''%%% +\footnote{%%% + Note that this terminology differs somewhat + from that used by Cactus in general; Cactus + would call these ``patch zones'' or ``symmetry + zones''. + }%%% +{} (2~grid points wide in +figure~\ref{AEIDevelopment_AHFinderDirect_fig-3patch}). Angular grid +function values in the ghost zone can be obtained by interpatch +interpolation%%% +\footnote{%%% + Due to the way the patch coordinates are defined, + adjacent patches always share a common ``perpendicular'' + angular coordinate, so only 1-D interpolation + is needed here. + }%%% +{} or by applying symmetry operations. Once this is done, then angular +finite differencing within the nominal extent of each patch can proceed +normally, ignoring the patch boundaries. \thorn{AHFinderDirect} can +be configured at compile time to use either 2nd~order or 4th~order +angular finite differencing (3~point or 5~point angular molecules); +the default is 4th~order (5~point). See \verb|src/include/config.h| +if you want to change this. + +By default \thorn{AHFinderDirect} will automagically choose a patch +system type for each apparent horizon searched for, based on the local +coordinate origin and the symmetries implicit in the Cactus grid type. +This generally works well, but if desired you can instead manually +specify the patch system type, the angular resolution, the width of +the ghost zones, etc. See the \verb|param.ccl| file for details. + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + \section{Using \thorn{AHFinderDirect}} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -664,7 +723,7 @@ parameters controlling if and how these should be written to data files: apparent horizon number \item the second \verb|%s| is the file name extension, which defaults to - \verb|gnuplot| (see the \verb|param.ccl| + \verb|"gp"| (see the \verb|param.ccl| file if you want to change this) \end{itemize} \item Comment lines begin with \verb|#|. @@ -711,9 +770,9 @@ parameters controlling if and how these should be written to data files: described above. \end{description} -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\subsection{Accuracy and Diagnostics} +\section{Accuracy and Diagnostics} The apparent horizon positions are typically computed very accurately; tests on Kerr spacetimes give typical errors of $10^{-4}m$ to $10^{-5}m$. @@ -754,8 +813,8 @@ tests. The \verb|Kerr-tiny.par| parameter file is close to a minimal # to the black hole, to make this a non-trivial test for the apparent # horizon finder. # -# This parameter is "tiny" in the sense that it sets only a small number -# of AHFinderDirect parameters. +# This parameter file is "tiny" in the sense that it sets only a +# small number of AHFinderDirect parameters. # next two lines is actually one long line # (Cactus doesn't seem to grok \-newline continuation here :( :( ) diff --git a/doc/makefile b/doc/makefile index 4592087..ab8680c 100644 --- a/doc/makefile +++ b/doc/makefile @@ -1,12 +1,18 @@ # Makefile for documentation -# $Header: /usr/local/svn/cvs-repositories/numrelcvs/AEIThorns/AHFinderDirect/doc/makefile,v 1.2 2002-09-16 18:08:06 jthorn Exp $ +# $Header: /usr/local/svn/cvs-repositories/numrelcvs/AEIThorns/AHFinderDirect/doc/makefile,v 1.3 2002-10-17 22:27:44 jthorn Exp $ # -# Targets: +# ***** main targets ***** # default latex & dvips everything # bib bibtex everything # clean remove all machine-generated files # +# ***** targets for rebuilding figures in the thorn guide ***** +# 3patch-cactus run Cactus to generate the data file for the +# 3patch figure +# 3patch-gnuplot run gnuplot to generate the postscript file +# for the 3patch figure +# # # Note this Makefile always re-TeXs things, even if the .tex files @@ -17,6 +23,10 @@ ############################################################################### +# +# ***** main targets ***** +# + .PHONY : default default : ps @@ -35,3 +45,17 @@ bib : .PHONY : clean clean : -rm *.dvi *.ps *.aux *.bbl *.log *.blg + +############################################################################### + +# +# ***** targets for rebuilding figures in the thorn guide ***** +# + +.PHONY : 3patch-cactus +3patch-cactus : + cactus_test-ahfinderdirect 3patch.par 2>&1 | tee 3patch.log + +.PHONY : 3patch-gnuplot +3patch-gnuplot : + gnuplot 3patch.gnuplot |