Clarify types of test suite files.

[ginac.git] / check / exam_archive.cpp

diff --git a/check/exam_archive.cpp b/check/exam_archive.cpp
index 97962d7a637bc8cd34b4ab7381f1b16b07a8e536..0abc752183e67c512dac67b2e577b3e837a1fca1 100644 (file)
--- a/check/exam_archive.cpp
+++ b/check/exam_archive.cpp
@@ -3,7 +3,7 @@
  *  Here we test GiNaC's archiving system. */
 
 /*
- *  GiNaC Copyright (C) 1999-2008 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2020 Johannes Gutenberg University Mainz, Germany
  *
  *  This program is free software; you can redistribute it and/or modify
  *  it under the terms of the GNU General Public License as published by
@@ -20,18 +20,20 @@
  *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
  */
 
-#include <iostream>
-#include <fstream>
 #include "ginac.h"
-using namespace std;
 using namespace GiNaC;
 
+#include <fstream>
+#include <iostream>
+using namespace std;
+
+#include <cln/cln.h>
+
+
 unsigned exam_archive()
 {
        unsigned result = 0;
        
-       cout << "examining archiving system" << flush;
-
        symbol x("x"), y("y"), mu("mu"), dim("dim", "\\Delta");
        ex e, f;
 
@@ -40,9 +42,9 @@ unsigned exam_archive()
        e = -42 * x * pow(y, sin(y*Catalan)) * dirac_ONE()
            * epsilon_tensor(idx(fail(), 3), idx(0, 3), idx(y/2, 3))
          + lorentz_g(
-             varidx(lst(x, -11*y, acos(2*x).series(x==3-5*I, 3)) * color_ONE()
+             varidx(lst{x, -11*y, acos(2*x).series(x==3-5*I, 3)} * color_ONE()
                * metric_tensor(varidx(log(cos(128.0/(x*y))), 5), varidx(2, 5)), zeta(3)),
-             varidx(diag_matrix(lst(-1, Euler, atan(x/y==-15*I/17)))
+             varidx(diag_matrix({-1, Euler, atan(x/y==-15*I/17)})
                * delta_tensor(idx(x, 2), idx(wild(7), 3)), zeta(3), true),
              true
            )
@@ -63,9 +65,10 @@ unsigned exam_archive()
                std::ifstream fin("exam.gar", std::ios_base::binary);
                fin >> ar;
        }
-       f = ar.unarchive_ex(lst(x, y, mu, dim), "expr 1");
+       f = ar.unarchive_ex(lst{x, y, mu, dim}, "expr 1");
 
-       if (!f.is_equal(e)) {
+       ex difference = (f - e).expand();
+       if (!difference.is_zero()) {
                clog << "archiving/unarchiving " << e << endl
                     << "erroneously returned " << f << endl;
                ++result;
@@ -74,7 +77,49 @@ unsigned exam_archive()
        return result;
 }
 
+/** numeric::archive used to fail if the real part of a complex number
+ *  is a rational number and the imaginary part is a floating point one. */
+unsigned numeric_complex_bug()
+{
+       using namespace cln;
+       struct archive_unarchive_check
+       {
+               unsigned operator()(const cl_N& n) const
+               {
+                       ex e = numeric(n);
+                       archive ar;
+                       ar.archive_ex(e, "test");
+                       ex check = ar.unarchive_ex(lst{}, "test");
+                       if (!check.is_equal(e)) {
+                               clog << __FILE__ << ':' << __LINE__ << ": expected: " << e << ", got " << check << endl;
+                               return 1;
+                       }
+                       return 0;
+               }
+       } checker;
+       unsigned result = 0;
+       const cl_I one(1);
+       const cl_R three_fp = cl_float(3.0);
+       std::vector<cl_N> numbers = {
+               complex(one, one),
+               complex(one, three_fp),
+               complex(three_fp, one),
+               complex(three_fp, three_fp)
+       };
+       for (auto & n : numbers) {
+               result += checker(n);
+       }
+       return result;
+}
+
 int main(int argc, char** argv)
 {
-       return exam_archive();
+       unsigned result = 0;
+
+       cout << "examining archiving system" << flush;
+
+       result += exam_archive();  cout << '.' << flush;
+       result += numeric_complex_bug();  cout << '.' << flush;
+
+       return result;
 }