Represent scalar equations as Vectorised calls

.github/workflows/TagBot.yml

@@ -1,9 +1,11 @@
 name: TagBot
 on:
-  schedule:
-    - cron: 0 * * * *
+  issue_comment:
+    types:
+      - created
 jobs:
   TagBot:
+    if: github.event_name == 'workflow_dispatch' || github.actor == 'JuliaTagBot'
     runs-on: ubuntu-latest
     steps:
       - uses: JuliaRegistries/TagBot@v1

Project.toml

@@ -7,8 +7,10 @@ version = "0.1.2"
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 PreallocationTools = "d236fae5-4411-538c-8e31-a6e3d9e00b46"
+SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 SymbolicUtils = "d1185830-fcd6-423d-90d6-eec64667417b"
+Symbolics = "0c5d862f-8b57-4792-8d23-62f2024744c7"
 
 [weakdeps]
 LinearSolve = "7ed4a6bd-45f5-4d41-b270-4a48e9bafcae"
@@ -24,6 +26,7 @@ LinearAlgebra = "1.11.0, 1.10"
 LinearSolve = "3.53.0"
 PreallocationTools = "0.4.34, 1"
 Rotations = "1.7.1"
+SparseArrays = "1"
 StaticArrays = "1.9.15, 1"
 SymbolicUtils = "4.1.0"
 Symbolics = "7.2"

src/SymbolicCompilerPasses.jl

@@ -9,7 +9,9 @@ import SymbolicUtils.Code: Code, OptimizationRule, substitute_in_ir, apply_optim
     Assignment, CSEState, lhs, rhs, apply_substitution_map, issym, isterm, toexpr,
     _is_array_of_symbolics, MakeArray, shape
 import SymbolicUtils: search_variables, search_variables!
+import SymbolicUtils as SU
 using StaticArrays
+using Symbolics
 
 using DataStructures
 
@@ -28,5 +30,6 @@ include("ldiv_opt.jl")
 include("la_opt.jl")
 
 include("mb_opt.jl")
+include("jacobian_vectorize.jl")
 
 end # module SymbolicCompilerPasses