1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
//------------------------------------------------------------------------------
// GB_aop: assign/subassign kernels with accum
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
// C(I,J)<M> += A
#include "GB_control.h"
#if defined (GxB_NO_INT64)
#define GB_TYPE_ENABLED 0
#else
#define GB_TYPE_ENABLED 1
#endif
#if GB_TYPE_ENABLED
#include "GB.h"
#include "FactoryKernels/GB_aop__include.h"
// accum operator
#define GB_ACCUM_OP(z,x,y) z = GB_pow_int64 (x, y)
#define GB_Z_TYPE int64_t
#define GB_X_TYPE int64_t
#define GB_Y_TYPE int64_t
#define GB_DECLAREY(ywork) int64_t ywork
#define GB_COPY_aij_to_ywork(ywork,Ax,pA,A_iso) ywork = Ax [(A_iso) ? 0 : (pA)]
// A and C matrices
#define GB_A_TYPE int64_t
#define GB_C_TYPE int64_t
#define GB_DECLAREC(cwork) int64_t cwork
#define GB_COPY_aij_to_cwork(cwork,Ax,pA,A_iso) cwork = Ax [A_iso ? 0 : (pA)]
#define GB_COPY_aij_to_C(Cx,pC,Ax,pA,A_iso,cwork,C_iso) Cx [pC] = (A_iso) ? cwork : Ax [pA]
#define GB_COPY_cwork_to_C(Cx,pC,cwork,C_iso) Cx [pC] = cwork
#define GB_AX_MASK(Ax,pA,asize) (Ax [pA] != 0)
// C(i,j) += ywork
#define GB_ACCUMULATE_scalar(Cx,pC,ywork,C_iso) \
GB_ACCUM_OP (Cx [pC], Cx [pC], ywork)
// C(i,j) += (ytype) A(i,j)
#define GB_ACCUMULATE_aij(Cx,pC,Ax,pA,A_iso,ywork,C_iso) \
{ \
if (A_iso) \
{ \
GB_ACCUMULATE_scalar (Cx, pC, ywork, C_iso) ; \
} \
else \
{ \
/* A and Y have the same type here */ \
GB_ACCUMULATE_scalar (Cx, pC, Ax [pA], C_iso) ; \
} \
}
// disable this operator and use the generic case if these conditions hold
#if (defined(GxB_NO_POW) || defined(GxB_NO_INT64) || defined(GxB_NO_POW_INT64))
#define GB_DISABLE 1
#else
#define GB_DISABLE 0
#endif
#include "assign/include/GB_assign_shared_definitions.h"
//------------------------------------------------------------------------------
// C += A, accumulate a sparse matrix into a dense matrix
//------------------------------------------------------------------------------
#undef GB_SCALAR_ASSIGN
#define GB_SCALAR_ASSIGN 0
GrB_Info GB (_subassign_23__pow_int64)
(
GrB_Matrix C,
const GrB_Matrix A,
GB_Werk Werk
)
{
#if GB_DISABLE
return (GrB_NO_VALUE) ;
#else
int nthreads_max = GB_Context_nthreads_max ( ) ;
double chunk = GB_Context_chunk ( ) ;
#include "assign/template/GB_subassign_23_template.c"
return (GrB_SUCCESS) ;
#endif
}
//------------------------------------------------------------------------------
// C += y, accumulate a scalar into a dense matrix
//------------------------------------------------------------------------------
#undef GB_SCALAR_ASSIGN
#define GB_SCALAR_ASSIGN 1
GrB_Info GB (_subassign_22__pow_int64)
(
GrB_Matrix C,
const GB_void *ywork_handle
)
{
#if GB_DISABLE
return (GrB_NO_VALUE) ;
#else
// get the scalar ywork for C += ywork, of type GB_Y_TYPE
GB_Y_TYPE ywork = (*((GB_Y_TYPE *) ywork_handle)) ;
int nthreads_max = GB_Context_nthreads_max ( ) ;
double chunk = GB_Context_chunk ( ) ;
#include "assign/template/GB_subassign_22_template.c"
return (GrB_SUCCESS) ;
#endif
}
#else
GB_EMPTY_PLACEHOLDER
#endif
|
