Random seed set to 2 while building models
Initial experimental model
==========================
beam s0 = (-0.0237, 0.0000, 0.8473)
sensor origin = (-126.7316, -144.5101, 200.1049)
sensor dir1 = (1.0000, 0.0007, -0.0012)
sensor dir2 = (-0.0007, 1.0000, 0.0003)
crystal unit cell = 14.7802, 13.3487, 12.9060, 90.1647, 90.2195, 89.4948
crystal orientation matrix U =
{{1.0, -0.0007, 0.0004},
 {0.0007, 1.0, 0.0041},
 {-0.0004, -0.0041, 1.0}}
Offsetting initial model
========================
beam s0 = (-0.0220, 0.0000, 0.8473)
sensor origin = (-124.4459, -142.7587, 202.0671)
sensor dir1 = (1.0000, 0.0027, -0.0032)
sensor dir2 = (-0.0027, 1.0000, 0.0023)
crystal unit cell = 14.7802, 13.3487, 12.9060, 90.1647, 90.2195, 89.4948
crystal orientation matrix U =
{{1.0, -0.0027, 0.0024},
 {0.0027, 1.0, 0.0021},
 {-0.0024, -0.0021, 1.0}}

Generating reflections
======================
Total number of reflections excited 1284
Total number of observations made 1284

Resetting to initial model
==========================
Initial values of parameters are
Parameters: 199.99987 -0.22230 0.03775 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 457.80561 561.25787 600.37861 -4.46370 1.99534 1.64939

Calculated gradients
====================
Target L           = 588118.460110
Gradients dL/dp    = -92177.264377 -250951.138479 -303829.717515 799.277958 11112.449592 -4869.690210 -68809.420135 56680.886687 -3136.023377 -1154.960708 -1204.210077 -6205.332469 -6870.086975 -6546.509675 -7243.154566 2035.237595 -2643.155958
Curvatures d2L/dp2 = 42457.517023 159004.867496 159004.867496 1698.298522 383.930599 161.385149 10694.405558 10633.430498 1544.232107 1637.030382 1589.660711 1015.177532 679.818385 613.181641 3319.089278 3146.308879 2370.417425

Finite difference gradients
===========================
Parameter 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

FD gradients dL[fd]/dp          = -92177.137267 -250951.137859 -303829.711629 799.560221 11112.580542 -4869.733239 -68809.394725 56680.864654 -3136.026207 -1154.956408 -1204.215223 -6205.311511 -6870.081415 -6546.569057 -7243.087748 2035.299549 -2643.212210
dL/dp - dL[fd]/dp               = -0.127110 -0.000619 -0.005886 -0.282262 -0.130950 0.043029 -0.025410 0.022033 0.002831 -0.004300 0.005146 -0.020959 -0.005559 0.059381 -0.066818 -0.061954 0.056252
Normalised differences:
(dL/dp - dL[fd]/dp) / dL[fd]/dp = 0.000001 0.000000 0.000000 -0.000353 -0.000012 -0.000009 0.000000 0.000000 -0.000001 0.000004 -0.000004 0.000003 0.000001 -0.000009 0.000009 -0.000030 -0.000021

Finite difference curvatures
============================
FD curvatures d2L[fd]/dp2             = 42457.516683 159004.867601 159004.867892 1716.735123 363.663494 150.963997 10657.323146 10600.638998 1544.232150 1646.877711 1598.682079 1019.734364 681.849597 614.901110 3324.500976 3153.338334 2372.854033
d2L/dp2 - d2L[fd]/dp2                 = 0.000340 -0.000105 -0.000396 -18.436601 20.267105 10.421153 37.082412 32.791501 -0.000043 -9.847328 -9.021368 -4.556831 -2.031211 -1.719469 -5.411699 -7.029456 -2.436608
Normalised differences:
(d2L/dp2 - d2L[fd]/dp2) / d2L[fd]/dp2 = 0.000000 -0.000000 -0.000000 -0.010739 0.055730 0.069031 0.003480 0.003093 -0.000000 -0.005979 -0.005643 -0.004469 -0.002979 -0.002796 -0.001628 -0.002229 -0.001027
OK