Changes sign of derivatives to match function changes.

diff --git a/source/base/polynomials_bdm.cc b/source/base/polynomials_bdm.cc
index 1700b16dcb2cc13742a594887ba683fe38d1b1d4..e4deadea47d77e73d81aa4c64b44b5a83cd1d12e 100644 (file)
--- a/source/base/polynomials_bdm.cc
+++ b/source/base/polynomials_bdm.cc
@@ -126,10 +126,10 @@ PolynomialsBDM<dim>::compute (const Point<dim>            &unit_point,
           grads[start][0][1] = 0.;
           grads[start][1][0] = -unit_point(1) * monovali[0][2];
           grads[start][1][1] = -monovali[0][1];
-          grads[start+1][0][0] = -monovali[1][1];
-          grads[start+1][0][1] = -unit_point(0) * monovali[1][2];
+          grads[start+1][0][0] = monovali[1][1];
+          grads[start+1][0][1] = unit_point(0) * monovali[1][2];
           grads[start+1][1][0] = 0.;
-          grads[start+1][1][1] = monovali[1][1];
+          grads[start+1][1][1] = -monovali[1][1];
         }
       if (grad_grads.size() != 0)
         {
@@ -142,13 +142,13 @@ PolynomialsBDM<dim>::compute (const Point<dim>            &unit_point,
           grad_grads[start][1][1][0] = -monovali[0][2];
           grad_grads[start][1][1][1] = 0.;
           grad_grads[start+1][0][0][0] = 0;
-          grad_grads[start+1][0][0][1] = -monovali[1][2];
-          grad_grads[start+1][0][1][0] = -monovali[1][2];
-          grad_grads[start+1][0][1][1] = -unit_point(0) * monovali[1][3];
+          grad_grads[start+1][0][0][1] = monovali[1][2];
+          grad_grads[start+1][0][1][0] = monovali[1][2];
+          grad_grads[start+1][0][1][1] = unit_point(0) * monovali[1][3];
           grad_grads[start+1][1][0][0] = 0.;
           grad_grads[start+1][1][0][1] = 0.;
           grad_grads[start+1][1][1][0] = 0.;
-          grad_grads[start+1][1][1][1] = monovali[1][2];
+          grad_grads[start+1][1][1][1] = -monovali[1][2];
         }
     }
   else // dim == 3