//------------------------------------------------------------------------------ // GxB_Matrix_export_HyperCSR: export a matrix in hypersparse CSR format //------------------------------------------------------------------------------ // SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2023, All Rights Reserved. // SPDX-License-Identifier: Apache-2.0 //------------------------------------------------------------------------------ #include "GB_export.h" #define GB_FREE_ALL ; GrB_Info GxB_Matrix_export_HyperCSR // export and free a hypersparse CSR matrix ( GrB_Matrix *A, // handle of matrix to export and free GrB_Type *type, // type of matrix exported GrB_Index *nrows, // number of rows of the matrix GrB_Index *ncols, // number of columns of the matrix GrB_Index **Ap, // row "pointers" GrB_Index **Ah, // row indices GrB_Index **Aj, // column indices void **Ax, // values GrB_Index *Ap_size, // size of Ap in bytes GrB_Index *Ah_size, // size of Ah in bytes GrB_Index *Aj_size, // size of Aj in bytes GrB_Index *Ax_size, // size of Ax in bytes bool *iso, // if true, A is iso GrB_Index *nvec, // number of rows that appear in Ah bool *jumbled, // if true, indices in each row may be unsorted const GrB_Descriptor desc ) { //-------------------------------------------------------------------------- // check inputs and get the descriptor //-------------------------------------------------------------------------- GB_WHERE1 ("GxB_Matrix_export_HyperCSR (&A, &type, &nrows, &ncols, " "&Ap, &Ah, &Aj, &Ax, &Ap_size, &Ah_size, &Aj_size, &Ax_size, " "&iso, &nvec, &jumbled, desc)") ; // GB_BURBLE_START ("GxB_Matrix_export_HyperCSR") ; GB_RETURN_IF_NULL (A) ; GB_RETURN_IF_NULL_OR_FAULTY (*A) ; GB_GET_DESCRIPTOR (info, desc, xx1, xx2, xx3, xx4, xx5, xx6, xx7) ; //-------------------------------------------------------------------------- // ensure the matrix is in by-row format //-------------------------------------------------------------------------- if ((*A)->is_csc) { // A = A', done in-place, to put A in by-row format GBURBLE ("(export transpose) ") ; GB_OK (GB_transpose_in_place (*A, false, Werk)) ; } //-------------------------------------------------------------------------- // finish any pending work //-------------------------------------------------------------------------- if (jumbled == NULL) { // the exported matrix cannot be jumbled GB_MATRIX_WAIT (*A) ; } else { // the exported matrix is allowed to be jumbled GB_MATRIX_WAIT_IF_PENDING_OR_ZOMBIES (*A) ; } //-------------------------------------------------------------------------- // ensure the matrix is hypersparse //-------------------------------------------------------------------------- GB_OK (GB_convert_any_to_hyper (*A, Werk)) ; //-------------------------------------------------------------------------- // export the matrix //-------------------------------------------------------------------------- ASSERT (GB_IS_HYPERSPARSE (*A)) ; ASSERT (!(*A)->is_csc) ; ASSERT (!GB_ZOMBIES (*A)) ; ASSERT (GB_IMPLIES (jumbled == NULL, !GB_JUMBLED (*A))) ; ASSERT (!GB_PENDING (*A)) ; int sparsity ; bool is_csc ; info = GB_export (false, A, type, ncols, nrows, false, Ap, Ap_size, // Ap Ah, Ah_size, // Ah NULL, NULL, // Ab Aj, Aj_size, // Aj Ax, Ax_size, // Ax NULL, jumbled, nvec, // jumbled or not &sparsity, &is_csc, // hypersparse by row iso, Werk) ; if (info == GrB_SUCCESS) { ASSERT (sparsity == GxB_HYPERSPARSE) ; ASSERT (!is_csc) ; } // GB_BURBLE_END ; return (info) ; }