matrix does not accept "-inf"

[Peter]

Is that a restriction created by yat or is it inherited from GSL?

If not inherited, why?

[Markus Ringnér]

You can test if your compiler supports infinity for doubles by checking the bool returned by:

std::numeric_limits<double>::has_infinity()

[Markus Ringnér]

Vhat is meant with does not accept -inf? If it is rad from a file to a matrix or something else?

If the problem is with reading -inf from a file, it could be a yat problem. If so could you check if one can read nan from a file?

[Peter]

Yes I'm using istream constructor reading from file.

The exception I guess is thrown from here

    if(utility::is_double(element)) {
 	          v.push_back(atof(element.c_str()));
 	        }
 	        else if (!element.size() || utility::is_nan(element)) {
 	          v.push_back(std::numeric_limits<double>::quiet_NaN());
 	        }
 	        else {
 	          std::stringstream ss("Warning: '");
 	          ss << element << "' is not accepted as a matrix element.";
 	          throw IO_error(ss.str());

[Peter]

Yes, I can read a nan from file, which is pretty obvious from the code above, but I tested anyway.

Should we modify is_double or shoule we add checks against -inf and inf?

[Markus Ringnér]

Maybe nan, inf and -inf should all be fixed in is_double ? After all they are in the (IEEE whatever) double specification. This means yat would be limited to systems with numeric_limits<double>::has_infinity and then (-)inf can be generated with (-)numeric_limits<double>::infinity(). This requirement is probably minor since yat already requires systems where doubles have quiet_nan. Would changing is_double influence any other yat functionality?

Perhaps we should start using std::isnan and std::isinf. They are not strictly speaking standard C++ as they were added to C99 (after the latest C++ standard C++98), but they are in the C++0X draft.

[Peter]

I removed is_double and re-compiled. Everything compiles except matrix and vector.

We have to document this change in the NEWS file.

[Peter]

I have started to use std::isnan in a couple of places (mainly assertions).

For this issue it does help though because is_doubel is a function checking whether a string corresponds to a double.

[Peter]

fixed in [917] and [918]