<div dir="ltr">Hi everyone,<div><br></div><div>I'm trying to use a timestepper to solve a nonlinear system with Petsc4py. This system includes a grid variables which I can manage using a DMDA, but it also includes a small number of globally coupled scalar values. This seems like what DMREDUNDANT was designed for: make DMDA and DMREDUNDANT objects, stick them into a composite, and use the global vector generated by that composite DM as the state variable in my nonlinear solver.</div><div><br></div><div>As best I can tell, Petsc4py doesn't have a method like DM().createRedundant()--perhaps it this DM just isn't implemented? Can anyone recommend an alternate approach to getting the same functionality?</div><div><br></div><div>Thanks,</div><div>Ben</div></div>