/* StarPU --- Runtime system for heterogeneous multicore architectures. * * Copyright (C) 2013, 2014, 2015, 2016, 2017 CNRS * * StarPU is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by * the Free Software Foundation; either version 2.1 of the License, or (at * your option) any later version. * * StarPU is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. * * See the GNU Lesser General Public License in COPYING.LGPL for more details. */ #include #include #include "helper.h" typedef void (*check_func)(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error); void check_void(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error) { FPRINTF_MPI(stderr, "Success with void value\n"); } void check_variable(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error) { float *v_s, *v_r; STARPU_ASSERT(starpu_variable_get_elemsize(handle_s) == starpu_variable_get_elemsize(handle_r)); v_s = (float *)starpu_variable_get_local_ptr(handle_s); v_r = (float *)starpu_variable_get_local_ptr(handle_r); if (*v_s == *v_r) { FPRINTF_MPI(stderr, "Success with variable value: %f == %f\n", *v_s, *v_r); } else { *error = 1; FPRINTF_MPI(stderr, "Error with variable value: %f != %f\n", *v_s, *v_r); } } void check_vector(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error) { int i; int nx; int *v_r, *v_s; STARPU_ASSERT(starpu_vector_get_elemsize(handle_s) == starpu_vector_get_elemsize(handle_r)); STARPU_ASSERT(starpu_vector_get_nx(handle_s) == starpu_vector_get_nx(handle_r)); nx = starpu_vector_get_nx(handle_r); v_r = (int *)starpu_vector_get_local_ptr(handle_r); v_s = (int *)starpu_vector_get_local_ptr(handle_s); for(i=0 ; i %d] value: %c == %c\n", x, y, index, matrix_s[index], matrix_r[index]); } else { *error = 1; FPRINTF_MPI(stderr, "Error with matrix[%d,%d --> %d] value: %c != %c\n", x, y, index, matrix_s[index], matrix_r[index]); } } } } void check_block(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error) { STARPU_ASSERT(starpu_block_get_elemsize(handle_s) == starpu_block_get_elemsize(handle_r)); STARPU_ASSERT(starpu_block_get_nx(handle_s) == starpu_block_get_nx(handle_r)); STARPU_ASSERT(starpu_block_get_ny(handle_s) == starpu_block_get_ny(handle_r)); STARPU_ASSERT(starpu_block_get_nz(handle_s) == starpu_block_get_nz(handle_r)); STARPU_ASSERT(starpu_block_get_local_ldy(handle_s) == starpu_block_get_local_ldy(handle_r)); STARPU_ASSERT(starpu_block_get_local_ldz(handle_s) == starpu_block_get_local_ldz(handle_r)); starpu_data_acquire(handle_s, STARPU_R); starpu_data_acquire(handle_r, STARPU_R); float *block_s = (float *)starpu_block_get_local_ptr(handle_s); float *block_r = (float *)starpu_block_get_local_ptr(handle_r); int nx = starpu_block_get_nx(handle_s); int ny = starpu_block_get_ny(handle_s); int nz = starpu_block_get_nz(handle_s); int ldy = starpu_block_get_local_ldy(handle_s); int ldz = starpu_block_get_local_ldz(handle_s); int x, y, z; for(z=0 ; z %d] value: %f == %f\n", x, y, z, index, block_s[index], block_r[index]); } else { *error = 1; FPRINTF_MPI(stderr, "Error with block[%d,%d,%d --> %d] value: %f != %f\n", x, y, z, index, block_s[index], block_r[index]); } } } starpu_data_release(handle_s); starpu_data_release(handle_r); } void check_bcsr(starpu_data_handle_t handle_s, starpu_data_handle_t handle_r, int *error) { STARPU_ASSERT(starpu_bcsr_get_elemsize(handle_s) == starpu_bcsr_get_elemsize(handle_r)); STARPU_ASSERT(starpu_bcsr_get_nnz(handle_s) == starpu_bcsr_get_nnz(handle_r)); STARPU_ASSERT(starpu_bcsr_get_nrow(handle_s) == starpu_bcsr_get_nrow(handle_r)); STARPU_ASSERT(starpu_bcsr_get_firstentry(handle_s) == starpu_bcsr_get_firstentry(handle_r)); STARPU_ASSERT(starpu_bcsr_get_r(handle_s) == starpu_bcsr_get_r(handle_r)); STARPU_ASSERT(starpu_bcsr_get_c(handle_s) == starpu_bcsr_get_c(handle_r)); // STARPU_ASSERT(starpu_bcsr_get_local_colind(handle_s) == starpu_bcsr_get_local_colind(handle_r)); // STARPU_ASSERT(starpu_bcsr_get_local_rowptr(handle_s) == starpu_bcsr_get_local_rowptr(handle_r)); starpu_data_acquire(handle_s, STARPU_R); starpu_data_acquire(handle_r, STARPU_R); int *bcsr_s = (int *)starpu_bcsr_get_local_nzval(handle_s); int *bcsr_r = (int *)starpu_bcsr_get_local_nzval(handle_r); int r = starpu_bcsr_get_r(handle_s); int c = starpu_bcsr_get_c(handle_s); int nnz = starpu_bcsr_get_nnz(handle_s); int x; for(x=0 ; x