From 26df7b041dc591529d5199d04e688a5f6a3d13cf Mon Sep 17 00:00:00 2001 From: Peter Schaefer Date: Thu, 29 Mar 2012 18:01:51 +0200 Subject: [PATCH] [src] eta wird anders berechnet [src] error = phi - phi_h/2 --- doc/doc.pdf | Bin 330452 -> 330452 bytes doc/doc.tex | 2 +- src/A_step.m | 12 ++++++++---- 3 files changed, 9 insertions(+), 5 deletions(-) diff --git a/doc/doc.pdf b/doc/doc.pdf index add9a02eabd437711c882856f64871fb358e9ade..cf94a25ed416e8393a25886ff49677bed02db60e 100644 GIT binary patch delta 117 zcmcaIRpiQ4k%kt=7N!>FEiCh%v0EBi7#NtCOrI^yCFEiCh%u~`@x7+DxjpC!#G%MM{1Z$I&jWsS14p}Co(iIIh+ nk%@(&g^8=HlYy~?fuWhRxvR6Qp|Ojhoq`P^CDXTavC0Af%{3ne diff --git a/doc/doc.tex b/doc/doc.tex index 3d663ab..1029755 100644 --- a/doc/doc.tex +++ b/doc/doc.tex @@ -631,7 +631,7 @@ Zum Plotten (\ref{exmplAA_2DQuad})werden noch folgende Schritte ausgeführt \item Berechne Galerkinlösung $\phi_{l}^{(i)} \in P^0(\T_l^{(i)})$ % \item $\enorm{\phi^{(i)}_{l/2}}$ % \item $\enorm{\phi^{(i)}_l}$ - \item $error = \sqrt{\enorm{\phi}^2 - \enorm{\phi_l^{(i)}}^2}$ + \item $error = \sqrt{\enorm{\phi}^2 - \enorm{\phi_{l}^{(i)}}^2}$ \item $\mu_{l,i} := \norm{\varrho^{1/2}(\phi_{l/2}^{(i)} - \phi_{l}^{(i)} )}$ \item $\eta_{l,i} = \enorm{\phi_{l/2}^{(i)} - \phi_{l}^{(i)}}$ diff --git a/src/A_step.m b/src/A_step.m index 038d47f..520155e 100644 --- a/src/A_step.m +++ b/src/A_step.m @@ -80,22 +80,26 @@ time = zeros(1,3); xd_fine = xo_fine'-x_fine; % |||h/2 -h||| - eta = xd_fine'*A_fine*xd_fine; +% eta = xd_fine'*A_fine*xd_fine; % \tilde \mu ( h/2 -h + L_2 ) mu = hmin.*b.*sum((x_fine(f2s)'-repmat(x',4,1)).^2)'/4; %Energienorm^2 Berechnen |||h||| & |||h/2||| - xe_fine = x_fine'*A_fine*x_fine; - xe = x'*A*x; +% xe_fine = x_fine'*A_fine*x_fine; + xe_fine = b_fine'*x_fine; +% xe = x'*A*x; + xe = b'*x; + + eta = xe_fine-xe; save_A{i} = A_fine; save_x{i} = x_fine; - data = [data type(i) sqrt(sum(tmu)) sqrt(eta) xe_fine sqrt(sum(mu))... + data = [data type(i) sqrt(sum(tmu)) sqrt(eta) xe sqrt(sum(mu))... min(hmin)/max(hmax) min(hmax)/max(hmax) min(hmin./hmax)]; end time(2) = toc; -- 2.47.3