Add vecjac operators

Author: Avik Pal

Project.toml

@@ -16,6 +16,7 @@ Requires = "ae029012-a4dd-5104-9daa-d747884805df"
 SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 VertexSafeGraphs = "19fa3120-7c27-5ec5-8db8-b0b0aa330d6f"
+Zygote = "e88e6eb3-aa80-5325-afca-941959d7151f"
 
 [compat]
 Adapt = "1, 2.0, 3.0"

src/SparseDiffTools.jl

@@ -8,6 +8,7 @@ using LightGraphs: SimpleGraph
 using Requires
 using VertexSafeGraphs
 using Adapt
+using Zygote
 
 using LinearAlgebra
 using SparseArrays, ArrayInterface
@@ -30,6 +31,8 @@ export  contract_color,
         ForwardColorJacCache,
         auto_jacvec,auto_jacvec!,
         num_jacvec,num_jacvec!,
+        auto_vecjac,auto_vecjac!,
+        num_vecjac,num_vecjac!,
         num_hesvec,num_hesvec!,
         numauto_hesvec,numauto_hesvec!,
         autonum_hesvec,autonum_hesvec!,
@@ -48,6 +51,7 @@ include("coloring/greedy_star2_coloring.jl")
 include("coloring/matrix2graph.jl")
 include("differentiation/compute_jacobian_ad.jl")
 include("differentiation/jaches_products.jl")
+include("differentiation/vecjac_products.jl")
 
 Base.@pure __parameterless_type(T) = Base.typename(T).wrapper
 parameterless_type(x) = parameterless_type(typeof(x))