exTS_odejacobian
Date: 2013-07-02

    function
       petsc-3.4.2 2013-07-02 [flg,err] = exTS_odejacobian(ts,time,x,xdot,shift,A,B,ctx)
    %
    %  Example of a nonlinear Jacobian needed by TS
    %  Use identity as approximation for Jacobian
    %
    err = 0;
    flg = PetscMat.SAME_NONZERO_PATTERN;
    for i=1:length(x(:))
      B.SetValues(i,i,1.0);
    end
    err = B.AssemblyBegin(PetscMat.FINAL_ASSEMBLY);
    err = B.AssemblyEnd(PetscMat.FINAL_ASSEMBLY);
    err = A.AssemblyBegin(PetscMat.FINAL_ASSEMBLY);
    err = A.AssemblyEnd(PetscMat.FINAL_ASSEMBLY);