raise ValueError for Gauss order out of the range (#220)

Co-authored-by: Matthew Scroggs <[email protected]>

bempp/api/integration/gauss.py

@@ -1045,6 +1045,10 @@ def rule(order):
 
     Orders between 1 and 30 are supported.
     """
+    if order < 1 or order > 30:
+        raise ValueError(
+            f"Gauss quadrature order must be between 1 and 30. Provided: {order}"
+        )
     npoints = order
     address = (npoints * (npoints - 1)) // 2
     return (

bempp/api/integration/triangle_gauss.py

@@ -2931,6 +2931,10 @@ def rule(order):
 
     The order must be between 1 and 20.
     """
+    if order < 1 or order > 20:
+        raise ValueError(
+            f"Symmetric Gauss quadrature order must be between 1 and 20. Provided: {order}"
+        )
     npoints = points_per_order[order - 1]
     address = points_address[npoints - 1]
     bary_coords = coords[3 * address : 3 * (address + npoints)].reshape(