Modified the library to use doubles instead of single precision floats
Changed initial() to take const double * from double *

Other trivial changes...