Merge branch 'hdivsurface' into 'master'

Hdivsurface

See merge request jschoeberl/ngsolve!282

comp/CMakeLists.txt

@@ -12,7 +12,7 @@ add_library( ngcomp ${NGS_LIB_TYPE}
         ${ngcomp_object_libs}
         bilinearform.cpp facetfespace.cpp fespace.cpp 
         gridfunction.cpp h1hofespace.cpp hcurlhdivfes.cpp hcurlhofespace.cpp 
-        hdivfes.cpp hdivhofespace.cpp hierarchicalee.cpp l2hofespace.cpp     
+        hdivfes.cpp hdivhofespace.cpp hdivhosurfacefespace.cpp hierarchicalee.cpp l2hofespace.cpp     
         linearform.cpp meshaccess.cpp ngsobject.cpp postproc.cpp	     
         preconditioner.cpp vectorfacetfespace.cpp numberfespace.cpp bddc.cpp 
         hypre_precond.cpp hdivdivfespace.cpp hdivdivsurfacespace.cpp tpfes.cpp 
@@ -34,7 +34,7 @@ endif(NOT WIN32)
 install( FILES
         bilinearform.hpp comp.hpp facetfespace.hpp	   
         fespace.hpp gridfunction.hpp h1hofespace.hpp hcurlhdivfes.hpp	   
-        hcurlhofespace.hpp hdivfes.hpp hdivhofespace.hpp		   
+        hcurlhofespace.hpp hdivfes.hpp hdivhofespace.hpp hdivhosurfacefespace.hpp		   	   
         l2hofespace.hpp hdivdivsurfacespace.hpp tpfes.hpp linearform.hpp meshaccess.hpp ngsobject.hpp	   
         postproc.hpp preconditioner.hpp vectorfacetfespace.hpp hypre_precond.hpp 
         pde.hpp numproc.hpp vtkoutput.hpp pmltrafo.hpp periodic.hpp  hypre_ams_precond.hpp facetsurffespace.hpp

comp/comp.hpp

@@ -67,7 +67,8 @@ namespace ngcomp
 #include "hdivfes.hpp"
 #include "h1hofespace.hpp"
 #include "l2hofespace.hpp"
-#include "hdivhofespace.hpp" 
+#include "hdivhofespace.hpp"
+#include "hdivhosurfacefespace.hpp" 
 #include "hcurlhofespace.hpp" 
 #include "facetfespace.hpp" 
 #include "vectorfacetfespace.hpp"

comp/facetsurffespace.cpp

@@ -302,7 +302,9 @@ namespace ngcomp
 	  break;
 	}
       case BBBND:
-        throw Exception("No BBBND GetFE implemented for FacetSurfaceFESpace");
+        //throw Exception("No BBBND GetFE implemented for FacetSurfaceFESpace");
+	return * new (lh) DummyFE<ET_POINT>();
+	break;
       }
   }
 
@@ -342,7 +344,13 @@ namespace ngcomp
 	    }
 	}
 	break;
-      case BBND: case BBBND:
+      case BBND:
+	{
+	  dnums += GetEdgeDofs(ma->GetElEdges(ei)[0]);
+	}
+	break;
+
+      case BBBND:
 	break;
       }
   }