From dc72d1d551562d76d07176876478bc2d49697e95 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@ulg.ac.be>
Date: Thu, 4 Feb 2016 22:30:33 +0000
Subject: [PATCH] patch from Wendy

---
 Numeric/polynomialBasis.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Numeric/polynomialBasis.cpp b/Numeric/polynomialBasis.cpp
index 3c3596b186..2e1c563dc9 100644
--- a/Numeric/polynomialBasis.cpp
+++ b/Numeric/polynomialBasis.cpp
@@ -141,7 +141,7 @@ void polynomialBasis::f(const fullMatrix<double> &coord, fullMatrix<double> &sf)
   double p[1256];
   sf.resize (coord.size1(), coefficients.size1());
   for (int iPoint = 0; iPoint < coord.size1(); iPoint++) {
-    evaluateMonomials(coord(iPoint, 0), coord(iPoint, 1),
+    evaluateMonomials(coord(iPoint, 0), coord.size2() > 1 ? coord(iPoint, 1) : 0,
                       coord.size2() > 2 ? coord(iPoint, 2) : 0, p);
     for (int i = 0; i < coefficients.size1(); i++) {
       sf(iPoint,i) = 0.;
-- 
GitLab