Hello,
is there any possibility to add a constant Matrix to the Bilinear Form?
Code:
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!