119     def __init__(self, ode):  argument
120         self.ode_ = ode
132     def __init__(self, ode):  argument
133         self.ode_ = ode
141     def __init__(self, ode):  argument
142         self.ode_ = ode
154     def __init__(self, ode):  argument
155         self.ode_ = ode
170 ode = VDP(mu_, mf_, imexform_)  variable
173     Jim = PETSc.Mat().createDense([ode.n, ode.n], comm=ode.comm)
175     JimP = PETSc.Mat().createDense([ode.n, 1], comm=ode.comm)
177     Jex = PETSc.Mat().createDense([ode.n, ode.n], comm=ode.comm)
179     JexP = PETSc.Mat().createDense([ode.n, 1], comm=ode.comm)
183     Jim.setSizes([ode.n, ode.n])
185     shell = IJacShell(ode)
190     JimP.setSizes([ode.n, 1])
192     shell = IJacPShell(ode)
197     Jex.setSizes([ode.n, ode.n])
199     shell = JacShell(ode)
204     JexP.setSizes([ode.n, 1])
206     shell = JacPShell(ode)
212 u = PETSc.Vec().createSeq(ode.n, comm=ode.comm)
215 adj_u.append(PETSc.Vec().createSeq(ode.n, comm=ode.comm))
216 adj_u.append(PETSc.Vec().createSeq(ode.n, comm=ode.comm))
218 adj_p.append(PETSc.Vec().createSeq(1, comm=ode.comm))
219 adj_p.append(PETSc.Vec().createSeq(1, comm=ode.comm))
221 ts = PETSc.TS().create(comm=ode.comm)
226     ts.setIFunction(ode.evalIFunction, f)
227     ts.setIJacobian(ode.evalIJacobian, Jim)
228     ts.setIJacobianP(ode.evalIJacobianP, JimP)
229     ts.setRHSFunction(ode.evalFunction, f)
230     ts.setRHSJacobian(ode.evalJacobian, Jex)
231     ts.setRHSJacobianP(ode.evalJacobianP, JexP)
235         ts.setIFunction(ode.evalIFunction, f)
236         ts.setIJacobian(ode.evalIJacobian, Jim)
237         ts.setIJacobianP(ode.evalIJacobianP, JimP)
240         ts.setRHSFunction(ode.evalFunction, f)
241         ts.setRHSJacobian(ode.evalJacobian, Jex)
242         ts.setRHSJacobianP(ode.evalJacobianP, JexP)
252 ode.initialCondition(u)
286 del ode, Jim, JimP, Jex, JexP, u, f, ts, adj_u, adj_p