projects/doodle: doodle/tk/geometry.d annotate

annotate doodle/tk/geometry.d @ 79:535bae7a7305

Checkpoint

author	"David Bryant <bagnose@gmail.com>"
date	Sun, 15 Aug 2010 23:18:05 +0930
parents	024a5608087f
children	a274d16ab6ce

rev	line source
28 1754cb773d41 Part-way through getting to compile with configure/builder. Graham St Jack <graham.stjack@internode.on.net> parents: 25 diff changeset	1 module doodle.tk.geometry;
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	2
16 9e63308b749c * Fix up public/private includes David Bryant <daveb@acres.com.au> parents: 15 diff changeset	3 private {
9e63308b749c * Fix up public/private includes David Bryant <daveb@acres.com.au> parents: 15 diff changeset	4 import std.stdio;
9e63308b749c * Fix up public/private includes David Bryant <daveb@acres.com.au> parents: 15 diff changeset	5 import std.math;
67 31d10176415d Checkpoint daveb parents: 66 diff changeset	6 import doodle.core.misc;
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	7 import doodle.core.logging;
16 9e63308b749c * Fix up public/private includes David Bryant <daveb@acres.com.au> parents: 15 diff changeset	8 }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	9
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	10 // In doodle x and y increase right/east and up/north respectively.
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	11
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	12 // TODO explain the strategy for ensuring numerical stability
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	13 // and the division of responsibility between users of these
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	14 // types and the types themselves.
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	15 //
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	16 // Explain how numerical instability is handled. The current policy
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	17 // is to correct bad user input (eg a gradient with miniscule length) and
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	18 // print warnings rather than have assertions that cause crashes.
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	19 //
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	20 // There are no mutating operations other than opAssign
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	21
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	22 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	23 // A location in 2D space
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	24 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	25
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	26 struct Point {
51 0eaf39fda206 First cut at palette class daveb parents: 48 diff changeset	27 static immutable Point DEFAULT = Point(0.0, 0.0);
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	28
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	29 this(in double x, in double y) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	30 _x = x;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	31 _y = y;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	32 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	33
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	34 Point opAdd(in Vector v) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	35 return Point(_x + v._x, _y + v._y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	36 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	37
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	38 Point opSub(in Vector v) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	39 return Point(_x - v._x, _y - v._y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	40 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	41
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	42 Vector opSub(in Point p) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	43 return Vector(_x - p._x, _y - p._y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	44 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	45
62 6c3993f4c3eb Checkpoint daveb parents: 51 diff changeset	46 string toString() {
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	47 return std.string.format("(%f, %f)", _x, _y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	48 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	49
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	50 double x() const { return _x; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	51 double y() const { return _y; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	52
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	53 private {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	54 double _x, _y;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	55 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	56 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	57
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	58 Point minExtents(in Point a, in Point b) {
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	59 return Point(min(a.x, b.x), min(a.y, b.y));
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	60 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	61
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	62 Point maxExtents(in Point a, in Point b) {
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	63 return Point(max(a.x, b.x), max(a.y, b.y));
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	64 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	65
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	66 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	67 // The displacement between two locations in 2D space
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	68 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	69
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	70 struct Vector {
51 0eaf39fda206 First cut at palette class daveb parents: 48 diff changeset	71 static Vector DEFAULT = Vector(0.0, 0.0);
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	72
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	73 this(in double x, in double y) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	74 _x = x;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	75 _y = y;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	76 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	77
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	78 Vector opAdd(in Vector v) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	79 return Vector(_x + v._x, _y + v._y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	80 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	81
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	82 Vector opSub(in Vector v) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	83 return Vector(_x - v._x, _y - v._y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	84 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	85
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	86 Vector opNeg() const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	87 return Vector(-_x, -_y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	88 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	89
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	90 Vector opMul_r(in double d) const {
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	91 assert(!isnan(d));
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	92 return Vector(d * _x, d * _y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	93 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	94
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	95 Vector opDiv(in double d) const {
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	96 assert(!isnan(d));
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	97 return Vector(_x / d, _y / d);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	98 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	99
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	100 double length() const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	101 return sqrt(_x * _x + _y * _y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	102 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	103
62 6c3993f4c3eb Checkpoint daveb parents: 51 diff changeset	104 string toString() {
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	105 return std.string.format("[%f, %f]", _x, _y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	106 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	107
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	108 double x() const { return _x; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	109 double y() const { return _y; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	110
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	111 private {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	112 double _x, _y;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	113 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	114 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	115
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	116 /*
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	117 Vector normal(in Vector v) {
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	118 double l = v.length;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	119
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	120 if (l < 1e-9) { // TODO consolidate numerical stability constants
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	121 writefln("Warning: normalising tiny vector. Length: %f", l);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	122 return Vector(1.0, 0.0);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	123 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	124 else {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	125 return v / l;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	126 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	127 }
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	128 */
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	129
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	130 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	131 // A rectangle in 2D space.
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	132 // Internally represented by:
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	133 // a point defining the bottom left corner
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	134 // a vector defining the displacement to the upper right corner
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	135 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	136
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	137 struct Rectangle {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	138 static Rectangle DEFAULT;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	139
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	140 static this() {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	141 DEFAULT = Rectangle(Point.DEFAULT, Vector.DEFAULT);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	142 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	143
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	144 /*
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	145 static Rectangle from_arbitrary_corners(in Point corner1, in Point corner) {
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	146 }
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	147 */
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	148
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	149 this(in Point position, in Vector size) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	150 this(position.x, position.y, size.x, size.y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	151 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	152
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	153 this(in Point corner1, in Point corner) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	154 this(corner1.x, corner1.y, corner.x - corner1.x, corner.y - corner1.y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	155 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	156
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	157 double x0() const { return _position.x; }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	158 double y0() const { return _position.y; }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	159 double w() const { return _size.x; }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	160 double h() const { return _size.y; }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	161 double x1() const { return x0 + w; }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	162 double y1() const { return y0 + h; }
70 0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	163
0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	164 alias position corner0;
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	165 Point position() const { return _position; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	166
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	167 Vector size() const { return _size; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	168
70 0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	169 Point corner1() const { return _position + _size; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	170
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	171 bool valid() const { return _size.x > 0.0 && _size.y > 0.0; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	172
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	173 bool invalid() const { return !valid(); }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	174
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	175 double area() const { return _size.x * _size.y; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	176
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	177 // Intersection
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	178 Rectangle opAnd(in Rectangle r) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	179 if (invalid() \|\| r.invalid()) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	180 return DEFAULT;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	181 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	182 else {
70 0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	183 Point max = minExtents(corner1(), r.corner1());
0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	184 Point min = maxExtents(corner0(), r.corner0());
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	185
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	186 if (max.x < min.x \|\| max.y < min.y) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	187 return DEFAULT;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	188 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	189 else {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	190 return Rectangle(min, max);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	191 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	192 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	193 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	194
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	195 // Union
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	196 Rectangle opOr(in Rectangle r) const {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	197 if (invalid()) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	198 return r;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	199 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	200 else if (r.invalid()) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	201 return this;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	202 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	203 else {
70 0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	204 return Rectangle(minExtents(corner0(), r.corner0()),
0e61702c6ea6 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 67 diff changeset	205 maxExtents(corner1(), r.corner1()));
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	206 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	207 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	208
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	209 //
15 2f79aab4d385 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 13 diff changeset	210
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	211 // FIXME this method is all about pixels. Not sure it belongs in
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	212 // this file, let alone this class.
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	213 void getQuantised(out int x, out int y, out int w, out int h) const {
17 c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	214 x = cast(int)floor(_position.x);
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	215 y = cast(int)floor(_position.y);
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	216 w = cast(int)ceil(_position.x + _size.x) - x;
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	217 h = cast(int)ceil(_position.y + _size.y) - y;
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	218 }
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	219
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	220 //
c643c04e3f5e Checkpoint David Bryant <daveb@acres.com.au> parents: 16 diff changeset	221
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	222 Point centre() const { return _position + _size / 2.0; }
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	223
62 6c3993f4c3eb Checkpoint daveb parents: 51 diff changeset	224 string toString() {
4 58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	225 return std.string.format("{%s, %s}", _position, _size);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	226 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	227
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	228 private {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	229 this(double x, double y, double w, double h) {
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	230 if (w < 0.0) { x += w; w = -w; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	231 if (h < 0.0) { y += h; h = -h; }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	232 _position = Point(x, y);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	233 _size = Vector(w, h);
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	234 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	235
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	236 Point _position;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	237 Vector _size;
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	238 }
58a8ad20b228 Ooops, this got left out of previous commit David Bryant <daveb@acres.com.au> parents: diff changeset	239 }
25 8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	240
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	241 Rectangle growCentre(in Rectangle r, in Vector amount) {
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	242 return Rectangle(r.x0 - amount.x / 2, r.y0 - amount.y / 2, r.w + amount.x, r.h + amount.y);
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	243 }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	244
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	245 Rectangle growCentre(in Rectangle r, in double amount) {
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	246 return Rectangle(r.x0 - amount / 2, r.y0 - amount / 2, r.w + amount, r.h + amount);
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	247 }
535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	248
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	249 // TODO review these functions.
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	250 // Want a clear and simple set.
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	251
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	252 /+
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	253 Rectangle move(in Rectangle r, in Vector displacement) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	254 return Rectangle(r.position + displacement, r.size);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	255 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	256
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	257 Rectangle reposition(in Rectangle r, in Point newPosition) {
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	258 return Rectangle(newPosition, r.size);
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	259 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	260
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	261 Rectangle resize(in Rectangle r, in Vector new_size) {
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	262 return Rectangle(r.position, new_size);
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	263 }
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	264
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	265 // Operations about the bottom left corner
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	266
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	267 Rectangle expand(in Rectangle r, in Vector expand_amount) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	268 return Rectangle(r.position, r.size + expand_amount);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	269 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	270
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	271 Rectangle shrink(in Rectangle r, in Vector shrink_amount) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	272 return Rectangle(r.position, r.size - shrink_amount);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	273 }
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	274 +/
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	275
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	276 // Operations about the centre
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	277
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	278 /+
78 024a5608087f Meh "David Bryant <bagnose@gmail.com>" parents: 70 diff changeset	279 Rectangle feather(in Rectangle r, double amount) { // feather isn't the right name
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	280 assert(amount >= 0.0);
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	281 assert(!isnan(amount));
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	282 return Rectangle(Point(r.position.x - amount, r.position.y - amount),
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	283 Vector(r.size.x + 2.0 * amount, r.size.y + 2.0 * amount));
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	284 }
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	285 +/
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	286
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	287 private {
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	288 // This function computes the intersection of two lines.
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	289 // The lines are defined by a start point and an end point, however they
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	290 // notionally extend infinitely in each direction.
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	291 // The out parameters specify the fraction along the line-segment at which
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	292 // intersection occurred.
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	293 //
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	294 // This is a commmon building block for computing intersection between lines, segments,
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	295 // rectangles, etc.
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	296 //
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	297 // The function returns false if the lines are parallel or nearly so.
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	298 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	299 // Influenced by http://ozviz.wasp.uwa.edu.au/~pbourke/geometry/lineline2d/
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	300 bool computeIntersection(in Point pa1, in Point pa2, out double ua,
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	301 in Point pb1, in Point pb2, out double ub) {
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	302 double den = (pb2.y - pb1.y) * (pa2.x - pa1.x) - (pb2.x - pb1.x) * (pa2.y - pa1.y);
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	303
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	304 if (abs(den) < 1e-9) { // TODO consolidate constants used for numerical stability
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	305 // Lines are parallel or nearly so
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	306 warning("Warning, parallel lines!");
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	307 return false;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	308 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	309 else {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	310 // It will be safe to divide by den
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	311 double numA = (pb2.x - pb1.x) * (pa1.y - pb1.y) - (pb2.y - pb1.y) * (pa1.x - pb1.x);
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	312 double numB = (pa2.x - pa1.x) * (pa1.y - pb1.y) - (pa2.y - pa1.y) * (pa1.x - pb1.x);
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	313
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	314 ua = numA / den;
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	315 ub = numB / den;
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	316
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	317 return true;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	318 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	319 }
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	320
f2e4e1d29b98 Bah daveb parents: 40 diff changeset	321 /+
66 43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	322 double compute_angle(in Point p1, in Point p2) {
43cc2135ced0 Some code cleanups "David Bryant <bagnose@gmail.com>" parents: 62 diff changeset	323 }
41 f2e4e1d29b98 Bah daveb parents: 40 diff changeset	324 +/
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	325 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	326
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	327 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	328 // A line (notionally infinitely extending in both directions) in 2D space.
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	329 // Internally represented by:
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	330 // a point at an arbitrary location along the line
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	331 // a vector defining the gradient of the line
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	332 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	333
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	334 struct Line {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	335 this(in Point p, in Vector g) {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	336 _point = p;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	337 _gradient = g;
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	338 // FIXME should we normalise (make unit length) the gradient?
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	339 assert(_gradient.length > 1e-6); // FIXME how to best deal with this
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	340 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	341
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	342 this(in Point a, in Point b) {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	343 _point = a;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	344 _gradient = b - a;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	345 assert(_gradient.length > 1e-6); // FIXME as above
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	346 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	347
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	348 Point point() const { return _point; }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	349 Vector gradient() const { return _gradient; }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	350
62 6c3993f4c3eb Checkpoint daveb parents: 51 diff changeset	351 string toString() {
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	352 return std.string.format("{%s %s}", _point, _gradient);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	353 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	354
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	355 private {
51 0eaf39fda206 First cut at palette class daveb parents: 48 diff changeset	356 Point _point; // Arbitrary point along line
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	357 Vector _gradient;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	358 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	359 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	360
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	361 // Calculate the point "p" where lines "a" and "b" intersect.
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	362 // Returns false if lines are parallel or too close for numerical stability
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	363 bool intersection(in Line a, in Line b, out Point p) {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	364 Point pa = a.point;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	365 Vector va = a.gradient;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	366 double ua;
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	367
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	368 Point pb = b.point;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	369 Vector vb = b.gradient;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	370 double ub;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	371
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	372 if (computeIntersection(pa, pa + va, ua, pb, pb + vb, ub)) {
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	373 // We could just have easily evaluated for line b...
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	374 p = pa + ua * va;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	375 // p = pb + ub * vb;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	376 return true;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	377 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	378 else {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	379 return false;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	380 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	381 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	382
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	383 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	384 // A line segment (has a beginning and an end) in 2D space.
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	385 //
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	386
25 8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	387 struct Segment {
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	388 this(in Point a, in Point b) {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	389 _begin = a;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	390 _end = b;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	391 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	392
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	393 Point begin() const { return _begin; }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	394 Point end() const { return _end; }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	395
62 6c3993f4c3eb Checkpoint daveb parents: 51 diff changeset	396 string toString() {
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	397 return std.string.format("{%s %s}", _begin, _end);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	398 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	399
25 8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	400 private {
8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	401 Point _begin, _end;
8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	402 }
8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	403 }
8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	404
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	405 /*
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	406 Segment reverse(in Segment s) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	407 return Segment(s.end, s.begin);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	408 }
79 535bae7a7305 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 78 diff changeset	409 */
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	410
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	411 bool intersection(in Segment a, in Segment b, out Point p) {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	412 Point pa1 = a.begin;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	413 Point pa2 = a.end;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	414 double ua;
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	415
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	416 Point pb1 = b.begin;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	417 Point pb2 = b.end;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	418 double ub;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	419
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	420 if (computeIntersection(pa1, pa2, ua, pb1, pb2, ub)) {
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	421 if (ua >= 0.0 && ua <= 1.0 && // inside of segment a
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	422 ub >= 0.0 && ub <= 1.0) { // inside of segment b
67 31d10176415d Checkpoint daveb parents: 66 diff changeset	423 // We could just as easily evaluated for line b...
33 157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	424 p = pa1 + ua * (pa2 - pa1);
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	425 // p = pa2 + ub * (pb2 - pb1);
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	426 return true;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	427 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	428 else {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	429 return false;
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	430 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	431 }
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	432 else {
157b4ad5615d Added intersection code for lines and segments. David Bryant <bagnose@gmail.com> parents: 28 diff changeset	433 return false;
25 8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	434 }
8f58a8f88735 Checkpoint "David Bryant <bagnose@gmail.com>" parents: 19 diff changeset	435 }
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	436
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	437 bool intersection(in Segment a, in Line b, out Point p) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	438 Point pa1 = a.begin;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	439 Point pa2 = a.end;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	440 double ua;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	441
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	442 Point pb = b.point;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	443 Vector vb = b.gradient;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	444 double ub;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	445
48 1b4c9ba58673 Stylistic overhaul. daveb parents: 41 diff changeset	446 if (computeIntersection(pa1, pa2, ua, pb, pb + vb, ub)) {
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	447 if (ua >= 0.0 && ua <= 1.0) { // inside of segment
67 31d10176415d Checkpoint daveb parents: 66 diff changeset	448 // We could just as easily evaluated for line b...
34 c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	449 p = pa1 + ua * (pa2 - pa1);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	450 // p = pb + ub * vb;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	451 return true;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	452 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	453 else {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	454 return false;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	455 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	456 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	457 else {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	458 return false;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	459 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	460 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	461
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	462 bool intersection(in Line a, in Segment b, out Point p) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	463 // Equivalent to intersection of segment and line. Just reverse the args.
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	464 return intersection(b, a, p);
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	465 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	466
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	467 /+
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	468 bool intersection(in Line l, in Rectangle r, out Point p1, out Point p2) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	469 // TODO
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	470 return false;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	471 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	472
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	473 bool intersection(in Segment s, in Rectangle r, out Point p1, out Point p2) {
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	474 // TODO
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	475 return false;
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	476 }
c2f11e1d7470 Geometry cleanup and checkpoint. David Bryant <bagnose@gmail.com> parents: 33 diff changeset	477 +/

Mercurial > projects > doodle

annotate doodle/tk/geometry.d @ 79:535bae7a7305