1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
#!/usr/bin/perl -w
# Copyright 2010, 2011, 2012, 2013 Kevin Ryde
# This file is part of Math-PlanePath.
#
# Math-PlanePath is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by the Free
# Software Foundation; either version 3, or (at your option) any later
# version.
#
# Math-PlanePath is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
# or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
# for more details.
#
# You should have received a copy of the GNU General Public License along
# with Math-PlanePath. If not, see <http://www.gnu.org/licenses/>.
use 5.004;
use strict;
use Math::PlanePath::Base::Digits 'round_down_pow';
{
require Image::Base::GD;
my $width = 810;
my $height = 810;
my $image = Image::Base::GD->new (-width => $width, -height => $height);
$image->rectangle (0,0, $width-1,$height-1, 'black');
my $foreground = 'white';
# *---------*
# | |
# *----* . |
# |
# *----* *----*
# | |
# * *
my $recurse;
$recurse = sub {
my ($x,$y, $dx,$dy, $level) = @_;
if (--$level < 0) {
$image->line($x,$y, $x+$dx,$y+$dy, $foreground);
$x += $dx;
$y += $dy;
($dx,$dy) = (-$dy,$dx); # rotate +90
$image->line($x,$y, $x+$dx,$y+$dy, $foreground);
} else {
$dx /= 2;
$dy /= 2;
$image->line($x,$y, $x+$dx,$y+$dy, $foreground);
$x += $dx;
$y += $dy;
($dx,$dy) = ($dy,-$dx); # rotate -90
$recurse->($x,$y, $dx,$dy, $level);
$x += $dx;
$y += $dy;
($dx,$dy) = (-$dy,$dx); # rotate +90
$x += $dx;
$y += $dy;
$recurse->($x,$y, $dx,$dy, $level);
$x += $dx;
$y += $dy;
($dx,$dy) = (-$dy,$dx); # rotate +90
$x += $dx;
$y += $dy;
$recurse->($x,$y, $dx,$dy, $level);
$x += $dx;
$y += $dy;
($dx,$dy) = (-$dy,$dx); # rotate +90
$x += $dx;
$y += $dy;
($dx,$dy) = ($dy,-$dx); # rotate -90
$image->line($x,$y, $x+$dx,$y+$dy, $foreground);
}
};
my $scale = 2;
my ($pow,$exp) = round_down_pow($height/$scale, 2);
foreach my $level (0 .. $exp) {
my $len = 2**$level * $scale;
$recurse->(0, $height-1 - $len, $len,0, $level);
}
$image->save('/tmp/x.png');
system('xzgv /tmp/x.png');
exit 0;
}
|
