//------------------------------------------------------------------------------ // GB_red: hard-coded functions for reductions //------------------------------------------------------------------------------ // SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2023, All Rights Reserved. // SPDX-License-Identifier: Apache-2.0 //------------------------------------------------------------------------------ #include "GB.h" #include "GB_control.h" #include "GB_red__include.h" // reduction operator and type: #define GB_UPDATE(z,a) if (a < z) { z = a ; } #define GB_ADD(z,zin,a) z = GB_IMIN (zin, a) #define GB_GETA_AND_UPDATE(z,Ax,p) if (Ax [p] < z) { z = Ax [p] ; } // A matrix (no typecasting to Z type here) #define GB_A_TYPE uint16_t #define GB_DECLAREA(aij) uint16_t aij #define GB_GETA(aij,Ax,pA,A_iso) aij = Ax [pA] // monoid properties: #define GB_Z_TYPE uint16_t #define GB_DECLARE_IDENTITY(z) uint16_t z = UINT16_MAX #define GB_DECLARE_IDENTITY_CONST(z) const uint16_t z = UINT16_MAX #define GB_MONOID_IS_TERMINAL 1 #define GB_TERMINAL_CONDITION(z,zterminal) (z == 0) #define GB_IF_TERMINAL_BREAK(z,zterminal) if (z == 0) { break ; } #define GB_DECLARE_TERMINAL_CONST(zterminal) const uint16_t zterminal = 0 // panel size #define GB_PANEL 16 // disable this operator and use the generic case if these conditions hold #if (defined(GxB_NO_MIN) || defined(GxB_NO_UINT16) || defined(GxB_NO_MIN_UINT16)) #define GB_DISABLE 1 #else #define GB_DISABLE 0 #endif #include "GB_monoid_shared_definitions.h" //------------------------------------------------------------------------------ // reduce to a non-iso matrix to scalar, for monoids only //------------------------------------------------------------------------------ GrB_Info GB (_red__min_uint16) ( GB_Z_TYPE *result, const GrB_Matrix A, GB_void *restrict W_space, bool *restrict F, int ntasks, int nthreads ) { #if GB_DISABLE return (GrB_NO_VALUE) ; #else GB_Z_TYPE z = (*result) ; GB_Z_TYPE *restrict W = (GB_Z_TYPE *) W_space ; if (A->nzombies > 0 || GB_IS_BITMAP (A)) { #include "GB_reduce_to_scalar_template.c" } else { #include "GB_reduce_panel.c" } (*result) = z ; return (GrB_SUCCESS) ; #endif }