I have a custom type. It's a small autodiff variable, templated on scalar type and size. Everything works great; however, things break down when I attempt to use the LU decomposition. The trouble is that the precision() and machine_epsilon() functions use template specializations to resolve the type, which for some reason, doesn't work.
Is it possible to use Scalar types that are themselves templated?