We show comparison of four different numerical methods for simulating Photonic-Crystal (PC) VCSELs. We present the theoretical basis behind each method and analyze the differences by studying a benchmark VCSEL structure, where the PC structure penetrates all VCSEL layers, the entire top-mirror DBR, a fraction of the top-mirror DBR or just the VCSEL cavity. The different models are evaluated by comparing the predicted resonance wavelengths and threshold gains for different hole diameters and pitches of the PC. The agreement between the models is relatively good, except for one model, which corresponds to the effective index method. The simulation results elucidate the strength and weaknesses of the analyzed methods; and outline the limits of applicability of the different models.