## The CompGeoJS Project

I’ve received a few emails from readers of this (now closed) blog regarding my current development direction. I’m actually taking some time off after working on a very long Flex gig to get involved in some new stuff. As I look back over the span of time from 2004 (when I started this blog) to nearly 2014, my two favorite open-source projects were Degrafa and Singularity. The Singularity project provided an Actionscript library for commonly used tools from the field of computational geometry. I wrote a series of free TechNotes to explain the mathematics behind the techniques. Unfortunately, the TechNotes did not follow any reasonable pedagogical order, and I kind of regret the approach.

I’m very interested in reviving this development for the new wave of interactive developers who are concentrating their efforts in HTML 5 and JS. The goal of the CompGeoJS project is to provide interactive developers with a comprehensive, open-source library for addressing problems arising from the fields of analytic and computational geometry. Unlike Singularity, the library will be developed first, along with a robust collection of demos and online documentation. The process will be much more organized this time around and will eventually enable developers to tackle a wide variety of applications such as architecture, floor planning, advanced programmatic animation, and gaming.

Once the library reaches a solid ‘version one’ state, I will write an eBook to accompany the online documentation that discusses the mathematics behind all the tools. I hope this will be of use to developers who wish to modify the library for very customized uses as well as instructors of numerical analysis and computational geometry.

It’s at a VERY early stage of development, but you can check out what is available right now and where the library is going at the CompGeoJS project home page.

Thanks and enjoy!

## Time To Say Goodbye

When I started this blog, I wasn’t sure how much time I would devote to blogging, so I took the free WordPress route. I’ve been incredibly busy over the last year and a half and the lack of posts reflects that situation. I’m now in a position to devote a small amount of time to blogging and I decided to take the leap into the 21st century and create a WordPress site at algorithmist.net (my business site).

Since I now have a blog at my site, I will no longer maintain this one. I will leave all the posts online since many of them have code samples from which readers may realize some benefit.

I will continue blogging about math and programming, just in a different place. I hope you enjoyed the content of this blog and will keep up with the new one.

thanks!

– jim armstrong

## FlashBuilder 4.6 Hangs on Startup

Sorry about the extreme lack of posts. I hope this tip makes up for it :) So, this morning I load FB 4.6 and it hangs. I tried all the usual suspects and still no luck. I searched online and read the amusing story of Adobe wanting to charge people $250 for a tech support call to fix what appears to be a bug. Clearly a one-trick pony company with Photoshop – they sure can’t make anything else work.

Don’t waste your time or money with DumbDobe. A bit more searching yielded this post and it worked for me without reinstall or any other hassle. I wanted to pass along some props to the original poster and give the fix some attention in case this happens to you.

And Adobe is leading the charge for HTML 5? Yikes! I knew I should have been a lawyer.

## Rotate A Point Collection About Another Point

Well, this will probably be the last, or nearly last post of a busy year that involved a lot of work on gigs and little focus on this blog. Sorry, maybe next year will be better :)

In the previous example, I showed how to rotate a box (rectangle) around an arbitrary point, but the algorithm and code never presumed anything about the geometric nature of the object. It’s possible to extend the exact same algorithm to a collection of points and that is the subject of the current post.

Instead of maintaining variables for the four vertices of a rectangle, the code was modified to work with an arbitrary point collection. A RotatablePoint class was used to hold data points and offload some of the computations. This greatly simplifies the actual demo and provides you with ample means for experimentation.

The online example starts with a small collection of points as shown below.

These points are rotated about the fixed point, drawn as a red dot. The drawing may be cleared after each rotation increment or continually updated from the prior drawing as shown below.

Have a Merry Christmas and I’ll try to post more next year!

## Trig Plus Algebra Plus Geometry Plus Vectors Equals Fun

Math teachers hear the most often-uttered student phrase in history a lot more than I do, but my ears have twitched to the infamous, “I’ll never use that” or its many variants over the years. Funny how people know right now what will and will not be needed or used for the rest of their lives :)

