Some implementation details of fully Lagrangian vortex loops method algorithm are presented. The vortex loops method can be applied to various engineering applications, where it is required to simulate the flow around the body of arbitrary shape and to estimate hydrodynamic loads. The original approach is implemented that permits to solve boundary integral equation with rather high accuracy at coarse and non-uniform body surface triangulation. Vorticity distribution in the flow is simulated by closed vortex loops. The presented algorithm allows to simulate movement of the vortex loops, which simulate vorticity distribution in the flow, their stretching, reconnection, etc. Note, that the viscosity influence is taken into account only through vorticity generation on the body surface.