İşlem sırası, sayısal algoritmalar

Bunu okudum

(1) Kötü şartlandırılmış operasyonlar iyi şartlandırılmış operasyonlardan önce yapılmalıdır.

Örnek olarak, çarpma işlemi yapılmadığı sürece çıkarma kötü koşullandırıldığı için $xz-yz$ olarak hesaplanmalıdır . $(x-y)z$

Ancak, her iki algoritmanın birinci dereceden hata analizi, sadece üç (*) faktörü ile farklılık gösterdiklerini ortaya çıkarır ve neden bunu (1) ifadesine genelleştirebildiğini anlamıyorum, ne de sezgisel olarak operasyonların sırası. İfadenin (1) kabul edilmiş bir kural olduğunu düşünüyor musunuz ve bunun için başka açıklamalarınız var mı?

*: daha spesifik olarak, ilk sürümün sınırlı bir bağıl hatası var

eps + 3 \frac{| x | + | y |}{| x | - | y |} eps

$\text{eps}+3\frac{|x|+|y|}{|x|-|y|}\text{eps}$ ise ikinci versiyonun göreceli hatası

3 eps + \frac{| x | + | y |}{| x | - | y |} eps

$3\text{eps}+\frac{|x|+|y|}{|x|-|y|}\text{eps}$

Burada $\text{eps}$ makine hassasiyetidir.

Bu analiz varsayımına dayanmaktadır $i$ -inci ara sonuç ile çarpılan $(1+\varepsilon_i)$ (yuvarlama nedeniyle hatalara), $\varepsilon_i$ tarafından sınırlanan iid rasgele değişkenlerdir $\text{eps}$ . "Birinci dereceden", gibi yüksek dereceli terimlerin $\epsilon_i\epsilon_jx$ ihmal edildiği anlamına gelir .

stability floating-point

— Bananach
kaynak

Bunu nereden okudun?

— David Ketcheson

ders notlarımda

— Bananach

Tarafından edelim anlamında olabildikleri kesin çarpım kayan nokta analoglarını ((ben bölünme operatörünün daire içine alınmış sürümünü almaya çalışırken tembel) ), ekleme ( () ve çıkarma sırasıyla). (IEEE-754) hepsine göre $\otimes,\oplus,\ominus$ $\times$ $+$ $-$ burada bir üst yüzünden yuvarlama göre hata bağlanmış veren makinesi epsilondur. Ayrıcakolayca kanıtlanabilecekaşağıdaki lemmayı da kullanacağız (tüm ve çok büyük değil),

[x \oplus y] = (x + y) (1 + δ_{\oplus}), | δ_{\oplus} | \leq ε_{m bir c h},

$[x\oplus y]=(x+ y)(1+\delta_\oplus),\quad |\delta_\oplus|\le\epsilon_\mathrm{mach},$

ϵ_{m a c h}

$\epsilon_\mathrm{mach}$

| δ_{i} | \leq ϵ_{m a c h}

$|\delta_i|\le\epsilon_\mathrm{mach}$

m

$m$

Π_{ben = 1}^{m} (1 + δ_{ben}) = 1 + θ (m), | θ (m) | \leq \frac{m ε_{m bir c h}}{1 - m ε_{m bir c h}}

$\prod\limits_{i=1}^{m}(1+\delta_i)=1+\theta(m),\quad |\theta(m)|\le\frac{m\epsilon_\mathrm{mach}}{1-m\epsilon_\mathrm{mach}}$

Let gerçek işlevi tanımlayan reel sayılar çalışır olarak $f$ $x,y,z$

f (x, y, z) = (x x z) - (y x z)

$f(x,y,z)=(x\times z)-(y\times z)$

ve işlev olarak, IEEE uyumlu kayan nokta aritmetik uygulanması iki sürümü ve kayan nokta üzerinde işlem temsilleri , aşağıdaki gibidir: $\tilde{f_1}$ $\tilde{f_2}$ $\tilde{x}=x(1+\delta_x),\tilde{y},\tilde{z}$

\tilde{f_{1}} (\tilde{x}, \tilde{y}, \tilde{z}) = (\tilde{x} \otimes \tilde{z}) ⊖ (\tilde{y} \otimes \tilde{z}),

$\tilde{f_1}(\tilde{x},\tilde{y},\tilde{z})=(\tilde{x}\otimes\tilde{z})\ominus(\tilde{y}\otimes\tilde{z}),$

\tilde{f_{2}} (\tilde{x}, \tilde{y}, \tilde{z}) = (\tilde{x} ⊖ \tilde{y}) \otimes \tilde{z} .

$\tilde{f_2}(\tilde{x},\tilde{y},\tilde{z})=(\tilde{x}\ominus\tilde{y})\otimes\tilde{z}.$

için hata analizi : $\tilde{f_1}$

\begin{aligned} \tilde{f_{1}} & = (\underset{(\tilde{x} \otimes \tilde{z})}{\underset{⏟}{(x (1 + δ_{x}) x z (1 + δ_{z})) (1 + δ_{\otimes_{x z}})}} - \underset{(\tilde{y} \otimes \tilde{z})}{\underset{⏟}{(y (1 + δ_{y}) x z (1 + δ_{z})) (1 + δ_{\otimes_{y z}})}}) (1 + δ_{⊖}) \\ = x z (1 + δ_{x}) (1 + δ_{z}) (1 + δ_{\otimes_{x z}}) (1 + δ_{⊖}) - y z (1 + δ_{y}) (1 + δ_{z}) (1 + δ_{\otimes_{y z}}) (1 + δ_{⊖}) \\ = x z (1 + θ_{x z, 1}) - y z (1 + θ_{y z, 1}) . \end{aligned}

