.\" Automatically generated by Pod::Man 2.28 (Pod::Simple 3.28) .\" .\" Standard preamble: .\" ======================================================================== .de Sp \" Vertical space (when we can't use .PP) .if t .sp .5v .if n .sp .. .de Vb \" Begin verbatim text .ft CW .nf .ne \\$1 .. .de Ve \" End verbatim text .ft R .fi .. .\" Set up some character translations and predefined strings. \*(-- will .\" give an unbreakable dash, \*(PI will give pi, \*(L" will give a left .\" double quote, and \*(R" will give a right double quote. \*(C+ will .\" give a nicer C++. Capital omega is used to do unbreakable dashes and .\" therefore won't be available. \*(C` and \*(C' expand to `' in nroff, .\" nothing in troff, for use with C<>. .tr \(*W- .ds C+ C\v'-.1v'\h'-1p'\s-2+\h'-1p'+\s0\v'.1v'\h'-1p' .ie n \{\ . ds -- \(*W- . ds PI pi . if (\n(.H=4u)&(1m=24u) .ds -- \(*W\h'-12u'\(*W\h'-12u'-\" diablo 10 pitch . if (\n(.H=4u)&(1m=20u) .ds -- \(*W\h'-12u'\(*W\h'-8u'-\" diablo 12 pitch . ds L" "" . ds R" "" . ds C` "" . ds C' "" 'br\} .el\{\ . ds -- \|\(em\| . ds PI \(*p . ds L" `` . ds R" '' . ds C` . ds C' 'br\} .\" .\" Escape single quotes in literal strings from groff's Unicode transform. .ie \n(.g .ds Aq \(aq .el .ds Aq ' .\" .\" If the F register is turned on, we'll generate index entries on stderr for .\" titles (.TH), headers (.SH), subsections (.SS), items (.Ip), and index .\" entries marked with X<> in POD. Of course, you'll have to process the .\" output yourself in some meaningful fashion. .\" .\" Avoid warning from groff about undefined register 'F'. .de IX .. .nr rF 0 .if \n(.g .if rF .nr rF 1 .if (\n(rF:(\n(.g==0)) \{ . if \nF \{ . de IX . tm Index:\\$1\t\\n%\t"\\$2" .. . if !\nF==2 \{ . nr % 0 . nr F 2 . \} . \} .\} .rr rF .\" .\" Accent mark definitions (@(#)ms.acc 1.5 88/02/08 SMI; from UCB 4.2). .\" Fear. Run. Save yourself. No user-serviceable parts. . \" fudge factors for nroff and troff .if n \{\ . ds #H 0 . ds #V .8m . ds #F .3m . ds #[ \f1 . ds #] \fP .\} .if t \{\ . ds #H ((1u-(\\\\n(.fu%2u))*.13m) . ds #V .6m . ds #F 0 . ds #[ \& . ds #] \& .\} . \" simple accents for nroff and troff .if n \{\ . ds ' \& . ds ` \& . ds ^ \& . ds , \& . ds ~ ~ . ds / .\} .if t \{\ . ds ' \\k:\h'-(\\n(.wu*8/10-\*(#H)'\'\h"|\\n:u" . ds ` \\k:\h'-(\\n(.wu*8/10-\*(#H)'\`\h'|\\n:u' . ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'^\h'|\\n:u' . ds , \\k:\h'-(\\n(.wu*8/10)',\h'|\\n:u' . ds ~ \\k:\h'-(\\n(.wu-\*(#H-.1m)'~\h'|\\n:u' . ds / \\k:\h'-(\\n(.wu*8/10-\*(#H)'\z\(sl\h'|\\n:u' .\} . \" troff and (daisy-wheel) nroff accents .ds : \\k:\h'-(\\n(.wu*8/10-\*(#H+.1m+\*(#F)'\v'-\*(#V'\z.\h'.2m+\*(#F'.\h'|\\n:u'\v'\*(#V' .ds 8 \h'\*(#H'\(*b\h'-\*(#H' .ds o \\k:\h'-(\\n(.wu+\w'\(de'u-\*(#H)/2u'\v'-.3n'\*(#[\z\(de\v'.3n'\h'|\\n:u'\*(#] .ds d- \h'\*(#H'\(pd\h'-\w'~'u'\v'-.25m'\f2\(hy\fP\v'.25m'\h'-\*(#H' .ds D- D\\k:\h'-\w'D'u'\v'-.11m'\z\(hy\v'.11m'\h'|\\n:u' .ds th \*(#[\v'.3m'\s+1I\s-1\v'-.3m'\h'-(\w'I'u*2/3)'\s-1o\s+1\*(#] .ds Th \*(#[\s+2I\s-2\h'-\w'I'u*3/5'\v'-.3m'o\v'.3m'\*(#] .ds ae a\h'-(\w'a'u*4/10)'e .ds Ae A\h'-(\w'A'u*4/10)'E . \" corrections for vroff .if v .ds ~ \\k:\h'-(\\n(.wu*9/10-\*(#H)'\s-2\u~\d\s+2\h'|\\n:u' .if v .ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'\v'-.4m'^\v'.4m'\h'|\\n:u' . \" for low resolution devices (crt and lpr) .if \n(.H>23 .if \n(.V>19 \ \{\ . ds : e . ds 8 ss . ds o a . ds d- d\h'-1'\(ga . ds D- D\h'-1'\(hy . ds th \o'bp' . ds Th \o'LP' . ds ae ae . ds Ae AE .\} .rm #[ #] #H #V #F C .\" ======================================================================== .\" .IX Title "Math::PlanePath::R5DragonMidpoint 3pm" .TH Math::PlanePath::R5DragonMidpoint 3pm "2014-08-26" "perl v5.20.1" "User Contributed Perl Documentation" .\" For nroff, turn off justification. Always turn off hyphenation; it makes .\" way too many mistakes in technical documents. .if n .ad l .nh .SH "NAME" Math::PlanePath::R5DragonMidpoint \-\- R5 dragon curve midpoints .SH "SYNOPSIS" .IX Header "SYNOPSIS" .Vb 3 \& use Math::PlanePath::R5DragonMidpoint; \& my $path = Math::PlanePath::R5DragonMidpoint\->new; \& my ($x, $y) = $path\->n_to_xy (123); .Ve .SH "DESCRIPTION" .IX Header "DESCRIPTION" This is midpoints of the R5 dragon curve by Jorg Arndt, .IX Xref "Arndt, Jorg" .PP .Vb 10 \& 31\-\-30 11 \& | | \& 32 29 10 \& | | \& 51\-\-50 35\-\-34\-\-33 28\-\-27\-\-26 9 \& | | | | \& 52 49 36\-\-37\-\-38 23\-\-24\-\-25 8 \& | | | | \& 55\-\-54\-\-53 48\-\-47\-\-46 41\-\-40\-\-39 22 7 \& | | | | \& 56\-\-57\-\-58 63\-\-64 45 42 19\-\-20\-\-21 6 \& | | | | | | \& 81\-\-80 59 62 65 44\-\-43 18\-\-17\-\-16 11\-\-10 5 \& | | | | | | | | \& 82 79 60\-\-61 66\-\-67\-\-68 15 12 9 4 \& | | | | | | \& ..\-83 78\-\-77\-\-76 71\-\-70\-\-69 14\-\-13 8\-\- 7\-\- 6 3 \& | | | \& 75 72 3\-\- 4\-\- 5 2 \& | | | \& 74\-\-73 2 1 \& | \& 0\-\- 1 <\- Y=0 \& \& ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ \& \-10 \-9 \-8 \-7 \-6 \-5 \-4 \-3 \-2 \-1 X=0 1 2 3 .Ve .PP The points are the middle of each edge of the \f(CW\*(C`R5DragonCurve\*(C'\fR, rotated \-45 degrees, shrunk by \fIsqrt\fR\|(2). and shifted to the origin. .PP .Vb 10 \& *\-\-11\-\-* *\-\-7\-\-* R5DragonCurve \& | | | | and its midpoints \& 12 10 8 6 \& | | | | \& *\-\-17\-\-*\-\-13\-\-*\-\-9\-\-*\-\-5\-\-* \& | | | | \& 18 16 14 4 \& | | | | \& ..\-* *\-\-15\-\-* *\-\-3\-\-* \& | \& 2 \& | \& +\-\-1\-\-* .Ve .SS "Arms" .IX Subsection "Arms" Multiple copies of the curve can be selected, each advancing successively. Like the main \f(CW\*(C`R5DragonCurve\*(C'\fR this midpoint curve covers 1/4 of the plane and 4 arms rotated by 0, 90, 180, 270 degrees mesh together perfectly. With 4 arms all integer X,Y points are visited. .PP \&\f(CW\*(C`arms => 4\*(C'\fR begins as follows. N=0,4,8,12,16,etc is the first arm (the same shape as the plain curve above), then N=1,5,9,13,17 the second, N=2,6,10,14 the third, etc. .PP .Vb 10 \& arms=>4 76\-\-80\-... 6 \& | \& 72\-\-68\-\-64 44\-\-40 5 \& | | | \& 25\-\-21 60 48 36 4 \& | | | | | \& 29 17 56\-\-52 32\-\-28\-\-24 75\-\-79 3 \& | | | | | \& 41\-\-37\-\-33 13\-\- 9\-\- 5 12\-\-16\-\-20 71 83 2 \& | | | | | \& 45\-\-49\-\-53 6\-\- 2 1 8 59\-\-63\-\-67 ... 1 \& | | | | \& ... 65\-\-61\-\-57 10 3 0\-\- 4 55\-\-51\-\-47 <\- Y=0 \& | | | | | \& 81 69 22\-\-18\-\-14 7\-\-11\-\-15 35\-\-39\-\-43 \-1 \& | | | | | \& 77\-\-73 26\-\-30\-\-34 54\-\-58 19 31 \-2 \& | | | | | \& 38 50 62 23\-\-27 \-3 \& | | | \& 42\-\-46 66\-\-70\-\-74 \-4 \& | \& ...\-82\-\-78 \-5 \& \& ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ \& \-6 \-5 \-4 \-3 \-2 \-1 X=0 1 2 3 4 5 .Ve .SH "FUNCTIONS" .IX Header "FUNCTIONS" See \*(L"\s-1FUNCTIONS\*(R"\s0 in Math::PlanePath for behaviour common to all path classes. .ie n .IP """$path = Math::PlanePath::R5DragonMidpoint\->new ()""" 4 .el .IP "\f(CW$path = Math::PlanePath::R5DragonMidpoint\->new ()\fR" 4 .IX Item "$path = Math::PlanePath::R5DragonMidpoint->new ()" Create and return a new path object. .ie n .IP """($x,$y) = $path\->n_to_xy ($n)""" 4 .el .IP "\f(CW($x,$y) = $path\->n_to_xy ($n)\fR" 4 .IX Item "($x,$y) = $path->n_to_xy ($n)" Return the X,Y coordinates of point number \f(CW$n\fR on the path. Points begin at 0 and if \f(CW\*(C`$n < 0\*(C'\fR then the return is an empty list. .Sp Fractional positions give an X,Y position along a straight line between the integer positions. .ie n .IP """$n = $path\->n_start()""" 4 .el .IP "\f(CW$n = $path\->n_start()\fR" 4 .IX Item "$n = $path->n_start()" Return 0, the first N in the path. .SS "Level Methods" .IX Subsection "Level Methods" .ie n .IP """($n_lo, $n_hi) = $path\->level_to_n_range($level)""" 4 .el .IP "\f(CW($n_lo, $n_hi) = $path\->level_to_n_range($level)\fR" 4 .IX Item "($n_lo, $n_hi) = $path->level_to_n_range($level)" Return \f(CW\*(C`(0, 5**$level \- 1)\*(C'\fR, or for multiple arms return \f(CW\*(C`(0, $arms * 5**$level \- 1)\*(C'\fR. .Sp There are 5^level segments comprising the curve, or arms*5^level when multiple arms, numbered starting from 0. .SH "FORMULAS" .IX Header "FORMULAS" .SS "X,Y to N" .IX Subsection "X,Y to N" An X,Y point can be turned into N by dividing out digits of a complex base 1+2i. At each step the low base\-5 digit is formed from X,Y and an adjustment applied to move X,Y to a multiple of 1+2i ready to divide out. .PP A 10x10 table is used for the digit and adjustments, indexed by Xmod10 and Ymod10. There's probably an a*X+b*Y mod 5 or mod 20 for a smaller table. But in any case once the adjustment is found the result is .PP .Vb 3 \& Ndigit = digit_table[X mod 10, Y mod 10] # low to high \& Xm = X + Xadj_table [X mod 10, Y mod 10] \& Ym = Y + Yadj_table [X mod 10, Y mod 10] \& \& new X,Y = (Xm,Ym) / (1+2i) \& = (Xm,Ym) * (1\-2i) / 5 \& = ((Xm+2*Ym)/5, (Ym\-2*Xm)/5) .Ve .PP These X,Y reductions eventually reach one of the starting points for the four arms .PP .Vb 6 \& X,Y endpoint Arm +\-\-\-+\-\-\-+ \& \-\-\-\-\-\-\-\-\-\-\-\- \-\-\- | 2 | 1 | Y=1 \& 0, 0 0 +\-\-\-+\-\-\-+ \& 0, 1 1 | 3 | 0 | Y=0 \& \-1, 1 2 +\-\-\-+\-\-\-+ \& \-1, 0 3 X=\-1 X=0 .Ve .PP For arms 1 and 3 the digits must be flipped 4\-digit, so 0,1,2,3,4 \-> 4,3,2,1,0. The arm number and hence whether this flip is needed is not known until reaching the endpoint. .PP .Vb 2 \& if arm odd \& then N = 5^numdigits \- 1 \- N .Ve .PP If only some of the arms are of interest then reaching one of the other arm numbers means the original X,Y was outside the desired curve. .SH "SEE ALSO" .IX Header "SEE ALSO" Math::PlanePath, Math::PlanePath::R5DragonCurve .PP Math::PlanePath::DragonMidpoint, Math::PlanePath::TerdragonMidpoint .SH "HOME PAGE" .IX Header "HOME PAGE" .SH "LICENSE" .IX Header "LICENSE" Copyright 2012, 2013, 2014 Kevin Ryde .PP This file is part of Math-PlanePath. .PP Math-PlanePath is free software; you can redistribute it and/or modify it under the terms of the \s-1GNU\s0 General Public License as published by the Free Software Foundation; either version 3, or (at your option) any later version. .PP Math-PlanePath is distributed in the hope that it will be useful, but \&\s-1WITHOUT ANY WARRANTY\s0; without even the implied warranty of \s-1MERCHANTABILITY\s0 or \s-1FITNESS FOR A PARTICULAR PURPOSE. \s0 See the \s-1GNU\s0 General Public License for more details. .PP You should have received a copy of the \s-1GNU\s0 General Public License along with Math-PlanePath. If not, see .