Fixed a bug in atan2. It gave a division by zero error for calls like

[ginac.git] / ginac / inifcns_trans.cpp

diff --git a/ginac/inifcns_trans.cpp b/ginac/inifcns_trans.cpp
index 39931dcd738cf5a72cebd371468726c88fa1809f..3f161d218c068e2a4d016fce5d5cfcb2d32f789e 100644 (file)
--- a/ginac/inifcns_trans.cpp
+++ b/ginac/inifcns_trans.cpp
@@ -842,7 +842,7 @@ static ex atan2_eval(const ex & y, const ex & x)
                        return _ex1_2*Pi;
 
                // atan(y, 0), y real and negative -> -Pi/2
                        return _ex1_2*Pi;
 
                // atan(y, 0), y real and negative -> -Pi/2

-               if (y.info(info_flags::negative))
+               if (y.info(info_flags::real) && !y.is_zero())

                        return _ex_1_2*Pi;
        }
 
                        return _ex_1_2*Pi;
        }