add multiplication operations

Author: abhinavmehndiratta

src/Operations/Multiplication.jl

@@ -0,0 +1,60 @@
+
+function GrB_mxm(              # C<Mask> = accum (C, A*B)
+    C::GrB_Matrix,             # input/output matrix for results
+    Mask::GrB_Matrix,          # optional mask for C, unused if NULL
+    accum::GrB_BinaryOp,       # optional accum for Z=accum(C,T)
+    semiring::GrB_Semiring,    # defines '+' and '*' for A*B
+    A::GrB_Matrix,             # first input:  matrix A
+    B::GrB_Matrix,             # second input: matrix B
+    desc::GrB_Descriptor       # descriptor for C, Mask, A, and B
+)
+
+    return GrB_Info(
+                ccall(
+                        dlsym(graphblas_lib, "GrB_mxm"),
+                        Cint,
+                        (Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}),
+                        C.p, Mask.p, accum.p, semiring.p, A.p, B.p, desc.p
+                    )
+                )
+end
+
+function GrB_vxm(              # w'<Mask> = accum (w, u'*A)
+    w::GrB_Vector,             # input/output vector for results
+    mask::GrB_Vector,          # optional mask for w, unused if NULL
+    accum::GrB_BinaryOp,       # optional accum for z=accum(w,t)
+    semiring::GrB_Semiring,    # defines '+' and '*' for u'*A
+    u::GrB_Vector,             # first input:  vector u
+    A::GrB_Matrix,             # second input: matrix A
+    desc::GrB_Descriptor       # descriptor for w, mask, and A
+)
+
+    return GrB_Info(
+                ccall(
+                        dlsym(graphblas_lib, "GrB_vxm"),
+                        Cint,
+                        (Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}),
+                        w.p, mask.p, accum.p, semiring.p, u.p, A.p, desc.p
+                    )
+                )
+end
+
+function GrB_mxv(               # w<Mask> = accum (w, A*u)
+    w::GrB_Vector,              # input/output vector for results
+    mask::GrB_Vector,           # optional mask for w, unused if NULL
+    accum::GrB_BinaryOp,        # optional accum for z=accum(w,t)
+    semiring::GrB_Semiring,     # defines '+' and '*' for A*B
+    A::GrB_Matrix,              # first input:  matrix A
+    u::GrB_Vector,              # second input: vector u
+    desc::GrB_Descriptor        # descriptor for w, mask, and A
+)
+
+    return GrB_Info(
+                ccall(
+                        dlsym(graphblas_lib, "GrB_mxv"),
+                        Cint,
+                        (Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}, Ptr{Cvoid}),
+                        w.p, mask.p, accum.p, semiring.p, A.p, u.p, desc.p
+                    )
+                )
+end

src/SuiteSparseGraphBLAS.jl

@@ -138,6 +138,7 @@ include("Object_Methods/Vector_Methods.jl")
 include("Object_Methods/Algebra_Methods.jl")
 include("Object_Methods/Descriptor_Methods.jl")
 include("Object_Methods/Print_Objects.jl")
+include("Operations/Multiplication.jl")
 
 export
 # Context Methods
@@ -160,7 +161,10 @@ GrB_UnaryOp_new, GrB_BinaryOp_new, GrB_Monoid_new, GrB_Semiring_new,
 
 # Print functions
 @GxB_UnaryOp_fprint, @GxB_BinaryOp_fprint, @GxB_Monoid_fprint, @GxB_Semiring_fprint,
-@GxB_Matrix_fprint, @GxB_Vector_fprint, @GxB_Descriptor_fprint, @GxB_fprint
+@GxB_Matrix_fprint, @GxB_Vector_fprint, @GxB_Descriptor_fprint, @GxB_fprint,
+
+# GraphBLAS Operations
+GrB_mxm, GrB_vxm, GrB_mxv
 
 # Export global variables