diff --git a/integration/VODE/vode_dvjac.H b/integration/VODE/vode_dvjac.H
index 848a31621..8c420e6f1 100644
--- a/integration/VODE/vode_dvjac.H
+++ b/integration/VODE/vode_dvjac.H
@@ -79,9 +79,6 @@ void dvjac (int& IERPJ, BurnT& state, DvodeT& vstate)
             // Indicate that the Jacobian is current for this solve.
             vstate.JCUR = 1;
 
-            // Initialize the Jacobian to zero
-            vstate.jac.zero();
-
             jac(vstate.tn, state, vstate, vstate.jac);
 
 #ifdef ALLOW_JACOBIAN_CACHING
diff --git a/sphinx_docs/source/integrators.rst b/sphinx_docs/source/integrators.rst
index 1e7b5fe54..cf6c87565 100644
--- a/sphinx_docs/source/integrators.rst
+++ b/sphinx_docs/source/integrators.rst
@@ -276,6 +276,11 @@ The form looks like:
    A network is not required to provide a Jacobian if a numerical
    Jacobian is used.
 
+.. important::
+
+   The integrator does not zero the Jacobian elements.  It is the responsibility
+   of the Jacobian implementation to zero the Jacobian array if necessary.
+
 
 Jacobian wrapper
 ^^^^^^^^^^^^^^^^