From: Bill Hoffman <bill.hoffman@ki...>  20020308 19:18:09

The following always aborts with a vnl_matrix<std::complex<float>>: //: Return false if any element of (*this) is inf or nan template <class T> bool vnl_matrix<T>::is_finite() const { for (unsigned int i = 0; i < this>rows(); ++i) for (unsigned int j = 0; j < this>columns(); ++j) if (!vnl_math_isfinite((*this)(i,j))) return false; return true; } There is no overload for vnl_math_isfinite for complex types.... Bill 