Reports:
Following the method set out in the interim report I built up the system in C++. I ran into some problems with the conjugate gradient solver but eventually was able to figure it out. After completing the mathematical sections of the program I set out to create a limited graphical representation using the GL_Points function in openGL and a pixel whose elements were based off of the voxel values. At this point I realized that there was some error in the math sections since for some reason the values of the velocities were becoming exponentially large and then turning into undefined values. Unfortunately I did not have enough time to determine the cause of this problem.
-Jason Poulos