Add as a template parameter & default to double precision. Might be faster/less memory intensive/more conducive to SIMD, depending on the situation?