- As advertised: we are calling the Gamma function tgamma() now!

[ginac.git] / ginac / normal.cpp

diff --git a/ginac/normal.cpp b/ginac/normal.cpp
index fffdc925b09fe5d4e83040f153ed97d258a12c26..a8f64be82a988a330c6cb70997b60b86ae858af2 100644 (file)
--- a/ginac/normal.cpp
+++ b/ginac/normal.cpp
@@ -1025,7 +1025,7 @@ static ex heur_gcd(const ex &a, const ex &b, ex *ca, ex *cb, sym_desc_vec::const
 {
 //clog << "heur_gcd(" << a << "," << b << ")\n";
 
-       // Trivial cases
+       // GCD of two numeric values -> CLN
     if (is_ex_exactly_of_type(a, numeric) && is_ex_exactly_of_type(b, numeric)) {
         numeric g = gcd(ex_to_numeric(a), ex_to_numeric(b));
         numeric rg;
@@ -1057,9 +1057,9 @@ static ex heur_gcd(const ex &a, const ex &b, ex *ca, ex *cb, sym_desc_vec::const
         xi = mp * _num2() + _num2();
 
     // 6 tries maximum
-    for (int t=0; t<6; t++) {
-        if (xi.int_length() * maxdeg > 100000) {
-//clog << "giving up heur_gcd, xi.int_length = " << xi.int_length() << ", maxdeg = " << maxdeg << endl;
+    for (int t=0; t<6; t++) {  // MAGIC
+        if (xi.int_length() * maxdeg > 100000) {  // MAGIC
+// clog << "giving up heur_gcd, xi.int_length = " << xi.int_length() << ", maxdeg = " << maxdeg << endl;
             throw gcdheu_failed();
                }