diff --git a/test/testGWCorrelation.cpp b/test/testGWCorrelation.cpp index 5c570b7f..be17a3e8 100644 --- a/test/testGWCorrelation.cpp +++ b/test/testGWCorrelation.cpp @@ -55,32 +55,22 @@ TEST_CASE("GWCorrelation: londonhp100") vec p = {0.0, 0.25, 0.5, 0.75, 1.0}; - // mat localcov_q0 = { - // {-59963.99642, 180641.0455, -8.952112965, -34.03759576}, - // {-55755.97844, 218271.5560, -5.884950904, -3.34033709}, - // {-53872.53438, 243590.8673, -3.185232937, 26.98199512}, - // {-46512.68133, 264445.5274, 2.355410295, 41.15969044}, - // {-36165.71804, 288239.6796, 9.330908683, 52.39072918}}; - - // mat localcov_q = quantile(algorithm.localCov(), p, 0); - // REQUIRE(approx_equal(localcov_q, localcov_q0, "absdiff", 1e-1)); - mat localcorr_q0 = { - {-0.3206001836, 0.2030111401, -0.12688297645, -0.08925682044}, - {-0.2968839387, 0.2466674205, -0.08536051649, -0.01004434385}, - {-0.2544518512, 0.2828306292, -0.04667167176, 0.08606678349}, - {-0.2409211543, 0.3241799401, 0.03032820334, 0.14221392650}, - {-0.2019074966, 0.3526574868, 0.10675955887, 0.16175140462}}; + {-0.3206001836, -0.12688297645, 0.2030111401,-0.08925682044}, + {-0.2974903966, -0.08555669604, 0.2465607269,-0.01080456730}, + {-0.2544518512, -0.04667167176, 0.2828306292, 0.08606678349}, + {-0.2408802856, 0.03056310093, 0.3243000920, 0.14224438342}, + {-0.2019074966, 0.10675955887, 0.3526574868, 0.16175140462}}; mat localcorr_q = quantile(algorithm.localCorr(), p, 0); REQUIRE(approx_equal(localcorr_q, localcorr_q0, "absdiff", 1e-3)); mat localscorr_q0 = { - {-0.3865373154, 0.2720981003, -0.1329130573468, -0.070690496147}, - {-0.3675980290, 0.2823359549, -0.1003033138009, 0.003969910873}, - {-0.3338697101, 0.3360144608, -0.0756778419096, 0.073938787835}, - {-0.3145486110, 0.3580575086, -0.0002331551565, 0.108199984250}, - {-0.2965442864, 0.3807852261, 0.0690739762835, 0.170298974147}}; + {-0.3865373154, -0.1329130573468, 0.2720981003, -0.070690496147}, + {-0.3675980290, -0.1003033138009, 0.2823359549, 0.003969910873}, + {-0.3338697101, -0.0756778419096, 0.3360144608, 0.073938787835}, + {-0.3145486110, -0.0002331551565, 0.3580575086, 0.108199984250}, + {-0.2965442864, 0.0690739762835, 0.3807852261, 0.170298974147}}; mat localscorr_q = quantile(algorithm.localSCorr(), p, 0); localscorr_q.print();