<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>

</head>

<body dir="ltr">

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Good afternoon,

<div><br>

</div>

<div>Thanks for all your help.</div>

<div>I got a<span style="background-color:rgb(255, 255, 255);display:inline !important"><span> </span>question about how to solve the DAE step by step using PETSC4PY.</span></div>

<div>I got two DAE systems, let's say f1(x1, u1) and f2(x2, u2). The simplified algorithm I need to implement is as follows:</div>

<div>Step 1: solve f1 for 1 step.</div>

<div>Step 2: use part of x1 as input for f2, which means u2 = x1[:n].</div>

<div>Step 3: solve f2 for one step with u2 = x1[:n].</div>

<div>Step 4: use part of x2 as input for f1, which means u1 = x2[:m].</div>

<div><br>

</div>

<div>I'm not able to find any examples of how to use PETSC4PY in such scenarios. If using the "scikits.odes.dae" package, it is like:</div>

<blockquote style="margin-top:0;margin-bottom:0">

<div>daesolver.init_step(timeInit, XInit, XpInit)</div>

<div>daesolver.step(time)</div>

</blockquote>

<div><br>

</div>

<div>Please let me know if there are any examples I can refer to. Thank you.</div>

<blockquote style="margin-top:0;margin-bottom:0">

<div><br>

</div>

</blockquote>

<span></span>Best,</div>

<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">

Jing</div>

</body>

</html>