  We test the distance--duality relation $\eta \equiv d_L / [ (1 + z)^2 d_A ] =
1$ between cosmological luminosity distance ($d_L$) from the JLA SNe Ia
compilation (arXiv:1401.4064) and angular-diameter distance ($d_A$) based on
Baryon Oscillation Spectroscopic Survey (BOSS; arXiv:1607.03155) and WiggleZ
baryon acoustic oscillation measurements (arXiv:1105.2862, arXiv:1204.3674).
The $d_L$ measurements are matched to $d_A$ redshift by a statistically
consistent compression procedure. With Monte Carlo methods, nontrivial and
correlated distributions of $\eta$ can be explored in a straightforward manner
without resorting to a particular evolution template $\eta(z)$. Assuming
independent constraints on cosmological parameters that are necessary to obtain
$d_L$ and $d_A$ values, we find 9% constraints consistent with $\eta = 1$ from
the analysis of SNIa + BOSS and an 18% bound results from SNIa + WiggleZ. These
results are contrary to previous claims that $\eta < 1$ has been found close to
or above the $1 \sigma$ level. We discuss the effect of different cosmological
parameter inputs and the use of the apparent deviation from distance--duality
as a proxy of systematic effects on cosmic distance measurements. The results
suggest possible systematic overestimation of SNIa luminosity distances
compared with $d_A$ data when a Planck {\Lambda}CDM cosmological parameter
inference (arXiv:1502.01589) is used to enhance the precision. If interpreted
as an extinction correction due to a gray dust component, the effect is broadly
consistent with independent observational constraints.
