adding a dot multiply for bfam and test
adding += versions of kron products. Also adding y = a.*b.*x.
Some more convenience macros
tmp: derivatives terms for velocity