Probe direction eta (np.random.seed=42, eta = np.random.randn(27) * 0.1): eta[ 0] = +0.0496714153 eta[ 1] = -0.0138264301 eta[ 2] = +0.0647688538 eta[ 3] = +0.1523029856 eta[ 4] = -0.0234153375 eta[ 5] = -0.0234136957 eta[ 6] = +0.1579212816 eta[ 7] = +0.0767434729 eta[ 8] = -0.0469474386 eta[ 9] = +0.0542560044 eta[10] = -0.0463417693 eta[11] = -0.0465729754 eta[12] = +0.0241962272 eta[13] = -0.1913280245 eta[14] = -0.1724917833 eta[15] = -0.0562287529 eta[16] = -0.1012831120 eta[17] = +0.0314247333 eta[18] = -0.0908024076 eta[19] = -0.1412303701 eta[20] = +0.1465648769 eta[21] = -0.0225776300 eta[22] = +0.0067528205 eta[23] = -0.1424748186 eta[24] = -0.0544382725 eta[25] = +0.0110922590 eta[26] = -0.1150993577 ==== Regge cubic along seed-42 direction ==== h cubic estimate 1e-02 -0.2392853326 1e-03 -0.2392893165 1e-04 -0.2131404429 1e-05 13.8581331095 Richardson extrapolations (assuming O(h^2) truncation): h=1e-02 -> h=1e-03: -0.2392893568 h=1e-03 -> h=1e-04: -0.2128763129 h=1e-04 -> h=1e-05: 14.0002671858 ==== J-cost cubic along seed-42 direction (must be ~0) ==== h cubic estimate 1e-02 0.0000000000 1e-03 0.0000000000 1e-04 0.0000000000 1e-05 0.0000000000 Richardson extrapolations (assuming O(h^2) truncation): h=1e-02 -> h=1e-03: 0.0000000000 h=1e-03 -> h=1e-04: 0.0000000000 h=1e-04 -> h=1e-05: 0.0000000000 ==== Regge cubic along constant direction (must be ~0) ==== h cubic estimate 1e-02 -0.0000001873 1e-03 0.0002941473 1e-04 0.1643626063 1e-05 -232.4834327233 Richardson extrapolations (assuming O(h^2) truncation): h=1e-02 -> h=1e-03: 0.0002971204 h=1e-03 -> h=1e-04: 0.1660198635 h=1e-04 -> h=1e-05: -234.8334104539 ==== J-cost cubic along constant direction (must be 0 identically) ==== h cubic estimate 1e-02 0.0000000000 1e-03 0.0000000000 1e-04 0.0000000000 1e-05 0.0000000000 Richardson extrapolations (assuming O(h^2) truncation): h=1e-02 -> h=1e-03: 0.0000000000 h=1e-03 -> h=1e-04: 0.0000000000 h=1e-04 -> h=1e-05: 0.0000000000