[src] eta wird anders berechnet

[src] error = phi - phi_h/2

diff --git a/doc/doc.pdf b/doc/doc.pdf
index add9a02eabd437711c882856f64871fb358e9ade..cf94a25ed416e8393a25886ff49677bed02db60e 100644 (file)
Binary files a/doc/doc.pdf and b/doc/doc.pdf differ

diff --git a/doc/doc.tex b/doc/doc.tex
index 3d663abc2b24ca1e23caf7c65d7c89a9a3efad3e..102975581c2495af0a5e5a3cb182a8c20bf2d07f 100644 (file)
--- 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 038d47f4de7927583e62bfdf3bc532169c4fef7f..520155e2e897d6e912ffbd2a19c0dfd521f16012 100644 (file)
--- 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;