def update_gradients_diag(self, dL_dKdiag, X):
        self.variance.gradient = np.sum(dL_dKdiag)
        self.lengthscale.gradient = 0.
Take a sample of data points which are drawn from a power law probability
dataset/ETHPy150Open SheffieldML/GPy/GPy/kern/src/stationary.py/Stationary.update_gradients_diag--crossover with--dataset/ETHPy150Open charanpald/APGL/apgl/util/Util.py/Util.fitPowerLaw