Convenience wrapper: builds the full canonical design via
.build_full_design and contracts it against a named
coefficient vector. Because the design is rebuilt in the model's slice
order, the positional contraction is correct even with nodal covariates.
.xbeta_newdata(beta, newdata, fitted_design, n, m)an n x m matrix of the linear predictor.
named numeric coefficient vector (a single draw or the posterior mean).
dyadic covariate array (see .build_full_design).
fitted design array (fit$X); held-fixed source
for intercept + nodal slices.
output dimensions.