patch from Wendy

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.;