$\begin{aligned} \tilde{f_1}&=\Big(\underbrace{\big(x(1+\delta_x)\times z(1+\delta_z)\big)\big(1+\delta_{\otimes_{xz}}\big)}_{(\tilde{x}\otimes\tilde{z})}-\underbrace{\big(y(1+\delta_y)\times z(1+\delta_z)\big)\big(1+\delta_{\otimes_{yz}}\big)}_{(\tilde{y}\otimes\tilde{z})}\Big)\Big(1+\delta_{\ominus}\Big)\\ &=xz(1+\delta_x)(1+\delta_z)(1+\delta_{\otimes_{xz}})(1+\delta_\ominus)-yz(1+\delta_y)(1+\delta_z)(1+\delta_{\otimes_{yz}})(1+\delta_\ominus)\\ &=xz(1+\theta_{xz,1})-yz(1+\theta_{yz,1}). \end{aligned}$

| θ_{x z, 1} |, | θ_{y z, 1} | \leq \frac{4 ϵ_{m a c h}}{1 - 4 ϵ_{m a c h}}

$|\theta_{xz,1}|,|\theta_{yz,1}|\le\frac{4\epsilon_\mathrm{mach}}{1-4\epsilon_\mathrm{mach}}$

$\tilde{f_2}$

\begin{aligned} \tilde{f_{2}} & = (((x (1 + δ_{x}) - y (1 + δ_{y}) (1 + δ_{⊖_{x y}})) x (z (1 + δ_{z}))) (1 + δ_{\otimes}) \\ = x z (1 + δ_{x}) (1 + δ_{z}) (1 + δ_{⊖_{x y}}) (1 + δ_{\otimes}) - y z (1 + δ_{y}) (1 + δ_{z}) (1 + δ_{⊖_{x y}}) (1 + δ_{\otimes}) \\ = x z (1 + θ_{x, 2}) - y z (1 + θ_{y, 2}) . \end{aligned}

$\begin{aligned} \tilde{f_2}&=\Bigg(\Big(\big( x(1+\delta_x)-y(1+\delta_y \big)\big(1+\delta_{\ominus_{xy}}\big)\Big)\times \Big(z(1+\delta_z)\Big)\Bigg)\Bigg(1+\delta_{\otimes}\Bigg)\\ &=xz(1+\delta_x)(1+\delta_z)(1+\delta_{\ominus_{xy}})(1+\delta_\otimes)-yz(1+\delta_y)(1+\delta_z)(1+\delta_{\ominus_{xy}})(1+\delta_\otimes)\\ &=xz(1+\theta_{x,2})-yz(1+\theta_{y,2}). \end{aligned}$

| θ_{x, 2} |, | θ_{y, 2} | \leq \frac{4 ϵ_{m a c h}}{1 - 4 ϵ_{m a c h}}

$|\theta_{x,2}|,|\theta_{y,2}|\le\frac{4\epsilon_\mathrm{mach}}{1-4\epsilon_\mathrm{mach}}$

$\tilde{f_1}$ $\tilde{f_2}$ $\tilde{f_2}$ $\tilde{f_1}$

$x$ $y$

\begin{aligned} \frac{| \tilde{f_{1}} - f |}{| f |} & = \frac{| x z + x z θ_{x z, 1} - y z - y z θ_{y z, 1} - (x z - y z) |}{| x z - y z |} = \frac{| x θ_{x z, 1} - y θ_{y z, 1} |}{| x - y |} \\ \leq \frac{| x | + | y |}{| x - y |} \frac{4 ε_{m bir c h}}{1 - 4 ε_{m bir c h}}, \end{aligned}

$\begin{aligned} \frac{|\tilde{f_1}-f|}{|f|}&=\frac{|xz+xz\theta_{xz,1}-yz-yz\theta_{yz,1}-(xz-yz)|}{|xz-yz|}=\frac{|x\theta_{xz,1}-y\theta_{yz,1}|}{|x-y|}\\ &\le\frac{|x|+|y|}{|x-y|}\frac{4\epsilon_\mathrm{mach}}{1-4\epsilon_\mathrm{mach}}, \end{aligned}$

\begin{aligned} \frac{| \tilde{f_{2}} - f |}{| f |} & = \frac{| x z + x z θ_{x, 2} - y z - y z θ_{y, 2} - (x z - y z) |}{| x z - y z |} = \frac{| x θ_{x, 2} - y θ_{y, 2} |}{| x - y |} \\ \leq \frac{| x | + | y |}{| x - y |} \frac{4 ε_{m bir c h}}{1 - 4 ε_{m bir c h}} . \end{aligned}

$\begin{aligned} \frac{|\tilde{f_2}-f|}{|f|}&=\frac{|xz+xz\theta_{x,2}-yz-yz\theta_{y,2}-(xz-yz)|}{|xz-yz|}=\frac{|x\theta_{x,2}-y\theta_{y,2}|}{|x-y|}\\ &\le\frac{|x|+|y|}{|x-y|}\frac{4\epsilon_\mathrm{mach}}{1-4\epsilon_\mathrm{mach}}. \end{aligned}$

$\theta$ $x,y,z$ $(x-y)$ $x$ $y$

$x,y,z,f(x,y,z)\in\mathbb F_0$ $\mathbb F_0$

— Anton Menshov
kaynak