added flow mtx checking script

Author: dperry17

experimental/test/LG_check_flow.c

@@ -0,0 +1,36 @@
+
+#include "LG_internal.h"
+#include <LG_Xtest.h>
+
+#undef LG_FREE_ALL
+#undef LG_FREE_WORK
+
+#define LG_FREE_WORK				\
+{						\
+ GrB_free(&flow_r);				\
+ GrB_free(&flow_c);				\
+ GrB_free(&result_vec);				\
+}
+
+
+#define LG_FREE_ALL				\
+{						\
+  LG_FREE_WORK;					\
+}
+
+
+int LG_check_flow(const GrB_Matrix *flow_mtx, char* msg) {
+  GrB_Vector flow_r=NULL, flow_c=NULL, result_vec=NULL ;
+  GrB_Index n ;
+  double net_flow = -1;
+  LG_TRY(GrB_Matrix_nrows(&n, *flow_mtx));
+  LG_TRY(GrB_Vector_new(&flow_r, GrB_FP64, n));
+  LG_TRY(GrB_Vector_new(&flow_c, GrB_FP64, n));
+  LG_TRY(GrB_Vector_new(&result_vec, GrB_FP64, n));
+  LG_TRY(GrB_reduce(flow_c, NULL, NULL, GrB_PLUS_MONOID_FP64, *flow_mtx, NULL));
+  LG_TRY(GrB_reduce(flow_r, NULL, NULL, GrB_PLUS_MONOID_FP64, *flow_mtx, GrB_DESC_T1));
+  LG_TRY(GrB_eWiseAdd(result_vec, NULL, NULL, GrB_MINUS_FP64, flow_r, flow_c, NULL));
+  LG_TRY(GrB_reduce(&net_flow, NULL, GrB_PLUS_MONOID_FP64, result_vec, NULL));
+  LG_ASSERT_MSG(net_flow == 0, GrB_INVALID_VALUE, "Flow conservation is not followed");
+  return GrB_SUCCESS;
+}

experimental/test/include/LG_Xtest.h

@@ -116,6 +116,12 @@ int LG_check_argminmax
     int dim,                    // dim=1: cols of A, dim=2: rows of A
     bool is_min,
     char *msg
+)
+
+int LG_check_flow
+(
+   const GrB_Matrix *flow_mtx,
+   char* msg
 ) ;
 
 #endif

experimental/test/test_MaxFlow.c

@@ -1,9 +1,9 @@
+
 #include <acutest.h>
 #include <LAGraphX.h>
 #include <LAGraph_test.h>
 #include <stdio.h>
-
-#include "LG_Xtest.h"
+#include <LG_Xtest.h>
 #include "LG_internal.h"
 
 
@@ -63,4 +63,45 @@ void test_MaxFlow(void) {
   LAGraph_Finalize(msg);
 }
 
-TEST_LIST = {{"MaxFlow", test_MaxFlow}, {NULL, NULL}};
+void test_MaxFlowMtx(void) {
+  LAGraph_Init(msg);
+//OK(LG_SET_BURBLE(1));
+  OK(LG_SET_BURBLE(0));
+  OK(GxB_Global_Option_set(GxB_JIT_C_CONTROL, 4));
+  for(uint8_t test = 0; test < NTESTS; test++){
+    GrB_Matrix A=NULL;
+    TEST_CASE(tests[test].filename);
+    snprintf(filename, LEN, LG_DATA_DIR "%s", tests[test].filename);
+    FILE* f = fopen(filename, "r");
+    TEST_CHECK(f != NULL);
+    OK(LAGraph_MMRead(&A, f, msg));
+
+    //create flow mtx
+    GrB_Matrix flow_mtx=NULL;
+    GrB_Index n;
+    OK(GrB_Matrix_nrows(&n, A));
+    OK(GrB_Matrix_new(&flow_mtx, GrB_FP64, n, n));
+    
+    OK(fclose(f));
+    OK(LAGraph_New(&G, &A, LAGraph_ADJACENCY_DIRECTED, msg));
+    OK(LAGraph_Cached_AT(G, msg));
+    OK(LAGraph_Cached_EMin(G, msg));
+
+    //begin test
+    double flow = 0;
+    OK(LAGr_MaxFlow(&flow, &flow_mtx, G, tests[test].S, tests[test].T, msg));
+    int status = LG_check_flow(&flow_mtx, msg);
+    printf("%d", status);
+    printf("%s\n", msg);
+    TEST_CHECK(flow == tests[test].F);
+    printf("flow is: %lf\n", flow);
+
+    //free work
+    GrB_free(&flow_mtx);
+    OK(LAGraph_Delete(&G, msg));
+  }
+  LAGraph_Finalize(msg);
+}
+
+
+TEST_LIST = {{"MaxFlow", test_MaxFlow}, {"MaxFlowMtx", test_MaxFlowMtx}, {NULL, NULL}};

include/LAGraphX.h

@@ -1492,6 +1492,7 @@ int LAGraph_coloring_MIS
     char *msg
 ) ;
 
+LAGRAPHX_PUBLIC
 int LAGr_MaxFlow(
     //outputs
     double* f,