Now, taken in strict isolation, a single trig or algebra formula might appear to have limited or no use in the myriad of life situations in someone’s future. It is, however, quite interesting how often problems can be solved by strategic use of a single formula from multiple areas of study. The problem discussed in this post is one such example.

In the prior post, I illustrated how some basic trig concepts could be used to solve a rotation and bounding-box problem without any presumptions or special considerations of the programming environment. This example is similar in that it illustrates how to rotate a box around an arbitrary point in its interior. The box is represented by a sequence of four coordinate pairs. It is not a Flash symbol or anything other than a sequence of coordinates. The drawing environment is simple; it can move the pen, draw lines, and draw filled circles. The programming environment contains the typical set of math functions.

Our math background includes a semester of trig, analytic geometry, and algebra. We know the very basics of vectors, but have not been introduced to matrices or matrix/coordinate transforms. As it happens, all we need to know to solve the stated problem is

1 – Polar coordinates, i.e. x = r cos(a) , y = r sin(a)

2 – Parametric equation of a straight line, i.e. P = (1-t)P0 + tP1

3 – How to add and subtract vectors

4 – How to compute the distance between two points

Refer to the diagram below.

The four points of the box are A (x1,y1), B (x2,y2), C (x3,y3), and D (x4,y4). The rotation point is R (rx,ry). The example code was written in Flex and the coordinate system in the Flash player is y-down. We are not given the coordinates of R. In fact, the box may be at an angle to the horizontal as shown in the lower, right part of the diagram. The only thing we know about R is that the component parallel to each side of the box is a fraction of that side’s length. This is where vector math can be useful.

The coordinates of R can be computed by adding the vectors **V1** and **V2**. The percentage inputs can be converted into numbers t1 and t2 in [0,1]. The exact coordinates of the terminal points of these vectors are obtained from the parametric equation of a line. Let the non-bolded notation V1 and V2 refer to vectors from the origin representing these terminal points. Then, V1 = (1-t1)A + t1B and V2 = (1-t2)A + t2D. Since the initial point of both **V1** and **V2** is A, the vector constituents (Δx and Δy) are immediately available. Add the vectors and add the result to A to obtain the coordinates for R. So, we can adjust sliders that change the percentage along each side and quickly compute R regardless of any prior rotation since we should always recompute the correct coordinates for A, B, C, and D.

When R is set, imagine lines drawn from R to each of the vectors A, B, C, and D. Let the distances be r1, r2, r3, and r4, respectively. Also compute the angle each line segment makes with the horizontal using the atan2() function. Since this returns a value in [-Π, Π], a small adjustment is made to always record the angles in [0, 2Π]. This makes studying the angle computations a bit easier if you want to trace out the results as part of deconstructing the code.

Also record the initial rotation angle when R is assigned. Subsequent changes to this value rotate the box around R by some delta. If the initial angle of the box is θ, and the delta is δ, then the four points are rotated around R by an angle θ+δ. Compute the coordinates of each new vector, A, B, C, and D by using the formula for polar coordinates and add the center, R.

These concepts are illustrated in an online demo.

The percentage along each side used to compute R is adjusted by sliders as well as the rotation angle. The red dot visually identifies the rotation point, R. You may optionally display the circle traced by each of the four corner points as the box is rotated. There is a *_clear* (Boolean) variable in the code that you may adjust to either clear the drawing each rotation update or choose to draw over each update of the box. If the value is false, you can produce some spirograph-style drawings.

The UI was created in Flex, but the actual code is very straightforward and does not rely on anything unique to Actionscript. It could be ported to many other environment quite easily. There is also some internal documentation on how to use trig identities to reduce the number of trig computations at each rotation update.

Even more importantly, there is nothing in the code or the algorithm that relies on a box being rotated other than computing R. The actual rotation is applied to a sequence of points using the simplest of concepts from a few math disciplines. I hope this provides some sense of appreciation of problems that can be solved using only the most fundamental techniques from a few branches of applied math.

## Trig 101

In the previous post, I promised to discuss the math behind rotating a rectangular sprite about its upper, left-hand registration point; that is, the origin of the sprite in its local coordinate space is the upper, left-hand corner. This is actually a segue into a more general discussion of rotation about an arbitrary point.

