pathrex-sys 0.1.0

//----------------------------------------------------------------------------
// LAGraph/src/test/test_SWrite.c: test cases for LAGraph_SWrite and SRead
// ----------------------------------------------------------------------------

// LAGraph, (c) 2019-2022 by The LAGraph Contributors, All Rights Reserved.
// SPDX-License-Identifier: BSD-2-Clause
//
// For additional details (including references to third party source code and
// other files) see the LICENSE file or contact permission@sei.cmu.edu. See
// Contributors.txt for a full list of contributors. Created, in part, with
// funding and support from the U.S. Government (see Acknowledgments.txt file).
// DM22-0790

// Contributed by Timothy A. Davis, Texas A&M University

//-----------------------------------------------------------------------------

#include <stdio.h>
#include <acutest.h>
#include <LAGraphX.h>
#include <LAGraph_test.h>
#include "LG_internal.h"

char msg [LAGRAPH_MSG_LEN] ;
LAGraph_Graph G = NULL ;
GrB_Matrix A = NULL ;
GrB_Matrix B = NULL ;
GrB_Matrix *S = NULL ;

#define LEN 512
char filename [LEN+1] ;

GrB_Type atype = NULL ;
char atypename [LAGRAPH_MAX_NAME_LEN] ;

#define NFILES 51
const char *files [ ] =
{
    "A.mtx",
    "cover.mtx",
    "cover_structure.mtx",
    "jagmesh7.mtx",
    "ldbc-cdlp-directed-example.mtx",
    "ldbc-cdlp-undirected-example.mtx",
    "ldbc-directed-example-bool.mtx",
    "ldbc-directed-example.mtx",
    "ldbc-directed-example-unweighted.mtx",
    "ldbc-undirected-example-bool.mtx",
    "ldbc-undirected-example.mtx",
    "ldbc-undirected-example-unweighted.mtx",
    "ldbc-wcc-example.mtx",
    "LFAT5.mtx",
    "msf1.mtx",
    "msf2.mtx",
    "msf3.mtx",
    "sample2.mtx",
    "sample.mtx",
    "sources_7.mtx",
    "olm1000.mtx",
    "bcsstk13.mtx",
    "cryg2500.mtx",
    "tree-example.mtx",
    "west0067.mtx",
    "lp_afiro.mtx",
    "lp_afiro_structure.mtx",
    "karate.mtx",
    "matrix_bool.mtx",
    "matrix_int8.mtx",
    "matrix_int16.mtx",
    "matrix_int32.mtx",
    "matrix_int64.mtx",
    "matrix_uint8.mtx",
    "matrix_uint16.mtx",
    "matrix_uint32.mtx",
    "matrix_uint64.mtx",
    "matrix_fp32.mtx",
    "matrix_fp32_structure.mtx",
    "matrix_fp64.mtx",
    "west0067_jumbled.mtx",
    "skew_fp32.mtx",
    "skew_fp64.mtx",
    "skew_int8.mtx",
    "skew_int16.mtx",
    "skew_int32.mtx",
    "skew_int64.mtx",
    "structure.mtx",
    "full.mtx",
    "full_symmetric.mtx",
    "empty.mtx",
    "",
} ;

//****************************************************************************

