Lawrence,
I added the following line:
rs_solver.setDM(a._dm)
And it gives this error:
Traceback (most recent call last):
File "sphere.py", line 128, in <module>
rs_solver.setDM(a._dm)
AttributeError: 'Form' object has no attribute '_dm'
Thanks,
Justin
On Tue, Nov 29, 2016 at 1:12 AM, Lawrence Mitchell <
[email protected]> wrote:
>
>
> > On 29 Nov 2016, at 02:02, Justin Chang <[email protected]> wrote:
> >
> > Why is this happening? It seems to me VINEWTONRSLS is not understanding
> that I have a two-field formulation.
>
> For the built in solvers, we set up a dmshell on the snes that knows how
> to do field splitting. Because you set your snes up by hand, it isn't
> available. You'll need to call snes.setDM. The mixed finctionspace provides
> the correct dm via a ._dm property I think.
>
> Lawrence
>
