New advances in the Gaussian-process approach to pulsar-timing data analysis