From 2bb98ef3ea9d6108b2b06475b3fd915d418aecad Mon Sep 17 00:00:00 2001
From: Ross Kynch <rosskynch@gmail.com>
Date: Wed, 7 Jun 2017 22:26:29 +0100
Subject: [PATCH] Added non-zero evaluation points to test

---
 .../polynomials_integrated_legendre_sz.cc     | 19 +++++++++++++-
 .../polynomials_integrated_legendre_sz.output | 25 +++++++++++++++++++
 2 files changed, 43 insertions(+), 1 deletion(-)

diff --git a/tests/base/polynomials_integrated_legendre_sz.cc b/tests/base/polynomials_integrated_legendre_sz.cc
index b027bfea72..61061da59f 100644
--- a/tests/base/polynomials_integrated_legendre_sz.cc
+++ b/tests/base/polynomials_integrated_legendre_sz.cc
@@ -48,13 +48,30 @@ void check_at_one (const std::vector<Polynomial<double> > &p)
 
 
 
+void check_at_half (const std::vector<Polynomial<double> > &p)
+{
+  // Ignore first two polynomials as the integrated Legendre polynomials are only defined
+  // for degree > 1, it is only added to maintain the recursive relation.
+
+  deallog << "Function value of polynomial at -0.5 | 0.5:" << std::endl;
+  for (unsigned int i=2; i<p.size(); ++i)
+    {
+      const double y = p[i].value(0.5);
+      const double z = p[i].value(-0.5);
+      deallog << y << " | " << z << std::endl;
+    }
+  deallog << std::endl;
+}
+
+
+
 void
 check_poly (const unsigned int n)
 {
   deallog << "Degree: " << n+1 << std::endl;
   std::vector<Polynomial<double> > p = IntegratedLegendreSZ::generate_complete_basis(n);
   check_at_one (p);
-  deallog << std::endl;
+  check_at_half (p);
 }
 
 
diff --git a/tests/base/polynomials_integrated_legendre_sz.output b/tests/base/polynomials_integrated_legendre_sz.output
index 8bc294c7bd..3b45c2bf59 100644
--- a/tests/base/polynomials_integrated_legendre_sz.output
+++ b/tests/base/polynomials_integrated_legendre_sz.output
@@ -1,4 +1,29 @@
 
 DEAL::Degree: 26
 DEAL::Function value of polynomial at right end point: ........................
+DEAL::Function value of polynomial at -0.5 | 0.5:
+DEAL::-0.375 | -0.375
+DEAL::-0.188 | 0.188
+DEAL::-0.0234 | -0.0234
+DEAL::0.0586 | -0.0586
+DEAL::0.0557 | 0.0557
+DEAL::0.0103 | -0.0103
+DEAL::-0.0265 | -0.0265
+DEAL::-0.0289 | 0.0289
+DEAL::-0.00603 | -0.00603
+DEAL::0.0158 | -0.0158
+DEAL::0.0183 | 0.0183
+DEAL::0.00408 | -0.00408
+DEAL::-0.0108 | -0.0108
+DEAL::-0.0130 | 0.0130
+DEAL::-0.00299 | -0.00299
+DEAL::0.00795 | -0.00795
+DEAL::0.00978 | 0.00978
+DEAL::0.00231 | -0.00231
+DEAL::-0.00617 | -0.00617
+DEAL::-0.00771 | 0.00771
+DEAL::-0.00186 | -0.00186
+DEAL::0.00497 | -0.00497
+DEAL::0.00628 | 0.00628
+DEAL::0.00153 | -0.00153
 DEAL::
-- 
2.39.5