fix minus in diffophesseboundary

NGSolve · Dec 18, 2017 · 01da27c · 01da27c
1 parent 3cc179e
commit 01da27c
Showing 1 changed file with 1 addition and 36 deletions.
diff --git a/fem/bdbequations.hpp b/fem/bdbequations.hpp
@@ -614,7 +614,7 @@ namespace ngfem
 		    {
 		      for( int l = 0; l < D-1; l++)
 			{
-			  mat(i*D+j,n) += jinv(k,i)*tmpmat(n,k*(D-1)+l)*jinv(l,j);
+			  mat(i*D+j,n) -= jinv(k,i)*tmpmat(n,k*(D-1)+l)*jinv(l,j);
 			}
 		    }
 		}		  
@@ -626,41 +626,6 @@ namespace ngfem
 
 
 
-  /*template <int D>
-  class DiffOpSurfaceHesse : public DiffOp<DiffOpSurfaceHesse<D>>
-  {
-  public:
-    enum { DIM = 1 };
-    enum { DIM_SPACE = D };
-    enum { DIM_ELEMENT = D-1 };
-    enum { DIM_DMAT = D*D };
-    enum { DIFFORDER = 2 };
-    
-    static string Name() { return "hesse"; }
-    
-    static auto & Cast (const FiniteElement & fel) 
-    { return static_cast<const ScalarFiniteElement<D-1>&> (fel); }
-    
-    template <typename MIP, typename MAT>
-    static void GenerateMatrix (const FiniteElement & fel, const MIP & mip,
-                                MAT && mat, LocalHeap & lh)
-    {
-      cout << "surface hesse, generate matrix!" << endl;
-      HeapReset hr(lh);
-      Cast(fel).CalcMappedDDShape(mip, Trans(mat));
-    }
-    };*/
-
-
-
-
-
-
-
-
-
-
-