is there any possibility to add a constant Matrix to the Bilinear Form?
a = BilinearForm(fes) a+= SymbolicBFI(weak_K(u,v)) a+= "const. Matrix from Numpy" #either here a.Assemble() a.mat += "const Matrix form Numpy" # or here (but does then a.AssembleLinearization() work later?)
The benefit would be that I can still use a.Apply, a.AssembleLinearization, etc… as I would with only considering the “analytical, continous” weak form of my problem.
Thanks for your help in advance!