On 6 July 2017 at 17:10, Alejandro D Otero wrote:
>
> # Set rhs function
> ts.setRHSFunction(self.evalRHSFunction, self.g)
> (where evalRHSFunction has the form f(self, ts, t, Vort, f))
>
Please double check that the RHS function "outputs" the result in the
passed-in vector
Hello,
I'trying to solve an ODE system using PETSc TS through its wrapper
petsc4py, and I'm having problems with the setting. The problem is of the
kind: u_t = G(t, u) with u being a vector with DoF from a FEM problem and G
is computed by a series of matrix operations and manipulations, which are