Whenever there is a NaN value inside an array, exp() will return an Inf instead of NaN.

Example bellow:


  #include <iostream>
  #include <limits>
  #include <Eigen/Dense>

  int main(){
    double nan = std::numeric_limits<double>::quiet_NaN();

    Eigen::MatrixXd m1(3,3);
    m1 << 10, 1.5, 3.2,
          nan, 10, 4.1,
          0.1, 4.1, 10;
         
    Eigen::MatrixXd m2(3,3);
    Eigen::ArrayXXd a1(3,3);

    //Has bug
    m2 = m1.array().exp().matrix();
    std::cout << m2 << std::endl << std::endl;
    
    //Has bug
    a1 = m1.array().exp();
    std::cout << a1 << std::endl << std::endl;

    //Output:
    //22026.5 4.48169 24.5325
    //inf 22026.5 60.3403
    //1.10517 60.3403 22026.5
    //
    //22026.5 4.48169 24.5325
    //inf 22026.5 60.3403
    //1.10517 60.3403 22026.5

    



    //Now, some other calls, just for comparison:
    a1 = m1.array();
    std::cout << a1 << std::endl << std::endl;  
      
    for (int i; i < m1.size(); i++)
      std::cout << std::exp(m1(i)) << std::endl;
    
    std::cout << std::endl;

    m2 = m1.array().log().matrix();
    std::cout << m2 << std::endl << std::endl;
    
    a1 = m1.array().log();
    std::cout << a1 << std::endl << std::endl;
  }