Added more methods to compute the SPH equations of motion.  Renamed
several variables. Added C wrappers for some thrust functions.