Hello,

i would like to set up an educational example where i assume not knowing that for magnetostatics one should use H(Curl) elements. Also neglecting all the stuff about gauging and stabilization.

I used the cmagnet.py as reference and did the following changes:

```
fes = H1(mesh, order=2, dirichlet="outer")
a+= SymbolicBFI(nu*div(u)*div(v))
```

NGsolve tells me now that it cannot form div.

Isn’t that the same structure as for linear elasticity problems?

Am i missing a mathematical trap or am i treating the program wrong?

Thank you for that great software and your strong support!

With best regards