This paper attempts to provide and explain a particle level set library to help complex simulations. The particle level set method was to solve the numerical dissipation problem. The dissipation problem led to using explicit particles to correct the dissipation of the implicit surface. Fast first order accurate semi-lagrangian advection method and velocity field was used for updating the level set. Particles were introduced within an interior band of the interface and used to correct volume losses. They were update using the velocity field using RK2 time integration. Error corrections using particles radius and position is done on the interface. Then fast marching method was used to re-initialize the grid. In the end, particle re-initialization is done if the distribution is too uneven and the gains of re-initialization outweigh the errors that it introduces.