Limit-cycle-preserving simulation of gene regulatory oscillators (Q1936031)

Summary: In order to simulate gene regulatory oscillators more effectively, Runge-Kutta (RK) integrators are adapted to the limit-cycle structure of the system. Taking into account the oscillatory features of gene regulatory oscillators, phase-fitted and amplification-fitted Runge-Kutta (FRK) methods are designed. New FRK methods with phase-fitted and amplification-fitted updates are also considered. The error coefficients and error constants for each of the new FRK methods are obtained. In the numerical simulations of the two-gene regulatory system, the new methods are shown to be more accurate and more efficient than their prototype RK methods in the long-term integration. It is a new discovery that the best fitting frequency not only depends on the problem to be solved, but also depends on the method.