Add implmentation of elementwise bounds of full matrices

src/Gradients.jl

@@ -284,6 +284,9 @@ function compute_standard_form(sys :: StaticSystem, ps; margin, debug = false, t
                             end
                             Abc.b[S_ixs[(k, :bound, bound_ix)]] += sign*value
                         end
+                        (value, Elementwise(), Full(nr,nc)) => begin
+                            [ set_bound(sign, k, bound_ix, value, r, c) for r in 1:nr for c in 1:nc ]
+                        end
                         (_, _, Full(nr, nc)) => @error "Full matrix not yet implemented"
                             
                         (_, _, _) => begin