fix bc

Author: yolhan83

src/1DModel/bc1d.jl

@@ -22,13 +22,17 @@ function boundary_condition_outflow(u_inner, orientation_or_normal,
     eq::BloodFlowEquations1D)
         # calculate the boundary flux
         if iseven(direction) # u_inner is "left" of boundary, u_boundary is "right" of boundary
-            flux = surface_flux_function(u_inner, u_inner, orientation_or_normal,
+            flux1 = surface_flux_function[1](u_inner, u_inner, orientation_or_normal,
+            eq)
+            flux2 = surface_flux_function[2](u_inner, u_inner, orientation_or_normal,
             eq)
             else # u_inner is "left" of boundary, u_inner is "right" of boundary
-            flux = surface_flux_function(u_inner, u_inner, orientation_or_normal,
+            flux1 = surface_flux_function[1](u_inner, u_inner, orientation_or_normal,
+            eq)
+            flux2 = surface_flux_function[2](u_inner, u_inner, orientation_or_normal,
             eq)
             end
-    return flux
+    return flux1,flux2
 end
 
 @doc raw"""
@@ -61,12 +65,16 @@ function boundary_condition_slip_wall(u_inner, orientation_or_normal,
 
     # calculate the boundary flux
     if iseven(direction) # u_inner is "left" of boundary, u_boundary is "right" of boundary
-    flux = surface_flux_function(u_inner, u_boundary, orientation_or_normal,
+    flux1 = surface_flux_function[1](u_inner, u_boundary, orientation_or_normal,
+    eq)
+    flux2 = surface_flux_function[2](u_inner, u_boundary, orientation_or_normal,
     eq)
     else # u_boundary is "left" of boundary, u_inner is "right" of boundary
-    flux = surface_flux_function(u_boundary, u_inner, orientation_or_normal,
+    flux1 = surface_flux_function[1](u_boundary, u_inner, orientation_or_normal,
+    eq)
+    flux2 = surface_flux_function[2](u_boundary, u_inner, orientation_or_normal,
     eq)
     end
 
-    return flux
+    return flux1,flux2
 end

src/2DModel/bc2d.jl

@@ -17,8 +17,9 @@ Boundary flux as an `SVector`.
 """
 function boundary_condition_outflow(u_inner, orientation_or_normal, direction, x, t, surface_flux_function, eq::BloodFlowEquations2D)
     # Calculate the boundary flux without reflection
-    flux = surface_flux_function(u_inner, u_inner, orientation_or_normal, eq)
-    return flux
+    flux1 = surface_flux_function[1](u_inner, u_inner, orientation_or_normal, eq)
+    flux2 = surface_flux_function[2](u_inner, u_inner, orientation_or_normal, eq)
+    return flux1,flux2
 end
 
 
@@ -40,8 +41,9 @@ Boundary flux as an `SVector`.
 """
 function boundary_condition_outflow(u_inner, orientation_or_normal, x, t, surface_flux_function, eq::BloodFlowEquations2D)
     # Calculate the boundary flux without reflection
-    flux = surface_flux_function(u_inner, u_inner, orientation_or_normal, eq)
-    return flux
+    flux1 = surface_flux_function[1](u_inner, u_inner, orientation_or_normal, eq)
+    flux2 = surface_flux_function[2](u_inner, u_inner, orientation_or_normal, eq)
+    return flux1,flux2
 end
 
 
@@ -68,9 +70,11 @@ function boundary_condition_slip_wall(u_inner, orientation_or_normal, direction,
 
     # Calculate the boundary flux based on direction
     if iseven(direction)
-        flux = surface_flux_function(u_inner, u_boundary, orientation_or_normal, eq)
+        flux1 = surface_flux_function[1](u_inner, u_boundary, orientation_or_normal, eq)
+        flux2 = surface_flux_function[2](u_inner, u_boundary, orientation_or_normal, eq)
     else
-        flux = surface_flux_function(u_boundary, u_inner, orientation_or_normal, eq)
+        flux1 = surface_flux_function[1](u_boundary, u_inner, orientation_or_normal, eq)
+        flux2 = surface_flux_function[2](u_boundary, u_inner, orientation_or_normal, eq)
     end
-    return flux
+    return flux1,flux2
 end