From e7220fadb6b3775e55afd9a07831e90136cfaf24 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Tue, 7 Nov 2017 10:43:35 -0500 Subject: many changes, including added references --- .gitignore | 3 +- fig_correlation.tex | 3547 +++++++++++++++++++++++++++++ fig_correlation_collapse-hL.tex | 1878 +++++++++++++++ figs/fig_correlation_data_collapse-hL.dat | 106 +- makefile | 2 +- monte-carlo.bib | 153 ++ monte-carlo.tex | 48 +- refs/alexandrowicz1989swendsen.pdf | Bin 0 -> 937546 bytes refs/coniglio1980clusters.pdf | Bin 0 -> 442501 bytes refs/coniglio1989exact.pdf | Bin 0 -> 463325 bytes refs/destri1992swendsen.pdf | Bin 0 -> 438528 bytes refs/dotsenko1991cluster.pdf | Bin 0 -> 172731 bytes refs/du2006dynamic.pdf | Bin 0 -> 602685 bytes refs/janke1998nonlocal.pdf | Bin 0 -> 269280 bytes refs/lauwers1989critical.pdf | Bin 0 -> 195270 bytes refs/liu2014dynamic.pdf | Bin 0 -> 2182262 bytes refs/ray1990metastability.pdf | Bin 0 -> 432932 bytes refs/swendsen1987nonuniversal.pdf | Bin 0 -> 136908 bytes refs/wang1989clusters.pdf | Bin 0 -> 929783 bytes refs/wang1990cluster.pdf | Bin 0 -> 652558 bytes refs/wolff1989collective.pdf | Bin 0 -> 201013 bytes refs/wolff1989comparison.pdf | Bin 0 -> 279147 bytes refs/wolff1990critical.pdf | Bin 0 -> 945386 bytes 23 files changed, 5674 insertions(+), 63 deletions(-) create mode 100644 fig_correlation.tex create mode 100644 fig_correlation_collapse-hL.tex create mode 100644 refs/alexandrowicz1989swendsen.pdf create mode 100644 refs/coniglio1980clusters.pdf create mode 100644 refs/coniglio1989exact.pdf create mode 100644 refs/destri1992swendsen.pdf create mode 100644 refs/dotsenko1991cluster.pdf create mode 100644 refs/du2006dynamic.pdf create mode 100644 refs/janke1998nonlocal.pdf create mode 100644 refs/lauwers1989critical.pdf create mode 100644 refs/liu2014dynamic.pdf create mode 100644 refs/ray1990metastability.pdf create mode 100644 refs/swendsen1987nonuniversal.pdf create mode 100644 refs/wang1989clusters.pdf create mode 100644 refs/wang1990cluster.pdf create mode 100644 refs/wolff1989collective.pdf create mode 100644 refs/wolff1989comparison.pdf create mode 100644 refs/wolff1990critical.pdf diff --git a/.gitignore b/.gitignore index 7156c49..9325c3c 100644 --- a/.gitignore +++ b/.gitignore @@ -4,6 +4,7 @@ *.blg *Notes.bib *.dvi -*.pdf +monte-carlo.pdf *.synctex.gz +fig-*.tex diff --git a/fig_correlation.tex b/fig_correlation.tex new file mode 100644 index 0000000..eb09910 --- /dev/null +++ b/fig_correlation.tex @@ -0,0 +1,3547 @@ +% GNUPLOT: LaTeX picture with Postscript +\begingroup% +\makeatletter% +\newcommand{\GNUPLOTspecial}{% + \@sanitize\catcode`\%=14\relax\special}% +\setlength{\unitlength}{0.0500bp}% +\begin{picture}(4920,9118)(0,0)% + {\GNUPLOTspecial{" +%!PS-Adobe-2.0 EPSF-2.0 +%%Creator: gnuplot 5.2 patchlevel 0 +%%CreationDate: Mon Nov 6 21:12:43 2017 +%%DocumentFonts: +%%BoundingBox: 0 0 246 455 +%%EndComments +%%BeginProlog +/gnudict 256 dict def +gnudict begin +% +% The following true/false flags may be edited by hand if desired. +% The unit line width and grayscale image gamma correction may also be changed. +% +/Color false def +/Blacktext true def +/Solid false def +/Dashlength 1 def +/Landscape false def +/Level1 false def +/Level3 false def +/Rounded false def +/ClipToBoundingBox false def +/SuppressPDFMark false def +/TransparentPatterns false def +/gnulinewidth 5.000 def +/userlinewidth gnulinewidth def +/Gamma 1.0 def +/BackgroundColor {-1.000 -1.000 -1.000} def +% +/vshift -66 def +/dl1 { + 10.0 Dashlength userlinewidth gnulinewidth div mul mul mul + Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if +} def +/dl2 { + 10.0 Dashlength userlinewidth gnulinewidth div mul mul mul + Rounded { currentlinewidth 0.75 mul add } if +} def +/hpt_ 31.5 def +/vpt_ 31.5 def +/hpt hpt_ def +/vpt vpt_ def +/doclip { + ClipToBoundingBox { + newpath 0 0 moveto 246 0 lineto 246 455 lineto 0 455 lineto closepath + clip + } if +} def +% +% Gnuplot Prolog Version 5.1 (Oct 2015) +% +%/SuppressPDFMark true def +% +/M {moveto} bind def +/L {lineto} bind def +/R {rmoveto} bind def +/V {rlineto} bind def +/N {newpath moveto} bind def +/Z {closepath} bind def +/C {setrgbcolor} bind def +/f {rlineto fill} bind def +/g {setgray} bind def +/Gshow {show} def % May be redefined later in the file to support UTF-8 +/vpt2 vpt 2 mul def +/hpt2 hpt 2 mul def +/Lshow {currentpoint stroke M 0 vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def +/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def +/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} 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 0 setgray Solid {pop []} if 0 setdash} ifelse} def +/BL {stroke userlinewidth 2 mul setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +/AL {stroke userlinewidth 2 div setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +/UL {dup gnulinewidth mul /userlinewidth exch def + dup 1 lt {pop 1} if 10 mul /udl exch def} def +/PL {stroke userlinewidth setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +3.8 setmiterlimit +% Classic Line colors (version 5.0) +/LCw {1 1 1} def +/LCb {0 0 0} def +/LCa {0 0 0} def +/LC0 {1 0 0} def +/LC1 {0 1 0} def +/LC2 {0 0 1} def +/LC3 {1 0 1} def +/LC4 {0 1 1} def +/LC5 {1 1 0} def +/LC6 {0 0 0} def +/LC7 {1 0.3 0} def +/LC8 {0.5 0.5 0.5} def +% Default dash patterns (version 5.0) +/LTB {BL [] LCb DL} def +/LTw {PL [] 1 setgray} def +/LTb {PL [] LCb DL} def +/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def +/LT0 {PL [] LC0 DL} def +/LT1 {PL [2 dl1 3 dl2] LC1 DL} def +/LT2 {PL [1 dl1 1.5 dl2] LC2 DL} def +/LT3 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC3 DL} def +/LT4 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def +/LT5 {PL [4 dl1 2 dl2] LC5 DL} def +/LT6 {PL [1.5 dl1 1.5 dl2 1.5 dl1 1.5 dl2 1.5 dl1 6 dl2] LC6 DL} def +/LT7 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC7 DL} def +/LT8 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC8 DL} def +/SL {[] 0 setdash} 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 +/Density { + /Fillden exch def + currentrgbcolor + /ColB exch def /ColG exch def /ColR exch def + /ColR ColR Fillden mul Fillden sub 1 add def + /ColG ColG Fillden mul Fillden sub 1 add def + /ColB ColB Fillden mul Fillden sub 1 add def + ColR ColG ColB setrgbcolor} def +/BoxColFill {gsave Rec PolyFill} def +/PolyFill {gsave Density fill grestore grestore} def +/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def +% +% PostScript Level 1 Pattern Fill routine for rectangles +% Usage: x y w h s a XX PatternFill +% x,y = lower left corner of box to be filled +% w,h = width and height of box +% a = angle in degrees between lines and x-axis +% XX = 0/1 for no/yes cross-hatch +% +/PatternFill {gsave /PFa [ 9 2 roll ] def + PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate + PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec + TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse + clip + currentlinewidth 0.5 mul setlinewidth + /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def + 0 0 M PFa 5 get rotate PFs -2 div dup translate + 0 1 PFs PFa 4 get div 1 add floor cvi + {PFa 4 get mul 0 M 0 PFs V} for + 0 PFa 6 get ne { + 0 1 PFs PFa 4 get div 1 add floor cvi + {PFa 4 get mul 0 2 1 roll M PFs 0 V} for + } if + stroke grestore} def +% +/languagelevel where + {pop languagelevel} {1} ifelse +dup 2 lt + {/InterpretLevel1 true def + /InterpretLevel3 false def} + {/InterpretLevel1 Level1 def + 2 gt + {/InterpretLevel3 Level3 def} + {/InterpretLevel3 false def} + ifelse } + ifelse +% +% PostScript level 2 pattern fill definitions +% +/Level2PatternFill { +/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8} + bind def +/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} +>> matrix makepattern +/Pat1 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke + 0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke} +>> matrix makepattern +/Pat2 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L + 8 8 L 8 0 L 0 0 L fill} +>> matrix makepattern +/Pat3 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L + 0 12 M 12 0 L stroke} +>> matrix makepattern +/Pat4 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L + 0 -4 M 12 8 L stroke} +>> matrix makepattern +/Pat5 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L + 0 12 M 8 -4 L 4 12 M 10 0 L stroke} +>> matrix makepattern +/Pat6 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L + 0 -4 M 8 12 L 4 -4 M 10 8 L stroke} +>> matrix makepattern +/Pat7 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L + 12 0 M -4 8 L 12 4 M 0 10 L stroke} +>> matrix makepattern +/Pat8 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L + -4 0 M 12 8 L -4 4 M 8 10 L stroke} +>> matrix makepattern +/Pat9 exch def +/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def +/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def +/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def +/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def +/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def +/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def +/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def +} def +% +% +%End of PostScript Level 2 code +% +/PatternBgnd { + TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse +} def +% +% Substitute for Level 2 pattern fill codes with +% grayscale if Level 2 support is not selected. +% +/Level1PatternFill { +/Pattern1 {0.250 Density} bind def +/Pattern2 {0.500 Density} bind def +/Pattern3 {0.750 Density} bind def +/Pattern4 {0.125 Density} bind def +/Pattern5 {0.375 Density} bind def +/Pattern6 {0.625 Density} bind def +/Pattern7 {0.875 Density} bind def +} def +% +% Now test for support of Level 2 code +% +Level1 {Level1PatternFill} {Level2PatternFill} ifelse +% +/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 +% +Level1 SuppressPDFMark or +{} { +/SDict 10 dict def +systemdict /pdfmark known not { + userdict /pdfmark systemdict /cleartomark get put +} if +SDict begin [ + /Title () + /Subject (gnuplot plot) + /Creator (gnuplot 5.2 patchlevel 0) +% /Producer (gnuplot) +% /Keywords () + /CreationDate (Mon Nov 6 21:12:43 2017) + /DOCINFO pdfmark +end +} ifelse +% +% Support for boxed text - Ethan A Merritt Sep 2016 +% +/InitTextBox { userdict /TBy2 3 -1 roll put userdict /TBx2 3 -1 roll put + userdict /TBy1 3 -1 roll put userdict /TBx1 3 -1 roll put + /Boxing true def } def +/ExtendTextBox { dup type /stringtype eq + { Boxing { gsave dup false charpath pathbbox + dup TBy2 gt {userdict /TBy2 3 -1 roll put} {pop} ifelse + dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse + dup TBy1 lt {userdict /TBy1 3 -1 roll put} {pop} ifelse + dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse + grestore } if } + {} ifelse} def +/PopTextBox { newpath TBx1 TBxmargin sub TBy1 TBymargin sub M + TBx1 TBxmargin sub TBy2 TBymargin add L + TBx2 TBxmargin add TBy2 TBymargin add L + TBx2 TBxmargin add TBy1 TBymargin sub L closepath } def +/DrawTextBox { PopTextBox stroke /Boxing false def} def +/FillTextBox { gsave PopTextBox fill grestore /Boxing false def} def +0 0 0 0 InitTextBox +/TBxmargin 20 def +/TBymargin 20 def +/Boxing false def +/textshow { ExtendTextBox Gshow } def +% +end +%%EndProlog +%%Page: 1 1 +gnudict begin +gsave +doclip +0 0 translate +0.050 0.050 scale +0 setgray +newpath +BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 4920.00 9118.00 BoxColFill} if +1.000 UL +LTb +LCb setrgbcolor +720 640 M +31 0 V +3928 0 R +-31 0 V +720 723 M +31 0 V +3928 0 R +-31 0 V +720 791 M +31 0 V +3928 0 R +-31 0 V +720 849 M +31 0 V +3928 0 R +-31 0 V +720 899 M +31 0 V +3928 0 R +-31 0 V +720 942 M +31 0 V +3928 0 R +-31 0 V +720 982 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 1240 M +31 0 V +3928 0 R +-31 0 V +720 1391 M +31 0 V +3928 0 R +-31 0 V +720 1499 M +31 0 V +3928 0 R +-31 0 V +720 1582 M +31 0 V +3928 0 R +-31 0 V +720 1650 M +31 0 V +3928 0 R +-31 0 V +720 1707 M +31 0 V +3928 0 R +-31 0 V +720 1757 M +31 0 V +3928 0 R +-31 0 V +720 1801 M +31 0 V +3928 0 R +-31 0 V +720 1841 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 2099 M +31 0 V +3928 0 R +-31 0 V +720 2250 M +31 0 V +3928 0 R +-31 0 V +720 2358 M +31 0 V +3928 0 R +-31 0 V +720 2441 M +31 0 V +3928 0 R +-31 0 V +720 2509 M +31 0 V +3928 0 R +-31 0 V +720 2566 M +31 0 V +3928 0 R +-31 0 V +720 2616 M +31 0 V +3928 0 R +-31 0 V +720 2660 M +31 0 V +3928 0 R +-31 0 V +720 2699 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 2958 M +31 0 V +3928 0 R +-31 0 V +720 3109 M +31 0 V +3928 0 R +-31 0 V +720 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +922 640 M +0 31 V +0 2438 R +0 -31 V +1040 640 M +0 31 V +0 2438 R +0 -31 V +1124 640 M +0 31 V +0 2438 R +0 -31 V +1189 640 M +0 31 V +0 2438 R +0 -31 V +1242 640 M +0 31 V +0 2438 R +0 -31 V +1287 640 M +0 31 V +0 2438 R +0 -31 V +1326 640 M +0 31 V +0 2438 R +0 -31 V +1360 640 M +0 31 V +0 2438 R +0 -31 V +1391 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1593 640 M +0 31 V +0 2438 R +0 -31 V +1711 640 M +0 31 V +0 2438 R +0 -31 V +1794 640 M +0 31 V +0 2438 R +0 -31 V +1859 640 M +0 31 V +0 2438 R +0 -31 V +1913 640 M +0 31 V +0 2438 R +0 -31 V +1957 640 M +0 31 V +0 2438 R +0 -31 V +1996 640 M +0 31 V +0 2438 R +0 -31 V +2031 640 M +0 31 V +0 2438 R +0 -31 V +2061 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2263 640 M +0 31 V +0 2438 R +0 -31 V +2381 640 M +0 31 V +0 2438 R +0 -31 V +2465 640 M +0 31 V +0 2438 R +0 -31 V +2530 640 M +0 31 V +0 2438 R +0 -31 V +2583 640 M +0 31 V +0 2438 R +0 -31 V +2628 640 M +0 31 V +0 2438 R +0 -31 V +2667 640 M +0 31 V +0 2438 R +0 -31 V +2701 640 M +0 31 V +0 2438 R +0 -31 V +2732 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2934 640 M +0 31 V +0 2438 R +0 -31 V +3052 640 M +0 31 V +0 2438 R +0 -31 V +3136 640 M +0 31 V +0 2438 R +0 -31 V +3201 640 M +0 31 V +0 2438 R +0 -31 V +3254 640 M +0 31 V +0 2438 R +0 -31 V +3299 640 M +0 31 V +0 2438 R +0 -31 V +3338 640 M +0 31 V +0 2438 R +0 -31 V +3372 640 M +0 31 V +0 2438 R +0 -31 V +3403 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +3605 640 M +0 31 V +0 2438 R +0 -31 V +3723 640 M +0 31 V +0 2438 R +0 -31 V +3806 640 M +0 31 V +0 2438 R +0 -31 V +3871 640 M +0 31 V +0 2438 R +0 -31 V +3925 640 M +0 31 V +0 2438 R +0 -31 V +3969 640 M +0 31 V +0 2438 R +0 -31 V +4008 640 M +0 31 V +0 2438 R +0 -31 V +4043 640 M +0 31 V +0 2438 R +0 -31 V +4073 640 M +0 63 V +0 2406 R +0 -63 V +stroke +LTb +LCb setrgbcolor +4275 640 M +0 31 V +0 2438 R +0 -31 V +4393 640 M +0 31 V +0 2438 R +0 -31 V +4477 640 M +0 31 V +0 2438 R +0 -31 V +4542 640 M +0 31 V +0 2438 R +0 -31 V +4595 640 M +0 31 V +0 2438 R +0 -31 V +4640 640 M +0 31 V +0 2438 R +0 -31 V +4679 640 M +0 31 V +0 2438 R +0 -31 V +stroke +LTb +LCb setrgbcolor +720 3109 N +720 640 L +3959 0 V +0 2469 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +LCb setrgbcolor +LTb +LCb setrgbcolor +LTb +1.000 UP +1.000 UL +LTb +0.37 0.51 0.71 C 720 1000 M +0 20 V +0 -20 R +31 0 V +-31 20 R +31 0 V +103 -2 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -31 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -23 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +104 -3 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +103 -7 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +103 -40 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 -6 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -13 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -26 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -16 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -17 R +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +104 -36 R +0 25 V +-31 -25 R +62 0 V +-62 25 R +62 0 V +103 -29 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 -7 R +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +103 -9 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -14 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 -16 R +0 19 V +-31 -19 R +62 0 V +stroke 3031 1000 M +-62 19 R +62 0 V +103 -33 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +104 -20 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +103 -27 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -14 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -10 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -34 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -7 R +0 29 V +-31 -29 R +62 0 V +-62 29 R +62 0 V +103 -25 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -25 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +104 -1 R +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +103 4 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +103 1 R +0 28 V +-31 -28 R +62 0 V +-62 28 R +62 0 V +720 1010 Pnt +854 1026 Pnt +988 1008 Pnt +1122 997 Pnt +1257 1002 Pnt +1391 1005 Pnt +1525 981 Pnt +1659 996 Pnt +1793 1004 Pnt +1927 996 Pnt +2061 996 Pnt +2195 1002 Pnt +2330 995 Pnt +2464 987 Pnt +2598 994 Pnt +2732 999 Pnt +2866 1005 Pnt +3000 1010 Pnt +3134 993 Pnt +3269 991 Pnt +3403 992 Pnt +3537 1001 Pnt +3671 1006 Pnt +3805 983 Pnt +3939 995 Pnt +4073 990 Pnt +4207 972 Pnt +4342 978 Pnt +4476 991 Pnt +4610 1009 Pnt +1.000 UP +1.000 UL +LTb +0.88 0.61 0.14 C 720 983 M +0 31 V +0 -31 R +31 0 V +-31 31 R +31 0 V +103 7 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -17 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -14 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +104 7 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 18 R +0 28 V +-31 -28 R +62 0 V +-62 28 R +62 0 V +103 32 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +103 28 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 115 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 126 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 228 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 309 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +104 465 R +0 22 V +-31 -22 R +62 0 V +-62 22 R +62 0 V +103 367 R +0 19 V +-31 -19 R +62 0 V +-62 19 R +62 0 V +103 -52 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 -177 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +103 -173 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -171 R +0 13 V +-31 -13 R +62 0 V +stroke 3031 2333 M +-62 13 R +62 0 V +103 -185 R +0 31 V +-31 -31 R +62 0 V +-62 31 R +62 0 V +104 -184 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +103 -188 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -195 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -126 R +0 19 V +-31 -19 R +62 0 V +-62 19 R +62 0 V +103 -146 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -71 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +103 -62 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +103 -77 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +104 -27 R +0 45 V +-31 -45 R +62 0 V +-62 45 R +62 0 V +103 -73 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -106 R +0 179 V +-31 -179 R +62 0 V +-62 179 R +62 0 V +720 999 Pnt +854 1023 Pnt +988 1014 Pnt +1122 1010 Pnt +1257 1023 Pnt +1391 1060 Pnt +1525 1109 Pnt +1659 1151 Pnt +1793 1285 Pnt +1927 1421 Pnt +2061 1660 Pnt +2195 1986 Pnt +2330 2469 Pnt +2464 2857 Pnt +2598 2823 Pnt +2732 2659 Pnt +2866 2497 Pnt +3000 2340 Pnt +3134 2177 Pnt +3269 2008 Pnt +3403 1839 Pnt +3537 1677 Pnt +3671 1578 Pnt +3805 1443 Pnt +3939 1374 Pnt +4073 1314 Pnt +4207 1248 Pnt +4342 1255 Pnt +4476 1216 Pnt +4610 1222 Pnt +1.000 UP +1.000 UL +LTb +0.56 0.69 0.20 C 720 1280 M +0 10 V +0 -10 R +31 0 V +-31 10 R +31 0 V +103 -11 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +103 4 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -29 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +104 -9 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -25 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -1 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -28 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -17 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 4 R +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +103 -18 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -26 R +0 30 V +-31 -30 R +62 0 V +-62 30 R +62 0 V +104 -20 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -9 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -18 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -8 R +0 30 V +-31 -30 R +62 0 V +-62 30 R +62 0 V +103 -22 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +103 -17 R +0 29 V +-31 -29 R +62 0 V +stroke 3031 1279 M +-62 29 R +62 0 V +103 -39 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +104 -3 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -33 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -28 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -12 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -15 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -6 R +0 35 V +-31 -35 R +62 0 V +-62 35 R +62 0 V +103 -54 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -23 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +104 -4 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +103 -32 R +0 86 V +-31 -86 R +62 0 V +-62 86 R +62 0 V +720 1285 Pnt +854 1280 Pnt +988 1294 Pnt +1122 1286 Pnt +1257 1296 Pnt +1391 1287 Pnt +1525 1294 Pnt +1659 1277 Pnt +1793 1275 Pnt +1927 1290 Pnt +2061 1283 Pnt +2195 1277 Pnt +2330 1284 Pnt +2464 1293 Pnt +2598 1286 Pnt +2732 1298 Pnt +2866 1294 Pnt +3000 1294 Pnt +3134 1282 Pnt +3269 1295 Pnt +3403 1275 Pnt +3537 1262 Pnt +3671 1261 Pnt +3805 1254 Pnt +3939 1271 Pnt +4073 1242 Pnt +4207 1235 Pnt +4342 1243 Pnt +4476 1260 Pnt +2.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +720 3109 N +720 640 L +3959 0 V +0 2469 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +720 3349 M +31 0 V +3928 0 R +-31 0 V +720 3439 M +31 0 V +3928 0 R +-31 0 V +720 3512 M +31 0 V +3928 0 R +-31 0 V +720 3574 M +31 0 V +3928 0 R +-31 0 V +720 3628 M +31 0 V +3928 0 R +-31 0 V +720 3675 M +31 0 V +3928 0 R +-31 0 V +720 3718 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 3997 M +31 0 V +3928 0 R +-31 0 V +720 4160 M +31 0 V +3928 0 R +-31 0 V +720 4276 M +31 0 V +3928 0 R +-31 0 V +720 4365 M +31 0 V +3928 0 R +-31 0 V +720 4439 M +31 0 V +3928 0 R +-31 0 V +720 4501 M +31 0 V +3928 0 R +-31 0 V +720 4555 M +31 0 V +3928 0 R +-31 0 V +720 4602 M +31 0 V +3928 0 R +-31 0 V +720 4644 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 4923 M +31 0 V +3928 0 R +-31 0 V +720 5086 M +31 0 V +3928 0 R +-31 0 V +720 5202 M +31 0 V +3928 0 R +-31 0 V +720 5292 M +31 0 V +3928 0 R +-31 0 V +720 5365 M +31 0 V +3928 0 R +-31 0 V +720 5427 M +31 0 V +3928 0 R +-31 0 V +720 5481 M +31 0 V +3928 0 R +-31 0 V +720 5529 M +31 0 V +3928 0 R +-31 0 V +720 5571 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 5850 M +31 0 V +3928 0 R +-31 0 V +720 6013 M +31 0 V +3928 0 R +-31 0 V +720 3349 M +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +922 3349 M +0 31 V +0 2633 R +0 -31 V +1040 3349 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1593 3349 M +0 31 V +0 2633 R +0 -31 V +1711 3349 M +0 31 V +0 2633 R +0 -31 V +83 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +54 -2633 R +0 31 V +0 2633 R +0 -31 V +44 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +35 -2633 R +0 31 V +0 2633 R +0 -31 V +30 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2263 3349 M +0 31 V +0 2633 R +0 -31 V +2381 3349 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2934 3349 M +0 31 V +0 2633 R +0 -31 V +3052 3349 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +3605 3349 M +0 31 V +0 2633 R +0 -31 V +3723 3349 M +0 31 V +0 2633 R +0 -31 V +83 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +54 -2633 R +0 31 V +0 2633 R +0 -31 V +44 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +35 -2633 R +0 31 V +0 2633 R +0 -31 V +30 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +4275 3349 M +0 31 V +0 2633 R +0 -31 V +4393 3349 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +stroke +LTb +LCb setrgbcolor +720 6013 N +0 -2664 V +3959 0 V +0 2664 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +LCb setrgbcolor +LTb +LCb setrgbcolor +LTb +1.000 UP +1.000 UL +LTb +0.37 0.51 0.71 C 854 4213 M +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -13 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -14 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +104 -20 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +103 -16 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -18 R +0 22 V +-31 -22 R +62 0 V +-62 22 R +62 0 V +103 -20 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +103 -20 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 6 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -19 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -43 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +104 -46 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +103 -49 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -51 R +0 38 V +-31 -38 R +62 0 V +-62 38 R +62 0 V +103 -79 R +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +103 -53 R +0 30 V +-31 -30 R +62 0 V +-62 30 R +62 0 V +103 -71 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 -45 R +0 21 V +-31 -21 R +62 0 V +stroke 3165 3951 M +-62 21 R +62 0 V +104 -68 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -41 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -74 R +0 41 V +-31 -41 R +62 0 V +-62 41 R +62 0 V +103 -55 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +103 -45 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -52 R +0 31 V +-31 -31 R +62 0 V +-62 31 R +62 0 V +103 -47 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -39 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +104 -12 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -2 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 12 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +854 4225 Pnt +988 4236 Pnt +1122 4241 Pnt +1257 4231 Pnt +1391 4228 Pnt +1525 4228 Pnt +1659 4226 Pnt +1793 4220 Pnt +1927 4240 Pnt +2061 4236 Pnt +2195 4209 Pnt +2330 4179 Pnt +2464 4145 Pnt +2598 4121 Pnt +2732 4067 Pnt +2866 4034 Pnt +3000 3987 Pnt +3134 3962 Pnt +3269 3913 Pnt +3403 3886 Pnt +3537 3839 Pnt +3671 3811 Pnt +3805 3781 Pnt +3939 3753 Pnt +4073 3727 Pnt +4207 3704 Pnt +4342 3708 Pnt +4476 3717 Pnt +4610 3752 Pnt +1.000 UP +1.000 UL +LTb +0.88 0.61 0.14 C 854 4211 M +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +103 -7 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -4 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +104 6 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +103 -1 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -27 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 -27 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 5 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 25 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 55 R +0 26 V +-31 -26 R +62 0 V +-62 26 R +62 0 V +103 153 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +104 230 R +0 19 V +-31 -19 R +62 0 V +-62 19 R +62 0 V +103 233 R +0 29 V +-31 -29 R +62 0 V +-62 29 R +62 0 V +103 -18 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -102 R +0 47 V +-31 -47 R +62 0 V +-62 47 R +62 0 V +103 -153 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 -134 R +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +103 -141 R +0 17 V +-31 -17 R +62 0 V +stroke 3165 4664 M +-62 17 R +62 0 V +104 -102 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +103 -141 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -132 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +103 -103 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -118 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +103 -59 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -84 R +0 19 V +-31 -19 R +62 0 V +-62 19 R +62 0 V +103 -37 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +104 -28 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -21 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -154 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +854 4227 Pnt +988 4238 Pnt +1122 4241 Pnt +1257 4258 Pnt +1391 4269 Pnt +1525 4261 Pnt +1659 4256 Pnt +1793 4282 Pnt +1927 4331 Pnt +2061 4411 Pnt +2195 4584 Pnt +2330 4831 Pnt +2464 5089 Pnt +2598 5095 Pnt +2732 5027 Pnt +2866 4902 Pnt +3000 4789 Pnt +3134 4673 Pnt +3269 4583 Pnt +3403 4455 Pnt +3537 4337 Pnt +3671 4256 Pnt +3805 4157 Pnt +3939 4108 Pnt +4073 4039 Pnt +4207 4018 Pnt +4342 4002 Pnt +4476 3997 Pnt +4610 3857 Pnt +1.000 UP +1.000 UL +LTb +0.56 0.69 0.20 C 854 5625 M +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -45 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 -19 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +104 9 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -12 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 -24 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -37 R +0 53 V +-31 -53 R +62 0 V +-62 53 R +62 0 V +103 -40 R +0 47 V +-31 -47 R +62 0 V +-62 47 R +62 0 V +103 -46 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +103 -12 R +0 41 V +-31 -41 R +62 0 V +-62 41 R +62 0 V +103 -59 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +104 -98 R +0 48 V +-31 -48 R +62 0 V +-62 48 R +62 0 V +103 -213 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -180 R +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +103 -182 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +103 -165 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -127 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -167 R +0 15 V +-31 -15 R +62 0 V +stroke 3165 4586 M +-62 15 R +62 0 V +104 -126 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -135 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -129 R +0 24 V +-31 -24 R +62 0 V +-62 24 R +62 0 V +103 -112 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +103 -46 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -85 R +0 28 V +-31 -28 R +62 0 V +-62 28 R +62 0 V +103 -54 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -43 R +0 18 V +-31 -18 R +62 0 V +-62 18 R +62 0 V +104 -27 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -23 R +0 26 V +-31 -26 R +62 0 V +-62 26 R +62 0 V +103 -59 R +0 63 V +-31 -63 R +62 0 V +-62 63 R +62 0 V +854 5632 Pnt +988 5606 Pnt +1122 5604 Pnt +1257 5620 Pnt +1391 5624 Pnt +1525 5620 Pnt +1659 5616 Pnt +1793 5626 Pnt +1927 5604 Pnt +2061 5615 Pnt +2195 5581 Pnt +2330 5513 Pnt +2464 5335 Pnt +2598 5183 Pnt +2732 5021 Pnt +2866 4865 Pnt +3000 4749 Pnt +3134 4593 Pnt +3269 4482 Pnt +3403 4365 Pnt +3537 4260 Pnt +3671 4169 Pnt +3805 4136 Pnt +3939 4070 Pnt +4073 4037 Pnt +4207 4011 Pnt +4342 4003 Pnt +4476 4003 Pnt +4610 3990 Pnt +2.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +720 6013 N +0 -2664 V +3959 0 V +0 2664 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +720 6253 M +31 0 V +3928 0 R +-31 0 V +720 6438 M +31 0 V +3928 0 R +-31 0 V +720 6589 M +31 0 V +3928 0 R +-31 0 V +720 6716 M +31 0 V +3928 0 R +-31 0 V +720 6827 M +31 0 V +3928 0 R +-31 0 V +720 6924 M +31 0 V +3928 0 R +-31 0 V +720 7011 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 7585 M +31 0 V +3928 0 R +-31 0 V +720 7921 M +31 0 V +3928 0 R +-31 0 V +720 8159 M +31 0 V +3928 0 R +-31 0 V +720 8343 M +31 0 V +3928 0 R +-31 0 V +720 8494 M +31 0 V +3928 0 R +-31 0 V +720 8622 M +31 0 V +3928 0 R +-31 0 V +720 8732 M +31 0 V +3928 0 R +-31 0 V +720 8830 M +31 0 V +3928 0 R +-31 0 V +720 8917 M +63 0 V +3896 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +720 6253 M +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +922 6253 M +0 31 V +0 2633 R +0 -31 V +1040 6253 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1593 6253 M +0 31 V +0 2633 R +0 -31 V +1711 6253 M +0 31 V +0 2633 R +0 -31 V +83 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +54 -2633 R +0 31 V +0 2633 R +0 -31 V +44 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +35 -2633 R +0 31 V +0 2633 R +0 -31 V +30 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2263 6253 M +0 31 V +0 2633 R +0 -31 V +2381 6253 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2934 6253 M +0 31 V +0 2633 R +0 -31 V +3052 6253 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +34 -2633 R +0 31 V +0 2633 R +0 -31 V +31 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +3605 6253 M +0 31 V +0 2633 R +0 -31 V +3723 6253 M +0 31 V +0 2633 R +0 -31 V +83 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +54 -2633 R +0 31 V +0 2633 R +0 -31 V +44 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +35 -2633 R +0 31 V +0 2633 R +0 -31 V +30 -2633 R +0 63 V +0 2601 R +0 -63 V +stroke +LTb +LCb setrgbcolor +4275 6253 M +0 31 V +0 2633 R +0 -31 V +4393 6253 M +0 31 V +0 2633 R +0 -31 V +84 -2633 R +0 31 V +0 2633 R +0 -31 V +65 -2633 R +0 31 V +0 2633 R +0 -31 V +53 -2633 R +0 31 V +0 2633 R +0 -31 V +45 -2633 R +0 31 V +0 2633 R +0 -31 V +39 -2633 R +0 31 V +0 2633 R +0 -31 V +stroke +LTb +LCb setrgbcolor +720 8917 N +0 -2664 V +3959 0 V +0 2664 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +LCb setrgbcolor +LTb +LCb setrgbcolor +LTb +1.000 UP +1.000 UL +LTb +0.37 0.51 0.71 C LCb setrgbcolor +1.000 UP +1.000 UL +LTb +0.37 0.51 0.71 C 3896 8754 M +543 0 V +-543 31 R +0 -62 V +543 62 R +0 -62 V +854 6695 M +0 19 V +-31 -19 R +62 0 V +-62 19 R +62 0 V +103 35 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +103 -50 R +0 50 V +-31 -50 R +62 0 V +-62 50 R +62 0 V +104 -77 R +0 29 V +-31 -29 R +62 0 V +-62 29 R +62 0 V +103 -29 R +0 53 V +-31 -53 R +62 0 V +-62 53 R +62 0 V +103 3 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 -34 R +0 35 V +-31 -35 R +62 0 V +-62 35 R +62 0 V +103 -63 R +0 66 V +-31 -66 R +62 0 V +-62 66 R +62 0 V +103 -61 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 -48 R +0 41 V +-31 -41 R +62 0 V +-62 41 R +62 0 V +103 -34 R +0 36 V +-31 -36 R +62 0 V +-62 36 R +62 0 V +104 -54 R +0 61 V +-31 -61 R +62 0 V +-62 61 R +62 0 V +103 -29 R +0 38 V +-31 -38 R +62 0 V +-62 38 R +62 0 V +103 -19 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 -80 R +0 35 V +-31 -35 R +62 0 V +-62 35 R +62 0 V +103 -11 R +0 70 V +-31 -70 R +62 0 V +-62 70 R +62 0 V +103 -89 R +0 70 V +-31 -70 R +62 0 V +stroke 3031 6674 M +-62 70 R +62 0 V +103 -93 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +104 -61 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -112 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 82 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 150 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 102 R +0 47 V +-31 -47 R +62 0 V +-62 47 R +62 0 V +103 60 R +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +103 -8 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -129 R +0 57 V +-31 -57 R +62 0 V +-62 57 R +62 0 V +104 -19 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +103 -8 R +0 49 V +-31 -49 R +62 0 V +-62 49 R +62 0 V +103 -81 R +0 68 V +-31 -68 R +62 0 V +-62 68 R +62 0 V +854 6704 Pnt +988 6754 Pnt +1122 6733 Pnt +1257 6696 Pnt +1391 6708 Pnt +1525 6749 Pnt +1659 6744 Pnt +1793 6731 Pnt +1927 6717 Pnt +2061 6703 Pnt +2195 6708 Pnt +2330 6702 Pnt +2464 6722 Pnt +2598 6736 Pnt +2732 6687 Pnt +2866 6729 Pnt +3000 6710 Pnt +3134 6655 Pnt +3269 6615 Pnt +3403 6522 Pnt +3537 6616 Pnt +3671 6790 Pnt +3805 6929 Pnt +3939 7028 Pnt +4073 7053 Pnt +4207 6970 Pnt +4342 6989 Pnt +4476 7017 Pnt +4610 6995 Pnt +4167 8754 Pnt +1.000 UP +1.000 UL +LTb +0.88 0.61 0.14 C LCb setrgbcolor +1.000 UP +1.000 UL +LTb +0.88 0.61 0.14 C 3896 8554 M +543 0 V +-543 31 R +0 -62 V +543 62 R +0 -62 V +854 6612 M +0 95 V +-31 -95 R +62 0 V +-62 95 R +62 0 V +103 -6 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +103 -54 R +0 39 V +-31 -39 R +62 0 V +-62 39 R +62 0 V +104 -38 R +0 104 V +-31 -104 R +62 0 V +-62 104 R +62 0 V +103 -38 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -22 R +0 79 V +-31 -79 R +62 0 V +-62 79 R +62 0 V +103 -159 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +103 47 R +0 46 V +-31 -46 R +62 0 V +-62 46 R +62 0 V +103 -63 R +0 36 V +-31 -36 R +62 0 V +-62 36 R +62 0 V +103 -2 R +0 25 V +-31 -25 R +62 0 V +-62 25 R +62 0 V +103 -130 R +0 89 V +-31 -89 R +62 0 V +-62 89 R +62 0 V +104 -48 R +0 73 V +-31 -73 R +62 0 V +-62 73 R +62 0 V +103 -23 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 -15 R +0 49 V +-31 -49 R +62 0 V +-62 49 R +62 0 V +103 -47 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +103 -24 R +0 35 V +-31 -35 R +62 0 V +-62 35 R +62 0 V +103 -5 R +0 19 V +-31 -19 R +62 0 V +stroke 3031 6754 M +-62 19 R +62 0 V +103 75 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +104 -1 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 145 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +103 86 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 153 R +0 61 V +-31 -61 R +62 0 V +-62 61 R +62 0 V +103 101 R +0 33 V +-31 -33 R +62 0 V +-62 33 R +62 0 V +103 27 R +0 17 V +-31 -17 R +62 0 V +-62 17 R +62 0 V +103 -11 R +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +103 19 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +104 -55 R +0 33 V +-31 -33 R +62 0 V +-62 33 R +62 0 V +103 -12 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +103 21 R +0 118 V +-31 -118 R +62 0 V +-62 118 R +62 0 V +854 6661 Pnt +988 6714 Pnt +1122 6693 Pnt +1257 6729 Pnt +1391 6749 Pnt +1525 6776 Pnt +1659 6666 Pnt +1793 6746 Pnt +1927 6724 Pnt +2061 6753 Pnt +2195 6681 Pnt +2330 6713 Pnt +2464 6738 Pnt +2598 6759 Pnt +2732 6742 Pnt +2866 6742 Pnt +3000 6764 Pnt +3134 6865 Pnt +3269 6889 Pnt +3403 7054 Pnt +3537 7159 Pnt +3671 7351 Pnt +3805 7499 Pnt +3939 7551 Pnt +4073 7554 Pnt +4207 7583 Pnt +4342 7549 Pnt +4476 7559 Pnt +4610 7646 Pnt +4167 8554 Pnt +1.000 UP +1.000 UL +LTb +0.56 0.69 0.20 C LCb setrgbcolor +1.000 UP +1.000 UL +LTb +0.56 0.69 0.20 C 3896 8354 M +543 0 V +-543 31 R +0 -62 V +543 62 R +0 -62 V +854 7085 M +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +103 -11 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -77 R +0 54 V +-31 -54 R +62 0 V +-62 54 R +62 0 V +104 -35 R +0 48 V +-31 -48 R +62 0 V +-62 48 R +62 0 V +103 -33 R +0 50 V +-31 -50 R +62 0 V +-62 50 R +62 0 V +103 -76 R +0 26 V +-31 -26 R +62 0 V +-62 26 R +62 0 V +103 35 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -56 R +0 25 V +-31 -25 R +62 0 V +-62 25 R +62 0 V +103 -53 R +0 28 V +-31 -28 R +62 0 V +-62 28 R +62 0 V +103 -78 R +0 75 V +-31 -75 R +62 0 V +-62 75 R +62 0 V +103 44 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +104 -56 R +0 42 V +-31 -42 R +62 0 V +-62 42 R +62 0 V +103 -35 R +0 52 V +-31 -52 R +62 0 V +-62 52 R +62 0 V +103 -20 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -68 R +0 37 V +-31 -37 R +62 0 V +-62 37 R +62 0 V +103 5 R +0 32 V +-31 -32 R +62 0 V +-62 32 R +62 0 V +103 86 R +0 38 V +-31 -38 R +62 0 V +stroke 3031 7222 M +-62 38 R +62 0 V +103 110 R +0 60 V +-31 -60 R +62 0 V +-62 60 R +62 0 V +104 78 R +0 95 V +-31 -95 R +62 0 V +-62 95 R +62 0 V +103 185 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +103 78 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +103 -43 R +0 29 V +-31 -29 R +62 0 V +-62 29 R +62 0 V +103 34 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +103 -237 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +103 -41 R +0 34 V +-31 -34 R +62 0 V +-62 34 R +62 0 V +103 -101 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +104 -57 R +0 21 V +-31 -21 R +62 0 V +-62 21 R +62 0 V +103 -7 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +103 -33 R +0 23 V +-31 -23 R +62 0 V +-62 23 R +62 0 V +854 7087 Pnt +988 7095 Pnt +1122 7063 Pnt +1257 7078 Pnt +1391 7094 Pnt +1525 7056 Pnt +1659 7114 Pnt +1793 7081 Pnt +1927 7054 Pnt +2061 7028 Pnt +2195 7111 Pnt +2330 7078 Pnt +2464 7091 Pnt +2598 7113 Pnt +2732 7080 Pnt +2866 7120 Pnt +3000 7242 Pnt +3134 7401 Pnt +3269 7557 Pnt +3403 7791 Pnt +3537 7879 Pnt +3671 7858 Pnt +3805 7917 Pnt +3939 7698 Pnt +4073 7682 Pnt +4207 7610 Pnt +4342 7574 Pnt +4476 7585 Pnt +4610 7569 Pnt +4167 8354 Pnt +2.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +720 8917 N +0 -2664 V +3959 0 V +0 2664 V +-3959 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +stroke +grestore +end +showpage + }}% + \put(3776,8354){\makebox(0,0)[r]{\strut{}Metropolis}}% + \put(3776,8554){\makebox(0,0)[r]{\strut{}Hybrid Wolff/Metropolis}}% + \put(3776,8754){\makebox(0,0)[r]{\strut{}Wolff}}% + \put(340,7585){% + \special{ps: gsave currentpoint currentpoint translate +630 rotate neg exch neg exch translate}% + \makebox(0,0){\strut{}$\tau/n$}% + \special{ps: currentpoint grestore moveto}% + }% + \put(600,8917){\makebox(0,0)[r]{\strut{}$10$}}% + \put(600,7011){\makebox(0,0)[r]{\strut{}$1$}}% + \put(220,4681){% + \special{ps: gsave currentpoint currentpoint translate +630 rotate neg exch neg exch translate}% + \makebox(0,0){\strut{}$\tau/n$}% + \special{ps: currentpoint grestore moveto}% + }% + \put(600,5571){\makebox(0,0)[r]{\strut{}$100$}}% + \put(600,4644){\makebox(0,0)[r]{\strut{}$10$}}% + \put(600,3718){\makebox(0,0)[r]{\strut{}$1$}}% + \put(2699,140){\makebox(0,0){\strut{}$\beta H$}}% + \put(220,1874){% + \special{ps: gsave currentpoint currentpoint translate +630 rotate neg exch neg exch translate}% + \makebox(0,0){\strut{}$\tau/n$}% + \special{ps: currentpoint grestore moveto}% + }% + \put(4073,440){\makebox(0,0){\strut{}$10^{0}$}}% + \put(3403,440){\makebox(0,0){\strut{}$10^{-1}$}}% + \put(2732,440){\makebox(0,0){\strut{}$10^{-2}$}}% + \put(2061,440){\makebox(0,0){\strut{}$10^{-3}$}}% + \put(1391,440){\makebox(0,0){\strut{}$10^{-4}$}}% + \put(720,440){\makebox(0,0){\strut{}$10^{-5}$}}% + \put(600,2699){\makebox(0,0)[r]{\strut{}$100$}}% + \put(600,1841){\makebox(0,0)[r]{\strut{}$10$}}% + \put(600,982){\makebox(0,0)[r]{\strut{}$1$}}% +\end{picture}% +\endgroup +\endinput diff --git a/fig_correlation_collapse-hL.tex b/fig_correlation_collapse-hL.tex new file mode 100644 index 0000000..9609378 --- /dev/null +++ b/fig_correlation_collapse-hL.tex @@ -0,0 +1,1878 @@ +% GNUPLOT: LaTeX picture with Postscript +\begingroup% +\makeatletter% +\newcommand{\GNUPLOTspecial}{% + \@sanitize\catcode`\%=14\relax\special}% +\setlength{\unitlength}{0.0500bp}% +\begin{picture}(4920,3038)(0,0)% + {\GNUPLOTspecial{" +%!PS-Adobe-2.0 EPSF-2.0 +%%Creator: gnuplot 5.2 patchlevel 0 +%%CreationDate: Mon Nov 6 21:12:43 2017 +%%DocumentFonts: +%%BoundingBox: 0 0 246 151 +%%EndComments +%%BeginProlog +/gnudict 256 dict def +gnudict begin +% +% The following true/false flags may be edited by hand if desired. +% The unit line width and grayscale image gamma correction may also be changed. +% +/Color false def +/Blacktext true def +/Solid false def +/Dashlength 1 def +/Landscape false def +/Level1 false def +/Level3 false def +/Rounded false def +/ClipToBoundingBox false def +/SuppressPDFMark false def +/TransparentPatterns false def +/gnulinewidth 5.000 def +/userlinewidth gnulinewidth def +/Gamma 1.0 def +/BackgroundColor {-1.000 -1.000 -1.000} def +% +/vshift -66 def +/dl1 { + 10.0 Dashlength userlinewidth gnulinewidth div mul mul mul + Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if +} def +/dl2 { + 10.0 Dashlength userlinewidth gnulinewidth div mul mul mul + Rounded { currentlinewidth 0.75 mul add } if +} def +/hpt_ 31.5 def +/vpt_ 31.5 def +/hpt hpt_ def +/vpt vpt_ def +/doclip { + ClipToBoundingBox { + newpath 0 0 moveto 246 0 lineto 246 151 lineto 0 151 lineto closepath + clip + } if +} def +% +% Gnuplot Prolog Version 5.1 (Oct 2015) +% +%/SuppressPDFMark true def +% +/M {moveto} bind def +/L {lineto} bind def +/R {rmoveto} bind def +/V {rlineto} bind def +/N {newpath moveto} bind def +/Z {closepath} bind def +/C {setrgbcolor} bind def +/f {rlineto fill} bind def +/g {setgray} bind def +/Gshow {show} def % May be redefined later in the file to support UTF-8 +/vpt2 vpt 2 mul def +/hpt2 hpt 2 mul def +/Lshow {currentpoint stroke M 0 vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def +/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} def +/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R + Blacktext {gsave 0 setgray textshow grestore} {textshow} ifelse} 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 0 setgray Solid {pop []} if 0 setdash} ifelse} def +/BL {stroke userlinewidth 2 mul setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +/AL {stroke userlinewidth 2 div setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +/UL {dup gnulinewidth mul /userlinewidth exch def + dup 1 lt {pop 1} if 10 mul /udl exch def} def +/PL {stroke userlinewidth setlinewidth + Rounded {1 setlinejoin 1 setlinecap} if} def +3.8 setmiterlimit +% Classic Line colors (version 5.0) +/LCw {1 1 1} def +/LCb {0 0 0} def +/LCa {0 0 0} def +/LC0 {1 0 0} def +/LC1 {0 1 0} def +/LC2 {0 0 1} def +/LC3 {1 0 1} def +/LC4 {0 1 1} def +/LC5 {1 1 0} def +/LC6 {0 0 0} def +/LC7 {1 0.3 0} def +/LC8 {0.5 0.5 0.5} def +% Default dash patterns (version 5.0) +/LTB {BL [] LCb DL} def +/LTw {PL [] 1 setgray} def +/LTb {PL [] LCb DL} def +/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def +/LT0 {PL [] LC0 DL} def +/LT1 {PL [2 dl1 3 dl2] LC1 DL} def +/LT2 {PL [1 dl1 1.5 dl2] LC2 DL} def +/LT3 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC3 DL} def +/LT4 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def +/LT5 {PL [4 dl1 2 dl2] LC5 DL} def +/LT6 {PL [1.5 dl1 1.5 dl2 1.5 dl1 1.5 dl2 1.5 dl1 6 dl2] LC6 DL} def +/LT7 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC7 DL} def +/LT8 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC8 DL} def +/SL {[] 0 setdash} 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 +/Density { + /Fillden exch def + currentrgbcolor + /ColB exch def /ColG exch def /ColR exch def + /ColR ColR Fillden mul Fillden sub 1 add def + /ColG ColG Fillden mul Fillden sub 1 add def + /ColB ColB Fillden mul Fillden sub 1 add def + ColR ColG ColB setrgbcolor} def +/BoxColFill {gsave Rec PolyFill} def +/PolyFill {gsave Density fill grestore grestore} def +/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def +% +% PostScript Level 1 Pattern Fill routine for rectangles +% Usage: x y w h s a XX PatternFill +% x,y = lower left corner of box to be filled +% w,h = width and height of box +% a = angle in degrees between lines and x-axis +% XX = 0/1 for no/yes cross-hatch +% +/PatternFill {gsave /PFa [ 9 2 roll ] def + PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate + PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec + TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse + clip + currentlinewidth 0.5 mul setlinewidth + /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def + 0 0 M PFa 5 get rotate PFs -2 div dup translate + 0 1 PFs PFa 4 get div 1 add floor cvi + {PFa 4 get mul 0 M 0 PFs V} for + 0 PFa 6 get ne { + 0 1 PFs PFa 4 get div 1 add floor cvi + {PFa 4 get mul 0 2 1 roll M PFs 0 V} for + } if + stroke grestore} def +% +/languagelevel where + {pop languagelevel} {1} ifelse +dup 2 lt + {/InterpretLevel1 true def + /InterpretLevel3 false def} + {/InterpretLevel1 Level1 def + 2 gt + {/InterpretLevel3 Level3 def} + {/InterpretLevel3 false def} + ifelse } + ifelse +% +% PostScript level 2 pattern fill definitions +% +/Level2PatternFill { +/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8} + bind def +/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} +>> matrix makepattern +/Pat1 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke + 0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke} +>> matrix makepattern +/Pat2 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L + 8 8 L 8 0 L 0 0 L fill} +>> matrix makepattern +/Pat3 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L + 0 12 M 12 0 L stroke} +>> matrix makepattern +/Pat4 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L + 0 -4 M 12 8 L stroke} +>> matrix makepattern +/Pat5 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L + 0 12 M 8 -4 L 4 12 M 10 0 L stroke} +>> matrix makepattern +/Pat6 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L + 0 -4 M 8 12 L 4 -4 M 10 8 L stroke} +>> matrix makepattern +/Pat7 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L + 12 0 M -4 8 L 12 4 M 0 10 L stroke} +>> matrix makepattern +/Pat8 exch def +<< Tile8x8 + /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L + -4 0 M 12 8 L -4 4 M 8 10 L stroke} +>> matrix makepattern +/Pat9 exch def +/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def +/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def +/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def +/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def +/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def +/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def +/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def +} def +% +% +%End of PostScript Level 2 code +% +/PatternBgnd { + TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse +} def +% +% Substitute for Level 2 pattern fill codes with +% grayscale if Level 2 support is not selected. +% +/Level1PatternFill { +/Pattern1 {0.250 Density} bind def +/Pattern2 {0.500 Density} bind def +/Pattern3 {0.750 Density} bind def +/Pattern4 {0.125 Density} bind def +/Pattern5 {0.375 Density} bind def +/Pattern6 {0.625 Density} bind def +/Pattern7 {0.875 Density} bind def +} def +% +% Now test for support of Level 2 code +% +Level1 {Level1PatternFill} {Level2PatternFill} ifelse +% +/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 +% +Level1 SuppressPDFMark or +{} { +/SDict 10 dict def +systemdict /pdfmark known not { + userdict /pdfmark systemdict /cleartomark get put +} if +SDict begin [ + /Title () + /Subject (gnuplot plot) + /Creator (gnuplot 5.2 patchlevel 0) +% /Producer (gnuplot) +% /Keywords () + /CreationDate (Mon Nov 6 21:12:43 2017) + /DOCINFO pdfmark +end +} ifelse +% +% Support for boxed text - Ethan A Merritt Sep 2016 +% +/InitTextBox { userdict /TBy2 3 -1 roll put userdict /TBx2 3 -1 roll put + userdict /TBy1 3 -1 roll put userdict /TBx1 3 -1 roll put + /Boxing true def } def +/ExtendTextBox { dup type /stringtype eq + { Boxing { gsave dup false charpath pathbbox + dup TBy2 gt {userdict /TBy2 3 -1 roll put} {pop} ifelse + dup TBx2 gt {userdict /TBx2 3 -1 roll put} {pop} ifelse + dup TBy1 lt {userdict /TBy1 3 -1 roll put} {pop} ifelse + dup TBx1 lt {userdict /TBx1 3 -1 roll put} {pop} ifelse + grestore } if } + {} ifelse} def +/PopTextBox { newpath TBx1 TBxmargin sub TBy1 TBymargin sub M + TBx1 TBxmargin sub TBy2 TBymargin add L + TBx2 TBxmargin add TBy2 TBymargin add L + TBx2 TBxmargin add TBy1 TBymargin sub L closepath } def +/DrawTextBox { PopTextBox stroke /Boxing false def} def +/FillTextBox { gsave PopTextBox fill grestore /Boxing false def} def +0 0 0 0 InitTextBox +/TBxmargin 20 def +/TBymargin 20 def +/Boxing false def +/textshow { ExtendTextBox Gshow } def +% +end +%%EndProlog +%%Page: 1 1 +gnudict begin +gsave +doclip +0 0 translate +0.050 0.050 scale +0 setgray +newpath +BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 4920.00 3038.00 BoxColFill} if +1.000 UL +LTb +LCb setrgbcolor +860 640 M +31 0 V +3668 0 R +-31 0 V +860 746 M +63 0 V +3636 0 R +-63 0 V +stroke +LTb +LCb setrgbcolor +860 1443 M +31 0 V +3668 0 R +-31 0 V +860 1851 M +31 0 V +3668 0 R +-31 0 V +860 2140 M +31 0 V +3668 0 R +-31 0 V +860 2364 M +31 0 V +3668 0 R +-31 0 V +860 2548 M +31 0 V +3668 0 R +-31 0 V +860 2703 M +31 0 V +3668 0 R +-31 0 V +860 2837 M +31 0 V +3668 0 R +-31 0 V +860 640 M +0 31 V +0 2166 R +0 -31 V +894 640 M +0 31 V +0 2166 R +0 -31 V +923 640 M +0 31 V +0 2166 R +0 -31 V +948 640 M +0 31 V +0 2166 R +0 -31 V +970 640 M +0 31 V +0 2166 R +0 -31 V +989 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1119 640 M +0 31 V +0 2166 R +0 -31 V +1195 640 M +0 31 V +0 2166 R +0 -31 V +1248 640 M +0 31 V +0 2166 R +0 -31 V +1290 640 M +0 31 V +0 2166 R +0 -31 V +1324 640 M +0 31 V +0 2166 R +0 -31 V +1353 640 M +0 31 V +0 2166 R +0 -31 V +1378 640 M +0 31 V +0 2166 R +0 -31 V +1400 640 M +0 31 V +0 2166 R +0 -31 V +1419 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1549 640 M +0 31 V +0 2166 R +0 -31 V +1625 640 M +0 31 V +0 2166 R +0 -31 V +1678 640 M +0 31 V +0 2166 R +0 -31 V +1720 640 M +0 31 V +0 2166 R +0 -31 V +1754 640 M +0 31 V +0 2166 R +0 -31 V +1783 640 M +0 31 V +0 2166 R +0 -31 V +1808 640 M +0 31 V +0 2166 R +0 -31 V +1830 640 M +0 31 V +0 2166 R +0 -31 V +1849 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +1979 640 M +0 31 V +0 2166 R +0 -31 V +2055 640 M +0 31 V +0 2166 R +0 -31 V +2108 640 M +0 31 V +0 2166 R +0 -31 V +2150 640 M +0 31 V +0 2166 R +0 -31 V +2184 640 M +0 31 V +0 2166 R +0 -31 V +2213 640 M +0 31 V +0 2166 R +0 -31 V +2238 640 M +0 31 V +0 2166 R +0 -31 V +2260 640 M +0 31 V +0 2166 R +0 -31 V +2279 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2409 640 M +0 31 V +0 2166 R +0 -31 V +2485 640 M +0 31 V +0 2166 R +0 -31 V +2538 640 M +0 31 V +0 2166 R +0 -31 V +2580 640 M +0 31 V +0 2166 R +0 -31 V +2614 640 M +0 31 V +0 2166 R +0 -31 V +2643 640 M +0 31 V +0 2166 R +0 -31 V +2668 640 M +0 31 V +0 2166 R +0 -31 V +2690 640 M +0 31 V +0 2166 R +0 -31 V +2710 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +2839 640 M +0 31 V +0 2166 R +0 -31 V +2915 640 M +0 31 V +0 2166 R +0 -31 V +2968 640 M +0 31 V +0 2166 R +0 -31 V +3010 640 M +0 31 V +0 2166 R +0 -31 V +3044 640 M +0 31 V +0 2166 R +0 -31 V +3073 640 M +0 31 V +0 2166 R +0 -31 V +3098 640 M +0 31 V +0 2166 R +0 -31 V +3120 640 M +0 31 V +0 2166 R +0 -31 V +3140 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +3269 640 M +0 31 V +0 2166 R +0 -31 V +3345 640 M +0 31 V +0 2166 R +0 -31 V +3398 640 M +0 31 V +0 2166 R +0 -31 V +3440 640 M +0 31 V +0 2166 R +0 -31 V +3474 640 M +0 31 V +0 2166 R +0 -31 V +3503 640 M +0 31 V +0 2166 R +0 -31 V +3528 640 M +0 31 V +0 2166 R +0 -31 V +3550 640 M +0 31 V +0 2166 R +0 -31 V +3570 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +3699 640 M +0 31 V +0 2166 R +0 -31 V +3775 640 M +0 31 V +0 2166 R +0 -31 V +3828 640 M +0 31 V +0 2166 R +0 -31 V +3870 640 M +0 31 V +0 2166 R +0 -31 V +3904 640 M +0 31 V +0 2166 R +0 -31 V +3933 640 M +0 31 V +0 2166 R +0 -31 V +3958 640 M +0 31 V +0 2166 R +0 -31 V +3980 640 M +0 31 V +0 2166 R +0 -31 V +4000 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +4129 640 M +0 31 V +0 2166 R +0 -31 V +4205 640 M +0 31 V +0 2166 R +0 -31 V +4258 640 M +0 31 V +0 2166 R +0 -31 V +4300 640 M +0 31 V +0 2166 R +0 -31 V +4334 640 M +0 31 V +0 2166 R +0 -31 V +4363 640 M +0 31 V +0 2166 R +0 -31 V +4388 640 M +0 31 V +0 2166 R +0 -31 V +4410 640 M +0 31 V +0 2166 R +0 -31 V +4430 640 M +0 63 V +0 2134 R +0 -63 V +stroke +LTb +LCb setrgbcolor +4559 640 M +0 31 V +0 2166 R +0 -31 V +stroke +LTb +LCb setrgbcolor +860 2837 N +860 640 L +3699 0 V +0 2197 V +-3699 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +LCb setrgbcolor +LTb +LCb setrgbcolor +LTb +1.000 UP +1.000 UL +LTb +0.58 0.00 0.83 C 1100 2588 M +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -2 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -5 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 3 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +55 -16 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -7 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -6 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -5 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -10 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -8 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -3 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 1 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -8 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +55 -22 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -35 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -80 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -112 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +55 -125 R +0 6 V +-31 -6 R +62 0 V +stroke 2593 2247 M +-62 6 R +62 0 V +55 -113 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -104 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +55 -101 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -98 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -88 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -84 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -76 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -50 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -34 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -15 R +0 989 V +-31 -989 R +62 0 V +-62 989 R +62 0 V +55 -989 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 0 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +1343 2612 M +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 2 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -8 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 0 R +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +55 -20 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -2 R +0 2 V +stroke 1773 2616 M +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 4 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -6 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -15 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -7 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -18 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +55 -37 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -74 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -107 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -111 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -100 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -99 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +55 -103 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -97 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -99 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -96 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -95 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -92 R +0 9 V +-31 -9 R +62 0 V +stroke 3266 1553 M +-62 9 R +62 0 V +55 -90 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -69 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -57 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -38 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -11 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -7 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -2 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +1586 2608 M +0 14 V +-31 -14 R +62 0 V +-62 14 R +62 0 V +55 -15 R +0 20 V +-31 -20 R +62 0 V +-62 20 R +62 0 V +55 -22 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 3 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -17 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -6 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 1 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -17 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -27 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +55 -70 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +stroke 2391 2515 M +55 -94 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -100 R +0 16 V +-31 -16 R +62 0 V +-62 16 R +62 0 V +55 -110 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -95 R +0 15 V +-31 -15 R +62 0 V +-62 15 R +62 0 V +55 -102 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -92 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +55 -101 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -94 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -101 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -96 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -100 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -92 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -92 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -85 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -69 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -49 R +0 27 V +-31 -27 R +62 0 V +-62 27 R +62 0 V +55 -58 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -10 R +0 1423 V +3877 1110 M +62 0 V +stroke 3939 1110 M +-62 1423 R +62 0 V +55 -1437 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 1 R +0 13 V +-31 -13 R +62 0 V +-62 13 R +62 0 V +1828 2585 M +0 11 V +-31 -11 R +62 0 V +-62 11 R +62 0 V +55 -11 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -18 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -6 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -11 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -27 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -56 R +0 1 V +-31 -1 R +62 0 V +-62 1 R +62 0 V +55 -88 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -97 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -82 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -86 R +0 12 V +-31 -12 R +62 0 V +-62 12 R +62 0 V +55 -81 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -83 R +0 9 V +-31 -9 R +62 0 V +-62 9 R +62 0 V +55 -90 R +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -87 R +0 6 V +-31 -6 R +62 0 V +-62 6 R +62 0 V +55 -99 R +0 9 V +stroke 3118 1765 M +-31 -9 R +62 0 V +-62 9 R +62 0 V +55 -89 R +0 10 V +-31 -10 R +62 0 V +-62 10 R +62 0 V +55 -99 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -97 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -94 R +0 3 V +-31 -3 R +62 0 V +-62 3 R +62 0 V +55 -98 R +0 2 V +-31 -2 R +62 0 V +-62 2 R +62 0 V +55 -87 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -87 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -81 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -68 R +0 8 V +-31 -8 R +62 0 V +-62 8 R +62 0 V +55 -55 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 72 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -119 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 -47 R +-31 0 R +62 0 V +-62 0 R +62 0 V +55 21 R +-31 0 R +62 0 V +-62 0 R +62 0 V +2071 2557 M +0 7 V +-31 -7 R +62 0 V +-62 7 R +62 0 V +55 -45 R +0 5 V +-31 -5 R +62 0 V +-62 5 R +62 0 V +55 -10 R +0 4 V +-31 -4 R +62 0 V +-62 4 R +62 0 V +55 -43 R +-31 0 R +62 0 V +-62 0 R +62 0 V +141 -156 R +-31 0 R +62 0 V +stroke 2532 2319 M +-62 0 R +62 0 V +1100 2591 Pnt +1186 2593 Pnt +1272 2591 Pnt +1358 2602 Pnt +1444 2595 Pnt +1530 2593 Pnt +1616 2593 Pnt +1702 2595 Pnt +1788 2590 Pnt +1874 2587 Pnt +1960 2589 Pnt +2046 2595 Pnt +2132 2594 Pnt +2218 2580 Pnt +2304 2549 Pnt +2390 2470 Pnt +2476 2365 Pnt +2562 2250 Pnt +2648 2140 Pnt +2734 2040 Pnt +2820 1945 Pnt +2906 1850 Pnt +2992 1766 Pnt +3078 1687 Pnt +3164 1617 Pnt +3250 1570 Pnt +3336 1539 Pnt +3422 2137 Pnt +3508 1527 Pnt +3594 1527 Pnt +1343 2613 Pnt +1429 2618 Pnt +1515 2614 Pnt +1601 2624 Pnt +1687 2614 Pnt +1773 2615 Pnt +1859 2623 Pnt +1945 2622 Pnt +2031 2613 Pnt +2117 2611 Pnt +2203 2600 Pnt +2289 2572 Pnt +2375 2506 Pnt +2461 2403 Pnt +2547 2296 Pnt +2633 2202 Pnt +2719 2113 Pnt +2805 2018 Pnt +2891 1924 Pnt +2977 1828 Pnt +3063 1735 Pnt +3149 1643 Pnt +3235 1557 Pnt +3321 1474 Pnt +3407 1410 Pnt +3493 1362 Pnt +3579 1331 Pnt +3665 1322 Pnt +3751 1315 Pnt +3837 1318 Pnt +1586 2615 Pnt +1672 2617 Pnt +1758 2608 Pnt +1844 2616 Pnt +1930 2605 Pnt +2016 2602 Pnt +2102 2605 Pnt +2188 2593 Pnt +2274 2576 Pnt +2360 2515 Pnt +2446 2424 Pnt +2532 2335 Pnt +2618 2237 Pnt +2704 2154 Pnt +2790 2063 Pnt +2876 1980 Pnt +2962 1890 Pnt +3048 1804 Pnt +3134 1709 Pnt +3220 1619 Pnt +3306 1523 Pnt +3392 1435 Pnt +3478 1347 Pnt +3564 1264 Pnt +3650 1197 Pnt +3736 1162 Pnt +3822 1118 Pnt +3908 2054 Pnt +3994 1096 Pnt +4080 1105 Pnt +1828 2591 Pnt +1914 2590 Pnt +2000 2580 Pnt +2086 2581 Pnt +2172 2577 Pnt +2258 2554 Pnt +2344 2501 Pnt +2430 2418 Pnt +2516 2327 Pnt +2602 2247 Pnt +2688 2168 Pnt +2774 2093 Pnt +2860 2015 Pnt +2946 1933 Pnt +3032 1852 Pnt +3118 1761 Pnt +3204 1681 Pnt +3290 1589 Pnt +3376 1496 Pnt +3462 1406 Pnt +3548 1310 Pnt +3634 1224 Pnt +3720 1137 Pnt +3806 1056 Pnt +3892 992 Pnt +3978 941 Pnt +4064 1013 Pnt +4150 894 Pnt +4236 847 Pnt +4322 868 Pnt +2071 2560 Pnt +2157 2522 Pnt +2243 2516 Pnt +2329 2475 Pnt +2501 2319 Pnt +1.000 UL +LTb +0.00 0.62 0.45 C 2041 2837 M +15 -12 V +37 -30 V +37 -30 V +38 -31 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -31 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -31 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +38 -30 V +37 -30 V +37 -30 V +38 -30 V +37 -30 V +stroke +2.000 UL +LTb +LCb setrgbcolor +1.000 UL +LTb +LCb setrgbcolor +860 2837 N +860 640 L +3699 0 V +0 2197 V +-3699 0 V +Z stroke +1.000 UP +1.000 UL +LTb +LCb setrgbcolor +stroke +grestore +end +showpage + }}% + \put(2709,140){\makebox(0,0){\strut{}$hL^{-\beta\delta/\nu}$}}% + \put(360,1738){% + \special{ps: gsave currentpoint currentpoint translate +630 rotate neg exch neg exch translate}% + \makebox(0,0){\strut{}$\tau L^{-z}$}% + \special{ps: currentpoint grestore moveto}% + }% + \put(4430,440){\makebox(0,0){\strut{}$10^{5}$}}% + \put(4000,440){\makebox(0,0){\strut{}$10^{4}$}}% + \put(3570,440){\makebox(0,0){\strut{}$10^{3}$}}% + \put(3140,440){\makebox(0,0){\strut{}$10^{2}$}}% + \put(2710,440){\makebox(0,0){\strut{}$10^{1}$}}% + \put(2279,440){\makebox(0,0){\strut{}$10^{0}$}}% + \put(1849,440){\makebox(0,0){\strut{}$10^{-1}$}}% + \put(1419,440){\makebox(0,0){\strut{}$10^{-2}$}}% + \put(989,440){\makebox(0,0){\strut{}$10^{-3}$}}% + \put(740,746){\makebox(0,0)[r]{\strut{}$0.1$}}% +\end{picture}% +\endgroup +\endinput diff --git a/figs/fig_correlation_data_collapse-hL.dat b/figs/fig_correlation_data_collapse-hL.dat index edeac8e..fd6cc85 100644 --- a/figs/fig_correlation_data_collapse-hL.dat +++ b/figs/fig_correlation_data_collapse-hL.dat @@ -12,18 +12,18 @@ 16 0.2868963133044901 0.6287962142376368 0.0015874742546196424 16 0.4547000138984317 0.6281007044870663 0.0028411030243642487 16 0.7206509566396645 0.6194074806843178 0.002284206653081905 -16 1.1421547953187856 0.6005833795290582 0.0003360816223467733 -16 1.8101933598375441 0.5552102153110481 0.00042902430497232524 -16 2.868963133044901 0.5004079283004712 0.003454317360817249 -16 4.547000138984715 0.44623957116763613 0.0011712122687413118 -16 7.2065095663966465 0.3998744161125095 0.00003466425520529114 -16 11.421547953187858 0.362269042972682 0.0015513027467385918 -16 18.101933598375602 0.3294778723097587 0.0002942282409112984 -16 28.68963133044861 0.29992664049688833 0.000836983839425126 -16 45.470001389847155 0.27576376985285456 0.0003864875755694997 -16 72.06509566396613 0.254868056433297 0.0009994565216572824 -16 114.21547953187881 0.23770773652299154 0.0003651557948488315 -16 181.01933598375615 0.22688948873452594 0.00037885049672886634 +16 1.1421547953187856 0.6006184183821923 0.00019717506651543908 +16 1.8101933598375441 0.5560556703286 0.0008809927003800477 +16 2.868963133044901 0.49975831305436635 0.0020974832802797415 +16 4.547000138984715 0.4458926608677604 0.000759995273434837 +16 7.2065095663966465 0.4001151038499891 0.000241518371632569 +16 11.421547953187858 0.3628701196069377 0.0010786441447148951 +16 18.101933598375602 0.329982708286965 0.0005326500882052419 +16 28.68963133044861 0.3007015251562208 0.0009132142233643782 +16 45.470001389847155 0.27599567606359066 0.00032182506521583575 +16 72.06509566396613 0.2543855219698998 0.0007522038428913224 +16 114.21547953187881 0.23751514034651802 0.0002855512890815055 +16 181.01933598375615 0.22708888224187773 0.0002959735411563314 16 286.89631330448634 0.2200196611705308 0.00037113377826057664 16 454.70001389847135 0.39891341844240713 0.18170524445924247 16 720.6509566396612 0.21741361507934745 0.00010145294308730456 @@ -43,18 +43,18 @@ 32 1.0523403171229262 0.6146304699498576 0.0029411128080925717 32 1.667847004760328 0.5757182738289438 0.0017891042515727899 32 2.643359363911546 0.519746074483992 0.00044408860844279895 -32 4.18944226109172 0.4670204427314064 0.001365894277698315 -32 6.639818519813105 0.42560291165019704 0.0015942680639849103 -32 10.523403171229258 0.38949748338986645 0.00240252809173605 -32 16.678470047604744 0.3542902285727013 0.0006252639838920104 -32 26.433593639115465 0.32260194466630127 0.00013593524267008217 -32 41.8944226109172 0.29316056163502646 0.0008711389856035645 -32 66.39818519813163 0.2673108060826328 0.0002561260604034554 -32 105.23403171229113 0.2439722389470263 0.00041353736622222465 -32 166.78470047604742 0.22408499741937016 0.0009748993257848032 -32 264.33593639115344 0.20622295754570164 0.0003192570981769976 -32 418.9442261091729 0.19350538611964102 0.0008098416927330293 -32 663.9818519813168 0.18447107210616082 0.0009386296608507316 +32 4.18944226109172 0.46667585126308353 0.0008605999879874023 +32 6.639818519813105 0.42446433865471184 0.0014640965425802752 +32 10.523403171229258 0.38868706435096145 0.0016064949722340582 +32 16.678470047604744 0.35368632421805146 0.0007035757382228955 +32 26.433593639115465 0.3222911790877719 0.0003205225548984619 +32 41.8944226109172 0.292848933548775 0.0005916697628071 +32 66.39818519813163 0.26788829254081953 0.0005961186646367341 +32 105.23403171229113 0.24369169553416772 0.00036838701241900546 +32 166.78470047604742 0.22389083902190107 0.000595404944929348 +32 264.33593639115344 0.20646837308888066 0.00030692640878987555 +32 418.9442261091729 0.19341427336960249 0.00047635706738658406 +32 663.9818519813168 0.18465992621572822 0.0005738824688994034 32 1052.3403171229122 0.1788427351021588 0.0001922835239056425 32 1667.8470047604735 0.17724655446618287 0.00014040830435756158 32 2643.359363911534 0.17612308355378462 0.00008859092338167617 @@ -73,18 +73,18 @@ 64 2.435496171525549 0.5305996387485918 0.0014959230318421645 64 3.8600013025160367 0.48557922109035934 0.0036751832028803612 64 6.117689787248085 0.44044962601219084 0.0018857447024898469 -64 9.695884897399262 0.4056083464205059 0.0030207911414076583 -64 15.36694196877451 0.37034496808040485 0.0011285047217174937 -64 24.35496171525549 0.3410597472226905 0.001984517358535637 -64 38.60001302516036 0.3119763635141043 0.0015139424073597196 -64 61.176897872486215 0.28648946213147886 0.0010216188795957037 -64 96.95884897399264 0.2606451714862178 0.0006902931938541912 -64 153.6694196877451 0.23837531396287037 0.0007131188988310638 -64 243.54961715255703 0.2166054845624757 0.00014709214343654187 -64 386.00013025159825 0.1984726135248987 0.0006330601894378868 -64 611.7689787248622 0.18185027101891627 0.0001649096826504964 -64 969.588489739922 0.16742136264578164 0.0002487946153808817 -64 1536.6941968774543 0.1565974673693383 0.00004431211808114676 +64 9.695884897399262 0.40563766234685533 0.0017443009479081834 +64 15.36694196877451 0.37085340478832524 0.0008264475310672823 +64 24.35496171525549 0.34124393616706267 0.001160472008568937 +64 38.60001302516036 0.3119031918247171 0.0008771324303633754 +64 61.176897872486215 0.28632943948682943 0.0006111537928509062 +64 96.95884897399264 0.2605250393511181 0.0004162530813228257 +64 153.6694196877451 0.23782958962548653 0.0006836138581365798 +64 243.54961715255703 0.21664193300132426 0.00009241494241417804 +64 386.00013025159825 0.19947059323368097 0.0010628037920008965 +64 611.7689787248622 0.1817052030857358 0.00017352167889887098 +64 969.588489739922 0.1674466007578578 0.0001458419778916332 +64 1536.6941968774543 0.15661382019143888 0.00003036340001161581 64 2435.496171525572 0.15118568965616075 0.001975073557340154 64 3860.001302515986 0.1448337200806049 0.0002595070942839617 64 6117.68978724862 0.367371300110232 0.22380440827218762 @@ -100,28 +100,28 @@ 128 1.415854720473395 0.5725642835937848 0.00014026941966266887 128 2.2439785079905805 0.5272959382383 0.0025756573485712577 128 3.5564662613452867 0.48174955931471536 0.0005001598678584231 -128 5.636619166821931 0.44491238376473485 0.0003661776347952949 -128 8.933439345993705 0.4113572807767111 0.0024356478537859434 -128 14.158547204730013 0.3817803172768068 0.00008389740373021536 -128 22.439785079913676 0.353218055274667 0.0015682350992511638 -128 35.56466261344893 0.3254983052094862 0.0011331832832819908 -128 56.36619166823111 0.3003941905165709 0.0008062274638504784 -128 89.33439345993705 0.27427695641426886 0.001109371145218918 -128 141.5854720473001 0.25338014351418314 0.0012608266730329276 -128 224.39785079915646 0.2312919586007066 0.0005096147290300223 -128 355.6466261344894 0.21077098134013394 0.00039002269998274214 -128 563.6619166823111 0.1927527400715816 0.00027052134638859415 -128 893.3439345993783 0.17525595209598482 0.00013628284241916366 +128 5.636619166821931 0.4438388073898844 0.0010941944918729596 +128 8.933439345993705 0.41092720925776854 0.0014705174828234917 +128 14.158547204730013 0.38121192764431805 0.0005704498509406503 +128 22.439785079913676 0.3529764733474239 0.0009370960123763205 +128 35.56466261344893 0.32488640337259517 0.0008958005595521174 +128 56.36619166823111 0.3006193023569792 0.0005170521397121607 +128 89.33439345993705 0.27456708026595844 0.0007031405468330241 +128 141.5854720473001 0.2536271681032865 0.0007687104661479951 +128 224.39785079915646 0.23100431012855005 0.000411473815548584 +128 355.6466261344894 0.21147268527569854 0.00073694933024878 +128 563.6619166823111 0.1931657687811568 0.0004415729247632602 +128 893.3439345993783 0.17567152403394282 0.00042295512764053896 128 1415.8547204729814 0.16083831181683386 0.000038393852729745597 128 2243.9785079915646 0.14755136241680297 0.000010427848840532938 128 3556.4662613448777 0.13613167291690653 0.00003824559764080415 128 5636.619166823122 0.12770211002072657 0.000480391007107131 128 8933.43934599379 0.12140474759127326 0.000021786704140921986 -128 14158.547204729826 0.1304089776503813 0. -128 22439.785079915637 0.11590410749694247 0. -128 35564.66261344877 0.11053344320283695 0. -128 56366.191668231244 0.11290716591905191 0. -128 89334.39345993791 42.3724331482014 0. +128 14158.547204729826 0.12450333905336947 0.005905638597011847 +128 22439.785079915637 0.11616358912822208 0.0002594816312796177 +128 35564.66261344877 0.11074819112465964 0.00021474792182268887 +128 56366.191668231244 0.11173511904556105 0.0011720468734908519 +128 89334.39345993791 43.80504345251086 1.432610304309456 256 0.32767999999811953 0.6075086695911571 0.002178147172579081 256 0.519337801301477 0.5847006386014801 0.0013928429465341339 256 0.8230949458839193 0.5812979760560055 0.0009786732427118916 diff --git a/makefile b/makefile index 2a9a11a..364b659 100644 --- a/makefile +++ b/makefile @@ -5,7 +5,7 @@ DATA=fig_correlation_data_collapse-hL all: ${DOC}.pdf -%.tex: figs/%.gplot ${DATA:%=data/%.dat} ${FUNCS:%=figs/%.dat} +%.tex: figs/%.gplot ${DATA:%=figs/%.dat} gnuplot $< > $@ ${DOC}.pdf: ${DOC}.tex ${DOC}.bib ${FIGS:%=%.tex} diff --git a/monte-carlo.bib b/monte-carlo.bib index 4413e2e..d527915 100644 --- a/monte-carlo.bib +++ b/monte-carlo.bib @@ -1,4 +1,59 @@ +@article{alexandrowicz1989swendsen, + title={Swendsen-Wang simulation of Ising spins and a precise definition of critical clusters}, + author={Alexandrowicz, Z}, + journal={Physica A: Statistical Mechanics and its Applications}, + volume={160}, + number={3}, + pages={310--320}, + year={1989}, + publisher={Elsevier} +} + +@article{coniglio1980clusters, + title={Clusters and Ising critical droplets: a renormalisation group approach}, + author={Coniglio, A and Klein, W}, + journal={Journal of Physics A: Mathematical and General}, + volume={13}, + number={8}, + pages={2775}, + year={1980}, + publisher={IOP Publishing} +} + +@article{coniglio1989exact, + title={Exact relations between droplets and thermal fluctuations in external field}, + author={Coniglio, A and de Liberto, F and Monroy, G and Peruggi, F}, + journal={Journal of Physics A: Mathematical and General}, + volume={22}, + number={17}, + pages={L837}, + year={1989}, + publisher={IOP Publishing} +} + +@article{destri1992swendsen, + title={Swendsen-Wang Monte Carlo study of the Ising model with external field}, + author={Destri, C and Di Renzo, F and Onofri, E and Rossi, P and Tecchiolli, GP}, + journal={Physics Letters B}, + volume={278}, + number={3}, + pages={311--316}, + year={1992}, + publisher={Elsevier} +} + +@article{dotsenko1991cluster, + title={Cluster Monte Carlo algorithms for random Ising models}, + author={Dotsenko, Vl S and Selke, W and Talapov, AL}, + journal={Physica A: Statistical Mechanics and its Applications}, + volume={170}, + number={2}, + pages={278--281}, + year={1991}, + publisher={Elsevier} +} + @article{du2006dynamic, title={Dynamic critical exponents for Swendsen--Wang and Wolff algorithms obtained by a nonequilibrium relaxation method}, author={Du, Jianqing and Zheng, Bo and Wang, Jian-Sheng}, @@ -10,6 +65,39 @@ publisher={IOP Publishing} } +@article{hastings1970monte, + title={Monte Carlo sampling methods using Markov chains and their applications}, + author={Hastings, W Keith}, + journal={Biometrika}, + volume={57}, + number={1}, + pages={97--109}, + year={1970}, + publisher={Biometrika Trust} +} + +@article{janke1998nonlocal, + title={Nonlocal Monte Carlo algorithms for statistical physics applications}, + author={Janke, Wolfhard}, + journal={Mathematics and computers in simulation}, + volume={47}, + number={2}, + pages={329--346}, + year={1998}, + publisher={Elsevier} +} + +@article{lauwers1989critical, + title={The critical 2D Ising model in a magnetic field. A Monte Carlo study using a Swendesen-Wang algorithm}, + author={Lauwers, Paul G and Rittenberg, Vladimir}, + journal={Physics Letters B}, + volume={233}, + number={1-2}, + pages={197--200}, + year={1989}, + publisher={Elsevier} +} + @article{liu2014dynamic, title={Dynamic scaling at classical phase transitions approached through nonequilibrium quenching}, author={Liu, Cheng-Wei and Polkovnikov, Anatoli and Sandvik, Anders W}, @@ -21,6 +109,61 @@ publisher={APS} } +@article{metropolis1953equation, + title={Equation of state calculations by fast computing machines}, + author={Metropolis, Nicholas and Rosenbluth, Arianna W and Rosenbluth, Marshall N and Teller, Augusta H and Teller, Edward}, + journal={The journal of chemical physics}, + volume={21}, + number={6}, + pages={1087--1092}, + year={1953}, + publisher={AIP} +} + +@article{ray1990metastability, + title={Metastability and nucleation in Ising models with Swendsen-Wang dynamics}, + author={Ray, TS and Wang, Jian-Sheng}, + journal={Physica A: Statistical Mechanics and its Applications}, + volume={167}, + number={3}, + pages={580--588}, + year={1990}, + publisher={Elsevier} +} + +@article{swendsen1987nonuniversal, + title={Nonuniversal critical dynamics in Monte Carlo simulations}, + author={Swendsen, Robert H and Wang, Jian-Sheng}, + journal={Physical review letters}, + volume={58}, + number={2}, + pages={86}, + year={1987}, + publisher={APS} +} + +@article{wang1989clusters, + title={Clusters in the three-dimensional Ising model with a magnetic field}, + author={Wang, Jian-Sheng}, + journal={Physica A: Statistical Mechanics and its Applications}, + volume={161}, + number={2}, + pages={249--268}, + year={1989}, + publisher={Elsevier} +} + +@article{wang1990cluster, + title={Cluster monte carlo algorithms}, + author={Wang, Jian-Sheng and Swendsen, Robert H}, + journal={Physica A: Statistical Mechanics and its Applications}, + volume={167}, + number={3}, + pages={565--579}, + year={1990}, + publisher={Elsevier} +} + @article{wolff1989collective, title={Collective Monte Carlo updating for spin systems}, author={Wolff, Ulli}, @@ -43,3 +186,13 @@ publisher={Elsevier} } +@article{wolff1990critical, + title={Critical slowing down}, + author={Wolff, Ulli}, + journal={Nuclear Physics B-Proceedings Supplements}, + volume={17}, + pages={93--102}, + year={1990}, + publisher={Elsevier} +} + diff --git a/monte-carlo.tex b/monte-carlo.tex index fe2a315..2490b31 100644 --- a/monte-carlo.tex +++ b/monte-carlo.tex @@ -82,7 +82,7 @@ \begin{document} -\title{An efficient cluster algorithm for the Ising model in an external field} +\title{Efficiently sampling Ising states in an external field} \author{Jaron Kent-Dobias} \author{James P.~Sethna} \affiliation{Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY, USA} @@ -90,15 +90,43 @@ \date\today \begin{abstract} - An abstract. + We introduce an extension of the Wolff algorithm that preforms efficiently + in an external magnetic field. Near the Ising critical point, the + correlation time of our algorithm has a conventional scaling form that + reduces to that of the Wolff algorithm at zero field. As an application, we + directly measure scaling functions of observables in the metastable state of + the 2D Ising model. \end{abstract} \maketitle -\section{Introduction} - The Ising model is a simple model of a magnet comprised of locally interacting -spins. +spins. Like most large thermal systems, computation of its properties cannot +be carried out explicitly and is preformed using Monte Carlo methods. Near its +continuous phase transition, divergent correlation length leads to divergent +correlation time in any locally-updating algorithm, hampering computation. +At zero external field, this was largely alleviated by cluster algorithms, +like the Wolff algorithm, whose dynamics are nonlocal and each step flips +groups of spins whose size diverges with the correlation length. However, the +Wolff algorithm only works at zero field. We describe an extension of this +algorithm that works in arbitrary external field while preserving the Wolff +algorithm's small dynamic exponent. + +The Wolff algorithm works by first choosing a random spin and adding it to an +empty cluster. Every neighbor of that spin pointed in the same direction as +the spin is added to the cluster with probability $1-e^{-2\beta J}$, where +$\beta=1/T$ and $J$ is the coupling between sites. This process is iterated +again for neighbors of every spin added to the cluster. When all sites +surrounding the cluster have been exhausted, the cluster is flipped. Our +algorithm is a simple extension of this. An extra spin is introduced (often +referred to as a ``ghost spin'') that couples to all others with coupling $H$. +The traditional Wolff algorithm is then preformed on this larger lattice exactly as described above, +with the extra spin treated no differently from any others. Observables in the +original system can be exactly estimated on the new one using a simple +mapping. As an application, we use our algorithm to measure critical scaling functions +of the 2D Ising model in its metastable phase. + +\section{Introduction} Consider an undirected graph $G=(V,E)$ describing a system of interacting spins. The set of vertices $V=\{1,\ldots,N\}$ enumerates the sites of the network, and the @@ -147,7 +175,8 @@ $S^n$ according to the Boltzmann distribution $e^{-\beta\H(s)}$, so that averages of observables made using their samples asymptotically approach the true expected value. -The Metropolis--Hastings algorithm is very popular for systems in statistical +The Metropolis--Hastings algorithm +\cite{metropolis1953equation,hastings1970monte} is very popular for systems in statistical physics. A starting state $s\in S^n$ is randomly perturbed to the state $s'$, usually by flipping one spin. The change in energy $\Delta\H=\H(s')-\H(s)$ due to the perturbation is then computed. If the change is negative the perturbed @@ -177,7 +206,7 @@ take many perturbations to move between in configuration space. \end{algorithm} \end{figure} -The Wolff algorithm solves many of these problems, but only at zero external +The Wolff algorithm \cite{wolff1989collective} solves many of these problems, but only at zero external field, $H=0$. This algorithm solves the problem of critical slowing-down by flipping carefully-constructed clusters of spins at once in a way that samples high-correlated states quickly while also always accepting prospective states. @@ -196,6 +225,9 @@ cluster acceptance or rejection across the spins in the cluster as they are added. In this version, the cluster is abandoned with probability $1-e^{-\beta H}$ every time a spin is added to it. +$z=0.29(1)$ \cite{wolff1989comparison,liu2014dynamic} $z=0.35(1)$ for +Swendsen--Wang \cite{swendsen1987nonuniversal} + \begin{figure} \begin{algorithm}[H] \begin{algorithmic} @@ -491,7 +523,7 @@ a given average if the associated state is in the reduced space of interest. Thanks! \end{acknowledgments} -%\bibliography{monte-carlo} +\bibliography{monte-carlo} \end{document} diff --git a/refs/alexandrowicz1989swendsen.pdf b/refs/alexandrowicz1989swendsen.pdf new file mode 100644 index 0000000..08c1c5e Binary files /dev/null and b/refs/alexandrowicz1989swendsen.pdf differ diff --git a/refs/coniglio1980clusters.pdf b/refs/coniglio1980clusters.pdf new file mode 100644 index 0000000..9ab9c4a Binary files /dev/null and b/refs/coniglio1980clusters.pdf differ diff --git a/refs/coniglio1989exact.pdf b/refs/coniglio1989exact.pdf new file mode 100644 index 0000000..1b145fa Binary files /dev/null and b/refs/coniglio1989exact.pdf differ diff --git a/refs/destri1992swendsen.pdf b/refs/destri1992swendsen.pdf new file mode 100644 index 0000000..591c8f6 Binary files /dev/null and b/refs/destri1992swendsen.pdf differ diff --git a/refs/dotsenko1991cluster.pdf b/refs/dotsenko1991cluster.pdf new file mode 100644 index 0000000..599984e Binary files /dev/null and b/refs/dotsenko1991cluster.pdf differ diff --git a/refs/du2006dynamic.pdf b/refs/du2006dynamic.pdf new file mode 100644 index 0000000..685ba66 Binary files /dev/null and b/refs/du2006dynamic.pdf differ diff --git a/refs/janke1998nonlocal.pdf b/refs/janke1998nonlocal.pdf new file mode 100644 index 0000000..f2666e0 Binary files /dev/null and b/refs/janke1998nonlocal.pdf differ diff --git a/refs/lauwers1989critical.pdf b/refs/lauwers1989critical.pdf new file mode 100644 index 0000000..09063c5 Binary files /dev/null and b/refs/lauwers1989critical.pdf differ diff --git a/refs/liu2014dynamic.pdf b/refs/liu2014dynamic.pdf new file mode 100644 index 0000000..af17e6a Binary files /dev/null and b/refs/liu2014dynamic.pdf differ diff --git a/refs/ray1990metastability.pdf b/refs/ray1990metastability.pdf new file mode 100644 index 0000000..f7f1e2f Binary files /dev/null and b/refs/ray1990metastability.pdf differ diff --git a/refs/swendsen1987nonuniversal.pdf b/refs/swendsen1987nonuniversal.pdf new file mode 100644 index 0000000..6c359eb Binary files /dev/null and b/refs/swendsen1987nonuniversal.pdf differ diff --git a/refs/wang1989clusters.pdf b/refs/wang1989clusters.pdf new file mode 100644 index 0000000..c06de0b Binary files /dev/null and b/refs/wang1989clusters.pdf differ diff --git a/refs/wang1990cluster.pdf b/refs/wang1990cluster.pdf new file mode 100644 index 0000000..d308d39 Binary files /dev/null and b/refs/wang1990cluster.pdf differ diff --git a/refs/wolff1989collective.pdf b/refs/wolff1989collective.pdf new file mode 100644 index 0000000..6407658 Binary files /dev/null and b/refs/wolff1989collective.pdf differ diff --git a/refs/wolff1989comparison.pdf b/refs/wolff1989comparison.pdf new file mode 100644 index 0000000..4465d2e Binary files /dev/null and b/refs/wolff1989comparison.pdf differ diff --git a/refs/wolff1990critical.pdf b/refs/wolff1990critical.pdf new file mode 100644 index 0000000..6526c0f Binary files /dev/null and b/refs/wolff1990critical.pdf differ -- cgit v1.2.3-70-g09d2