Morley: Bugfix that fixes the Bug!

dune/functions/functionspacebases/morleybasis.hh

@@ -12,6 +12,9 @@
 #include <dune/istl/bcrsmatrix.hh>
 #include <dune/localfunctions/morley.hh>
 
+// TODO remvoe
+// #include <dune/istl/io.hh>
+
 namespace Dune
 {
   namespace Functions
@@ -97,15 +100,15 @@ namespace Dune
           {
             B_11[i] = -referenceTangents[i][1]
                         * (-globalTangents[i][1] * jacobianTransposed[0][0]
-                           + globalTangents[i][0] * jacobianTransposed[1][0])
+                           + globalTangents[i][0] * jacobianTransposed[0][1])
                     + referenceTangents[i][0]
-                          * (-globalTangents[i][1] * jacobianTransposed[0][1]
+                          * (-globalTangents[i][1] * jacobianTransposed[1][0]
                              + globalTangents[i][0] * jacobianTransposed[1][1]);
             B_12[i] = -referenceTangents[i][1]
                         * (globalTangents[i][0] * jacobianTransposed[0][0]
-                           + globalTangents[i][1] * jacobianTransposed[1][0])
+                           + globalTangents[i][1] * jacobianTransposed[0][1])
                     + referenceTangents[i][0]
-                          * (globalTangents[i][0] * jacobianTransposed[0][1]
+                          * (globalTangents[i][0] * jacobianTransposed[1][0]
                              + globalTangents[i][1] * jacobianTransposed[1][1]);
           }
 
@@ -117,11 +120,14 @@ namespace Dune
             for (std::size_t j = 0; j < 3; ++j)
               if (j != i)
               {
-                mat_[i][3 + j] = sign * B_12[i] * l_inv[i];
+                mat_[j][3 + i] = sign * B_12[i] * l_inv[i];
                 sign *= -1;
               }
             mat_[3 + i][3 + i] = B_11[i];
           }
+
+          // TODO remvoe
+          // printmatrix(std::cout, mat_, "morley transformation matrix", "");
         }
         BCRSMatrix<R> mat_;