#704 added option to finite element boundary_value_or_flux

Author: Scottmar93

pybamm/spatial_methods/scikit_finite_element.py

@@ -322,7 +322,7 @@ def integral_form(v, dv, w):
 
         return pybamm.Matrix(integration_vector[np.newaxis, :])
 
-    def boundary_value_or_flux(self, symbol, discretised_child):
+    def boundary_value_or_flux(self, symbol, discretised_child, extrapolation="linear"):
         """
         Returns the average value of the symbol over the negative tab ("negative tab")
         or the positive tab ("positive tab") in the Finite Element Method.