Alternatively, it also compute the product of an mfd object with a numeric scalar.
Let coefficient arrays have dimensions \((nbasis, nobs, nvar)\).
The function:
requires both inputs to be mfd objects;
requires identical basis systems (checked with identical());
requires the same number of variables;
for observations: if both \(nobs_1\) and \(nobs_2\) are greater
than one, they must be equal; otherwise, the object with
\(nobs = 1\) is replicated to match the other.
Internally, coefficient arrays are converted to fd objects
and multiplied via times.fd, with basisobj set to
the common basis so that the result is re-expanded on the same basis.