Continuity constrained least-squares interpolation for SFO suppression in immersed boundary methods