inline static Real dummy_precision() { Real r; r.fill(NumTraits<T>::dummy_precision()); return r; }
inline static Real epsilon() { Real r; r.fill(NumTraits<T>::epsilon()); return r; }