obsolete
full matrix multilpication
#rhs from linear term
depfun deprecated, but still way faster, disable warnings