Make matrix::solve() work with non-normal zeros.

[ginac.git] / check / time_gammaseries.cpp

diff --git a/check/time_gammaseries.cpp b/check/time_gammaseries.cpp
index db3d442a27c59e5cddfe6b475e403631109535b5..0c8c6575c0b0ba28b5c2e1a66cb9c36bca8d5d0b 100644 (file)
--- a/check/time_gammaseries.cpp
+++ b/check/time_gammaseries.cpp
@@ -3,7 +3,7 @@
  *  Some timings on series expansion of the Gamma function around a pole. */
 
 /*
- *  GiNaC Copyright (C) 1999-2016 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2018 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
@@ -55,15 +55,10 @@ unsigned time_gammaseries()
 
        cout << "timing Laurent series expansion of Gamma function" << flush;
 
-       vector<unsigned> sizes;
+       vector<unsigned> sizes = {20, 25, 30, 35};
        vector<double> times;
        timer omega;
 
-       sizes.push_back(20);
-       sizes.push_back(25);
-       sizes.push_back(30);
-       sizes.push_back(35);
-
        for (vector<unsigned>::iterator i=sizes.begin(); i!=sizes.end(); ++i) {
                omega.start();
                result += tgammaseries(*i);