初期値と終了基準

初期値は x0=[x^0y^0s^0τ^0]=[0001],y0=0 x_0= \left[ \begin{matrix} \hat x_0 \\ \hat y_0 \\ \hat s_0 \\ \hat \tau_0 \end{matrix} \right] = \left[ \begin{matrix} 0 \\ 0 \\ 0 \\ 1 \end{matrix} \right] ,\qquad y_0=0 とする。

Pock/Chambolleの一次法の反復において、xkx_kproxGτ\mathbf{prox}^\tau_{G} による射影の後のため必ず xkRn×K×K×R+x_k \in \mathbb{R}^n \times \mathcal{K}^\ast \times \mathcal{K} \times \mathbb{R}_+ となっており、 主・双対の錐の制約条件を満たしている。 この状況のもと、終了基準は参考文献[1]と同等のものを適用する。

τ^k>εzero\hat \tau^k > \varepsilon_\mathrm{zero} の場合

収束判定

Ax^k/τ^k+s^k/τ^kb1+bεaccATy^k/τ^k+c1+cεacccTx^k/τ^k+bTy^k/τ^k1+cTx^k/τ^k+bTy^k/τ^kεacc \begin{array}{rl} & \frac{\|A \hat x_k / \hat\tau_k + \hat s_k / \hat\tau_k - b\|}{1 + \|b\|} \le \varepsilon_\mathrm{acc} \\ \land& \frac{\|A^T \hat y_k / \hat\tau_k + c\|}{1 + \|c\|} \le \varepsilon_\mathrm{acc} \\ \land& \frac{\|c^T \hat x_k / \hat\tau_k + b^T \hat y_k / \hat\tau_k\|}{1 + |c^T \hat x_k / \hat\tau_k| + |b^T \hat y_k / \hat\tau_k|} \le \varepsilon_\mathrm{acc} \end{array} を満たすとき、解 x=x^k/τ^k, y=y^k/τ^kx^\star=\hat x_k / \hat\tau_k,\ y^\star=\hat y_k / \hat\tau_k に収束したと判定し、終了する。

τ^kεzero\hat \tau^k \le \varepsilon_\mathrm{zero} の場合

下限なし判定

cTx^k>εzeroAx^k+s^kccTx^kεinf -c^T \hat x_k > \varepsilon_\mathrm{zero} \quad \land \quad \|A \hat x_k + \hat s_k\|\frac{\|c\|}{-c^T \hat x_k} \le \varepsilon_\mathrm{inf} を満たすとき、主問題の下限なし(双対問題の実行可能な解なし)と判定し、終了する。

実行不可能判定

bTy^k>εzeroATy^kbbTy^kεinf -b^T \hat y_k > \varepsilon_\mathrm{zero} \quad \land \quad \|A^T \hat y_k\|\frac{\|b\|}{-b^T \hat y_k} \le \varepsilon_\mathrm{inf} を満たすとき、主問題の実行可能な解なしと判定し、終了する。

results matching ""

    No results matching ""