FIM method implemented for 2D GPU and FIM-FSM implemented for 2D CPU (parallel).

FIM method implemented. Neighbours are being found on CPU. 3D parallel method disabled because of Array changes.

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.

tnlDirectEikonalMethodsBase_impl.h: 2D method repaired, 3D method ready for debug
tnlDirectEikonalMethod2D_impl.h: Added save of aux-* as results from computation
tnlDirectEikonalMethod3D_impl.h: Implemented all sweaps and saving results from computation