So I have a finite element space like this:
fes_V = H1(Om_mesh, order=1, dim=5)
where the mesh is a 3-D cube, but the number of components of the functions in H^1 is 5.
Eventually, I try to do this:
a_form = BilinearForm (fes_V, symmetric=False)
and I get this error:
netgen.libngpy._meshing.NgException: Could not create BilinearForm, space-dimension is 5
either define MAX_SYS_DIM with higher value and recompile,
or set flag 'nonassemble'
I understand what it is saying, but where do I set this parameter? Also, am I going to run into issues when I try to do a solve with multi-grid? I'm trying to solve an elliptic *system*, which has 5 coupled equations, and is "very strongly elliptic".
Any advice would be helpful.