While some programming environments offer more direct capability than others, low-level math serves as a common denominator between them. If you understand the math, then it’s easy and efficient to port capability from one environment to another. There are also times where inlining direct computations offers opportunities for performance optimization in mobile apps or games. So, let’s start with that age-old question of what will I ever use trig for?

The programming example in the prior post dealt with computing the axis-aligned bounding box of a rectangular sprite after rotation about its origin. Refer to the diagram below.

The original sprite in light blue is rotated through an angle, *a*, to the orientation shown in light red. The origin (upper, left-hand corner) remains the same while the points A, B, and C are rotated to new positions A’, B’, and C’. Once the coordinates of these new points are known, it is trivial to compute the AABB (axis-aligned bounding box) shown in red.

From the outset drawing to the lower, right, let the width and height of the sprite be represented by *w* and *h*, respectively. The distance, r, is simply sqrt(w^{2} + h^{2}) and the angle, *b*, may be precomputed using the atan2 function available in most any programming language. Using the base definitions of sine and cosine, the coordinates of point A’ are ( *w*cos(a), w*sin(a)* ). The coordinates of B’ are ( *r*cos(a+b), r*sin(a+b)* ). If a+b = c, then the coordinates of C’ are computed as ( *h*cos(c+Π/2), h*sin(c+Π/2)* ) .

There is no need for an additional trig computation for C’ since we can use the well-known identity, *sin(x+y) = sin(x)*cos(y) + sin(y)*cos(x)*. cos(x+y) = *cos(x)*cos(y) – sin(x)*sin(y)*. *cos(Π/2) = 1 *and *sin(Π/2) = 0. *Plug these into the identities and you should see that we can re-use the computation of *cos(c)* and *sin(c)*. Refer to the code (the *rotation()* mutator function) in the prior post to check your computations. It should be easy to follow whether Actionscript, Javascript, C++, or Obj C is your language of choice.

Once we have the coordinates, the AABB is trivially computed and the computations can be easily ported to any programming environment or further optimized inside the application. All with some very basic trig and your friendly, neighborhood pythagorean theorem :) The operations could be optimized beyond what is illustrated in the code and that is left as an exercise.

Notice that we computed the rotated points in the sprite’s local coordinate space. The AABB is often computed in parent space, so it’s necessary to translate to the origin of that space. This is illustrated in the *boundingBox()* accessor function.

There is, of course, nothing special about rotation about the upper, left-hand corner. The next example will cover rotation about the centroid of the sprite, which is a special case of rotation about an arbitrary point in the sprite’s coordinate space. No matter what application you are working on or language you are working with, understanding what is happening ‘under the hood’ can go a long way. I hope you found something useful in this post.

## Dragging Bounds After Rotation

Ah, finally! Some code to give away! I was recently asked to develop a prototype for a friend. I don’t mind helping someone out, but my general policy is that if I do something free for one person, then I do it for everyone. All such developments are performed with the understanding that I reserve the right to post results to my blog.

This particular prototype dealt with simple dragging inside bounds with the twist that the rectangular item (most often an image) could be rotated about its upper, left-hand registration point yet still dragged only within the prescribed bounds. This was accomplished with a DraggableItem class that extends MovieClip. The visual symbol is associated with an instance of this class. Rectangular bounds are assigned and the DraggableItem class dispatches events on start drag, end drag, and collision with a boundary.

Bounds are in parent space and the DraggableItem rotation method recomputes the axis-aligned bounding box after rotation. The demo draws the AABB during rotation.

The computations are a simple application of trig and the entire demo could be ported to a number of environments. I may convert this to html/js and deconstruct the math in a subsequent series of posts. In the mean time, realize that like most demos, this one was created hastily and only lightly tested. Documentation is sparse, so you are on your own for the deconstruction. There may be an issue or two that requires polishing (for example, you can rotate the visual symbol beyond bounds) and I invite people to make this demo better. That’s what open source is all about :)