1 //------------------------------------------------------------------------------
2 // GB_convert_bitmap_to_sparse: convert a matrix from bitmap to sparse
3 //------------------------------------------------------------------------------
4
5 // SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2021, All Rights Reserved.
6 // SPDX-License-Identifier: Apache-2.0
7
8 //------------------------------------------------------------------------------
9
10 #include "GB.h"
11
12 #define GB_FREE_ALL \
13 { \
14 GB_FREE (&Ap, Ap_size) ; \
15 GB_FREE (&Ai, Ai_size) ; \
16 GB_FREE (&Ax, Ax_size) ; \
17 }
18
GB_convert_bitmap_to_sparse(GrB_Matrix A,GB_Context Context)19 GrB_Info GB_convert_bitmap_to_sparse // convert matrix from bitmap to sparse
20 (
21 GrB_Matrix A, // matrix to convert from bitmap to sparse
22 GB_Context Context
23 )
24 {
25
26 //--------------------------------------------------------------------------
27 // check inputs
28 //--------------------------------------------------------------------------
29
30 GrB_Info info ;
31 ASSERT_MATRIX_OK (A, "A converting bitmap to sparse", GB0) ;
32 ASSERT (!GB_IS_FULL (A)) ;
33 ASSERT (GB_IS_BITMAP (A)) ;
34 ASSERT (!GB_IS_SPARSE (A)) ;
35 ASSERT (!GB_IS_HYPERSPARSE (A)) ;
36 ASSERT (!GB_PENDING (A)) ; // bitmap never has pending tuples
37 ASSERT (!GB_JUMBLED (A)) ; // bitmap is never jumbled
38 ASSERT (!GB_ZOMBIES (A)) ; // bitmap never has zomies
39 GBURBLE ("(bitmap to sparse) ") ;
40
41 //--------------------------------------------------------------------------
42 // allocate Ap, Ai, and Ax
43 //--------------------------------------------------------------------------
44
45 const int64_t anz = GB_NNZ (A) ;
46 const int64_t anzmax = GB_IMAX (anz, 1) ;
47 int64_t anvec_nonempty ;
48 const int64_t avdim = A->vdim ;
49 const size_t asize = A->type->size ;
50 int64_t *restrict Ap = NULL ; size_t Ap_size = 0 ;
51 int64_t *restrict Ai = NULL ; size_t Ai_size = 0 ;
52 GB_void *restrict Ax = NULL ; size_t Ax_size = 0 ;
53 Ap = GB_MALLOC (avdim+1, int64_t, &Ap_size) ;
54 Ai = GB_MALLOC (anzmax, int64_t, &Ai_size) ;
55 Ax = GB_MALLOC (anzmax * asize, GB_void, &Ax_size) ;
56 if (Ap == NULL || Ai == NULL || Ax == NULL)
57 {
58 // out of memory
59 GB_FREE_ALL ;
60 return (GrB_OUT_OF_MEMORY) ;
61 }
62
63 //--------------------------------------------------------------------------
64 // convert to sparse format (Ap, Ai, and Ax)
65 //--------------------------------------------------------------------------
66
67 GB_OK (GB_convert_bitmap_worker (Ap, Ai, NULL, Ax, &anvec_nonempty, A,
68 Context)) ;
69
70 //--------------------------------------------------------------------------
71 // free prior content of A and transplant the new content
72 //--------------------------------------------------------------------------
73
74 GB_phbix_free (A) ;
75
76 A->p = Ap ; A->p_size = Ap_size ; A->p_shallow = false ;
77 A->i = Ai ; A->i_size = Ai_size ; A->i_shallow = false ;
78 A->x = Ax ; A->x_size = Ax_size ; A->x_shallow = false ;
79
80 A->nzmax = anzmax ;
81 A->nvals = 0 ; // only used when A is bitmap
82
83 A->plen = avdim ;
84 A->nvec = avdim ;
85 A->nvec_nonempty = anvec_nonempty ;
86
87 A->magic = GB_MAGIC ;
88
89 //--------------------------------------------------------------------------
90 // return result
91 //--------------------------------------------------------------------------
92
93 ASSERT_MATRIX_OK (A, "A converted from to bitmap to sparse", GB0) ;
94 ASSERT (GB_IS_SPARSE (A)) ;
95 ASSERT (!GB_ZOMBIES (A)) ;
96 ASSERT (!GB_JUMBLED (A)) ;
97 ASSERT (!GB_PENDING (A)) ;
98 return (GrB_SUCCESS) ;
99 }
100
101