Check for fixed bug in multiple polylogs.
authorJens Vollinga <vollinga@thep.physik.uni-mainz.de>
Wed, 13 Dec 2006 20:04:23 +0000 (20:04 +0000)
committerJens Vollinga <jensv@nikhef.nl>
Fri, 4 Apr 2008 12:51:35 +0000 (14:51 +0200)
check/exam_inifcns_nstdsums.cpp

index 247fa80d8ebb2edd12b6ac442c1047fe800f7237..d0f6869effc033164f24f9d7a40e1764c4f2fe3e 100644 (file)
@@ -290,6 +290,30 @@ static unsigned inifcns_test_LiG()
 }
 
 
+////////////////////////////////////////////////////////////////////////////////
+////////////////////////////////////////////////////////////////////////////////
+//  legacy exam - checking for historical bugs
+////////////////////////////////////////////////////////////////////////////////
+////////////////////////////////////////////////////////////////////////////////
+
+
+static unsigned inifcns_test_legacy()
+{
+       Digits = 17;
+       ex prec = 5 * pow(10, -(int)Digits);
+
+       unsigned result = 0;
+
+       ex r1 = zeta(lst(1,1,1,1,1,1),lst(-1,-1,-1,1,1,1));
+       if ((r1.evalf() - numeric("-0.0012588769028204890704")) > prec) {
+               clog << "zeta({1,1,1,1,1,1},{-1,-1,-1,1,1,1}) seems to be wrong." << endl;
+               result++;
+       }
+
+       return result;
+}
+
+
 unsigned exam_inifcns_nstdsums(void)
 {
        unsigned result = 0;
@@ -300,6 +324,7 @@ unsigned exam_inifcns_nstdsums(void)
        result += inifcns_test_S();
        result += inifcns_test_HLi();
        result += inifcns_test_LiG();
+       result += inifcns_test_legacy();
        
        return result;
 }