You are right: ```fabs``` is limited to certain float or double. According to the [docu](http://en.cppreference.com/w/cpp/numeric/math/abs) ```abs``` is limited to integers. I don't know if there is a templated abs function. But in general I suspect the matrix entries are floating point values.