void test_SWrite (void)
{
    #if LAGRAPH_SUITESPARSE
    LAGraph_Init (msg) ;

    for (int k = 0 ; k < NFILES ; k++)
    {

        // load the matrix as A
        const char *aname = files [k] ;
        if (strlen (aname) == 0) break;
        printf ("\n================================== %d %s:\n", k, aname) ;
        TEST_CASE (aname) ;
        snprintf (filename, LEN, LG_DATA_DIR "%s", aname) ;
        FILE *f = fopen (filename, "r") ;
        TEST_CHECK (f != NULL) ;
        OK (LAGraph_MMRead (&A, f, msg)) ;
        fclose (f) ;

        // get the name of the C typedef for the matrix
        OK (LAGraph_Matrix_TypeName (atypename, A, msg)) ;
        OK (LAGraph_TypeFromName (&atype, atypename, msg)) ;

        for (int scon = 1 ; scon <= 8 ; scon = 2*scon)
        {
            // for SuiteSparse only: test all sparsity formats
            OK (LG_SET_FORMAT_HINT (A, scon)) ;

            // open a temporary *.lagraph file to hold the matrix
            f = tmpfile ( ) ;
            // snprintf (filename, LEN, "%s.lagraph", aname) ;
            // f = fopen (filename, "w") ;
            TEST_CHECK (f != NULL) ;

            // serialize the matrix
            void *blob = NULL ;
            GrB_Index blob_size = 0 ;

            if (k % 2 == 0)
            {
                // for SuiteSparse
                OK (GxB_Matrix_serialize (&blob, &blob_size, A, NULL)) ;
            }
            else
            {
                // try GrB version
                OK (GrB_Matrix_serializeSize (&blob_size, A)) ;
                GrB_Index blob_size_old = blob_size ;
                OK (LAGraph_Malloc ((void **) &blob, blob_size,
                    sizeof (uint8_t), msg)) ;
                TEST_CHECK (blob != NULL) ;
                OK (GrB_Matrix_serialize (blob, &blob_size, A)) ;
                OK (LAGraph_Realloc ((void **) &blob, blob_size,
                    blob_size_old, sizeof (uint8_t), msg)) ;
            }

            // deserialize the matrix
            int rr = (GrB_Matrix_deserialize (&B, atype, blob, blob_size)) ;
            printf ("deserialize result: %d\n", rr) ;
            GxB_print (B, 2) ;
            OK (rr) ;

            // ensure the matrices A and B are the same
            bool ok = false ;
            OK (LAGraph_Matrix_IsEqual (&ok, A, B, msg)) ;
            TEST_CHECK (ok) ;
            OK (GrB_free (&B)) ;

            // write the header for a single matrix
            OK (LAGraph_SWrite_HeaderStart (f, "lagraph_test", msg)) ;
            OK (LAGraph_SWrite_HeaderItem (f, LAGraph_matrix_kind, "A",
                atypename, 0, blob_size, msg)) ;
            OK (LAGraph_SWrite_HeaderEnd (f, msg)) ;

            // write the binary blob to the file then free the blob
            OK (LAGraph_SWrite_Item (f, blob, blob_size, msg)) ;
            LAGraph_Free (&blob, NULL) ;

            // open the file and load back the contents
            rewind (f) ;
            // fclose (f) ;
            // f = fopen (filename, "r") ;

            char *collection = NULL ;
            LAGraph_Contents *Contents = NULL ;
            GrB_Index ncontents ;
            OK (LAGraph_SRead (f, &collection, &Contents, &ncontents, msg)) ;
            TEST_CHECK (collection != NULL) ;
            if (collection == NULL) abort ( ) ;
            // printf ("collection %s\n", collection) ;
            TEST_CHECK (strcmp (collection, "lagraph_test") == 0) ;
            TEST_CHECK (ncontents == 1) ;
            fclose (f) ;

            // convert the contents to a matrix B
            void *blob2 = Contents [0].blob ;
            size_t blob_size2 = Contents [0].blob_size ;
            // printf ("blob_size2 %lu\n", blob_size2) ;
            TEST_CHECK (blob_size == blob_size2) ;

            OK (GrB_Matrix_deserialize (&B, atype, blob2, blob_size2)) ;

            // ensure the matrices A and B are the same
            OK (LAGraph_Matrix_IsEqual (&ok, A, B, msg)) ;
            TEST_CHECK (ok) ;
            OK (GrB_free (&B)) ;

            // free the contents: todo make this a utility function
            LAGraph_Free ((void **) &collection, NULL) ;
            for (int i = 0 ; i < ncontents ; i++)
            {
                LAGraph_Contents *Item = &(Contents [i]) ;
                LAGraph_Free ((void **) &(Item->blob), NULL) ;
            }
            LAGraph_Free ((void **) &Contents, NULL) ;
        }

        OK (GrB_free (&A)) ;
    }

    LAGraph_Finalize (msg) ;
    #endif
}

//------------------------------------------------------------------------------

