In immersed interface methods, solids in a fluid are represented by singular forces in the Navier-Stokes equations, and flow jump conditions induced by the singular forces directly enter into numerical schemes. This paper focuses on the implementation of an immersed interface method for simulating fluid-solid interaction in 3D. The method employs the MAC scheme for the spatial discretization, the RK4 scheme for the time integration, and an FFT-based Poisson solver for the pressure Poisson equation. A fluid-solid interface is tracked by Lagrangian markers. Intersections of the interface with MAC grid lines identify finite difference stencils on which jump contributions to finite difference schemes are needed. To find the intersections and to interpolate jump conditions from the Lagrangian markers to the intersections, parametric triangulation of the interface is used. The velocity of the Lagrangian markers is interpolated directly from surrounding MAC grid nodes with interpolation schemes accounting for jump conditions. Numerical examples demonstrate that (1) the method has near second-order accuracy in the infinity norm for velocity, and the accuracy for pressure is between first and second order; (2) the method conserves the volume enclosed by a no-penetration boundary; and (3) the method can efficiently handle multiple moving solids with ease.