aboutsummaryrefslogtreecommitdiff
path: root/src/ReductionNorm3.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/ReductionNorm3.c')
-rw-r--r--src/ReductionNorm3.c224
1 files changed, 112 insertions, 112 deletions
diff --git a/src/ReductionNorm3.c b/src/ReductionNorm3.c
index c45b8b6..c7b55be 100644
--- a/src/ReductionNorm3.c
+++ b/src/ReductionNorm3.c
@@ -158,16 +158,16 @@ static int ReductionL3 (int N_dims, int operator_handle,
CCTK_REAL weight_sum;
/* temporary complex variables*/
- CCTK_COMPLEX cmplx_tmp;
- #ifdef HAVE_CCTK_COMPLEX8
- CCTK_COMPLEX8 cmplx_tmp8;
- #endif
- #ifdef HAVE_CCTK_COMPLEX16
- CCTK_COMPLEX16 cmplx_tmp16;
- #endif
- #ifdef HAVE_CCTK_COMPLEX32
- CCTK_COMPLEX32 cmplx_tmp32;
- #endif
+ /* CCTK_COMPLEX cmplx_tmp; */
+ /* #ifdef HAVE_CCTK_COMPLEX8 */
+ /* CCTK_COMPLEX8 cmplx_tmp8; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX16 */
+ /* CCTK_COMPLEX16 cmplx_tmp16; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX32 */
+ /* CCTK_COMPLEX32 cmplx_tmp32; */
+ /* #endif */
/* prevent warnings for unused vars */
(void)(operator_handle + 0);
@@ -303,24 +303,24 @@ static int ReductionL3 (int N_dims, int operator_handle,
LocalReduce_Norm3_REAL16(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle);
break;
#endif
- case CCTK_VARIABLE_COMPLEX:
- LocalReduce_Norm3_COMPLEX(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle);
- break;
- #ifdef HAVE_CCTK_COMPLEX8
- case CCTK_VARIABLE_COMPLEX8:
- LocalReduce_Norm3_COMPLEX8(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle);
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX16
- case CCTK_VARIABLE_COMPLEX16:
- LocalReduce_Norm3_COMPLEX16(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle);
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX32
- case CCTK_VARIABLE_COMPLEX32:
- LocalReduce_Norm3_COMPLEX32(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle);
- break;
- #endif
+ /* case CCTK_VARIABLE_COMPLEX: */
+ /* LocalReduce_Norm3_COMPLEX(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle); */
+ /* break; */
+ /* #ifdef HAVE_CCTK_COMPLEX8 */
+ /* case CCTK_VARIABLE_COMPLEX8: */
+ /* LocalReduce_Norm3_COMPLEX8(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle); */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX16 */
+ /* case CCTK_VARIABLE_COMPLEX16: */
+ /* LocalReduce_Norm3_COMPLEX16(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle); */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX32 */
+ /* case CCTK_VARIABLE_COMPLEX32: */
+ /* LocalReduce_Norm3_COMPLEX32(i, weight_on, weight, input_array_offsets, indices, max_iter, actual_indices, input_array_strides, input_array_min_subscripts, input_array_dims, num_points, actual_iters_per_dim, iters_per_dim, N_dims, input_arrays, output_number_type_codes, output_numbers, param_table_handle); */
+ /* break; */
+ /* #endif */
}
}
@@ -403,48 +403,48 @@ static int ReductionL3 (int N_dims, int operator_handle,
*( (CCTK_REAL16 *) output_numbers[i]) = pow(*( (CCTK_REAL16 *) output_numbers[i]), 1.0/3.0);
break;
#endif
- case CCTK_VARIABLE_COMPLEX:
- (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / num_points;
- cmplx_tmp.Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re;
- cmplx_tmp.Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX *) output_numbers[i])).Re = cmplx_tmp.Re;
- (*( (CCTK_COMPLEX *) output_numbers[i])).Im = cmplx_tmp.Im;
- break;
- #ifdef HAVE_CCTK_COMPLEX8
- case CCTK_VARIABLE_COMPLEX8:
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / num_points;
- cmplx_tmp8.Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re;
- cmplx_tmp8.Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = cmplx_tmp8.Re;
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = cmplx_tmp8.Im;
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX16
- case CCTK_VARIABLE_COMPLEX16:
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / num_points;
- cmplx_tmp16.Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re;
- cmplx_tmp16.Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = cmplx_tmp16.Re;
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = cmplx_tmp16.Im;
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX32
- case CCTK_VARIABLE_COMPLEX32:
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / num_points;
- cmplx_tmp32.Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re;
- cmplx_tmp32.Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re;
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im;
- break;
- #endif
+ /* case CCTK_VARIABLE_COMPLEX: */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp.Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re; */
+ /* cmplx_tmp.Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Re = cmplx_tmp.Re; */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Im = cmplx_tmp.Im; */
+ /* break; */
+ /* #ifdef HAVE_CCTK_COMPLEX8 */
+ /* case CCTK_VARIABLE_COMPLEX8: */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp8.Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re; */
+ /* cmplx_tmp8.Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = cmplx_tmp8.Re; */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = cmplx_tmp8.Im; */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX16 */
+ /* case CCTK_VARIABLE_COMPLEX16: */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp16.Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re; */
+ /* cmplx_tmp16.Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = cmplx_tmp16.Re; */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = cmplx_tmp16.Im; */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX32 */
+ /* case CCTK_VARIABLE_COMPLEX32: */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp32.Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re; */
+ /* cmplx_tmp32.Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re; */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im; */
+ /* break; */
+ /* #endif */
}
}
}
@@ -515,48 +515,48 @@ static int ReductionL3 (int N_dims, int operator_handle,
*( (CCTK_REAL16 *) output_numbers[i]) = pow(*( (CCTK_REAL16 *) output_numbers[i]), 1.0/3.0);
break;
#endif
- case CCTK_VARIABLE_COMPLEX:
- (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / num_points;
- cmplx_tmp.Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re;
- cmplx_tmp.Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX *) output_numbers[i])).Re = cmplx_tmp.Re;
- (*( (CCTK_COMPLEX *) output_numbers[i])).Im = cmplx_tmp.Im;
- break;
- #ifdef HAVE_CCTK_COMPLEX8
- case CCTK_VARIABLE_COMPLEX8:
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / num_points;
- cmplx_tmp8.Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re;
- cmplx_tmp8.Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = cmplx_tmp8.Re;
- (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = cmplx_tmp8.Im;
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX16
- case CCTK_VARIABLE_COMPLEX16:
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / num_points;
- cmplx_tmp16.Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re;
- cmplx_tmp16.Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = cmplx_tmp16.Re;
- (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = cmplx_tmp16.Im;
- break;
- #endif
- #ifdef HAVE_CCTK_COMPLEX32
- case CCTK_VARIABLE_COMPLEX32:
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re / num_points;
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / num_points;
- cmplx_tmp32.Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re;
- cmplx_tmp32.Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im;
- cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0);
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re;
- (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im;
- break;
- #endif
+ /* case CCTK_VARIABLE_COMPLEX: */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp.Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re; */
+ /* cmplx_tmp.Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Re = cmplx_tmp.Re; */
+ /* (*( (CCTK_COMPLEX *) output_numbers[i])).Im = cmplx_tmp.Im; */
+ /* break; */
+ /* #ifdef HAVE_CCTK_COMPLEX8 */
+ /* case CCTK_VARIABLE_COMPLEX8: */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp8.Re = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re; */
+ /* cmplx_tmp8.Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Re = cmplx_tmp8.Re; */
+ /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = cmplx_tmp8.Im; */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX16 */
+ /* case CCTK_VARIABLE_COMPLEX16: */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp16.Re = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re; */
+ /* cmplx_tmp16.Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Re = cmplx_tmp16.Re; */
+ /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = cmplx_tmp16.Im; */
+ /* break; */
+ /* #endif */
+ /* #ifdef HAVE_CCTK_COMPLEX32 */
+ /* case CCTK_VARIABLE_COMPLEX32: */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re / num_points; */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / num_points; */
+ /* cmplx_tmp32.Re = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re; */
+ /* cmplx_tmp32.Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im; */
+ /* cmplx_tmp = CCTK_CmplxPow(cmplx_tmp, 1.0/3.0); */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re; */
+ /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im; */
+ /* break; */
+ /* #endif */
}
}
}
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-18 21:35:37 +0000