Why do problems have this method anyway? It was obviously never used...

When type traits for the "double-float" or "double-double" types are
needed, they should be added as specializations of the
std::numeric_limits template in some of the TNL header files.

Implementation of TimeIndependentPDESolver.