Awesome

This code is translated from the code provided in Jos Stam's GDC2003 paper