[petsc-users] is TS_EQ_DAE_SEMI_EXPLICIT_INDEX functional

2019-09-02 Thread Huck, Moritz via petsc-users
Hi, TS_EQ_DAE_SEMI_EXPLICIT_INDEX(?) are defined in TSEquationType but not mentioned in the manual. Is this feature functional ? If yes how do I have to define the RHSFunction? (I am asking since the ODE variant has it defined as G= M^-1 g, which cannot work for a DAE) Best Regards, Moritz

Re: [petsc-users] is TS_EQ_DAE_SEMI_EXPLICIT_INDEX functional

2019-09-02 Thread Jed Brown via petsc-users
I believe this is intended to work with most any implicit solver, *provided* the initial conditions are compatible. It was added by Emil, but I don't see it explicitly tested in PETSc. "Huck, Moritz via petsc-users" writes: > Hi, > TS_EQ_DAE_SEMI_EXPLICIT_INDEX(?) are defined in TSEquationType

Re: [petsc-users] is TS_EQ_DAE_SEMI_EXPLICIT_INDEX functional

2019-09-02 Thread Constantinescu, Emil M. via petsc-users
Indeed, various time steppers can take advantage of the differential form provided and also can serve as a sanity check (e.g., warn users before they use an explicit solver on an index-2 DAE). To my knowledge, we do not have solvers that take advantage of semi-explicit DAEs, but it's good pract

Re: [petsc-users] is TS_EQ_DAE_SEMI_EXPLICIT_INDEX functional

2019-09-05 Thread Zhang, Hong via petsc-users
You do not need to worry about these equation flags in the beginning. To solve a DAE, you need to define IFunction and IJacobian instead of RHSFunction. ts/examples/tutorials/ex19.c is the simplest DAE example that you can refer to. Hong (Mr.) > On Sep 2, 2019, at 4:51 AM, Huck, Moritz via pets