# Multi-output BBL mass constraint

- ID: 117
- Canonical URL: https://synthcore.org/idea/117/multi-output-bbl-mass-constraint
- API JSON: https://synthcore.org/api/idea/117.json
- API Markdown: https://synthcore.org/api/idea/117.md
- Verification status: mechanism_failed
- Source: [arXiv:2608.23963](https://arxiv.org/abs/2608.23963)
- Category: regularization
- Solves: stability, generalization
- ML areas: loss, regularization, attention
- Math tags: inequalities, convex-analysis, measure-theory
- Ratings: usefulness 6/10; difficulty 5/10; novelty 8/10

## Idea description

Represent each of m neural branches by a positive input field f_i and a positive output field g_i, then penalize violations of the paper's multi-output Borell-Brascamp-Lieb bound at weighted barycenters. The constraint couples branches through both local normalized ratios and global mass ratios, encouraging calibrated multi-view predictions without requiring all output functions to be identical.

## Mathematical statement

Let lambda_i>0 with sum_i lambda_i=1, and define the barycenter z_lambda(x)=sum_i lambda_i x_i. For p not equal to zero, the weighted power mean is M_p^lambda(a_1,...,a_m)=(sum_i lambda_i a_i^p)^(1/p); for p=0 it is M_0^lambda(a)=product_i a_i^(lambda_i), and for p=-infinity it is min_i a_i. The main inequality states that for nonnegative integrable functions f_i,g_i on R^d with 0<||f_i||_1<infinity and p<=1/d, ess inf_x M_p^lambda(g_1(z_lambda(x))/f_1(x_1),...,g_m(z_lambda(x))/f_m(x_m)) is at most M_{Q_d(p)}^lambda(||g_1||_1/||f_1||_1,...,||g_m||_1/||f_m||_1), where Q_d(p)=p/(1-dp). Define r_i^(n)=g_i(z^(n))/f_i(x_i^(n)) and A_i=||g_i||_1/||f_i||_1. Replace the essential infimum by the soft minimum smin_tau(u_1,...,u_N)=-tau log sum_n exp(-u_n/tau). The practical penalty is max(0, smin_tau(M_p^lambda(r^(n)))-M_{Q_d(p)}^lambda(A))^2. Positivity and nonzero mass must be enforced separately because otherwise g_i=0 is a trivial solution.

## Key formulas

- $$z_{\lambda}(x)=\sum_{i=1}^{m}\lambda_i x_i,\qquad \lambda_i>0,\qquad \sum_{i=1}^{m}\lambda_i=1.$$
- $$M_p^{\lambda}(a_1,\ldots,a_m)=\left(\sum_{i=1}^{m}\lambda_i a_i^p\right)^{1/p}\quad(p\neq 0),\qquad M_0^{\lambda}(a_1,\ldots,a_m)=\prod_{i=1}^{m}a_i^{\lambda_i}.$$
- $$\operatorname*{ess\,inf}_{x}M_p^{\lambda}\left(\frac{g_1(z_{\lambda}(x))}{f_1(x_1)},\ldots,\frac{g_m(z_{\lambda}(x))}{f_m(x_m)}\right)\le M_{Q_d(p)}^{\lambda}\left(\frac{\|g_1\|_1}{\|f_1\|_1},\ldots,\frac{\|g_m\|_1}{\|f_m\|_1}\right),\qquad Q_d(p)=\frac{p}{1-dp}.$$
- $$\mathcal L_{\mathrm{BBL}}=\left[\max\left\{0,\operatorname{smin}_{\tau,n}\left(M_p^{\lambda}(r_1^{(n)},\ldots,r_m^{(n)})\right)-M_{Q_d(p)}^{\lambda}(A_1,\ldots,A_m)\right\}\right]^2.$$

## Implementation notes

Integrate the constraint after the final positive prediction heads of a multi-view or multi-sensor model. For branch i, produce f_i(u)=softplus(F_i(u))+epsilon from its input representation and g_i(v)=softplus(G_i(v))+epsilon from its output representation, where epsilon is 1e-6. Sample N tuples (x_1^(n),...,x_m^(n)) from the branch coordinate grids, set z^(n)=sum_i lambda_i x_i^(n), and evaluate r_i^(n)=g_i(z^(n))/f_i(x_i^(n)). Estimate masses with Monte Carlo: ||f_i||_1 approximately equals B^(-1) sum_b f_i(x_i^b), and ||g_i||_1 approximately equals B^(-1) sum_b g_i(v_i^b). Compute q=Q_d(p)=p/(1-dp), then evaluate the right-hand power mean. Use logarithmic formulas for p=0 and clamp values before negative powers. Since the theorem uses an essential infimum, use smin_tau(u)=-tau log sum_n exp(-u_n/tau), with tau initially 0.2 times the median u and annealed toward 0.02 times the median. Add eta times L_BBL to the task loss. Prevent collapse by normalizing each predicted field to a fixed mass, or by adding a mass-matching penalty against a target estimated from labels. Use a detached sample pool or periodically refreshed hard-min samples to prevent the network from exploiting favorable samples. The theorem determines the inequality, valid p range, and exponent q; lambda, p, eta, tau, and sampling are empirical choices. First test a two-view MNIST or CIFAR-10 classifier with two positive confidence heads, comparing ordinary training against the penalty at equal FLOPs. Measure clean accuracy, accuracy under one-view corruption, calibration error, cross-view disagreement, and the lower quantile of normalized ratios. Success is improved corruption calibration or lower disagreement at unchanged clean accuracy without predicted-mass collapse.

## Verification

- Status: mechanism_failed
- Mechanism evidence: yes
- Mechanism confirmed: no
- Verdict: Built bbl_experiment.py with stable weighted power means, Q_d(p), soft-min BBL penalty, positive heads, Monte Carlo mass ratios, Gaussian equality verification, and a fixed-seed two-view classifier comparison. The Gaussian sanity check reproduced the claimed equality (minimum ratio 1.0 and RHS 1.0), but the practical penalty had a negative gap on every training step, zero penalty, and zero gradient. Baseline and idea were bit-identical: 99.50% clean accuracy, 88.06% corrupted accuracy, ECE 0.00518, and disagreement 0.0991; therefore no observed regularization or robustness win.

### Mechanism check

- Verdict: Built bbl_experiment.py with stable weighted power means, Q_d(p), soft-min BBL penalty, positive heads, Monte Carlo mass ratios, Gaussian equality verification, and a fixed-seed two-view classifier comparison. The Gaussian sanity check reproduced the claimed equality (minimum ratio 1.0 and RHS 1.0), but the practical penalty had a negative gap on every training step, zero penalty, and zero gradient. Baseline and idea were bit-identical: 99.50% clean accuracy, 88.06% corrupted accuracy, ECE 0.00518, and disagreement 0.0991; therefore no observed regularization or robustness win.
- Confidence: 9/10
- Limitations: Only a small synthetic two-view binary classifier was tested, with one p=0 setting, fixed equal weights, approximate shared-sample mass estimates, and a single seed. No MNIST/CIFAR benchmark, hyperparameter sweep, multi-branch study, or comparison of alternative mass-normalization schemes was run.

## Artifacts

- [bbl_experiment.py](https://synthcore.org/code/1181/bbl_experiment.py)
- [report.md](https://synthcore.org/code/1181/report.md)
- [results.json](https://synthcore.org/code/1181/results.json)
- [Download all files as ZIP](https://synthcore.org/download/1181)

## Disclaimer

AI-generated research hypothesis, automatically tested. Not peer-reviewed.
