Some people (including me) have been adding double versions of some of the math and geometry routines into Blender’s internal library. These are only used right now for intermediate calculations, for greater precision as you say. When results get stored back into the permanent geometry in Blender, it is still always as floats.
For some discussion on why going to all doubles would be problematic in terms of memory and performance, see Arithmetic types in Blender