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
|
//------------------------------------------------------------------------------
// GB_subref_sparse_jit: C=A(I,J) when C and A are sparse/hypersparse
//------------------------------------------------------------------------------
// 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_SUBREF_SPARSE_PROTO ((*GB_jit_dl_function)) ;
GrB_Info GB_subref_sparse_jit
(
// output matrix
GrB_Matrix C, // same type as A
// from phase1:
const GB_task_struct *restrict TaskList, // list of tasks
const int ntasks, // # of tasks
const int nthreads, // # of threads to use
const bool post_sort, // true if post-sort needed
const void *Ihead, // for I inverse buckets, size A->vlen
const void *Inext, // for I inverse buckets, size nI
const bool Ihead_is_32, // if true, Ihead/Inext 32-bit; else 64
const bool I_has_duplicates, // true if I has duplicates
// from phase0:
const void *Ap_start,
const void *Ap_end,
const bool need_qsort,
const int Ikind,
const int64_t nI,
const int64_t Icolon [3],
// original input:
const GrB_Matrix A,
const void *I,
const bool I_is_32
)
{
//--------------------------------------------------------------------------
// encodify the problem
//--------------------------------------------------------------------------
GB_jit_encoding encoding ;
char *suffix ;
uint64_t hash = GB_encodify_subref (&encoding, &suffix,
GB_JIT_KERNEL_SUBREF_SPARSE, C, I_is_32, false, Ikind, 0,
need_qsort, Ihead_is_32, I_has_duplicates, A) ;
//--------------------------------------------------------------------------
// get the kernel function pointer, loading or compiling it if needed
//--------------------------------------------------------------------------
void *dl_function ;
GrB_Info info = GB_jitifyer_load (&dl_function,
GB_jit_subref_family, "subref_sparse",
hash, &encoding, suffix, NULL, NULL,
NULL, C->type, NULL, NULL) ;
if (info != GrB_SUCCESS) return (info) ;
//--------------------------------------------------------------------------
// call the jit kernel and return result
//--------------------------------------------------------------------------
#include "include/GB_pedantic_disable.h"
GB_jit_dl_function GB_jit_kernel = (GB_jit_dl_function) dl_function ;
return (GB_jit_kernel (C, TaskList, ntasks, nthreads, post_sort, Ihead,
Inext, Ap_start, Ap_end, nI, Icolon, A, I, &GB_callback)) ;
}
|
