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
|
//------------------------------------------------------------------------------
// GB_union_jit: C=A+B, C<#M>=A+B eWiseUnion method, via the JIT
//------------------------------------------------------------------------------
// 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_UNION_PROTO ((*GB_jit_dl_function)) ;
GrB_Info GB_union_jit // C=A+B, C<#M>=A+B, eWiseUnion, via the JIT
(
// input/output:
GrB_Matrix C,
// input:
const int C_sparsity,
const GrB_Matrix M,
const bool Mask_struct,
const bool Mask_comp,
const GrB_BinaryOp binaryop,
const bool flipij,
const GrB_Matrix A,
const GrB_Matrix B,
const GB_void *alpha_scalar_in,
const GB_void *beta_scalar_in,
const bool Ch_is_Mh,
const int64_t *restrict C_to_M,
const int64_t *restrict C_to_A,
const int64_t *restrict C_to_B,
const GB_task_struct *restrict TaskList,
const int C_ntasks,
const int C_nthreads,
const int64_t *restrict M_ek_slicing,
const int M_nthreads,
const int M_ntasks,
const int64_t *restrict A_ek_slicing,
const int A_nthreads,
const int A_ntasks,
const int64_t *restrict B_ek_slicing,
const int B_nthreads,
const int B_ntasks
)
{
//--------------------------------------------------------------------------
// encodify the problem
//--------------------------------------------------------------------------
GB_jit_encoding encoding ;
char *suffix ;
uint64_t hash = GB_encodify_ewise (&encoding, &suffix,
GB_JIT_KERNEL_UNION, false,
false, false, C_sparsity, C->type, C->p_is_32, C->j_is_32, C->i_is_32,
M, Mask_struct, Mask_comp, binaryop, flipij, false, A, B) ;
//--------------------------------------------------------------------------
// get the kernel function pointer, loading or compiling it if needed
//--------------------------------------------------------------------------
void *dl_function ;
GrB_Info info = GB_jitifyer_load (&dl_function,
GB_jit_ewise_family, "union",
hash, &encoding, suffix, NULL, NULL,
(GB_Operator) binaryop, C->type, A->type, B->type) ;
if (info != GrB_SUCCESS) return (info) ;
//--------------------------------------------------------------------------
// call the jit kernel and return result
//--------------------------------------------------------------------------
// for the "easy mask" condition:
bool M_is_A = GB_all_aliased (M, A) ;
bool M_is_B = GB_all_aliased (M, B) ;
#include "include/GB_pedantic_disable.h"
GB_jit_dl_function GB_jit_kernel = (GB_jit_dl_function) dl_function ;
return (GB_jit_kernel (C, M, A, B, alpha_scalar_in, beta_scalar_in,
Ch_is_Mh, C_to_M, C_to_A, C_to_B, TaskList, C_ntasks, C_nthreads,
M_ek_slicing, M_nthreads, M_ntasks, A_ek_slicing, A_nthreads, A_ntasks,
B_ek_slicing, B_nthreads, B_ntasks, M_is_A, M_is_B,
binaryop->theta, &GB_callback)) ;
}
|
