Remove compute_flow_rate function and use the one in VaMPy instead

src/fsipy/automatedPreprocessing/automated_preprocessing.py

@@ -14,13 +14,13 @@
 from vampy.automatedPreprocessing import ToolRepairSTL
 from vampy.automatedPreprocessing.preprocessing_common import read_polydata, get_centers_for_meshing, \
     dist_sphere_diam, dist_sphere_curvature, dist_sphere_constant, get_regions_to_refine, add_flow_extension, \
-    write_mesh, mesh_alternative, find_boundaries, setup_model_network, \
+    write_mesh, mesh_alternative, find_boundaries, compute_flow_rate, setup_model_network, \
     radiusArrayName, scale_surface, get_furtest_surface_point, check_if_closed_surface
 from vampy.automatedPreprocessing.simulate import run_simulation
 from vampy.automatedPreprocessing.visualize import visualize_model
 
 from fsipy.automatedPreprocessing.preprocessing_common import generate_mesh, distance_to_spheres_solid_thickness, \
-    dist_sphere_spheres, convert_xml_mesh_to_hdf5, convert_vtu_mesh_to_xdmf, compute_flow_rate
+    dist_sphere_spheres, convert_xml_mesh_to_hdf5, convert_vtu_mesh_to_xdmf
 
 
 def run_pre_processing(input_model, verbose_print, smoothing_method, smoothing_factor, smoothing_iterations,

src/fsipy/automatedPreprocessing/preprocessing_common.py

@@ -1,4 +1,3 @@
-from typing import List, Dict
 from pathlib import Path
 
 import meshio
@@ -251,32 +250,3 @@ def convert_vtu_mesh_to_xdmf(file_name_vtu_mesh: str, file_name_xdmf_mesh: str)
 
     print(f"Tetra mesh XDMF file written to: {tetra_xdmf_path}")
     print(f"Triangle mesh XDMF file written to: {triangle_xdmf_path}\n")
-
-
-def compute_flow_rate(is_atrium: bool, inlet: List[float], parameters: Dict[str, float],
-                      flow_rate_factor: float) -> float:
-    """
-    Computes mean flow rate used as a boundary condition for pulsatile flow condition
-
-    Args:
-        is_atrium (bool): Determines if model is atrium or artery
-        inlet (list): List of points representing midpoint of boundaries
-        parameters (dict): Dictionary containing model parameters
-        flow_rate_factor (float): Factor to adjust flow rate calculation
-
-    Returns:
-        mean_inflow_rate (float): Mean inflow rate
-    """
-    if is_atrium:
-        # Calculate total inlet area for atrium case
-        total_inlet_area = 0.0
-        num_inlets = len(inlet)
-        for i in range(num_inlets):
-            inlet_area_key = "inlet{}_area".format(i)
-            total_inlet_area += parameters.get(inlet_area_key, 0)
-        mean_inflow_rate = flow_rate_factor * total_inlet_area
-    else:
-        # Use single inlet area for artery case
-        mean_inflow_rate = flow_rate_factor * parameters.get("inlet_area", 0)
-
-    return mean_inflow_rate