Identifying the dynamic gene regulatory network during latent HIV-1 reactivation using high-dimensional ordinary differential equations