Trying to compute sigma_sq via MuyGPS.get_sigma_optim and get_analytic_sigma (for BenchmarkGP) appears to not be working, and will produce wrong results when trying to find a known sigma_sq value. See the test cases GPSigmaSqBaselineTest and GPSigmaSqOptimTest inside of tests/optimize.py. Using "l2" rather than "F2" to create the pairwise distances appears to solve the issue in the latter case, but I am at a loss as to why.
bwpriest
commented
2 years ago
Trying to compute
sigma_sqviaMuyGPS.get_sigma_optimandget_analytic_sigma(forBenchmarkGP) appears to not be working, and will produce wrong results when trying to find a knownsigma_sqvalue. See the test casesGPSigmaSqBaselineTestandGPSigmaSqOptimTestinside oftests/optimize.py. Using"l2"rather than"F2"to create the pairwise distances appears to solve the issue in the latter case, but I am at a loss as to why.