diff options
Diffstat (limited to 'src/ReductionNorm2.c')
-rw-r--r-- | src/ReductionNorm2.c | 224 |
1 files changed, 112 insertions, 112 deletions
diff --git a/src/ReductionNorm2.c b/src/ReductionNorm2.c index 0e25b00..59200fa 100644 --- a/src/ReductionNorm2.c +++ b/src/ReductionNorm2.c @@ -158,16 +158,16 @@ static int ReductionL2 (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); @@ -306,24 +306,24 @@ static int ReductionL2 (int N_dims, int operator_handle, num_points = LocalReduce_Norm2_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: - num_points = LocalReduce_Norm2_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: - num_points = LocalReduce_Norm2_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: - num_points = LocalReduce_Norm2_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: - num_points = LocalReduce_Norm2_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: */ + /* num_points = LocalReduce_Norm2_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: */ + /* num_points = LocalReduce_Norm2_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: */ + /* num_points = LocalReduce_Norm2_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: */ + /* num_points = LocalReduce_Norm2_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 */ } } @@ -406,48 +406,48 @@ static int ReductionL2 (int N_dims, int operator_handle, *( (CCTK_REAL16 *) output_numbers[i]) = sqrt(*( (CCTK_REAL16 *) output_numbers[i])); break; #endif - case CCTK_VARIABLE_COMPLEX: - (*( (CCTK_COMPLEX *) output_numbers[i])).Re = (*( (CCTK_COMPLEX *) output_numbers[i])).Re / weight_sum; - (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / weight_sum; - 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/2.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 / weight_sum; - (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / weight_sum; - 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/2.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 / weight_sum; - (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / weight_sum; - 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/2.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 / weight_sum; - (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / weight_sum; - 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/2.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 / weight_sum; */ + /* (*( (CCTK_COMPLEX *) output_numbers[i])).Im = (*( (CCTK_COMPLEX *) output_numbers[i])).Im / weight_sum; */ + /* 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/2.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 / weight_sum; */ + /* (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX8 *) output_numbers[i])).Im / weight_sum; */ + /* 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/2.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 / weight_sum; */ + /* (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX16 *) output_numbers[i])).Im / weight_sum; */ + /* 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/2.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 / weight_sum; */ + /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im / weight_sum; */ + /* 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/2.0); */ + /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re; */ + /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im; */ + /* break; */ + /* #endif */ } } } @@ -518,48 +518,48 @@ static int ReductionL2 (int N_dims, int operator_handle, *( (CCTK_REAL16 *) output_numbers[i]) = sqrt(*( (CCTK_REAL16 *) output_numbers[i])); 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/2.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/2.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/2.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/2.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/2.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/2.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/2.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/2.0); */ + /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Re = cmplx_tmp32.Re; */ + /* (*( (CCTK_COMPLEX32 *) output_numbers[i])).Im = cmplx_tmp32.Im; */ + /* break; */ + /* #endif */ } } } |