Factorise dr/da part of gradient calculation for PR and prediction

1 files changed, 27 insertions, 27 deletions

diff --git a/tests/prediction_gradient_check.c b/tests/prediction_gradient_check.c
index 486e4ce1..21c4088e 100644
--- a/tests/prediction_gradient_check.c
+++ b/tests/prediction_gradient_check.c
@@ -110,15 +110,15 @@ static UnitCell *new_shifted_cell(UnitCell *input, int k, double shift)
 	                    &csx, &csy, &csz);
 	switch ( k )
 	{
-		case REF_ASX :  asx += shift;  break;
-		case REF_ASY :  asy += shift;  break;
-		case REF_ASZ :  asz += shift;  break;
-		case REF_BSX :  bsx += shift;  break;
-		case REF_BSY :  bsy += shift;  break;
-		case REF_BSZ :  bsz += shift;  break;
-		case REF_CSX :  csx += shift;  break;
-		case REF_CSY :  csy += shift;  break;
-		case REF_CSZ :  csz += shift;  break;
+		case GPARAM_ASX :  asx += shift;  break;
+		case GPARAM_ASY :  asy += shift;  break;
+		case GPARAM_ASZ :  asz += shift;  break;
+		case GPARAM_BSX :  bsx += shift;  break;
+		case GPARAM_BSY :  bsy += shift;  break;
+		case GPARAM_BSZ :  bsz += shift;  break;
+		case GPARAM_CSX :  csx += shift;  break;
+		case GPARAM_CSY :  csy += shift;  break;
+		case GPARAM_CSZ :  csz += shift;  break;
 	}
 	cell_set_reciprocal(cell, asx, asy, asz, bsx, bsy, bsz, csx, csy, csz);
 
@@ -141,15 +141,15 @@ static Crystal *new_shifted_crystal(Crystal *cr, int refine, double incr_val)
 
 	switch ( refine ) {
 
-		case REF_ASX :
-		case REF_ASY :
-		case REF_ASZ :
-		case REF_BSX :
-		case REF_BSY :
-		case REF_BSZ :
-		case REF_CSX :
-		case REF_CSY :
-		case REF_CSZ :
+		case GPARAM_ASX :
+		case GPARAM_ASY :
+		case GPARAM_ASZ :
+		case GPARAM_BSX :
+		case GPARAM_BSY :
+		case GPARAM_BSZ :
+		case GPARAM_CSX :
+		case GPARAM_CSY :
+		case GPARAM_CSZ :
 		cell = new_shifted_cell(crystal_get_cell(cr), refine,
 		                        incr_val);
 		crystal_set_cell(cr_new, cell);
@@ -468,41 +468,41 @@ int main(int argc, char *argv[])
 		                    &bx, &by, &bz, &cx, &cy, &cz);
 
 		incr_val = incr_frac * ax;
-		val = test_gradients(cr, incr_val, REF_ASX, "ax*", "ax",
+		val = test_gradients(cr, incr_val, GPARAM_ASX, "ax*", "ax",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * bx;
-		val = test_gradients(cr, incr_val, REF_BSX, "bx*", "bx",
+		val = test_gradients(cr, incr_val, GPARAM_BSX, "bx*", "bx",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * cx;
-		val = test_gradients(cr, incr_val, REF_CSX, "cx*", "cx",
+		val = test_gradients(cr, incr_val, GPARAM_CSX, "cx*", "cx",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 
 		incr_val = incr_frac * ay;
-		val = test_gradients(cr, incr_val, REF_ASY, "ay*", "ay",
+		val = test_gradients(cr, incr_val, GPARAM_ASY, "ay*", "ay",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * by;
-		val = test_gradients(cr, incr_val, REF_BSY, "by*", "by",
+		val = test_gradients(cr, incr_val, GPARAM_BSY, "by*", "by",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * cy;
-		val = test_gradients(cr, incr_val, REF_CSY, "cy*", "cy",
+		val = test_gradients(cr, incr_val, GPARAM_CSY, "cy*", "cy",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 
 		incr_val = incr_frac * az;
-		val = test_gradients(cr, incr_val, REF_ASZ, "az*", "az",
+		val = test_gradients(cr, incr_val, GPARAM_ASZ, "az*", "az",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * bz;
-		val = test_gradients(cr, incr_val, REF_BSZ, "bz*", "bz",
+		val = test_gradients(cr, incr_val, GPARAM_BSZ, "bz*", "bz",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;
 		incr_val = incr_frac * cz;
-		val = test_gradients(cr, incr_val, REF_CSZ, "cz*", "cz",
+		val = test_gradients(cr, incr_val, GPARAM_CSZ, "cz*", "cz",
 		                     quiet, plot);
 		if ( val < 0.99 ) fail = 1;