> I was thinking of basing other such operations if needed on:
> bool numbers::is_finite (const std::complex< long double > x)
Upon second thought, I think this is also not a bad idea, and might make
the code better readable. If you like this better, feel free to add
namespace numbers {
template <typename T>
bool is_zero (const T &t) { return t == T(); }
}
or its reverse, is_nonzero, or both, and use them wherever useful.
Best
Wolfgang
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
