A Gaussian process regression using a Gaussian kernel for both one-sample and two-sample cases. Includes non-stationary Gaussian kernel (exponential decay function) and several likelihood ratio tests for differential testing along target points.