The Monte Carlo method is the most commonly used dose calculation method in the field of boron neutron capture therapy (BNCT). General-purpose Monte Carlo (MC) code (e.g., MCNP) has been used in most treatment planning systems (TPS) to calculate dose distribution, which takes overmuch time in radiotherapy planning. Based on this, we developed COMPASS (COMpact PArticle Simulation System), an MC engine specifically for BNCT dose calculation. Several optimization algorithms are used in COMPASS to make it faster than general-purpose MC code. The parallel computation of COMPASS is performed by the message passing interface (MPI) library and OpenMP commands, which allows the user to increase computational speed by increasing the computer configurations. The physical dose of each voxel is calculated for developing a treatment plan. Comparison results show that the computed dose distribution of COMPASS is in good agreement with MCNP, and the computational efficiency is better than MCNP. These results validate that COMPASS has better performance than MCNP in BNCT dose calculation.