[en] The collision less reconnection process in Tokamaks due to the nonlinear development of m=1 (poloidal mode number) and n=1 (toroidal mode number) kinetic internal kink mode is simulated by the gyro-particle code (GYR3D), the gyro-reduced MHD code (GRM3D-2F), and the particle-fluid hybrid code (Hybrid3D). These codes are based on the nonlinear gyrokinetic Vlasov-Poisson-Ampere system and have exact energy invariance. GYR3D is a three-dimensional gyrokinetic magneto-inductive particle code with δf method. GRM3D-2F is a two-field and two-fluid model including the effects of electron inertia and the perturbed electron pressure gradients along the magnetic field. In Hybrid3D, electrons are treated as fluid, while ions are treated as particles. The results of these three codes agreed very well. We believe that the better understanding of the physics associated with the kinetic MHD phenomena in Tokamaks will be achieved by executing simultaneously these codes. (author)