Subdivision Surfaces with Loop's algorithm

My first assignment on Geometric Modeling was about subdivision surfaces. I implemented Loop's algorithm, and work on VC7 using CGAL library. Most of the time on this assignment i was dealing with how to work with the library, since i found it very confusing = =" (not the library's fault though). I also use some of the code from PBRT for counting the vertices' weights.
Here's the result of my program:

original shape



















after 3 levels of subdivision