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
|
//------------------------------------------------------------------------------
// GB_mex_random: construct a random matrix, double or Complex
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2022, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
#include "GB_mex.h"
#define USAGE "A = GB_mex_random (nrows, ncols, ntuples," \
" complex, seed, make_symmetric, no_self_edges, method)"
#define GET_DEEP_COPY ;
#define FREE_DEEP_COPY ;
#define FREE_ALL GB_mx_put_global (true) ;
void mexFunction
(
int nargout,
mxArray *pargout [ ],
int nargin,
const mxArray *pargin [ ]
)
{
bool malloc_debug = GB_mx_get_global (true) ;
GrB_Matrix A = NULL ;
// check inputs
if (nargout > 1 || nargin == 0 || nargin > 8)
{
mexErrMsgTxt ("Usage: " USAGE) ;
}
int64_t GET_SCALAR (0, int64_t , nrows, 4) ;
int64_t GET_SCALAR (1, int64_t , ncols, 4) ;
int64_t GET_SCALAR (2, int64_t , ntuples, 0) ;
bool GET_SCALAR (3, bool , A_complex, false) ;
int64_t GET_SCALAR (4, uint64_t, seed, 1) ;
bool GET_SCALAR (5, bool , make_symmetric, false) ;
bool GET_SCALAR (6, bool , no_self_edges, false) ;
int GET_SCALAR (7, int , method, 0) ;
// construct the random matrix
simple_rand_seed (seed) ;
if (method == 3)
{
// test out-of-memory condition
METHOD (GB_mx_random_matrix (&A, make_symmetric, no_self_edges,
nrows, ncols, ntuples, method, A_complex)) ;
}
else
{
GrB_Info info = GB_mx_random_matrix (&A, make_symmetric, no_self_edges,
nrows, ncols, ntuples, method, A_complex) ;
if (info != GrB_SUCCESS)
{
mexErrMsgTxt ("GB_mx_random_matrix failed") ;
}
}
// return A
pargout [0] = GB_mx_Matrix_to_mxArray (&A, "A final", false) ;
FREE_ALL ;
}
|