diff --git a/dune/gfe/orthogonalmatrix.hh b/dune/gfe/orthogonalmatrix.hh
index 62c294805a088c27c7deb0c3a45b59f5c8676750..b565c4e09604139f47562407b38f7e997d6f05cb 100644
--- a/dune/gfe/orthogonalmatrix.hh
+++ b/dune/gfe/orthogonalmatrix.hh
@@ -45,7 +45,7 @@ public:
         
         for (int i=0; i<N; i++)
             for (int j=0; j<N; j++) {
-                IdMinusXXT = (i==j);
+                IdMinusXXT[i][j] = (i==j);
                 for (int k=0; k<N; k++)
                     IdMinusXXT[i][j] -= X[i][k] * X[j][k];
             }