This thesis falls within the framework of advancing numerical computation techniques for reactor physics. Specifically, it focuses on the implementation of methods that incorporate higher-order spatial expansions for neutron flux and cross-sections. The primary objective is to accelerate both existing algorithms and those that will be developed through GPU programming. By harnessing the computational power of GPUs, this research aims to enhance the efficiency and accuracy of simulations in reactor physics, thereby contributing to the broader field of nuclear engineering and safety.