Cleanup code and do some bugfixes

bin/analyzegrid/test_2d_structured.ini

-__exec_suffix = deg2_symdiff_nonquadvec_nongradvec
-__name = /home/rene/phd/dune-les/build-optimized-clang60-python3/dune-perftool/test/sumfact/poisson/sumfact_poisson_dg_2d_deg2_symdiff_nonquadvec_nongradvec
 cells = 16 16
 deg_suffix = deg2
 diff_suffix = symdiff

bin/analyzegrid/test_2d_unstructured.ini

-__exec_suffix = deg2_symdiff_nonquadvec_nongradvec
-__name = /home/rene/phd/dune-les/build-optimized-clang60-python3/dune-perftool/test/sumfact/poisson/sumfact_poisson_dg_2d_unstructured_deg2_symdiff_nonquadvec_nongradvec
 deg_suffix = deg2
 diff_suffix = symdiff
 elementType = quadrilateral

bin/analyzegrid/test_3d_structured.ini

-__exec_suffix = deg2_symdiff_nonquadvec_nongradvec
-__name = /home/rene/phd/dune-les/build-optimized-clang60-python3/dune-perftool/test/sumfact/poisson/sumfact_poisson_dg_3d_deg2_symdiff_nonquadvec_nongradvec
 cells = 8 8 8
 deg_suffix = deg2
 diff_suffix = symdiff

bin/analyzegrid/test_3d_unstructured.ini

-__exec_suffix = deg2_symdiff_nonquadvec_nongradvec
-__name = /home/rene/phd/dune-les/build-optimized-clang60-python3/dune-perftool/test/sumfact/poisson/sumfact_poisson_dg_3d_unstructured_deg2_symdiff_nonquadvec_nongradvec
 deg_suffix = deg2
 diff_suffix = symdiff
 elementType = quadrilateral

dune/perftool/sumfact/invertgeometry.hh

@@ -5,7 +5,7 @@
 template<typename T>
 inline T invert_and_return_determinant(const T a00, const T a10, const T a01, const T a11, T inverse[4]){
   T det = a00 * a11 - a10 * a01;
-  assert (det != 0.0);
+  assert (std::abs(det) > 1e-12);
 
   inverse[0] = a11 / det;
   inverse[1] = -a10 / det;
@@ -35,7 +35,7 @@ inline T invert_and_return_determinant(const T a00, const T a10, const T a20,
   det += t12 * a12;
   det -= t14 * a11;
 
-  assert (det != 0.0);
+  assert (std::abs(det) > 1e-12);
 
   T t17 = 1.0/det;
 

python/dune/perftool/sumfact/switch.py

@@ -62,7 +62,11 @@ def decide_if_kernel_is_necessary(facedir_s, facemod_s, facedir_n, facemod_n):
             with open(filename) as csv_file:
                 csv_reader = csv.reader(csv_file, delimiter=" ")
                 for row in csv_reader:
-                    if (facedir_s == int(row[1])) and (facemod_s == int(row[2])) and (facedir_n == int(row[3])) and (facemod_n == int(row[4])):
+                    if (row[0] == 'skeleton' and
+                            facedir_s == int(row[1])) and \
+                            (facemod_s == int(row[2])) and \
+                            (facedir_n == int(row[3])) and \
+                            (facemod_n == int(row[4])):
                         return True
                 return False
         else:
@@ -110,7 +114,7 @@ def generate_exterior_facet_switch():
                                                                               ),
                                                               args))
 
-    block.append("    default: assert(false);")
+    block.append("    default: throw(std::exception());")
     block.append("  }")
     block.append("}")
 
@@ -143,7 +147,7 @@ def generate_interior_facet_switch():
                                                                                           ),
                                                                           args))
 
-    block.append("    default: assert(false);")
+    block.append("    default: throw(std::exception());")
     block.append("  }")
     block.append("}")