diff options
| author | Thomas White <taw@physics.org> | 2023-06-01 14:23:22 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2023-07-28 13:22:05 +0200 |
| commit | 236c5eff565f13d44dfc747a916e78cdb9f11798 (patch) | |
| tree | 6eb137d10cb14e51537c5bf1547d68add8cfc82b /libcrystfel/src | |
| parent | 535957ea56cb7c2624452ce64da4c7020903572e (diff) | |
Add clen refinement
Diffstat (limited to 'libcrystfel/src')
| -rw-r--r-- | libcrystfel/src/predict-refine.c | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/libcrystfel/src/predict-refine.c b/libcrystfel/src/predict-refine.c index af58b5f7..88dfbf36 100644 --- a/libcrystfel/src/predict-refine.c +++ b/libcrystfel/src/predict-refine.c @@ -580,8 +580,8 @@ static void write_mille(Mille *mille, int n, UnitCell *cell, #ifdef HAVE_MILLEPEDE int i; float local_gradients[9]; - float global_gradients[2]; - int labels[2]; + float global_gradients[6]; + int labels[6]; profile_start("mille-calc"); @@ -618,12 +618,15 @@ static void write_mille(Mille *mille, int n, UnitCell *cell, cell, rps[i].panel); global_gradients[1] = x_gradient(GPARAM_DETY, rps[i].refl, cell, rps[i].panel); + global_gradients[2] = x_gradient(GPARAM_CLEN, rps[i].refl, + cell, rps[i].panel); labels[0] = 1; labels[1] = 2; + labels[2] = 3; mille_add_measurement(mille, 9, local_gradients, - 2, global_gradients, labels, + 3, global_gradients, labels, x_dev(&rps[i], image->detgeom, dx, dy), 0.65*rps[i].panel->pixel_pitch); } @@ -641,12 +644,15 @@ static void write_mille(Mille *mille, int n, UnitCell *cell, cell, rps[i].panel); global_gradients[1] = y_gradient(GPARAM_DETY, rps[i].refl, cell, rps[i].panel); + global_gradients[2] = y_gradient(GPARAM_CLEN, rps[i].refl, + cell, rps[i].panel); labels[0] = 1; labels[1] = 2; + labels[2] = 3; mille_add_measurement(mille, 9, local_gradients, - 2, global_gradients, labels, + 3, global_gradients, labels, y_dev(&rps[i], image->detgeom, dx, dy), 0.65*rps[i].panel->pixel_pitch); } |