void test_SWrite_errors (void)
{
    #if LAGRAPH_SUITESPARSE
    LAGraph_Init (msg) ;

    // create a simple test matrix
    GrB_Index n = 5 ;
    OK (GrB_Matrix_new (&A, GrB_FP32, n, n)) ;
    OK (GrB_assign (A, NULL, NULL, 0, GrB_ALL, n, GrB_ALL, n, NULL)) ;
    OK (GrB_apply (A, NULL, NULL, GrB_ROWINDEX_INT64, A, 0, NULL)) ;
    printf ("\nTest matrix:\n") ;
    OK (LAGraph_Matrix_Print (A, LAGraph_COMPLETE, stdout, msg)) ;

    // serialize the matrix
    bool ok ;
    void *blob = NULL ;
    GrB_Index blob_size = 0 ;

    // for SuiteSparse
    OK (GxB_Matrix_serialize (&blob, &blob_size, A, NULL)) ;

    FILE *f = tmpfile ( )  ;
    TEST_CHECK (f != NULL) ;

    int result = LAGraph_SWrite_HeaderItem (f, -2, "A",
        "float", 0, blob_size, msg) ;
    printf ("result: %d [%s]\n", result, msg) ;
    TEST_CHECK (result == GrB_INVALID_VALUE) ;
    fclose (f) ;

    f = fopen ("error.lagraph", "wb") ;
    TEST_CHECK (f != NULL) ;

    result = LAGraph_SWrite_HeaderStart (f, NULL, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    result = LAGraph_SWrite_HeaderStart (NULL, "stuff", msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    OK (LAGraph_SWrite_HeaderStart (f, "lagraph_test", msg)) ;

    result = LAGraph_SWrite_HeaderItem (NULL, LAGraph_matrix_kind, "A",
        "float", 0, blob_size, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    result = LAGraph_SWrite_HeaderItem (NULL, -2, "A",
        "float", 0, blob_size, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    OK (LAGraph_SWrite_HeaderItem (f, LAGraph_matrix_kind, "A",
        "float", 0, blob_size, msg)) ;

    result = LAGraph_SWrite_HeaderEnd (NULL, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    OK (LAGraph_SWrite_HeaderEnd (f, msg)) ;

    // write the binary blob to the file then free the blob
    OK (LAGraph_SWrite_Item (f, blob, blob_size, msg)) ;
    LAGraph_Free (&blob, NULL) ;

    result = LAGraph_SWrite_Item (NULL, blob, blob_size, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    result = LAGraph_SWrite_Item (f, NULL, blob_size, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    // close the file and reopen it
    fclose (f) ;
    f = fopen ("error.lagraph", "r") ;
    TEST_CHECK (f != NULL) ;

    // load in the matrix
    GrB_Matrix *Set = NULL ;
    GrB_Index nmatrices = 0 ;
    char *collection = NULL ;

    result = LAGraph_SLoadSet ("error.lagraph", NULL, &nmatrices,
        &collection, msg) ;
    TEST_CHECK (result == GrB_NULL_POINTER) ;

    result = LAGraph_SLoadSet ("unknown.lagraph", &Set, &nmatrices,
        &collection, msg) ;
    printf ("\nresult %d, [%s]\n", result, msg) ;
    TEST_CHECK (result == -1002) ;

    // mangled file
    result = LAGraph_SLoadSet (LG_DATA_DIR "garbage.lagraph",
        &Set, &nmatrices, &collection, msg) ;
    printf ("\nresult %d, [%s]\n", result, msg) ;
    TEST_CHECK (result == LAGRAPH_IO_ERROR) ;

    // finally works
    OK (LAGraph_SLoadSet ("error.lagraph", &Set, &nmatrices, &collection,
        msg)) ;
    TEST_CHECK (Set != NULL) ;
    TEST_CHECK (collection != NULL) ;
    TEST_CHECK (nmatrices == 1) ;

    ok = false ;
    OK (LAGraph_Matrix_IsEqual (&ok, A, Set [0], msg)) ;
    TEST_CHECK (ok) ;

    // free everything
    LAGraph_SFreeSet (&Set, nmatrices) ;
    LAGraph_Free ((void **) &collection, NULL) ;
    fclose (f) ;

    // read garbage with LAGraph_SRead
    f = fopen (LG_DATA_DIR "garbage.lagraph", "r") ;
    TEST_CHECK (f != NULL) ;
    LAGraph_Contents *Contents = NULL ;
    GrB_Index ncontents = 0 ;
    result = LAGraph_SRead (f, &collection, &Contents, &ncontents, msg) ;
    TEST_CHECK (result = -1001) ;
    TEST_CHECK (collection == NULL) ;
    TEST_CHECK (Contents == NULL) ;
    TEST_CHECK (ncontents == 0) ;
    fclose (f) ;

    OK (GrB_free (&A)) ;
    LAGraph_Finalize (msg) ;
    #endif
}

//****************************************************************************

TEST_LIST = {
    {"SWrite", test_SWrite},
    {"SWrite_errors", test_SWrite_errors},
    {NULL, NULL}
};