I’ve been looking for an efficient mean-preserving calculation for quadratic interpolation. It doesn’t seem to exist. So I tried a bare-bones basic system of linear equations in python and it’s numerically unstable even if the input array is completely flat.