This paper proposes a new algorithm to speed-up the GPR computation:
https://arxiv.org/pdf/1703.09710.pdf
And they have a python code that do that for us:
http://celerite.readthedocs.io/en/latest/
I will try to implement that in PyCS and see how it compares to the current implementation, in terms of speed and efficiency