int
prod_to_all(int me, int npes)
{
int i, pass=0;
int float_rounding_err=0;
int double_rounding_err=0;
int ldouble_rounding_err=0;
memset(ok,0,sizeof(ok));
for (i = 0; i < N; i++) {
src0[i] = src1[i] = src2[i] = src3[i] = src4[i] = src5[i] = src6[i]
= me + 1;
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst3[i] = -9;
dst4[i] = -9;
dst5[i] = -9;
dst6[i] = -9;
}
expected_result0 = expected_result1 = expected_result2 =
expected_result6 = 1;
expected_result3 = expected_result4 = expected_result5 = 1.0;
for(i=1; i <= npes; i++) {
expected_result0 *= i;
expected_result1 *= i;
expected_result2 *= i;
expected_result3 *= (float)i;
expected_result4 *= (double)i;
if ((double)expected_result3 != expected_result4) {
if (!float_rounding_err && Verbose > 2 && me == 0)
printf("float_err @ npes %d\n",i);
float_rounding_err = 1;
}
expected_result5 *= (long double)i;
if ((long double)expected_result4 != expected_result5) {
if (!double_rounding_err && Verbose > 2 && me == 0)
printf("double_err @ npes %d\n",i);
ldouble_rounding_err = double_rounding_err = 1;
}
expected_result6 *= i;
}
shmem_barrier_all();
shmem_short_prod_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_prod_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_prod_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_float_prod_to_all(dst3, src3, N, 0, 0, npes, pWrk3, pSync1);
shmem_double_prod_to_all(dst4, src4, N, 0, 0, npes, pWrk4, pSync);
shmem_longdouble_prod_to_all(dst5, src5, N, 0, 0, npes, pWrk5, pSync1);
shmem_longlong_prod_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync);
if(me == 0) {
for (i = 0; i < N; i++) {
if(dst0[i] != expected_result0) ok[0] = 1;
if(dst1[i] != expected_result1) ok[1] = 1;
if(dst2[i] != expected_result2) ok[2] = 1;
/* check for overflow */
if(!float_rounding_err && dst3[i] != expected_result3) { ok[3] = 1;
printf("dst3[%d]: %f, expected val: %f\n",i, dst3[i], expected_result3);
}
if(!double_rounding_err && dst4[i] != expected_result4) {ok[4] = 1;
printf("dst4[%d]: %f, expected val: %f\n",i, dst4[i], expected_result4);
}
if(!ldouble_rounding_err && dst5[i] != expected_result5) {ok[5] = 1;
printf("dst5[%d]: %Lf, expected val: %Lf T4 %f\n",i, dst5[i], expected_result5,dst4[i]);
}
if(dst6[i] != expected_result6) ok[6] = 1;
}
if(ok[0]==1)
printf("Reduction operation shmem_short_prod_to_all: Failed\n");
else {
Vprintf("Reduction operation shmem_short_prod_to_all: Passed\n");
pass++;
}
if(ok[1]==1)
printf("Reduction operation shmem_int_prod_to_all: Failed\n");
else {
Vprintf("Reduction operation shmem_int_prod_to_all: Passed\n");
pass++;
}
if(ok[2]==1)
printf("Reduction operation shmem_long_prod_to_all: Failed\n");
else {
Vprintf("Reduction operation shmem_long_prod_to_all: Passed\n");
pass++;
}
if(ok[3]==1)
printf("Reduction operation shmem_float_prod_to_all: Failed\n");
else {
if (float_rounding_err) {
Vprintf("Reduction operation shmem_float_prod_to_all: skipped due to float rounding error\n");
}
else {
Vprintf("Reduction operation shmem_float_prod_to_all: Passed\n");
}
pass++;
}
if(ok[4]==1)
printf("Reduction operation shmem_double_prod_to_all: Failed\n");
else {
if (double_rounding_err) {
Vprintf("Reduction operation shmem_double_prod_to_all: skipped due to double rounding error\n");
}
else {
Vprintf("Reduction operation shmem_double_prod_to_all: Passed\n");
}
pass++;
}
if(ok[5]==1)
printf("Reduction operation shmem_longdouble_prod_to_all: Failed\n");
else {
if (double_rounding_err) {
Vprintf("Reduction operation shmem_longdouble_prod_to_all: skipped due to long double rounding error\n");
}
else {
Vprintf("Reduction operation shmem_longdouble_prod_to_all: Passed\n");
}
pass++;
}
if(ok[6]==1)
printf("Reduction operation shmem_longlong_prod_to_all: Failed\n");
else {
Vprintf("Reduction operation shmem_longlong_prod_to_all: Passed\n");
pass++;
}
Vprintf("\n");
}
if (Serialize) shmem_barrier_all();
return (pass == 7 ? 1 : 0);
}
int
main()
{
int i,j;
int me, npes;
int success0, success1, success2, success3, success4, success5, success6;
success0 = success1 = success2 = success3 = success4 = success5 = success6 = 0;
start_pes(0);
me = _my_pe();
npes = _num_pes();
for (i = 0; i < _SHMEM_REDUCE_SYNC_SIZE; i += 1) {
pSync[i] = _SHMEM_SYNC_VALUE;
pSync1[i] = _SHMEM_SYNC_VALUE;
}
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src3[i] = src4[i] = src5[i] = src6[i] = me + i;
}
/*Test MAX: shmem_double_max_to_all, shmem_float_max_to_all, shmem_int_max_to_all, shmem_long_max_to_all, shmem_longdouble_max_to_all, shmem_longlong_max_to_all, shmem_short_max_to_all */
shmem_barrier_all();
shmem_short_max_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_max_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_max_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_float_max_to_all(dst3, src3, N, 0, 0, npes, pWrk3, pSync1);
shmem_double_max_to_all(dst4, src4, N, 0, 0, npes, pWrk4, pSync);
shmem_longdouble_max_to_all(dst5, src5, N, 0, 0, npes, pWrk5, pSync1);
shmem_longlong_max_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync);
if(me == 0){
for (i = 0,j=-1; i < N; i++,j++) {
if(dst0[i] != npes+j)
success0 =1;
if(dst1[i] != npes+j)
success1 =1;
if(dst2[i] != npes+j)
success2 =1;
if(dst3[i] != npes+j)
success3 =1;
if(dst4[i] != npes+j)
success4 =1;
if(dst5[i] != npes+j)
success5 =1;
if(dst6[i] != npes+j)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_max_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_max_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_max_to_all: Passed\n");
}
if(success3==1){
printf("Reduction operation shmem_float_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_float_max_to_all: Passed\n");
}
if(success4==1){
printf("Reduction operation shmem_double_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_double_max_to_all: Passed\n");
}
if(success5==1){
printf("Reduction operation shmem_longdouble_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longdouble_max_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_max_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_max_to_all: Passed\n");
}
}
/*Test MIN: shmem_double_min_to_all, shmem_float_min_to_all, shmem_int_min_to_all, shmem_long_min_to_all, shmem_longdouble_min_to_all, shmem_longlong_min_to_all, shmem_short_min_to_all*/
success0 = success1 = success2 = success3 = success4 = success5 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src3[i] = src4[i] = src5[i] = src6[i] = me + i;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst3[i] = -9;
dst4[i] = -9;
dst5[i] = -9;
dst6[i] = -9;
}
shmem_barrier_all();
shmem_short_min_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_min_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_min_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_float_min_to_all(dst3, src3, N, 0, 0, npes, pWrk3, pSync1);
shmem_double_min_to_all(dst4, src4, N, 0, 0, npes, pWrk4, pSync);
shmem_longdouble_min_to_all(dst5, src5, N, 0, 0, npes, pWrk5, pSync1);
shmem_longlong_min_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync);
if(me == 0){
for (i = 0; i < N; i++) {
if(dst0[i] != i)
success0 =1;
if(dst1[i] != i)
success1 =1;
if(dst2[i] != i)
success2 =1;
if(dst3[i] != i)
success3 =1;
if(dst4[i] != i)
success4 =1;
if(dst5[i] != i)
success5 =1;
if(dst6[i] != i)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_min_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_min_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_min_to_all: Passed\n");
}
if(success3==1){
printf("Reduction operation shmem_float_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_float_min_to_all: Passed\n");
}
if(success4==1){
printf("Reduction operation shmem_double_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_double_min_to_all: Passed\n");
}
if(success5==1){
printf("Reduction operation shmem_longdouble_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longdouble_min_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_min_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_min_to_all: Passed\n");
}
}
/*Test SUM: shmem_double_sum_to_all, shmem_float_sum_to_all, shmem_int_sum_to_all, shmem_long_sum_to_all, shmem_longdouble_sum_to_all, shmem_longlong_sum_to_all, shmem_short_sum_to_all*/
success0 = success1 = success2 = success3 = success4 = success5 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src3[i] = src4[i] = src5[i] = src6[i] = me;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst3[i] = -9;
dst4[i] = -9;
dst5[i] = -9;
dst6[i] = -9;
}
shmem_barrier_all();
shmem_short_sum_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_sum_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_sum_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_float_sum_to_all(dst3, src3, N, 0, 0, npes, pWrk3, pSync1);
shmem_double_sum_to_all(dst4, src4, N, 0, 0, npes, pWrk4, pSync);
shmem_longdouble_sum_to_all(dst5, src5, N, 0, 0, npes, pWrk5, pSync1);
shmem_longlong_sum_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync);
if(me == 0){
for (i = 0; i < N; i++) {
if(dst0[i] != (npes * (npes-1)/2))
success0 =1;
if(dst1[i] != (npes * (npes-1)/2))
success1 =1;
if(dst2[i] != (npes * (npes-1)/2))
success2 =1;
if(dst3[i] != (npes * (npes-1)/2))
success3 =1;
if(dst4[i] != (npes * (npes-1)/2))
success4 =1;
if(dst5[i] != (npes * (npes-1)/2))
success5 =1;
if(dst6[i] != (npes * (npes-1)/2))
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_sum_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_sum_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_sum_to_all: Passed\n");
}
if(success3==1){
printf("Reduction operation shmem_float_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_float_sum_to_all: Passed\n");
}
if(success4==1){
printf("Reduction operation shmem_double_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_double_sum_to_all: Passed\n");
}
if(success5==1){
printf("Reduction operation shmem_longdouble_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longdouble_sum_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_sum_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_sum_to_all: Passed\n");
}
}
/*Test AND: shmem_int_and_to_all, shmem_long_and_to_all, shmem_longlong_and_to_all, shmem_short_and_to_all,*/
success0 = success1 = success2 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src6[i] = me;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst6[i] = -9;
}
shmem_barrier_all();
shmem_short_and_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_and_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_and_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_longlong_and_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync1);
if(me==0){
for (i = 0; i < N; i++) {
if(dst0[i] != 0)
success0 =1;
if(dst1[i] != 0)
success1 =1;
if(dst2[i] != 0)
success2 =1;
if(dst6[i] != 0)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_and_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_and_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_and_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_and_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_and_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_and_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_and_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_and_to_all: Passed\n");
}
}
/*Test PROD: shmem_double_prod_to_all, shmem_float_prod_to_all, shmem_int_prod_to_all, shmem_long_prod_to_all, shmem_longdouble_prod_to_all, shmem_longlong_prod_to_all, shmem_short_prod_to_all, */
success0 = success1 = success2 = success3 = success4 = success5 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src3[i] = src4[i] = src5[i] = src6[i] = me + 1;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst3[i] = -9;
dst4[i] = -9;
dst5[i] = -9;
dst6[i] = -9;
}
expected_result0 = expected_result1 = expected_result2 = expected_result3 = expected_result4 = expected_result5 = expected_result6 =1;
for(i=1;i<=npes;i++){
expected_result0 = expected_result0 * i;
expected_result1 = expected_result1 * i;
expected_result2 = expected_result2 * i;
expected_result3 = expected_result3 * i;
expected_result4 = expected_result4 * i;
expected_result5 = expected_result5 * i;
expected_result6 = expected_result6 * i;
}
shmem_barrier_all();
shmem_short_prod_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_prod_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_prod_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_float_prod_to_all(dst3, src3, N, 0, 0, npes, pWrk3, pSync1);
shmem_double_prod_to_all(dst4, src4, N, 0, 0, npes, pWrk4, pSync);
shmem_longdouble_prod_to_all(dst5, src5, N, 0, 0, npes, pWrk5, pSync1);
shmem_longlong_prod_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync);
if(me == 0){
for (i = 0; i < N; i++) {
/*printf("dst2[%d]: %ld, expected val: %ld\n",i, dst2[i], (long)expected_result2);*/
if(dst0[i] != expected_result0)
success0 =1;
if(dst1[i] != expected_result1)
success1 =1;
if(dst2[i] != expected_result2)
success2 =1;
if(dst3[i] != expected_result3)
success3 =1;
if(dst4[i] != expected_result4)
success4 =1;
if(dst5[i] != expected_result5)
success5 =1;
if(dst6[i] != expected_result6)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_prod_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_prod_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_prod_to_all: Passed\n");
}
if(success3==1){
printf("Reduction operation shmem_float_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_float_prod_to_all: Passed\n");
}
if(success4==1){
printf("Reduction operation shmem_double_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_double_prod_to_all: Passed\n");
}
if(success5==1){
printf("Reduction operation shmem_longdouble_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longdouble_prod_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_prod_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_prod_to_all: Passed\n");
}
}
/*Test OR: shmem_int_or_to_all, shmem_long_or_to_all, shmem_longlong_or_to_all, shmem_short_or_to_all,*/
success0 = success1 = success2 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src6[i] = (me + 1)%4;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst6[i] = -9;
}
shmem_barrier_all();
shmem_short_or_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_or_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_or_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_longlong_or_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync1);
if(me==0){
for (i = 0; i < N; i++) {
if(dst0[i] != 3)
success0 =1;
if(dst1[i] != 3)
success1 =1;
if(dst2[i] != 3)
success2 =1;
if(dst6[i] != 3)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_or_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_or_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_or_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_or_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_or_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_or_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_or_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_or_to_all: Passed\n");
}
}
/*Test XOR: shmem_int_xor_to_all, shmem_long_xor_to_all, shmem_longlong_xor_to_all, shmem_short_xor_to_all*/
success0 = success1 = success2 = success6 = 0;
for (i = 0; i < N; i += 1) {
src0[i] = src1[i] = src2[i] = src6[i] = me%2;
}
for (i = 0; i < N; i += 1) {
dst0[i] = -9;
dst1[i] = -9;
dst2[i] = -9;
dst6[i] = -9;
}
int expected_result = ((int)(npes/2) % 2);
shmem_barrier_all();
shmem_short_xor_to_all(dst0, src0, N, 0, 0, npes, pWrk0, pSync);
shmem_int_xor_to_all(dst1, src1, N, 0, 0, npes, pWrk1, pSync1);
shmem_long_xor_to_all(dst2, src2, N, 0, 0, npes, pWrk2, pSync);
shmem_longlong_xor_to_all(dst6, src6, N, 0, 0, npes, pWrk6, pSync1);
if(me==0){
for (i = 0; i < N; i++) {
if(dst0[i] != expected_result)
success0 =1;
if(dst1[i] != expected_result)
success1 =1;
if(dst2[i] != expected_result)
success2 =1;
if(dst6[i] != expected_result)
success6 =1;
}
if(success0==1){
printf("Reduction operation shmem_short_xor_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_short_xor_to_all: Passed\n");
}
if(success1==1){
printf("Reduction operation shmem_int_xor_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_int_xor_to_all: Passed\n");
}
if(success2==1){
printf("Reduction operation shmem_long_xor_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_long_xor_to_all: Passed\n");
}
if(success6==1){
printf("Reduction operation shmem_longlong_xor_to_all: Failed\n");
}
else{
printf("Reduction operation shmem_longlong_xor_to_all: Passed\n");
}
}
return 0;
}