As the reactor physics computation develops,the neutron transport equation would be solved on three dimensional complex geometry accurately,where the MOC can meet all these needs.However,the slow convergence and long computation time limit its application.Both numerical acceleration method and parallel computing technique are used in this paper to accelerate the MOC.Two-level GCMFD and angular parallelization are implemented into the three dimensional characteristics method code TCM to carry out test calculations in some benchmark problems.The results show a good performance on reducing the number of iterations and speeding up the code while keeping the accuracy of the computation.