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
|
//------------------------------------------------------------------------------
// GB_conform_hyper: conform a sparse matrix to its desired hypersparse format
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
// The input matrix must be sparse or hypersparse, and it may be left as-is,
// or converted to sparse/hypersparse.
// The input matrix can have shallow A->p and/or A->h components. If the
// hypersparsity is changed, these components are no longer shallow.
#include "GB.h"
#define GB_FREE_ALL ;
GrB_Info GB_conform_hyper // conform a matrix to sparse/hypersparse
(
GrB_Matrix A, // matrix to conform
GB_Werk Werk
)
{
//--------------------------------------------------------------------------
// check inputs
//--------------------------------------------------------------------------
GrB_Info info = GrB_SUCCESS ;
ASSERT_MATRIX_OK (A, "A to conform_hyper", GB0) ;
ASSERT (!GB_IS_FULL (A)) ;
ASSERT (!GB_IS_BITMAP (A)) ;
ASSERT (GB_ZOMBIES_OK (A)) ;
ASSERT (GB_JUMBLED_OK (A)) ;
ASSERT (GB_PENDING_OK (A)) ;
//--------------------------------------------------------------------------
// convert to sparse or hypersparse
//--------------------------------------------------------------------------
// A->nvec_nonempty is used to select sparse vs hypersparse
int64_t nvec_nonempty = GB_nvec_nonempty_update (A) ;
if (A->h == NULL && GB_convert_sparse_to_hyper_test (A->hyper_switch,
nvec_nonempty, A->vdim))
{
// A is sparse but should be converted to hypersparse
GB_OK (GB_convert_sparse_to_hyper (A, Werk)) ;
}
else if (A->h != NULL && GB_convert_hyper_to_sparse_test (A->hyper_switch,
nvec_nonempty, A->vdim))
{
// A is hypersparse but should be converted to sparse
GB_OK (GB_convert_hyper_to_sparse (A, true)) ;
}
else
{
// leave the matrix as-is
;
}
//--------------------------------------------------------------------------
// return result
//--------------------------------------------------------------------------
ASSERT_MATRIX_OK (A, "A conform_hyper result", GB0) ;
return (GrB_SUCCESS) ;
}
|
