If mfdobj2 is missing, the function returns mfdobj1 unchanged
(i.e., unary plus).
Let the coefficient arrays have dimensions \((nbasis, nobs, nvar)\).
The following checks/rules are enforced:
Both inputs must be mfd objects; otherwise an error is thrown.
The basis systems must be identical (checked via identical()).
The number of variables must match.
For the number of 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.