New paper out!

Haven’t updated this blog in a long time, but I figure that since it is summer now maybe I will finally be able to keep this thing regularly updated (see how long that lasts…)  Anyway, some good news is that my latest paper has been accepted to the SIAM/SIGGRAPH conference on solid and physical modeling.  Here is a link to the download:

There’s a bunch of neat ideas in here, but the big idea here is the Minkowski product.  The motivation for this comes from the basic Minkowski sum.  If we recall, for two sets A, B \subseteq \mathbb R^n, their Minkowski sum is defined as:

A \oplus B = \{ a + b | a \in A, b \in B \}

Or alternatively,

A \oplus B = \bigcup \limits_{a \in A} a + B

The idea behind the Minkowski product is to replace \mathbb R^n with an arbitrary group, G.  If we do this, then we can define the Minkowski product over G in the following way:

A \stackrel{G}{\otimes} B = \bigcup \limits_{a \in A} a B

Much like how the Minkowski sum is useful for collision detection of translating objects, the generalized Minkowski product can be used for collision detection between translating and rotating bodies.  Similarly, we can also define a Minkowski quotient which is analogous to the Minkowski difference:

A \stackrel{G}{\oslash} B = \bigcap \limits_{a \in A} a B

And even better yet, we show how to compute these quantities using convolution algebras!  These operators turn out to be very useful in understanding things like partial symmetries of solids, mechanism workspaces and robotics.  All of the gory details and more are in the paper!

BSP Trees and Collision Detection

Recently I published a paper in CAD on BSP tree merging via linear programming.  At a very high level, the idea behind this paper was to replace the tree partitioning step with a linear programming feasibility test.  Doing this leads to a bunch of nice things, such as the removal of the partition function and all 9 of its special cases.  It also has the benefit that that for some suitable linear programming algorithms, the cost of evaluating feasibility can be amortized over the tree traversal giving an optimal linear-time output sensitive algorithm.

Of course as it always turns out with these things, I recently discovered that I was not the first to think of using linear programming for solving this problem.  Just last week Prof. Alberto Paoluzzi visited our lab, and I was lucky enough to get a chance to talk to him about BSP trees.  As it just so happens, he proposed using linear programming for tree partitioning in his book Geometric Programming for Computer-Aided Design (which is actually a very good read.)  Of course, this shouldn’t diminish the contribution of our paper, as we do present a much simpler merge algorithm along with an analysis of both time complexity and robustness.  In any future works, I will definitely add a reference to his book.

Anyway, I wanted to talk about one idea which didn’t make it into that paper, but I still happen to think is rather interesting.  As an accidental side effect of our merging algorithm, we discovered a straightforward method for solving exact continuous collision detection between arbitrary moving non-convex objects.  The key idea to this process was initially suggested to me by my colleague Nicholas Smolinske (who sadly is no longer in graduate school.)  Here’s how it goes:

Suppose you have two convex shapes moving at a constant velocity.  For each shape, the object Minkowski summed with its trajectory traces out a 4D convex set space and time.  If you take these two sets for both shapes and intersect them, then the result is also a convex shape in 4D.  Moreover, the lowest point in time within the overlap region must be the point at which the two objects initially made contact.  This idea is illustrated in the following diagram:


Because all of the sets we are dealing with are convex, it is possible to exactly solve for the point of intersection using linear programming with an objective function c(x,t) = -t.  Now, the next question is how to do this with BSP trees?  Well, it turns out that one side effect from the algorithm proposed in our paper is that it ends up evaluating linear programming on each convex cell contained in the final BSP tree.  We can make use of this by simply picking the same objective function as above, and storing the minimal point over all convex cells as our result.  Alternatively, if we only want to get the point of intersection, then we don’t even need to compute the result of the merge so we may instead use it as a kind of driver for solving linear programming over the cell-complex hierarchically.

The idea of BSP trees representing shapes in space/time is not new, and to the best of my knowledge was first explored by Agarwal et al. (they give some novel tree construction strategies, but don’t really deal with merging or related issues).  Another interesting way to look at this idea is that one could understand the BSP tree as generalizing some of the ideas of polytope hierarchies (such as Dobkin-Kirkpatrick) to non-convex sets.  Of course this still doesn’t answer the big question which is how to go about building decent BSP trees (or hierarchies for that matter) and also how to pick which order to merge them.  Someday I might go back and try refining these ideas to the point where they could make for a nice short conference paper, but for now it is just something fun to think about.

Hello World

Hi.  My name is Mikola Lysenko and this is my first post on 0FPS. I am currently a first year graduate student at the University of Wisconsin-Madison, working in the spatial automation lab.  I used to run the website, but I ended up taking it down to save money on bandwidth expenses.  Anyway, I figure that lately I’ve been too lazy in my writing process and that a blog is a pretty good way to bring discipline back into my life.  I intend to focus on research topics in solid modeling and computer graphics, as well as other ideas I find interesting.