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
|
//------------------------------------------------------------------------------
// GB_subassign_jit: interface to JIT kernels for all assign/subassign methods
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
#include "GB.h"
#include "jitifyer/GB_stringify.h"
typedef GB_JIT_KERNEL_SUBASSIGN_PROTO ((*GB_jit_dl_function)) ;
GrB_Info GB_subassign_jit
(
// input/output:
GrB_Matrix C,
// input:
const bool C_replace,
// I:
const void *I,
const bool I_is_32,
const int64_t ni,
const int64_t nI,
const int Ikind,
const int64_t Icolon [3],
// J:
const void *J,
const bool J_is_32,
const int64_t nj,
const int64_t nJ,
const int Jkind,
const int64_t Jcolon [3],
// mask M:
const GrB_Matrix M,
const bool Mask_comp,
const bool Mask_struct,
// accum, if present:
const GrB_BinaryOp accum, // may be NULL
// A matrix or scalar:
const GrB_Matrix A, // NULL for scalar assignment
const void *scalar,
const GrB_Type scalar_type,
// S matrix:
const GrB_Matrix S, // NULL if not constructed
// kind and kernel:
const int assign_kind, // row assign, col assign, assign, or subassign
const int assign_kernel, // GB_JIT_KERNEL_SUBASSIGN_01, ... etc
const char *kname, // kernel name
GB_Werk Werk
)
{
//--------------------------------------------------------------------------
// encodify the problem
//--------------------------------------------------------------------------
GB_jit_encoding encoding ;
char *suffix ;
uint64_t hash = GB_encodify_assign (&encoding, &suffix, assign_kernel,
C, C_replace, I_is_32, J_is_32, Ikind, Jkind, M, Mask_comp,
Mask_struct, accum, A, scalar_type, S, assign_kind) ;
//--------------------------------------------------------------------------
// get the kernel function pointer, loading or compiling it if needed
//--------------------------------------------------------------------------
GrB_Type atype = (A == NULL) ? scalar_type : A->type ;
void *dl_function ;
GrB_Info info = GB_jitifyer_load (&dl_function,
GB_jit_assign_family, kname,
hash, &encoding, suffix, NULL, NULL,
(GB_Operator) accum, C->type, atype, NULL) ;
if (info != GrB_SUCCESS) return (info) ;
//--------------------------------------------------------------------------
// call the jit kernel and return result
//--------------------------------------------------------------------------
double chunk = GB_Context_chunk ( ) ;
int nthreads_max = GB_Context_nthreads_max ( ) ;
#include "include/GB_pedantic_disable.h"
GB_jit_dl_function GB_jit_kernel = (GB_jit_dl_function) dl_function ;
return (GB_jit_kernel (C, C_replace,
I, ni, nI, Ikind, Icolon,
J, nj, nJ, Jkind, Jcolon,
M, Mask_comp, Mask_struct,
accum, A, scalar, scalar_type, S, assign_kind, Werk,
nthreads_max, chunk, &GB_callback)) ;
}
|
