diff options
Diffstat (limited to 'src/template.c')
-rw-r--r-- | src/template.c | 293 |
1 files changed, 138 insertions, 155 deletions
diff --git a/src/template.c b/src/template.c index 9c012c4..a22401e 100644 --- a/src/template.c +++ b/src/template.c @@ -217,43 +217,43 @@ #ifdef HAVE_OP_I static - void compute_coeffs_I(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_I); + void compute_coeffs_I(FP_XYZ, struct COEFFS_STRUCT *coeffs_I); #endif #ifdef HAVE_OP_DX static - void compute_coeffs_dx(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dx); + void compute_coeffs_dx(FP_XYZ, struct COEFFS_STRUCT *coeffs_dx); #endif #ifdef HAVE_OP_DY static - void compute_coeffs_dy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dy); + void compute_coeffs_dy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dy); #endif #ifdef HAVE_OP_DZ static - void compute_coeffs_dz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dz); + void compute_coeffs_dz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dz); #endif #ifdef HAVE_OP_DXX static - void compute_coeffs_dxx(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxx); + void compute_coeffs_dxx(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxx); #endif #ifdef HAVE_OP_DXY static - void compute_coeffs_dxy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxy); + void compute_coeffs_dxy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxy); #endif #ifdef HAVE_OP_DXZ static - void compute_coeffs_dxz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxz); + void compute_coeffs_dxz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxz); #endif #ifdef HAVE_OP_DYY static - void compute_coeffs_dyy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dyy); + void compute_coeffs_dyy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dyy); #endif #ifdef HAVE_OP_DYZ static - void compute_coeffs_dyz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dyz); + void compute_coeffs_dyz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dyz); #endif #ifdef HAVE_OP_DZZ static - void compute_coeffs_dzz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dzz); + void compute_coeffs_dzz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dzz); #endif /******************************************************************************/ @@ -335,38 +335,38 @@ @@*/ int FUNCTION_NAME(/***** coordinate system *****/ - const CCTK_REAL *restrict coord_origin, - const CCTK_REAL *restrict coord_delta, + const CCTK_REAL coord_origin[], + const CCTK_REAL coord_delta[], /***** interpolation points *****/ int N_interp_points, int interp_coords_type_code, - const void *restrict const *restrict interp_coords, - const CCTK_INT *restrict N_boundary_points_to_omit, - const CCTK_REAL *restrict boundary_off_centering_tolerance, - const CCTK_REAL *restrict boundary_extrapolation_tolerance, + const void* const interp_coords[], + const CCTK_INT N_boundary_points_to_omit[], + const CCTK_REAL boundary_off_centering_tolerance[], + const CCTK_REAL boundary_extrapolation_tolerance[], /***** input arrays *****/ int N_input_arrays, - const CCTK_INT *restrict input_array_offsets, - const CCTK_INT *restrict input_array_strides, - const CCTK_INT *restrict input_array_min_subscripts, - const CCTK_INT *restrict input_array_max_subscripts, - const CCTK_INT *restrict input_array_type_codes, - const void *restrict const *restrict input_arrays, + const CCTK_INT input_array_offsets[], + const CCTK_INT input_array_strides[], + const CCTK_INT input_array_min_subscripts[], + const CCTK_INT input_array_max_subscripts[], + const CCTK_INT input_array_type_codes[], + const void* const input_arrays[], /***** output arrays *****/ int N_output_arrays, - const CCTK_INT *restrict output_array_type_codes, - void* const *restrict output_arrays, + const CCTK_INT output_array_type_codes[], + void* const output_arrays[], /***** operation info *****/ - const CCTK_INT *restrict operand_indices, - const CCTK_INT *restrict operation_codes, + const CCTK_INT operand_indices[], + const CCTK_INT operation_codes[], /***** debugging *****/ int debug, FILE* log_fp, /***** other return results *****/ struct error_info* error_info, struct molecule_structure_flags* molecule_structure_flags, struct molecule_min_max_m_info* molecule_min_max_m_info, - CCTK_INT *restrict const *restrict molecule_positions, - struct Jacobian_info *restrict Jacobian_info) + CCTK_INT* const molecule_positions[], + struct Jacobian_info* Jacobian_info) { /* * ***** Naming conventions ***** @@ -824,7 +824,6 @@ if (debug > 0) fflush(stdout); } -#pragma omp parallel for for (pt = 0 ; pt < N_interp_points ; ++pt) { /* this struct holds a molecule-sized piece of a single */ @@ -872,7 +871,7 @@ if (debug > 0) * from the interp_coords[] arrays * FIXME: Maybe it would be better (faster) to do this * with N_DIMS open-coded calls on a function? - * But then we'd have to have a sentinel value + * But then we'd have to have a sentinal value * return for the unknown-type-code error case. * Yuk! :( :( */ @@ -886,11 +885,16 @@ if (debug > 0) switch (interp_coords_type_code) { + case CCTK_VARIABLE_REAL: + { + const CCTK_REAL *const interp_coords_ptr_real + = (const CCTK_REAL *) interp_coords_ptr; + interp_coords_fp[axis] = interp_coords_ptr_real[pt]; + break; + } + #ifdef HAVE_CCTK_REAL4 case CCTK_VARIABLE_REAL4: - #ifdef CCTK_REAL_PRECISION_4 - case CCTK_VARIABLE_REAL: - #endif { const CCTK_REAL4 *const interp_coords_ptr_real4 = (const CCTK_REAL4 *) interp_coords_ptr; @@ -901,9 +905,6 @@ if (debug > 0) #ifdef HAVE_CCTK_REAL8 case CCTK_VARIABLE_REAL8: - #ifdef CCTK_REAL_PRECISION_8 - case CCTK_VARIABLE_REAL: - #endif { const CCTK_REAL8 *const interp_coords_ptr_real8 = (const CCTK_REAL8 *) interp_coords_ptr; @@ -914,9 +915,6 @@ if (debug > 0) #ifdef HAVE_CCTK_REAL16 case CCTK_VARIABLE_REAL16: - #ifdef CCTK_REAL_PRECISION_16 - case CCTK_VARIABLE_REAL: - #endif { /* FIXME: maybe we should warn (once per cactus run) */ /* that we may be doing arithmetic in lower */ @@ -937,9 +935,7 @@ if (debug > 0) , interp_coords_type_code); /*NOTREACHED*/ -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; /*** ERROR RETURN ***/ /* end of switch (interp_coords_type_code) */ } @@ -1156,12 +1152,8 @@ if (debug > 0) #endif } - if (this_point_status < 0) - then { -#pragma omp critical - if (this_point_status < return_status) - then return_status = this_point_status; - } + if (this_point_status < return_status) + then return_status = this_point_status; if (debug >= 6) then { @@ -1301,13 +1293,13 @@ if (debug > 0) * execute the ***load*** the first time in the test at the * top of the part loop below */ - const void *restrict input_array_ptr__last_load = NULL; + const void* input_array_ptr__last_load = NULL; int part__last_load = -1; for (out = 0 ; out < N_output_arrays ; ++out) { const int in = operand_indices[out]; - const void *restrict const input_array_ptr = input_arrays[in]; + const void* const input_array_ptr = input_arrays[in]; /* * ***decode*** the output array datatype @@ -1331,9 +1323,7 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, , out, (int) output_array_type_codes[out], N_output_parts); /*NOTREACHED*/ -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; /*** ERROR RETURN ***/ } { @@ -1371,9 +1361,8 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, in, (int) input_array_type_codes[in], N_input_parts, out, (int) output_array_type_codes[out], N_output_parts); /*NOTREACHED*/ -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; + /*** ERROR RETURN ***/ } /* @@ -1386,32 +1375,13 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, + input_array_offsets[in]; switch (input_array_type_codes[in]) { -#ifdef HAVE_CCTK_REAL4 -case CCTK_VARIABLE_REAL4: -#ifdef CCTK_REAL_PRECISION_4 -case CCTK_VARIABLE_REAL: -#endif - { - const CCTK_REAL4 *restrict const input_array_ptr_real4 - = ((const CCTK_REAL4 *) input_array_ptr) + input_posn; - LOAD_DATA_REAL4(input_array_ptr_real4, - STRIDE_IJK, - &data); - break; - } -#endif - -#ifdef HAVE_CCTK_REAL8 -case CCTK_VARIABLE_REAL8: -#ifdef CCTK_REAL_PRECISION_8 case CCTK_VARIABLE_REAL: -#endif { - const CCTK_REAL8 *restrict const input_array_ptr_real8 - = ((const CCTK_REAL8 *) input_array_ptr) + input_posn; - LOAD_DATA_REAL8(input_array_ptr_real8, - STRIDE_IJK, - &data); + const CCTK_REAL *const input_array_ptr_real + = ((const CCTK_REAL *) input_array_ptr) + input_posn; + LOAD_DATA_REAL(input_array_ptr_real, + STRIDE_IJK, + &data); #if (N_DIMS == 2) && (MOLECULE_SIZE == 4) if (debug >= 10) then { @@ -1438,17 +1408,37 @@ case CCTK_VARIABLE_REAL: #endif break; } + +#ifdef HAVE_CCTK_REAL4 +case CCTK_VARIABLE_REAL4: + { + const CCTK_REAL4 *const input_array_ptr_real4 + = ((const CCTK_REAL4 *) input_array_ptr) + input_posn; + LOAD_DATA_REAL4(input_array_ptr_real4, + STRIDE_IJK, + &data); + break; + } +#endif + +#ifdef HAVE_CCTK_REAL8 +case CCTK_VARIABLE_REAL8: + { + const CCTK_REAL8 *const input_array_ptr_real8 + = ((const CCTK_REAL8 *) input_array_ptr) + input_posn; + LOAD_DATA_REAL8(input_array_ptr_real8, + STRIDE_IJK, + &data); + break; + } #endif #ifdef HAVE_CCTK_REAL16 case CCTK_VARIABLE_REAL16: -#ifdef CCTK_REAL_PRECISION_16 -case CCTK_VARIABLE_REAL: -#endif { /* FIXME: maybe we should warn (once per cactus run) that we may be */ /* doing arithmetic in lower precision than the data type? */ - const CCTK_REAL16 *restrict const input_array_ptr_real16 + const CCTK_REAL16 *const input_array_ptr_real16 = ((const CCTK_REAL16 *) input_array_ptr) + input_posn; LOAD_DATA_REAL16(input_array_ptr_real16, STRIDE_IJK, @@ -1457,13 +1447,21 @@ case CCTK_VARIABLE_REAL: } #endif +case CCTK_VARIABLE_COMPLEX: + { + const CCTK_REAL (*const input_array_ptr_complex)[COMPLEX_N_PARTS] + = ((const CCTK_REAL (*)[COMPLEX_N_PARTS]) input_array_ptr) + + input_posn; + LOAD_DATA_COMPLEX(input_array_ptr_complex, + STRIDE_IJK, part, + &data); + break; + } + #ifdef HAVE_CCTK_COMPLEX8 case CCTK_VARIABLE_COMPLEX8: -#ifdef CCTK_COMPLEX_PRECISION_8 -case CCTK_VARIABLE_COMPLEX: -#endif { - const CCTK_REAL4 (*restrict const input_array_ptr_complex8)[COMPLEX_N_PARTS] + const CCTK_REAL4 (*const input_array_ptr_complex8)[COMPLEX_N_PARTS] = ((const CCTK_REAL4 (*)[COMPLEX_N_PARTS]) input_array_ptr) + input_posn; LOAD_DATA_COMPLEX8(input_array_ptr_complex8, @@ -1475,11 +1473,8 @@ case CCTK_VARIABLE_COMPLEX: #ifdef HAVE_CCTK_COMPLEX16 case CCTK_VARIABLE_COMPLEX16: -#ifdef CCTK_COMPLEX_PRECISION_16 -case CCTK_VARIABLE_COMPLEX: -#endif { - const CCTK_REAL8 (*restrict const input_array_ptr_complex16)[COMPLEX_N_PARTS] + const CCTK_REAL8 (*const input_array_ptr_complex16)[COMPLEX_N_PARTS] = ((const CCTK_REAL8 (*)[COMPLEX_N_PARTS]) input_array_ptr) + input_posn; LOAD_DATA_COMPLEX16(input_array_ptr_complex16, @@ -1492,13 +1487,10 @@ case CCTK_VARIABLE_COMPLEX: #ifdef HAVE_CCTK_COMPLEX32 case CCTK_VARIABLE_COMPLEX32: -#ifdef CCTK_COMPLEX_PRECISION_32 -case CCTK_VARIABLE_COMPLEX: -#endif { /* FIXME: maybe we should warn (once per cactus run) that we may be */ /* doing arithmetic in lower precision than the data type? */ - const CCTK_REAL16 (*restrict const input_array_ptr_complex32)[COMPLEX_N_PARTS] + const CCTK_REAL16 (*const input_array_ptr_complex32)[COMPLEX_N_PARTS] = ((const CCTK_REAL16 (*)[COMPLEX_N_PARTS]) input_array_ptr) + input_posn; LOAD_DATA_COMPLEX32(input_array_ptr_complex32, @@ -1516,9 +1508,7 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, , (int) input_array_type_codes[in], in); /*NOTREACHED*/ -#pragma omp critical -return_status = UTIL_ERROR_BAD_INPUT; -goto next_iteration; +return UTIL_ERROR_BAD_INPUT; /*** ERROR RETURN ***/ /* end of switch (input_array_type_codes[in]) */ } } @@ -1618,9 +1608,8 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, , (int) operation_codes[out], out); -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; + /*** ERROR RETURN ***/ /* end of switch (operation_codes[out]) */ } @@ -1640,13 +1629,29 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, switch (output_array_type_codes[out]) { +case CCTK_VARIABLE_REAL: + { + CCTK_REAL *const output_array_ptr_real + = (CCTK_REAL *) output_arrays[out]; + if (debug >= 10) + then { + if ((pt & (pt-1)) == 0) /* pt is 0 or a power of 2 */ + then { + printf(" result addr is %p\n", + (void *) &output_array_ptr_real[pt]); + printf(" previous value there was %g\n", + output_array_ptr_real[pt]); + fflush(stdout); + } + } + output_array_ptr_real[pt] = (CCTK_REAL) result; + break; + } + #ifdef HAVE_CCTK_REAL4 case CCTK_VARIABLE_REAL4: -#ifdef CCTK_REAL_PRECISION_4 -case CCTK_VARIABLE_REAL: -#endif { - CCTK_REAL4 *restrict const output_array_ptr_real4 + CCTK_REAL4 *const output_array_ptr_real4 = (CCTK_REAL4 *) output_arrays[out]; output_array_ptr_real4[pt] = (CCTK_REAL4) result; break; @@ -1655,48 +1660,36 @@ case CCTK_VARIABLE_REAL: #ifdef HAVE_CCTK_REAL8 case CCTK_VARIABLE_REAL8: -#ifdef CCTK_REAL_PRECISION_8 -case CCTK_VARIABLE_REAL: -#endif { - CCTK_REAL8 *restrict const output_array_ptr_real8 + CCTK_REAL8 *const output_array_ptr_real8 = (CCTK_REAL8 *) output_arrays[out]; output_array_ptr_real8[pt] = (CCTK_REAL8) result; - if (debug >= 10) - then { - if ((pt & (pt-1)) == 0) /* pt is 0 or a power of 2 */ - then { - printf(" result addr is %p\n", - (void *) &output_array_ptr_real8[pt]); - printf(" previous value there was %.17g\n", - (double) output_array_ptr_real8[pt]); - fflush(stdout); - } - } break; } #endif #ifdef HAVE_CCTK_REAL16 case CCTK_VARIABLE_REAL16: -#ifdef CCTK_REAL_PRECISION_16 -case CCTK_VARIABLE_REAL: -#endif { - CCTK_REAL16 *restrict const output_array_ptr_real16 + CCTK_REAL16 *const output_array_ptr_real16 = (CCTK_REAL16 *) output_arrays[out]; output_array_ptr_real16[pt] = (CCTK_REAL16) result; break; } #endif +case CCTK_VARIABLE_COMPLEX: + { + CCTK_REAL (*const output_array_ptr_complex)[COMPLEX_N_PARTS] + = (CCTK_REAL (*)[COMPLEX_N_PARTS]) output_arrays[out]; + output_array_ptr_complex[pt][part] = (CCTK_REAL) result; + break; + } + #ifdef HAVE_CCTK_COMPLEX8 case CCTK_VARIABLE_COMPLEX8: -#ifdef CCTK_COMPLEX_PRECISION_8 -case CCTK_VARIABLE_COMPLEX: -#endif { - CCTK_REAL4 (*restrict const output_array_ptr_complex8)[COMPLEX_N_PARTS] + CCTK_REAL4 (*const output_array_ptr_complex8)[COMPLEX_N_PARTS] = (CCTK_REAL4 (*)[COMPLEX_N_PARTS]) output_arrays[out]; output_array_ptr_complex8[pt][part] = (CCTK_REAL4) result; break; @@ -1705,11 +1698,8 @@ case CCTK_VARIABLE_COMPLEX: #ifdef HAVE_CCTK_COMPLEX16 case CCTK_VARIABLE_COMPLEX16: -#ifdef CCTK_COMPLEX_PRECISION_16 -case CCTK_VARIABLE_COMPLEX: -#endif { - CCTK_REAL8 (*restrict const output_array_ptr_complex16)[COMPLEX_N_PARTS] + CCTK_REAL8 (*const output_array_ptr_complex16)[COMPLEX_N_PARTS] = (CCTK_REAL8 (*)[COMPLEX_N_PARTS]) output_arrays[out]; output_array_ptr_complex16[pt][part] = (CCTK_REAL8) result; break; @@ -1718,11 +1708,8 @@ case CCTK_VARIABLE_COMPLEX: #ifdef HAVE_CCTK_COMPLEX32 case CCTK_VARIABLE_COMPLEX32: -#ifdef CCTK_COMPLEX_PRECISION_32 -case CCTK_VARIABLE_COMPLEX: -#endif { - CCTK_REAL16 (*restrict const output_array_ptr_complex32)[COMPLEX_N_PARTS] + CCTK_REAL16 (*const output_array_ptr_complex32)[COMPLEX_N_PARTS] = (CCTK_REAL16 (*)[COMPLEX_N_PARTS]) output_arrays[out]; output_array_ptr_complex32[pt][part] = (CCTK_REAL16) result; break; @@ -1737,9 +1724,7 @@ default: , (int) output_array_type_codes[out], out); -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; /*** ERROR RETURN ***/ /* end of switch (output type code) */ } @@ -1753,7 +1738,7 @@ default: if ( (Jacobian_info != NULL) && (Jacobian_info->Jacobian_pointer[out] != NULL)) then { - CCTK_REAL *restrict const Jacobian_ptr + CCTK_REAL *const Jacobian_ptr = Jacobian_info->Jacobian_pointer[out] + Jacobian_info->Jacobian_offset[out] + Jacobian_info->Jacobian_interp_point_stride*pt @@ -1843,9 +1828,8 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, , (int) operation_codes[out], out); -#pragma omp critical - return_status = UTIL_ERROR_BAD_INPUT; - goto next_iteration; + return UTIL_ERROR_BAD_INPUT; + /*** ERROR RETURN ***/ /* end of switch(operation_codes[out])*/ } /* end of Jacobian-query code */ @@ -1861,7 +1845,6 @@ CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING, } } - next_iteration:; /* end of for (pt = ...) loop */ } @@ -1926,7 +1909,7 @@ return return_status; #ifdef HAVE_OP_I static - void compute_coeffs_I(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_I) + void compute_coeffs_I(FP_XYZ, struct COEFFS_STRUCT *coeffs_I) { #include COEFFS_I_COMPUTE_FILE_NAME } @@ -1934,7 +1917,7 @@ return return_status; #ifdef HAVE_OP_DX static - void compute_coeffs_dx(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dx) + void compute_coeffs_dx(FP_XYZ, struct COEFFS_STRUCT *coeffs_dx) { #include COEFFS_DX_COMPUTE_FILE_NAME } @@ -1942,7 +1925,7 @@ return return_status; #ifdef HAVE_OP_DY static - void compute_coeffs_dy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dy) + void compute_coeffs_dy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dy) { #include COEFFS_DY_COMPUTE_FILE_NAME } @@ -1950,7 +1933,7 @@ return return_status; #ifdef HAVE_OP_DZ static - void compute_coeffs_dz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dz) + void compute_coeffs_dz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dz) { #include COEFFS_DZ_COMPUTE_FILE_NAME } @@ -1958,7 +1941,7 @@ return return_status; #ifdef HAVE_OP_DXX static - void compute_coeffs_dxx(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxx) + void compute_coeffs_dxx(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxx) { #include COEFFS_DXX_COMPUTE_FILE_NAME } @@ -1966,7 +1949,7 @@ return return_status; #ifdef HAVE_OP_DXY static - void compute_coeffs_dxy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxy) + void compute_coeffs_dxy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxy) { #include COEFFS_DXY_COMPUTE_FILE_NAME } @@ -1974,7 +1957,7 @@ return return_status; #ifdef HAVE_OP_DXZ static - void compute_coeffs_dxz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dxz) + void compute_coeffs_dxz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dxz) { #include COEFFS_DXZ_COMPUTE_FILE_NAME } @@ -1982,7 +1965,7 @@ return return_status; #ifdef HAVE_OP_DYY static - void compute_coeffs_dyy(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dyy) + void compute_coeffs_dyy(FP_XYZ, struct COEFFS_STRUCT *coeffs_dyy) { #include COEFFS_DYY_COMPUTE_FILE_NAME } @@ -1990,7 +1973,7 @@ return return_status; #ifdef HAVE_OP_DYZ static - void compute_coeffs_dyz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dyz) + void compute_coeffs_dyz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dyz) { #include COEFFS_DYZ_COMPUTE_FILE_NAME } @@ -1998,7 +1981,7 @@ return return_status; #ifdef HAVE_OP_DZZ static - void compute_coeffs_dzz(FP_XYZ, struct COEFFS_STRUCT *restrict coeffs_dzz) + void compute_coeffs_dzz(FP_XYZ, struct COEFFS_STRUCT *coeffs_dzz) { #include COEFFS_DZZ_COMPUTE_FILE_NAME } |