[Blitz-support] isinf?

[Michael W. <th...@gm...>]

Hi all

I started on using Blitz++ and now am left a bit wondering as to why 
there are implementations of finite() and blitz_isnan(), but not of 
isinf(). Of course, one could construct it from the others, possibly like

template<class T, int rank>
blitz::Array<bool,rank>
isinf( const blitz::Array<T,rank>& a ) {
   blitz::Array<bool,rank> result(a.shape());
   result = !(finite(a) || blitz_isnan(a));
   return result;
}

but i think it would be much better to define it like the others using 
the implementation provided by the system library. So, is there some 
specific reason as to why it isn't there, or has it just been forgotten? 
I couldn't find any reference in the list archives.



Regards

Michael