From 236c5eff565f13d44dfc747a916e78cdb9f11798 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Thu, 1 Jun 2023 14:23:22 +0200 Subject: Add clen refinement --- libcrystfel/src/predict-refine.c | 14 ++++++++++---- 1 file changed, 10 insertions(+), 4 deletions(-) (limited to 'libcrystfel/src') 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); } -- cgit v1.2.3