Manifold-constrained Gaussian process inference for ODE parameters, by profiling the states out and tempering the parameters onto them