orientations and position of the LA and SA sections are optimized using an optimization algorithm based on the assumption of a fixed apex.

5.  **Spatial Interpolation**\
    The available optimized cross-sections are interpolated using a bivariate method (depending on the available sections) to generate surface points at the chosen spatial resolution.

6.  **Temporal Smoothing**\
    To make the geometry suitable for CFD applications, the volume flux is smoothed in time using the Fourier curve fitting method considering the periodic nature of the dependent variable.

7.  **Mesh Generation and Output**\
    The final surface points are turned into a triangular mesh as needed by the immersed boundary (IB) node classification algorithm in our in-house the CFD solver[@borazjani2008curvilinear].

The following optimizations are performed on the previous reconstruction method [@rajan2016automated] to improve the accuracy and as well as incorporating physiological assumptions:

Data smoothing {#Data smoothing}
--------------

Rajan et al. [@rajan2016automated] used a b-spline curve fitting method which resulted in small non-physical ringing near the basal section in final reconstructed LV’s surface. To get rid of this ringing in the reconstructed geometry, the weighted moving average (WMA), and finally, the variable span weighted moving average (VSWMA) have been employed. Using WMA method, data at each node is 