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
|
//------------------------------------------------------------------------------
// GxB_Matrix_build_Scalar: build a sparse GraphBLAS matrix
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
// GxB_Matrix_build_Scalar builds a matrix C whose values in its sparsity
// pattern are all equal to a value given by a GrB_Scalar. Unlike the
// GrB_Matrix_build_* methods, there is no binary dup operator. Instead, any
// duplicate indices are ignored, which is not an error condition. The I and J
// arrays are of size nvals, just like GrB_Matrix_build_*.
#include "builder/GB_build.h"
#define GB_FREE_ALL ;
GrB_Info GxB_Matrix_build_Scalar
(
GrB_Matrix C, // matrix to build
const uint64_t *I, // array of row indices of tuples
const uint64_t *J, // array of column indices of tuples
GrB_Scalar scalar, // value for all tuples
uint64_t nvals // number of tuples
)
{
//--------------------------------------------------------------------------
// check inputs
//--------------------------------------------------------------------------
GB_RETURN_IF_NULL (C) ;
GB_RETURN_IF_NULL (scalar) ;
GB_RETURN_IF_OUTPUT_IS_READONLY (C) ;
GB_WHERE2 (C, scalar, "GxB_Matrix_build_Scalar (C, I, J, scalar, nvals)") ;
GB_BURBLE_START ("GxB_Matrix_build_Scalar") ;
GB_MATRIX_WAIT (scalar) ;
if (GB_nnz ((GrB_Matrix) scalar) != 1)
{
GB_ERROR (GrB_EMPTY_OBJECT, "Scalar value is %s", "missing") ;
}
//--------------------------------------------------------------------------
// build the matrix, ignoring duplicates
//--------------------------------------------------------------------------
info = GB_build (C, I, J, scalar->x, nvals, GxB_IGNORE_DUP, scalar->type,
/* is_matrix: */ true, /* X_iso: */ true,
/* I,J is 32: */ false, false, Werk) ;
GB_BURBLE_END ;
return (info) ;
}
|
