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