Interpolation with varying means

[e-eight]

Hi,

I am trying to use this package to calculate the truncation errors for some observables calculated with chiral EFT. I am following the examples given in the documentation.

My issue is that the mean varies from one order to the next. If I use center = 0 in the TruncationGP then while interpolating, it predicts the mean to be zero at the interpolation, which does not make sense. I think this means that I need to set the center to some nonzero value. But I am not sure how to ensure that it has the appropriate value for each order.

Any help with this is appreciated.