[src] nutzlose Funktionen entfernt

author Peter Schaefer <peter.schaefer@tuwien.ac.at>

Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)

committer Peter Schaefer <peter.schaefer@tuwien.ac.at>

Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)
author Peter Schaefer <peter.schaefer@tuwien.ac.at>
Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)
committer Peter Schaefer <peter.schaefer@tuwien.ac.at>
Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)
diff --git a/doc/doc.pdf b/doc/doc.pdf

index 8597ce2bb69b53a29b1b282f29e53f680935b08f..add9a02eabd437711c882856f64871fb358e9ade 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 5f20f08ac083d89729166bdd25527b1f7a4457d9..3d663abc2b24ca1e23caf7c65d7c89a9a3efad3e 100644 (file)
--- a/doc/doc.tex
+++ b/doc/doc.tex
@@ -249,13 +249,12 @@ wobei $\Pi_h$ die $L_2$ Projektion auf $P^0(\T_h)$ ist.
   \mu_h &= \norm{\varrho^{1/2}(\phi_{h/2} - \phi_h)}_{L^2(\Gamma)}\\
   \tilde\mu_h &= \norm{\varrho^{1/2}(\phi_{h/2} - \Pi_h\phi_{h/2})}_{L^2(\Gamma)}
  \end{align}
-Dann gilt:
+Dann gilt auf isotropen Netzen:
  \begin{itemize}
- \item Schätzer sind equivalent\\ $\eta_h \sim \tilde \eta_h \sim \mu_h \sim \tilde \mu_h$
+ \item Schätzer sind equivalent\\ $\tilde \mu_h \leq \mu_h \leq \sqrt{2}C_3 \eta_h$
   \item sie sind effizient
   \item sie sind unter Saturationsannahme auch zuverlässig
  \end{itemize}
-
  \end{sat}
  \begin{bew}
   Siehe S.F. Paper $\mapsto$ THM 3.2 \& 3.4
diff --git a/src/A_plot.m b/src/A_plot.m

deleted file mode 100644 (file)

index dffc3a9..0000000
--- a/src/A_plot.m
+++ /dev/null
@@ -1,35 +0,0 @@
-function A_plot()
-
-type2str = ['Analytisch  ' ; 'Quad Element' ; 'Quad Achse  '; 'Quad Seite  '];
-
-%%G_D
-global G_D;
-
-[m n] = size(G_D);
-
-step = round(n/3);
-
-
-if step<1
-    disp ('Error: No Data to show.')
-else
-
-figure(4)
-loglog(G_D(:,1),[G_D(:,[3+(0:step-1)*3]),20*G_D(:,1).^(-1/2),6*G_D(:,1).^(-1/4),5*G_D(:,1).^(-3/4)],...
-    G_D(1:end-1,1), sqrt(G_D(2:end,[4+(0:step-1)*3])-G_D(1:end-1,[4+(0:step-1)*3])),'--')
-
-title('Fehler')
-xlabel('Elemente');
-ylabel('Schaetzer');
-legend([['m';'m';'m'] type2str(G_D(1,[2+(0:step-1)*3])',:); 'N^{-1/2}     '; 'N^{-1/4}     '; 'N^{-3/4}     ';['g';'g';'g'] type2str(G_D(1,[2+(0:step-1)*3])',:)] ,'location','best');
-
-
-figure(5)
-loglog(G_D(:,1),G_D(:,[4+(0:step-1)*3]))
-title('Energie Norm')
-xlabel('Elemente');
-ylabel('eNorm');
-legend(type2str(G_D(1,[2+(0:step-1)*3])',:));
-end
-
-end
-\ No newline at end of file
diff --git a/src/A_plots.m b/src/A_plots.m

index 5c694fe8e78fd4d45ca29b05b084b95ec5633320..0b603b726371a5312df282a6512933a6be4ff125 100644 (file)
--- a/src/A_plots.m
+++ b/src/A_plots.m
@@ -53,7 +53,7 @@ for i = 1:length(files)
        leg2 = {leg2{:}...
            ['min hmin/max hmax ' l0 l1{i}]...
            ['min hmax/max hmax ' l0 l1{i}]...
-          ['hmin/hmax ' l0 l1{i}]...
+          ['min hmin/hmax ' l0 l1{i}]...
            }';
        sym = {sym{:} type2sym{data(1,[2+(i-1)*rows])}}'
     end
@@ -71,7 +71,7 @@ else
      
  %     sol = interp1(1./X((round(1)):(end),4)',G_D((round(1)):(end),4)',0,'spline')
  %     sol = 8.28466; % LShape
-    sol = 4.60957; % Quad
+    sol = 4.609193; % Quad
  
  % G_D
  
diff --git a/src/A_run.m b/src/A_run.m

index be5943d54c26752ea1d4d65fb9624d30ed1a540b..1476eef66ec098966108753edda155722c766ff3 100644 (file)
--- a/src/A_run.m
+++ b/src/A_run.m
@@ -1,8 +1,8 @@
-function A_run(file,times,mu,type,eta,eps,out)
+function A_run(file,times,zeta,type,eta,eps,out)
  % A_AnIso(file,times,mu,type,eta,eps,out)
  % file - starting mesh
  % times - how often
-% mu - Zulaessigkeitsbed
+% zeta - Zulaessigkeitsbed
  % type - art des Tests
  % eta - alle verfeinern oder nur wichtige? 
  % eps - wie verfeinern iso or aniso
@@ -11,7 +11,7 @@ function A_run(file,times,mu,type,eta,eps,out)
  % P. Schaefer
  
  %Parameter testen
-assert(length(mu)==max(type)-1||length(mu)==1,...
+assert(length(zeta)==max(type)-1||length(zeta)==1,...
    'Pleas set right type and mu parameters');
  
  %Netz laden
@@ -36,7 +36,7 @@ for i = 1:times
   end
   
   % Ein kompletter Verfeinerungschritt
- A_step(mu,type,eta,eps,out);
+ A_step(zeta,type,eta,eps,out);
   
   %Zeit Speichern
   usedTime = G_T(end,:);
diff --git a/src/A_step.m b/src/A_step.m

index 27ef6d5c18f9826919dfb1cc2fb38fbd5d182a47..71657c138f150460a57d2bcfa09bcef976e1c65d 100644 (file)
--- a/src/A_step.m
+++ b/src/A_step.m
@@ -38,7 +38,7 @@ time = zeros(1,3);
  
    %Flaecheninhalte Berechnen (rhs)
  %   b = sqrt(sum(quadNorm(coo_fine,ele_fine,'w').^2,2));
-  b_fine = area(coo_fine,ele_fine);
+  b_fine = areaQuad(coo_fine,ele_fine);
  
    tic
    %data -> ErgebnisMatrix
@@ -50,7 +50,7 @@ time = zeros(1,3);
    %Alle MatrixBrechenungsArten mit dem selben Netz berechnen
    for i = 1:length(type)
      %Matrix aufbauen -> MEX
-    A_fine = mex_build_AU(coo_fine,ele_fine,zeta,type(i));
+    A_fine = mex_build_V(coo_fine,ele_fine,zeta,type(i));
      
      %Testet auf Fehlerhafte Einträge (NaN +/-Inf)
      [r c] = find(isnan(A_fine)~=isinf(A_fine));
@@ -64,9 +64,7 @@ time = zeros(1,3);
      %Lösung Berechnen
      x_fine = A_fine\b_fine;
  
-
-    
-    [b,tnormal,s] = area(G_C,G_E);
+    [b,tnormal,s] = areaQuad(G_C,G_E);
      hmin = min(s,[],2);
      hmax = max(s,[],2);
      
@@ -75,17 +73,17 @@ time = zeros(1,3);
  
      
      %Fehlerschätzer 2 aufbauen
-    A = mex_build_AU(G_C,G_E,zeta,type(i));
+    A = mex_build_V(G_C,G_E,zeta,type(i));
      x = A\b;
      
      xo_fine(f2s) = repmat(x,1,4);
      xd_fine = xo_fine'-x_fine;
      
      % |||h/2 -h|||
-    tmu2 = xd_fine'*A_fine*xd_fine;
+    eta = xd_fine'*A_fine*xd_fine;
      
      % \tilde \mu ( h/2 -h + L_2 )
-    tmu3 = hmin.*b.*sum((x_fine(f2s)'-repmat(x',4,1)).^2)'/4;
+    mu = hmin.*b.*sum((x_fine(f2s)'-repmat(x',4,1)).^2)'/4;
  
  
      %Energienorm^2 Berechnen |||h||| & |||h/2|||
@@ -97,7 +95,7 @@ time = zeros(1,3);
      
      
      
-    data = [data type(i) sqrt(sum(tmu)) sqrt(tmu2) xe_fine sqrt(sum(tmu3))...
+    data = [data type(i) sqrt(sum(tmu)) sqrt(eta) xe_fine sqrt(sum(mu))...
          min(hmin)/max(hmax) min(hmax)/max(hmax) sqrt(sum((hmin./hmax).^2))];
    end
    time(2) = toc;
diff --git a/src/Aelements.m b/src/Aelements.m

deleted file mode 100644 (file)

index aad35f5..0000000
--- a/src/Aelements.m
+++ /dev/null
@@ -1,19 +0,0 @@
-function A = Aelements(coordinates, elements)
-% This function calculates the Sum of all Surfaces
-%
-% A=Aelements(coordinates,elements)
-%
-% P.Schaefer
-
-elles = size(elements,1);
-A=0;
-for i = 1:elles
-    tri = elements(i,:);
-    a = (coordinates(tri(1),:)-coordinates(tri(2),:));
-    b = (coordinates(tri(1),:)-coordinates(tri(4),:));
-    A1 = norm(cross(a',b'));
-    A = A+ A1;
-end
-
-end
-
diff --git a/src/area.m b/src/area.m

deleted file mode 100644 (file)

index f4c657a..0000000
--- a/src/area.m
+++ /dev/null
@@ -1,36 +0,0 @@
-function [area, vector, sites] = area(coordinates, elements,varargin)
-%
-% [area, vector, sites] = area(coordinates, elements)
-%
-% Diese Funktion Berechnet den Flächeninhalt zu jedem Element (area) und 
-% zu jedem Element die Orthogonalen mit Laenge 1, des Weiteren werden auch
-% die Seitenlaengen zu jedem Element gespeichert
-%
-% P.Schaefer
-
-    vector = zeros(size(elements,1),3);
-    area = zeros(size(elements,1),1);
-    sites = zeros(size(elements,1),2);
-    
-%%  Flaecheninhalt und Normalen berechnen
-
-    for i = 1:size(elements,1)
-        % normalized Vector on every triangle
-        tri = elements(i,:);
-        a = (coordinates(tri(2),:)-coordinates(tri(1),:));
-        b = (coordinates(tri(4),:)-coordinates(tri(1),:));
-        N = cross(a',b');
-        
-%         area(i) = norm(N);
-        
-        N = N/norm(N);
-        vector(i,:) = N;
-        area(i) = norm(a) * norm(b);
-        sites(i,:) = [norm(a) norm(b)];
-        
-    end
-
-end
-
-
-
diff --git a/src/areaQuad.m b/src/areaQuad.m

new file mode 100644 (file)

index 0000000..8473339
--- /dev/null
+++ b/src/areaQuad.m
@@ -0,0 +1,36 @@
+function [area, vector, sites] = areaQuad(coordinates, elements,varargin)
+%
+% [area, vector, sites] = area(coordinates, elements)
+%
+% Diese Funktion Berechnet den Flächeninhalt zu jedem Element (area) und 
+% zu jedem Element die Orthogonalen mit Laenge 1, des Weiteren werden auch
+% die Seitenlaengen zu jedem Element gespeichert
+%
+% P.Schaefer
+
+    vector = zeros(size(elements,1),3);
+    area = zeros(size(elements,1),1);
+    sites = zeros(size(elements,1),2);
+    
+%%  Flaecheninhalt und Normalen berechnen
+
+    for i = 1:size(elements,1)
+        % normalized Vector on every triangle
+        tri = elements(i,:);
+        a = (coordinates(tri(2),:)-coordinates(tri(1),:));
+        b = (coordinates(tri(4),:)-coordinates(tri(1),:));
+        N = cross(a',b');
+        
+%         area(i) = norm(N);
+        
+        N = N/norm(N);
+        vector(i,:) = N;
+        area(i) = norm(a) * norm(b);
+        sites(i,:) = [norm(a) norm(b)];
+        
+    end
+
+end
+
+
+
diff --git a/src/computeEstSlpMuTilde.m b/src/computeEstSlpMuTilde.m

deleted file mode 100644 (file)

index dce7564..0000000
--- a/src/computeEstSlpMuTilde.m
+++ /dev/null
@@ -1,34 +0,0 @@
-function ind=computeEstSlpMuTilde(x_fine,coo,ele,f2s)
-% ind=computeEstSlpMuTilde(x_fine,coo,ele,f2s)
-% Berechnet den Fehlerschätzer über die Elemente
-%
-% x_fine - Phi
-% coo - coordinaten
-% ele - elemente
-% f2s - Vater Sohn beziehungen
-%
-% P. Schaefer
-
-xF2S = x_fine(f2s)';
-
-if(size(xF2S,1)==1)
-  xF2S = xF2S';
-end
-
-tmp = abs([sum(coo(ele(:,2),:)-coo(ele(:,1),:),2)...
-  sum(coo(ele(:,4),:)-coo(ele(:,1),:),2)]);
-hmin = min(tmp,[],2);
-% hmax = max(tmp,[],2);
-
-% len = min(abs([sum(coo(ele(:,2),:)-coo(ele(:,1),:),2)...
-%   sum(coo(ele(:,4),:)-coo(ele(:,1),:),2)]),[],2);
-
-sizes = sqrt(sum(quadNorm(coo,ele,'w').^2,2));
-
-% max(sizes ~= hmin.*hmax)
-
-x_var = sum((xF2S-repmat(sum(xF2S,1)/4,4,1)).^2)';
-
-ind = hmin.* sizes .* x_var/4;
-
-end
-\ No newline at end of file
diff --git a/src/expmlAA_2DQuad_error.eps b/src/expmlAA_2DQuad_error.eps

deleted file mode 100644 (file)

index 7a86d86..0000000
--- a/src/expmlAA_2DQuad_error.eps
+++ /dev/null
@@ -1,824 +0,0 @@
-%!PS-Adobe-2.0 EPSF-1.2
-%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64.
-%%Title: ./expmlAA_2DQuad_error.eps
-%%CreationDate: 03/26/2012  23:43:43
-%%DocumentNeededFonts: Helvetica
-%%DocumentProcessColors: Cyan Magenta Yellow Black
-%%Extensions: CMYK
-%%Pages: 1
-%%BoundingBox:    53   195   550   602
-%%EndComments
-
-%%BeginProlog
-% MathWorks dictionary
-/MathWorks 160 dict begin
-% definition operators
-/bdef {bind def} bind def
-/ldef {load def} bind def
-/xdef {exch def} bdef
-/xstore {exch store} bdef
-% operator abbreviations
-/c  /clip ldef
-/cc /concat ldef
-/cp /closepath ldef
-/gr /grestore ldef
-/gs /gsave ldef
-/mt /moveto ldef
-/np /newpath ldef
-/cm /currentmatrix ldef
-/sm /setmatrix ldef
-/rm /rmoveto ldef
-/rl /rlineto ldef
-/s {show newpath} bdef
-/sc {setcmykcolor} bdef
-/sr /setrgbcolor ldef
-/sg /setgray ldef
-/w /setlinewidth ldef
-/j /setlinejoin ldef
-/cap /setlinecap ldef
-/rc {rectclip} bdef
-/rf {rectfill} bdef
-% page state control
-/pgsv () def
-/bpage {/pgsv save def} bdef
-/epage {pgsv restore} bdef
-/bplot /gsave ldef
-/eplot {stroke grestore} bdef
-% orientation switch
-/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
-% coordinate system mappings
-/dpi2point 0 def
-% font control
-/FontSize 0 def
-/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
-  makefont setfont} bdef
-/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
-  exch dup 3 1 roll findfont dup length dict begin
-  { 1 index /FID ne {def}{pop pop} ifelse } forall
-  /Encoding exch def currentdict end definefont pop} bdef
-/isroman {findfont /CharStrings get /Agrave known} bdef
-/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
-  exch FMS} bdef
-/csm {1 dpi2point div -1 dpi2point div scale neg translate
- dup landscapeMode eq {pop -90 rotate}
-  {rotateMode eq {90 rotate} if} ifelse} bdef
-% line types: solid, dotted, dashed, dotdash
-/SO { [] 0 setdash } bdef
-/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
-/DA { [6 dpi2point mul] 0 setdash } bdef
-/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
-  dpi2point mul] 0 setdash } bdef
-% macros for lines and objects
-/L {lineto stroke} bdef
-/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
-/AP {{rlineto} repeat} bdef
-/PDlw -1 def
-/W {/PDlw currentlinewidth def setlinewidth} def
-/PP {closepath eofill} bdef
-/DP {closepath stroke} bdef
-/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
-  neg 0 exch rlineto closepath} bdef
-/FR {MR stroke} bdef
-/PR {MR fill} bdef
-/L1i {{currentfile picstr readhexstring pop} image} bdef
-/tMatrix matrix def
-/MakeOval {newpath tMatrix currentmatrix pop translate scale
-0 0 1 0 360 arc tMatrix setmatrix} bdef
-/FO {MakeOval stroke} bdef
-/PO {MakeOval fill} bdef
-/PD {currentlinewidth 2 div 0 360 arc fill
-   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
-/FA {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
-/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
-/FAn {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
-/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
-/vradius 0 def /hradius 0 def /lry 0 def
-/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
-/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
-  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
-  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
-  tMatrix setmatrix lrx hradius sub uly vradius add translate
-  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
-  lrx hradius sub lry vradius sub translate hradius vradius scale
-  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
-  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
-  closepath} bdef
-/FRR {MRR stroke } bdef
-/PRR {MRR fill } bdef
-/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
-  closepath} bdef
-/FlrRR {MlrRR stroke } bdef
-/PlrRR {MlrRR fill } bdef
-/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
-  closepath} bdef
-/FtbRR {MtbRR stroke } bdef
-/PtbRR {MtbRR fill } bdef
-/stri 6 array def /dtri 6 array def
-/smat 6 array def /dmat 6 array def
-/tmat1 6 array def /tmat2 6 array def /dif 3 array def
-/asub {/ind2 exch def /ind1 exch def dup dup
-  ind1 get exch ind2 get sub exch } bdef
-/tri_to_matrix {
-  2 0 asub 3 1 asub 4 0 asub 5 1 asub
-  dup 0 get exch 1 get 7 -1 roll astore } bdef
-/compute_transform {
-  dmat dtri tri_to_matrix tmat1 invertmatrix 
-  smat stri tri_to_matrix tmat2 concatmatrix } bdef
-/ds {stri astore pop} bdef
-/dt {dtri astore pop} bdef
-/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
-  currentfile exch readhexstring pop
-  dup 0 3 index getinterval /rbmap xdef
-  dup 2 index dup getinterval /gbmap xdef
-  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
-/it {gs np dtri aload pop moveto lineto lineto cp c
-  cols rows 8 compute_transform 
-  rbmap gbmap bbmap true 3 colorimage gr}bdef
-/il {newpath moveto lineto stroke}bdef
-currentdict end def
-%%EndProlog
-
-%%BeginSetup
-MathWorks begin
-
-0 cap
-
-end
-%%EndSetup
-
-%%Page: 1 1
-%%BeginPageSetup
-%%PageBoundingBox:    53   195   550   602
-MathWorks begin
-bpage
-%%EndPageSetup
-
-%%BeginObject: obj1
-bplot
-
-/dpi2point 8.33333 def
-portraitMode 0150 5100 csm
-
-  294    79  4147  3390 MR c np
-76 dict begin %Colortable dictionary
-/c0 { 0.000000 0.000000 0.000000 sr} bdef
-/c1 { 1.000000 1.000000 1.000000 sr} bdef
-/c2 { 0.900000 0.000000 0.000000 sr} bdef
-/c3 { 0.000000 0.820000 0.000000 sr} bdef
-/c4 { 0.000000 0.000000 0.800000 sr} bdef
-/c5 { 0.910000 0.820000 0.320000 sr} bdef
-/c6 { 1.000000 0.260000 0.820000 sr} bdef
-/c7 { 0.000000 0.820000 0.820000 sr} bdef
-c0
-1 j
-1 sg
-   0    0 4803 3605 PR
-4.16667 w
-0 2938 3722 0 0 -2938 624 3208 4 MP
-PP
--3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
- 624 3208 mt 4346 3208 L
- 624 3208 mt  624  270 L
- 624 3208 mt  624 3189 L
- 624  270 mt  624  288 L
- 624 3208 mt  624 3170 L
- 624  270 mt  624  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 563 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 655 3283 mt 
-(0) s
- 904 3208 mt  904 3189 L
- 904  270 mt  904  288 L
-1067 3208 mt 1067 3189 L
-1067  270 mt 1067  288 L
-1184 3208 mt 1184 3189 L
-1184  270 mt 1184  288 L
-1274 3208 mt 1274 3189 L
-1274  270 mt 1274  288 L
-1348 3208 mt 1348 3189 L
-1348  270 mt 1348  288 L
-1410 3208 mt 1410 3189 L
-1410  270 mt 1410  288 L
-1464 3208 mt 1464 3189 L
-1464  270 mt 1464  288 L
-1511 3208 mt 1511 3189 L
-1511  270 mt 1511  288 L
-1554 3208 mt 1554 3189 L
-1554  270 mt 1554  288 L
-1554 3208 mt 1554 3170 L
-1554  270 mt 1554  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1493 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-1585 3283 mt 
-(1) s
-1834 3208 mt 1834 3189 L
-1834  270 mt 1834  288 L
-1998 3208 mt 1998 3189 L
-1998  270 mt 1998  288 L
-2114 3208 mt 2114 3189 L
-2114  270 mt 2114  288 L
-2204 3208 mt 2204 3189 L
-2204  270 mt 2204  288 L
-2278 3208 mt 2278 3189 L
-2278  270 mt 2278  288 L
-2340 3208 mt 2340 3189 L
-2340  270 mt 2340  288 L
-2394 3208 mt 2394 3189 L
-2394  270 mt 2394  288 L
-2442 3208 mt 2442 3189 L
-2442  270 mt 2442  288 L
-2485 3208 mt 2485 3189 L
-2485  270 mt 2485  288 L
-2485 3208 mt 2485 3170 L
-2485  270 mt 2485  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2424 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-2516 3283 mt 
-(2) s
-2765 3208 mt 2765 3189 L
-2765  270 mt 2765  288 L
-2928 3208 mt 2928 3189 L
-2928  270 mt 2928  288 L
-3045 3208 mt 3045 3189 L
-3045  270 mt 3045  288 L
-3135 3208 mt 3135 3189 L
-3135  270 mt 3135  288 L
-3209 3208 mt 3209 3189 L
-3209  270 mt 3209  288 L
-3271 3208 mt 3271 3189 L
-3271  270 mt 3271  288 L
-3325 3208 mt 3325 3189 L
-3325  270 mt 3325  288 L
-3372 3208 mt 3372 3189 L
-3372  270 mt 3372  288 L
-3415 3208 mt 3415 3189 L
-3415  270 mt 3415  288 L
-3415 3208 mt 3415 3170 L
-3415  270 mt 3415  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-3354 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-3446 3283 mt 
-(3) s
-3695 3208 mt 3695 3189 L
-3695  270 mt 3695  288 L
-3859 3208 mt 3859 3189 L
-3859  270 mt 3859  288 L
-3975 3208 mt 3975 3189 L
-3975  270 mt 3975  288 L
-4065 3208 mt 4065 3189 L
-4065  270 mt 4065  288 L
-4139 3208 mt 4139 3189 L
-4139  270 mt 4139  288 L
-4201 3208 mt 4201 3189 L
-4201  270 mt 4201  288 L
-4255 3208 mt 4255 3189 L
-4255  270 mt 4255  288 L
-4303 3208 mt 4303 3189 L
-4303  270 mt 4303  288 L
-4346 3208 mt 4346 3189 L
-4346  270 mt 4346  288 L
-4346 3208 mt 4346 3170 L
-4346  270 mt 4346  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-4285 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-4377 3283 mt 
-(4) s
- 624 3208 mt  642 3208 L
-4346 3208 mt 4327 3208 L
- 624 3208 mt  661 3208 L
-4346 3208 mt 4308 3208 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 3238 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 3187 mt 
-(-3) s
- 624 2986 mt  642 2986 L
-4346 2986 mt 4327 2986 L
- 624 2857 mt  642 2857 L
-4346 2857 mt 4327 2857 L
- 624 2765 mt  642 2765 L
-4346 2765 mt 4327 2765 L
- 624 2694 mt  642 2694 L
-4346 2694 mt 4327 2694 L
- 624 2636 mt  642 2636 L
-4346 2636 mt 4327 2636 L
- 624 2587 mt  642 2587 L
-4346 2587 mt 4327 2587 L
- 624 2544 mt  642 2544 L
-4346 2544 mt 4327 2544 L
- 624 2507 mt  642 2507 L
-4346 2507 mt 4327 2507 L
- 624 2473 mt  642 2473 L
-4346 2473 mt 4327 2473 L
- 624 2473 mt  661 2473 L
-4346 2473 mt 4308 2473 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 2503 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 2452 mt 
-(-2) s
- 624 2252 mt  642 2252 L
-4346 2252 mt 4327 2252 L
- 624 2123 mt  642 2123 L
-4346 2123 mt 4327 2123 L
- 624 2031 mt  642 2031 L
-4346 2031 mt 4327 2031 L
- 624 1960 mt  642 1960 L
-4346 1960 mt 4327 1960 L
- 624 1901 mt  642 1901 L
-4346 1901 mt 4327 1901 L
- 624 1852 mt  642 1852 L
-4346 1852 mt 4327 1852 L
- 624 1810 mt  642 1810 L
-4346 1810 mt 4327 1810 L
- 624 1772 mt  642 1772 L
-4346 1772 mt 4327 1772 L
- 624 1739 mt  642 1739 L
-4346 1739 mt 4327 1739 L
- 624 1739 mt  661 1739 L
-4346 1739 mt 4308 1739 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 1769 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 1718 mt 
-(-1) s
- 624 1517 mt  642 1517 L
-4346 1517 mt 4327 1517 L
- 624 1388 mt  642 1388 L
-4346 1388 mt 4327 1388 L
- 624 1296 mt  642 1296 L
-4346 1296 mt 4327 1296 L
- 624 1225 mt  642 1225 L
-4346 1225 mt 4327 1225 L
- 624 1167 mt  642 1167 L
-4346 1167 mt 4327 1167 L
- 624 1118 mt  642 1118 L
-4346 1118 mt 4327 1118 L
- 624 1075 mt  642 1075 L
-4346 1075 mt 4327 1075 L
- 624 1038 mt  642 1038 L
-4346 1038 mt 4327 1038 L
- 624 1004 mt  642 1004 L
-4346 1004 mt 4327 1004 L
- 624 1004 mt  661 1004 L
-4346 1004 mt 4308 1004 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 1034 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537  983 mt 
-(0) s
- 624  783 mt  642  783 L
-4346  783 mt 4327  783 L
- 624  654 mt  642  654 L
-4346  654 mt 4327  654 L
- 624  562 mt  642  562 L
-4346  562 mt 4327  562 L
- 624  491 mt  642  491 L
-4346  491 mt 4327  491 L
- 624  432 mt  642  432 L
-4346  432 mt 4327  432 L
- 624  383 mt  642  383 L
-4346  383 mt 4327  383 L
- 624  341 mt  642  341 L
-4346  341 mt 4327  341 L
- 624  303 mt  642  303 L
-4346  303 mt 4327  303 L
- 624  270 mt  642  270 L
-4346  270 mt 4327  270 L
- 624  270 mt  661  270 L
-4346  270 mt 4308  270 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445  300 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537  249 mt 
-(1) s
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
-gs 624 270 3723 2939 MR c np
-/c8 { 0.000000 0.000000 1.000000 sr} bdef
-c8
-68 -537 100 -204 65 7 95 -61 107 14 100 72 97 70 93 67 
-87 68 72 63 55 60 84 58 55 56 64 55 60 55 60 57 
-62 48 45 49 62 50 48 50 51 53 71 51 51 52 90 53 
-71 56 104 77 122 51 190 65 254 54 560 71 1184 933 31 MP stroke
-gs 1133 882 3146 1578 MR c np
-  25   25 1184  933 FO
-  25   25 1744 1004 FO
-  25   25 1998 1058 FO
-  25   25 2188 1123 FO
-  25   25 2310 1174 FO
-  25   25 2414 1251 FO
-  25   25 2485 1307 FO
-  25   25 2575 1360 FO
-  25   25 2626 1412 FO
-  25   25 2697 1463 FO
-  25   25 2748 1516 FO
-  25   25 2796 1566 FO
-  25   25 2858 1616 FO
-  25   25 2903 1665 FO
-  25   25 2965 1713 FO
-  25   25 3025 1770 FO
-  25   25 3085 1825 FO
-  25   25 3149 1880 FO
-  25   25 3204 1936 FO
-  25   25 3288 1994 FO
-  25   25 3343 2054 FO
-  25   25 3415 2117 FO
-  25   25 3502 2185 FO
-  25   25 3595 2252 FO
-  25   25 3692 2322 FO
-  25   25 3792 2394 FO
-  25   25 3899 2408 FO
-  25   25 3994 2347 FO
-  25   25 4059 2354 FO
-  25   25 4159 2150 FO
-  25   25 4227 1613 FO
-gr
-
-/c9 { 0.000000 0.500000 0.000000 sr} bdef
-c9
-68 5 100 11 65 11 95 56 107 64 100 52 97 36 93 52 
-87 44 72 29 55 33 84 52 55 40 64 43 60 46 60 50 
-62 48 45 50 62 46 48 52 51 47 71 54 51 54 90 46 
-71 61 104 69 122 59 190 59 254 59 560 92 1184 1285 31 MP stroke
-gs 1133 1234 3146 1523 MR c np
-  25   25 1184 1285 FO
-  25   25 1744 1377 FO
-  25   25 1998 1436 FO
-  25   25 2188 1495 FO
-  25   25 2310 1554 FO
-  25   25 2414 1623 FO
-  25   25 2485 1684 FO
-  25   25 2575 1730 FO
-  25   25 2626 1784 FO
-  25   25 2697 1838 FO
-  25   25 2748 1885 FO
-  25   25 2796 1937 FO
-  25   25 2858 1983 FO
-  25   25 2903 2033 FO
-  25   25 2965 2081 FO
-  25   25 3025 2131 FO
-  25   25 3085 2177 FO
-  25   25 3149 2220 FO
-  25   25 3204 2260 FO
-  25   25 3288 2312 FO
-  25   25 3343 2345 FO
-  25   25 3415 2374 FO
-  25   25 3502 2418 FO
-  25   25 3595 2470 FO
-  25   25 3692 2506 FO
-  25   25 3792 2558 FO
-  25   25 3899 2622 FO
-  25   25 3994 2678 FO
-  25   25 4059 2689 FO
-  25   25 4159 2700 FO
-  25   25 4227 2705 FO
-gr
-
-/c10 { 1.000000 0.000000 0.000000 sr} bdef
-c10
-68 1 100 0 65 0 95 2 107 2 100 2 97 3 93 5 
-87 6 72 7 55 12 84 19 55 23 64 28 60 33 60 38 
-62 49 45 49 62 51 48 51 51 51 71 53 51 52 90 53 
-71 57 104 76 122 56 190 68 254 64 560 97 1184 1252 31 MP stroke
-gs 1133 1201 3146 1111 MR c np
-  25   25 1184 1252 FO
-  25   25 1744 1349 FO
-  25   25 1998 1413 FO
-  25   25 2188 1481 FO
-  25   25 2310 1537 FO
-  25   25 2414 1613 FO
-  25   25 2485 1670 FO
-  25   25 2575 1723 FO
-  25   25 2626 1775 FO
-  25   25 2697 1828 FO
-  25   25 2748 1879 FO
-  25   25 2796 1930 FO
-  25   25 2858 1981 FO
-  25   25 2903 2030 FO
-  25   25 2965 2079 FO
-  25   25 3025 2117 FO
-  25   25 3085 2150 FO
-  25   25 3149 2178 FO
-  25   25 3204 2201 FO
-  25   25 3288 2220 FO
-  25   25 3343 2232 FO
-  25   25 3415 2239 FO
-  25   25 3502 2245 FO
-  25   25 3595 2250 FO
-  25   25 3692 2253 FO
-  25   25 3792 2255 FO
-  25   25 3899 2257 FO
-  25   25 3994 2259 FO
-  25   25 4059 2259 FO
-  25   25 4159 2259 FO
-  25   25 4227 2260 FO
-gr
-
-/c11 { 0.000000 0.750000 0.750000 sr} bdef
-c11
-68 -525 100 118 65 -379 95 -65 107 14 100 70 97 70 93 63 
-87 68 72 62 55 59 84 59 55 56 64 55 60 55 60 57 
-62 49 45 48 62 51 48 50 51 54 71 51 51 51 90 54 
-71 55 104 80 122 48 190 66 254 51 560 70 1184 931 31 MP stroke
-gs 1133 880 3146 1569 MR c np
-  25   25 1184  931 FO
-  25   25 1744 1001 FO
-  25   25 1998 1052 FO
-  25   25 2188 1118 FO
-  25   25 2310 1166 FO
-  25   25 2414 1246 FO
-  25   25 2485 1301 FO
-  25   25 2575 1355 FO
-  25   25 2626 1406 FO
-  25   25 2697 1457 FO
-  25   25 2748 1511 FO
-  25   25 2796 1561 FO
-  25   25 2858 1612 FO
-  25   25 2903 1660 FO
-  25   25 2965 1709 FO
-  25   25 3025 1766 FO
-  25   25 3085 1821 FO
-  25   25 3149 1876 FO
-  25   25 3204 1932 FO
-  25   25 3288 1991 FO
-  25   25 3343 2050 FO
-  25   25 3415 2112 FO
-  25   25 3502 2180 FO
-  25   25 3595 2243 FO
-  25   25 3692 2313 FO
-  25   25 3792 2383 FO
-  25   25 3899 2397 FO
-  25   25 3994 2332 FO
-  25   25 4059 1953 FO
-  25   25 4159 2071 FO
-  25   25 4227 1546 FO
-gr
-
-DD
-c8
-68 27 100 40 65 26 95 37 107 42 100 39 97 39 93 37 
-87 34 72 29 55 21 84 33 55 22 64 25 60 24 60 24 
-62 24 45 18 62 24 48 19 51 21 71 27 51 21 90 35 
-71 28 104 41 122 48 190 75 254 101 560 221 1184 604 31 MP stroke
-c9
-68 14 100 20 65 12 95 19 107 21 100 20 97 19 93 18 
-87 18 72 14 55 11 84 16 55 11 64 13 60 12 60 11 
-62 13 45 8 62 13 48 9 51 10 71 14 51 10 90 18 
-71 14 104 21 122 24 190 37 254 50 560 111 1184 764 31 MP stroke
-c10
-68 40 100 59 65 39 95 56 107 64 100 58 97 58 93 55 
-87 52 72 42 55 33 84 50 55 32 64 38 60 35 60 36 
-62 36 45 27 62 37 48 28 51 31 71 41 51 31 90 53 
-71 42 104 61 122 73 190 112 254 151 560 331 1184 1115 31 MP stroke
-gr
-
-c10
-DD
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2367  206 mt 
-(Fehler) s
-2308 3434 mt 
-(Elemente) s
- 388 1927 mt  -90 rotate
-(Schaetzer) s
-90 rotate
- 612 3239 mt 
-( ) s
-4335  300 mt 
-( ) s
-SO
-1 sg
-0 854 1238 0 0 -854 665 3168 4 MP
-PP
--1238 0 0 854 1238 0 0 -854 665 3168 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
- 665 3168 mt 1903 3168 L
- 665 3168 mt  665 2314 L
- 998 2414 mt 
-(tilde ) s
-%%IncludeResource: font Symbol
-/Symbol /ISOLatin1Encoding 83.3333 FMSR
-
-1174 2414 mt 
-(m) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1222 2414 mt 
-( testAA Analytisch) s
-gs 665 2314 1239 855 MR c np
-c8
-256 0 716 2383 2 MP stroke
-gs 793 2332 103 103 MR c np
-  25   25  844 2383 FO
-gr
-
-gr
-
-c8
-0 sg
-%%IncludeResource: font Symbol
-/Symbol /ISOLatin1Encoding 83.3333 FMSR
-
- 998 2527 mt 
-(h) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1048 2527 mt 
-( testAA Analytisch) s
-gs 665 2314 1239 855 MR c np
-c9
-256 0 716 2496 2 MP stroke
-gs 793 2445 103 103 MR c np
-  25   25  844 2496 FO
-gr
-
-gr
-
-c9
-0 sg
- 998 2635 mt 
-(error testAA Analytisch) s
-gs 665 2314 1239 855 MR c np
-c10
-256 0 716 2605 2 MP stroke
-gs 793 2554 103 103 MR c np
-  25   25  844 2605 FO
-gr
-
-gr
-
-c10
-0 sg
-%%IncludeResource: font Symbol
-/Symbol /ISOLatin1Encoding 83.3333 FMSR
-
- 998 2744 mt 
-(m) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1046 2744 mt 
-( testAA Analytisch) s
-gs 665 2314 1239 855 MR c np
-c11
-256 0 716 2713 2 MP stroke
-gs 793 2662 103 103 MR c np
-  25   25  844 2713 FO
-gr
-
-gr
-
-c11
-0 sg
- 998 2875 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1058 2834 mt 
-(-1/2) s
-gs 665 2314 1239 855 MR c np
-DD
-c8
-256 0 716 2834 2 MP stroke
-SO
-gr
-
-c8
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 998 3003 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1058 2962 mt 
-(-1/4) s
-gs 665 2314 1239 855 MR c np
-DD
-c9
-256 0 716 2962 2 MP stroke
-SO
-gr
-
-c9
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 998 3131 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1058 3090 mt 
-(-3/4) s
-gs 665 2314 1239 855 MR c np
-DD
-c10
-256 0 716 3090 2 MP stroke
-SO
-gr
-
-c10
-
-end %%Color Dict
-
-eplot
-%%EndObject
-
-epage
-end
-
-showpage
-
-%%Trailer
-%%EOF
diff --git a/src/expmlAA_2DQuad_hminmax.eps b/src/expmlAA_2DQuad_hminmax.eps

deleted file mode 100644 (file)

index 892747d..0000000
--- a/src/expmlAA_2DQuad_hminmax.eps
+++ /dev/null
@@ -1,865 +0,0 @@
-%!PS-Adobe-2.0 EPSF-1.2
-%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64.
-%%Title: ./expmlAA_2DQuad_hminmax.eps
-%%CreationDate: 03/26/2012  23:43:44
-%%DocumentNeededFonts: Helvetica
-%%DocumentProcessColors: Cyan Magenta Yellow Black
-%%Extensions: CMYK
-%%Pages: 1
-%%BoundingBox:    53   195   550   602
-%%EndComments
-
-%%BeginProlog
-% MathWorks dictionary
-/MathWorks 160 dict begin
-% definition operators
-/bdef {bind def} bind def
-/ldef {load def} bind def
-/xdef {exch def} bdef
-/xstore {exch store} bdef
-% operator abbreviations
-/c  /clip ldef
-/cc /concat ldef
-/cp /closepath ldef
-/gr /grestore ldef
-/gs /gsave ldef
-/mt /moveto ldef
-/np /newpath ldef
-/cm /currentmatrix ldef
-/sm /setmatrix ldef
-/rm /rmoveto ldef
-/rl /rlineto ldef
-/s {show newpath} bdef
-/sc {setcmykcolor} bdef
-/sr /setrgbcolor ldef
-/sg /setgray ldef
-/w /setlinewidth ldef
-/j /setlinejoin ldef
-/cap /setlinecap ldef
-/rc {rectclip} bdef
-/rf {rectfill} bdef
-% page state control
-/pgsv () def
-/bpage {/pgsv save def} bdef
-/epage {pgsv restore} bdef
-/bplot /gsave ldef
-/eplot {stroke grestore} bdef
-% orientation switch
-/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
-% coordinate system mappings
-/dpi2point 0 def
-% font control
-/FontSize 0 def
-/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
-  makefont setfont} bdef
-/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
-  exch dup 3 1 roll findfont dup length dict begin
-  { 1 index /FID ne {def}{pop pop} ifelse } forall
-  /Encoding exch def currentdict end definefont pop} bdef
-/isroman {findfont /CharStrings get /Agrave known} bdef
-/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
-  exch FMS} bdef
-/csm {1 dpi2point div -1 dpi2point div scale neg translate
- dup landscapeMode eq {pop -90 rotate}
-  {rotateMode eq {90 rotate} if} ifelse} bdef
-% line types: solid, dotted, dashed, dotdash
-/SO { [] 0 setdash } bdef
-/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
-/DA { [6 dpi2point mul] 0 setdash } bdef
-/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
-  dpi2point mul] 0 setdash } bdef
-% macros for lines and objects
-/L {lineto stroke} bdef
-/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
-/AP {{rlineto} repeat} bdef
-/PDlw -1 def
-/W {/PDlw currentlinewidth def setlinewidth} def
-/PP {closepath eofill} bdef
-/DP {closepath stroke} bdef
-/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
-  neg 0 exch rlineto closepath} bdef
-/FR {MR stroke} bdef
-/PR {MR fill} bdef
-/L1i {{currentfile picstr readhexstring pop} image} bdef
-/tMatrix matrix def
-/MakeOval {newpath tMatrix currentmatrix pop translate scale
-0 0 1 0 360 arc tMatrix setmatrix} bdef
-/FO {MakeOval stroke} bdef
-/PO {MakeOval fill} bdef
-/PD {currentlinewidth 2 div 0 360 arc fill
-   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
-/FA {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
-/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
-/FAn {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
-/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
-/vradius 0 def /hradius 0 def /lry 0 def
-/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
-/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
-  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
-  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
-  tMatrix setmatrix lrx hradius sub uly vradius add translate
-  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
-  lrx hradius sub lry vradius sub translate hradius vradius scale
-  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
-  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
-  closepath} bdef
-/FRR {MRR stroke } bdef
-/PRR {MRR fill } bdef
-/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
-  closepath} bdef
-/FlrRR {MlrRR stroke } bdef
-/PlrRR {MlrRR fill } bdef
-/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
-  closepath} bdef
-/FtbRR {MtbRR stroke } bdef
-/PtbRR {MtbRR fill } bdef
-/stri 6 array def /dtri 6 array def
-/smat 6 array def /dmat 6 array def
-/tmat1 6 array def /tmat2 6 array def /dif 3 array def
-/asub {/ind2 exch def /ind1 exch def dup dup
-  ind1 get exch ind2 get sub exch } bdef
-/tri_to_matrix {
-  2 0 asub 3 1 asub 4 0 asub 5 1 asub
-  dup 0 get exch 1 get 7 -1 roll astore } bdef
-/compute_transform {
-  dmat dtri tri_to_matrix tmat1 invertmatrix 
-  smat stri tri_to_matrix tmat2 concatmatrix } bdef
-/ds {stri astore pop} bdef
-/dt {dtri astore pop} bdef
-/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
-  currentfile exch readhexstring pop
-  dup 0 3 index getinterval /rbmap xdef
-  dup 2 index dup getinterval /gbmap xdef
-  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
-/it {gs np dtri aload pop moveto lineto lineto cp c
-  cols rows 8 compute_transform 
-  rbmap gbmap bbmap true 3 colorimage gr}bdef
-/il {newpath moveto lineto stroke}bdef
-currentdict end def
-%%EndProlog
-
-%%BeginSetup
-MathWorks begin
-
-0 cap
-
-end
-%%EndSetup
-
-%%Page: 1 1
-%%BeginPageSetup
-%%PageBoundingBox:    53   195   550   602
-MathWorks begin
-bpage
-%%EndPageSetup
-
-%%BeginObject: obj1
-bplot
-
-/dpi2point 8.33333 def
-portraitMode 0150 5100 csm
-
-  294    78  4147  3391 MR c np
-76 dict begin %Colortable dictionary
-/c0 { 0.000000 0.000000 0.000000 sr} bdef
-/c1 { 1.000000 1.000000 1.000000 sr} bdef
-/c2 { 0.900000 0.000000 0.000000 sr} bdef
-/c3 { 0.000000 0.820000 0.000000 sr} bdef
-/c4 { 0.000000 0.000000 0.800000 sr} bdef
-/c5 { 0.910000 0.820000 0.320000 sr} bdef
-/c6 { 1.000000 0.260000 0.820000 sr} bdef
-/c7 { 0.000000 0.820000 0.820000 sr} bdef
-c0
-1 j
-1 sg
-   0    0 4803 3605 PR
-4.16667 w
-0 2938 3722 0 0 -2938 624 3208 4 MP
-PP
--3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
- 624 3208 mt 4346 3208 L
- 624 3208 mt  624  270 L
- 624 3208 mt  624 3189 L
- 624  270 mt  624  288 L
- 624 3208 mt  624 3170 L
- 624  270 mt  624  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 563 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 655 3283 mt 
-(0) s
- 904 3208 mt  904 3189 L
- 904  270 mt  904  288 L
-1067 3208 mt 1067 3189 L
-1067  270 mt 1067  288 L
-1184 3208 mt 1184 3189 L
-1184  270 mt 1184  288 L
-1274 3208 mt 1274 3189 L
-1274  270 mt 1274  288 L
-1348 3208 mt 1348 3189 L
-1348  270 mt 1348  288 L
-1410 3208 mt 1410 3189 L
-1410  270 mt 1410  288 L
-1464 3208 mt 1464 3189 L
-1464  270 mt 1464  288 L
-1511 3208 mt 1511 3189 L
-1511  270 mt 1511  288 L
-1554 3208 mt 1554 3189 L
-1554  270 mt 1554  288 L
-1554 3208 mt 1554 3170 L
-1554  270 mt 1554  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1493 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-1585 3283 mt 
-(1) s
-1834 3208 mt 1834 3189 L
-1834  270 mt 1834  288 L
-1998 3208 mt 1998 3189 L
-1998  270 mt 1998  288 L
-2114 3208 mt 2114 3189 L
-2114  270 mt 2114  288 L
-2204 3208 mt 2204 3189 L
-2204  270 mt 2204  288 L
-2278 3208 mt 2278 3189 L
-2278  270 mt 2278  288 L
-2340 3208 mt 2340 3189 L
-2340  270 mt 2340  288 L
-2394 3208 mt 2394 3189 L
-2394  270 mt 2394  288 L
-2442 3208 mt 2442 3189 L
-2442  270 mt 2442  288 L
-2485 3208 mt 2485 3189 L
-2485  270 mt 2485  288 L
-2485 3208 mt 2485 3170 L
-2485  270 mt 2485  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2424 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-2516 3283 mt 
-(2) s
-2765 3208 mt 2765 3189 L
-2765  270 mt 2765  288 L
-2928 3208 mt 2928 3189 L
-2928  270 mt 2928  288 L
-3045 3208 mt 3045 3189 L
-3045  270 mt 3045  288 L
-3135 3208 mt 3135 3189 L
-3135  270 mt 3135  288 L
-3209 3208 mt 3209 3189 L
-3209  270 mt 3209  288 L
-3271 3208 mt 3271 3189 L
-3271  270 mt 3271  288 L
-3325 3208 mt 3325 3189 L
-3325  270 mt 3325  288 L
-3372 3208 mt 3372 3189 L
-3372  270 mt 3372  288 L
-3415 3208 mt 3415 3189 L
-3415  270 mt 3415  288 L
-3415 3208 mt 3415 3170 L
-3415  270 mt 3415  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-3354 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-3446 3283 mt 
-(3) s
-3695 3208 mt 3695 3189 L
-3695  270 mt 3695  288 L
-3859 3208 mt 3859 3189 L
-3859  270 mt 3859  288 L
-3975 3208 mt 3975 3189 L
-3975  270 mt 3975  288 L
-4065 3208 mt 4065 3189 L
-4065  270 mt 4065  288 L
-4139 3208 mt 4139 3189 L
-4139  270 mt 4139  288 L
-4201 3208 mt 4201 3189 L
-4201  270 mt 4201  288 L
-4255 3208 mt 4255 3189 L
-4255  270 mt 4255  288 L
-4303 3208 mt 4303 3189 L
-4303  270 mt 4303  288 L
-4346 3208 mt 4346 3189 L
-4346  270 mt 4346  288 L
-4346 3208 mt 4346 3170 L
-4346  270 mt 4346  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-4285 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-4377 3283 mt 
-(4) s
- 624 3208 mt  642 3208 L
-4346 3208 mt 4327 3208 L
- 624 3208 mt  661 3208 L
-4346 3208 mt 4308 3208 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 3238 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 3187 mt 
-(-7) s
- 624 3097 mt  642 3097 L
-4346 3097 mt 4327 3097 L
- 624 3032 mt  642 3032 L
-4346 3032 mt 4327 3032 L
- 624 2986 mt  642 2986 L
-4346 2986 mt 4327 2986 L
- 624 2951 mt  642 2951 L
-4346 2951 mt 4327 2951 L
- 624 2922 mt  642 2922 L
-4346 2922 mt 4327 2922 L
- 624 2897 mt  642 2897 L
-4346 2897 mt 4327 2897 L
- 624 2876 mt  642 2876 L
-4346 2876 mt 4327 2876 L
- 624 2857 mt  642 2857 L
-4346 2857 mt 4327 2857 L
- 624 2840 mt  642 2840 L
-4346 2840 mt 4327 2840 L
- 624 2840 mt  661 2840 L
-4346 2840 mt 4308 2840 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 2870 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 2819 mt 
-(-6) s
- 624 2730 mt  642 2730 L
-4346 2730 mt 4327 2730 L
- 624 2665 mt  642 2665 L
-4346 2665 mt 4327 2665 L
- 624 2619 mt  642 2619 L
-4346 2619 mt 4327 2619 L
- 624 2584 mt  642 2584 L
-4346 2584 mt 4327 2584 L
- 624 2554 mt  642 2554 L
-4346 2554 mt 4327 2554 L
- 624 2530 mt  642 2530 L
-4346 2530 mt 4327 2530 L
- 624 2509 mt  642 2509 L
-4346 2509 mt 4327 2509 L
- 624 2490 mt  642 2490 L
-4346 2490 mt 4327 2490 L
- 624 2473 mt  642 2473 L
-4346 2473 mt 4327 2473 L
- 624 2473 mt  661 2473 L
-4346 2473 mt 4308 2473 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 2503 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 2452 mt 
-(-5) s
- 624 2362 mt  642 2362 L
-4346 2362 mt 4327 2362 L
- 624 2298 mt  642 2298 L
-4346 2298 mt 4327 2298 L
- 624 2252 mt  642 2252 L
-4346 2252 mt 4327 2252 L
- 624 2216 mt  642 2216 L
-4346 2216 mt 4327 2216 L
- 624 2187 mt  642 2187 L
-4346 2187 mt 4327 2187 L
- 624 2163 mt  642 2163 L
-4346 2163 mt 4327 2163 L
- 624 2141 mt  642 2141 L
-4346 2141 mt 4327 2141 L
- 624 2123 mt  642 2123 L
-4346 2123 mt 4327 2123 L
- 624 2106 mt  642 2106 L
-4346 2106 mt 4327 2106 L
- 624 2106 mt  661 2106 L
-4346 2106 mt 4308 2106 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 2136 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 2085 mt 
-(-4) s
- 624 1995 mt  642 1995 L
-4346 1995 mt 4327 1995 L
- 624 1931 mt  642 1931 L
-4346 1931 mt 4327 1931 L
- 624 1885 mt  642 1885 L
-4346 1885 mt 4327 1885 L
- 624 1849 mt  642 1849 L
-4346 1849 mt 4327 1849 L
- 624 1820 mt  642 1820 L
-4346 1820 mt 4327 1820 L
- 624 1795 mt  642 1795 L
-4346 1795 mt 4327 1795 L
- 624 1774 mt  642 1774 L
-4346 1774 mt 4327 1774 L
- 624 1755 mt  642 1755 L
-4346 1755 mt 4327 1755 L
- 624 1739 mt  642 1739 L
-4346 1739 mt 4327 1739 L
- 624 1739 mt  661 1739 L
-4346 1739 mt 4308 1739 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 1769 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 1718 mt 
-(-3) s
- 624 1628 mt  642 1628 L
-4346 1628 mt 4327 1628 L
- 624 1563 mt  642 1563 L
-4346 1563 mt 4327 1563 L
- 624 1517 mt  642 1517 L
-4346 1517 mt 4327 1517 L
- 624 1482 mt  642 1482 L
-4346 1482 mt 4327 1482 L
- 624 1453 mt  642 1453 L
-4346 1453 mt 4327 1453 L
- 624 1428 mt  642 1428 L
-4346 1428 mt 4327 1428 L
- 624 1407 mt  642 1407 L
-4346 1407 mt 4327 1407 L
- 624 1388 mt  642 1388 L
-4346 1388 mt 4327 1388 L
- 624 1371 mt  642 1371 L
-4346 1371 mt 4327 1371 L
- 624 1371 mt  661 1371 L
-4346 1371 mt 4308 1371 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 1401 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537 1350 mt 
-(-2) s
- 624 1261 mt  642 1261 L
-4346 1261 mt 4327 1261 L
- 624 1196 mt  642 1196 L
-4346 1196 mt 4327 1196 L
- 624 1150 mt  642 1150 L
-4346 1150 mt 4327 1150 L
- 624 1115 mt  642 1115 L
-4346 1115 mt 4327 1115 L
- 624 1085 mt  642 1085 L
-4346 1085 mt 4327 1085 L
- 624 1061 mt  642 1061 L
-4346 1061 mt 4327 1061 L
- 624 1040 mt  642 1040 L
-4346 1040 mt 4327 1040 L
- 624 1021 mt  642 1021 L
-4346 1021 mt 4327 1021 L
- 624 1004 mt  642 1004 L
-4346 1004 mt 4327 1004 L
- 624 1004 mt  661 1004 L
-4346 1004 mt 4308 1004 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445 1034 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537  983 mt 
-(-1) s
- 624  893 mt  642  893 L
-4346  893 mt 4327  893 L
- 624  829 mt  642  829 L
-4346  829 mt 4327  829 L
- 624  783 mt  642  783 L
-4346  783 mt 4327  783 L
- 624  747 mt  642  747 L
-4346  747 mt 4327  747 L
- 624  718 mt  642  718 L
-4346  718 mt 4327  718 L
- 624  694 mt  642  694 L
-4346  694 mt 4327  694 L
- 624  672 mt  642  672 L
-4346  672 mt 4327  672 L
- 624  654 mt  642  654 L
-4346  654 mt 4327  654 L
- 624  637 mt  642  637 L
-4346  637 mt 4327  637 L
- 624  637 mt  661  637 L
-4346  637 mt 4308  637 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445  667 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537  616 mt 
-(0) s
- 624  526 mt  642  526 L
-4346  526 mt 4327  526 L
- 624  462 mt  642  462 L
-4346  462 mt 4327  462 L
- 624  416 mt  642  416 L
-4346  416 mt 4327  416 L
- 624  380 mt  642  380 L
-4346  380 mt 4327  380 L
- 624  351 mt  642  351 L
-4346  351 mt 4327  351 L
- 624  326 mt  642  326 L
-4346  326 mt 4327  326 L
- 624  305 mt  642  305 L
-4346  305 mt 4327  305 L
- 624  286 mt  642  286 L
-4346  286 mt 4327  286 L
- 624  270 mt  642  270 L
-4346  270 mt 4327  270 L
- 624  270 mt  661  270 L
-4346  270 mt 4308  270 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 445  300 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 537  249 mt 
-(1) s
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
-gs 624 269 3723 2939 MR c np
-/c8 { 0.000000 0.000000 1.000000 sr} bdef
-c8
-68 111 100 0 65 110 95 111 107 110 100 111 97 110 93 111 
-87 111 72 0 55 110 84 111 55 110 64 111 60 0 60 110 
-62 111 45 110 62 0 48 111 51 0 71 110 51 111 90 0 
-71 111 104 0 122 110 190 111 254 110 560 0 1184 637 31 MP stroke
-gs 1133 586 3146 2535 MR c np
-  25   25 1184  637 FO
-  25   25 1744  637 FO
-  25   25 1998  747 FO
-  25   25 2188  858 FO
-  25   25 2310  968 FO
-  25   25 2414  968 FO
-  25   25 2485 1079 FO
-  25   25 2575 1079 FO
-  25   25 2626 1190 FO
-  25   25 2697 1300 FO
-  25   25 2748 1300 FO
-  25   25 2796 1411 FO
-  25   25 2858 1411 FO
-  25   25 2903 1521 FO
-  25   25 2965 1632 FO
-  25   25 3025 1742 FO
-  25   25 3085 1742 FO
-  25   25 3149 1853 FO
-  25   25 3204 1963 FO
-  25   25 3288 2074 FO
-  25   25 3343 2184 FO
-  25   25 3415 2184 FO
-  25   25 3502 2295 FO
-  25   25 3595 2406 FO
-  25   25 3692 2516 FO
-  25   25 3792 2627 FO
-  25   25 3899 2737 FO
-  25   25 3994 2848 FO
-  25   25 4059 2958 FO
-  25   25 4159 2958 FO
-  25   25 4227 3069 FO
-gr
-
-/c9 { 0.000000 0.500000 0.000000 sr} bdef
-c9
-68 0 100 0 65 111 95 110 107 111 100 0 97 110 93 0 
-87 111 72 110 55 0 84 111 55 0 64 0 60 110 60 0 
-62 111 45 0 62 110 48 0 51 0 71 111 51 0 90 111 
-71 0 104 0 122 110 190 111 254 110 560 0 1184 637 31 MP stroke
-gs 1133 586 3146 1761 MR c np
-  25   25 1184  637 FO
-  25   25 1744  637 FO
-  25   25 1998  747 FO
-  25   25 2188  858 FO
-  25   25 2310  968 FO
-  25   25 2414  968 FO
-  25   25 2485  968 FO
-  25   25 2575 1079 FO
-  25   25 2626 1079 FO
-  25   25 2697 1190 FO
-  25   25 2748 1190 FO
-  25   25 2796 1190 FO
-  25   25 2858 1300 FO
-  25   25 2903 1300 FO
-  25   25 2965 1411 FO
-  25   25 3025 1411 FO
-  25   25 3085 1521 FO
-  25   25 3149 1521 FO
-  25   25 3204 1521 FO
-  25   25 3288 1632 FO
-  25   25 3343 1632 FO
-  25   25 3415 1742 FO
-  25   25 3502 1853 FO
-  25   25 3595 1853 FO
-  25   25 3692 1963 FO
-  25   25 3792 1963 FO
-  25   25 3899 2074 FO
-  25   25 3994 2184 FO
-  25   25 4059 2295 FO
-  25   25 4159 2295 FO
-  25   25 4227 2295 FO
-gr
-
-/c10 { 1.000000 0.000000 0.000000 sr} bdef
-c10
-68 110 100 0 65 111 95 0 107 110 100 111 97 110 93 111 
-87 111 72 0 55 110 84 111 55 110 64 111 60 0 60 110 
-62 111 45 110 62 0 48 111 51 0 71 110 51 111 90 0 
-71 111 104 110 122 111 190 0 254 110 560 0 1184 637 31 MP stroke
-gs 1133 586 3146 2424 MR c np
-  25   25 1184  637 FO
-  25   25 1744  637 FO
-  25   25 1998  747 FO
-  25   25 2188  747 FO
-  25   25 2310  858 FO
-  25   25 2414  968 FO
-  25   25 2485 1079 FO
-  25   25 2575 1079 FO
-  25   25 2626 1190 FO
-  25   25 2697 1300 FO
-  25   25 2748 1300 FO
-  25   25 2796 1411 FO
-  25   25 2858 1411 FO
-  25   25 2903 1521 FO
-  25   25 2965 1632 FO
-  25   25 3025 1742 FO
-  25   25 3085 1742 FO
-  25   25 3149 1853 FO
-  25   25 3204 1963 FO
-  25   25 3288 2074 FO
-  25   25 3343 2184 FO
-  25   25 3415 2184 FO
-  25   25 3502 2295 FO
-  25   25 3595 2406 FO
-  25   25 3692 2516 FO
-  25   25 3792 2627 FO
-  25   25 3899 2737 FO
-  25   25 3994 2737 FO
-  25   25 4059 2848 FO
-  25   25 4159 2848 FO
-  25   25 4227 2958 FO
-gr
-
-DD
-c8
-68 14 100 20 65 13 95 18 107 21 100 20 97 19 93 19 
-87 17 72 14 55 11 84 16 55 11 64 13 60 12 60 12 
-62 12 45 9 62 12 48 9 51 11 71 13 51 11 90 17 
-71 14 104 21 122 24 190 37 254 51 560 110 1184 437 31 MP stroke
-c9
-68 7 100 10 65 6 95 9 107 11 100 10 97 9 93 9 
-87 9 72 7 55 6 84 8 55 5 64 7 60 6 60 5 
-62 7 45 4 62 6 48 5 51 5 71 7 51 5 90 9 
-71 7 104 10 122 12 190 19 254 25 560 55 1184 517 31 MP stroke
-c10
-68 20 100 30 65 19 95 28 107 32 100 29 97 29 93 28 
-87 26 72 21 55 16 84 25 55 16 64 19 60 18 60 18 
-62 18 45 13 62 19 48 14 51 15 71 21 51 15 90 27 
-71 21 104 30 122 37 190 56 254 75 560 166 1184 692 31 MP stroke
-gr
-
-c10
-DD
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2280  205 mt 
-(hmin hmax) s
-2308 3434 mt 
-(Elemente) s
- 388 1927 mt  -90 rotate
-(Schaetzer) s
-90 rotate
- 612 3238 mt 
-( ) s
-4335  299 mt 
-( ) s
-SO
-1 sg
-0 714 1792 0 0 -714 665 3167 4 MP
-PP
--1792 0 0 714 1792 0 0 -714 665 3167 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
- 665 3167 mt 2457 3167 L
- 665 3167 mt  665 2453 L
- 989 2548 mt 
-(min hmin/max hmax testAA Analytisch) s
-gs 665 2453 1793 715 MR c np
-c8
-249 0 714 2517 2 MP stroke
-gs 788 2466 103 103 MR c np
-  25   25  839 2517 FO
-gr
-
-gr
-
-c8
-0 sg
- 989 2651 mt 
-(min hmax/max hmax testAA Analytisch) s
-gs 665 2453 1793 715 MR c np
-c9
-249 0 714 2621 2 MP stroke
-gs 788 2570 103 103 MR c np
-  25   25  839 2621 FO
-gr
-
-gr
-
-c9
-0 sg
- 989 2754 mt 
-(hmin/hmax testAA Analytisch) s
-gs 665 2453 1793 715 MR c np
-c10
-249 0 714 2724 2 MP stroke
-gs 788 2673 103 103 MR c np
-  25   25  839 2724 FO
-gr
-
-gr
-
-c10
-0 sg
- 989 2879 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1049 2838 mt 
-(-1/2) s
-gs 665 2453 1793 715 MR c np
-DD
-c8
-249 0 714 2838 2 MP stroke
-SO
-gr
-
-c8
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 989 3005 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1049 2964 mt 
-(-1/4) s
-gs 665 2453 1793 715 MR c np
-DD
-c9
-249 0 714 2964 2 MP stroke
-SO
-gr
-
-c9
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 989 3132 mt 
-(N) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
-1049 3091 mt 
-(-3/4) s
-gs 665 2453 1793 715 MR c np
-DD
-c10
-249 0 714 3090 2 MP stroke
-SO
-gr
-
-c10
-
-end %%Color Dict
-
-eplot
-%%EndObject
-
-epage
-end
-
-showpage
-
-%%Trailer
-%%EOF
diff --git a/src/expmlAA_2DQuad_norm.eps b/src/expmlAA_2DQuad_norm.eps

deleted file mode 100644 (file)

index d50751f..0000000
--- a/src/expmlAA_2DQuad_norm.eps
+++ /dev/null
@@ -1,552 +0,0 @@
-%!PS-Adobe-2.0 EPSF-1.2
-%%Creator: MATLAB, The MathWorks, Inc. Version 7.12.0.635 (R2011a). Operating System: Linux 3.0.0-16-server #28-Ubuntu SMP Fri Jan 27 18:03:45 UTC 2012 x86_64.
-%%Title: ./expmlAA_2DQuad_norm.eps
-%%CreationDate: 03/26/2012  23:43:44
-%%DocumentNeededFonts: Helvetica
-%%DocumentProcessColors: Cyan Magenta Yellow Black
-%%Extensions: CMYK
-%%Pages: 1
-%%BoundingBox:    44   195   550   602
-%%EndComments
-
-%%BeginProlog
-% MathWorks dictionary
-/MathWorks 160 dict begin
-% definition operators
-/bdef {bind def} bind def
-/ldef {load def} bind def
-/xdef {exch def} bdef
-/xstore {exch store} bdef
-% operator abbreviations
-/c  /clip ldef
-/cc /concat ldef
-/cp /closepath ldef
-/gr /grestore ldef
-/gs /gsave ldef
-/mt /moveto ldef
-/np /newpath ldef
-/cm /currentmatrix ldef
-/sm /setmatrix ldef
-/rm /rmoveto ldef
-/rl /rlineto ldef
-/s {show newpath} bdef
-/sc {setcmykcolor} bdef
-/sr /setrgbcolor ldef
-/sg /setgray ldef
-/w /setlinewidth ldef
-/j /setlinejoin ldef
-/cap /setlinecap ldef
-/rc {rectclip} bdef
-/rf {rectfill} bdef
-% page state control
-/pgsv () def
-/bpage {/pgsv save def} bdef
-/epage {pgsv restore} bdef
-/bplot /gsave ldef
-/eplot {stroke grestore} bdef
-% orientation switch
-/portraitMode 0 def /landscapeMode 1 def /rotateMode 2 def
-% coordinate system mappings
-/dpi2point 0 def
-% font control
-/FontSize 0 def
-/FMS {/FontSize xstore findfont [FontSize 0 0 FontSize neg 0 0]
-  makefont setfont} bdef
-/reencode {exch dup where {pop load} {pop StandardEncoding} ifelse
-  exch dup 3 1 roll findfont dup length dict begin
-  { 1 index /FID ne {def}{pop pop} ifelse } forall
-  /Encoding exch def currentdict end definefont pop} bdef
-/isroman {findfont /CharStrings get /Agrave known} bdef
-/FMSR {3 1 roll 1 index dup isroman {reencode} {pop pop} ifelse
-  exch FMS} bdef
-/csm {1 dpi2point div -1 dpi2point div scale neg translate
- dup landscapeMode eq {pop -90 rotate}
-  {rotateMode eq {90 rotate} if} ifelse} bdef
-% line types: solid, dotted, dashed, dotdash
-/SO { [] 0 setdash } bdef
-/DO { [.5 dpi2point mul 4 dpi2point mul] 0 setdash } bdef
-/DA { [6 dpi2point mul] 0 setdash } bdef
-/DD { [.5 dpi2point mul 4 dpi2point mul 6 dpi2point mul 4
-  dpi2point mul] 0 setdash } bdef
-% macros for lines and objects
-/L {lineto stroke} bdef
-/MP {3 1 roll moveto 1 sub {rlineto} repeat} bdef
-/AP {{rlineto} repeat} bdef
-/PDlw -1 def
-/W {/PDlw currentlinewidth def setlinewidth} def
-/PP {closepath eofill} bdef
-/DP {closepath stroke} bdef
-/MR {4 -2 roll moveto dup  0 exch rlineto exch 0 rlineto
-  neg 0 exch rlineto closepath} bdef
-/FR {MR stroke} bdef
-/PR {MR fill} bdef
-/L1i {{currentfile picstr readhexstring pop} image} bdef
-/tMatrix matrix def
-/MakeOval {newpath tMatrix currentmatrix pop translate scale
-0 0 1 0 360 arc tMatrix setmatrix} bdef
-/FO {MakeOval stroke} bdef
-/PO {MakeOval fill} bdef
-/PD {currentlinewidth 2 div 0 360 arc fill
-   PDlw -1 eq not {PDlw w /PDlw -1 def} if} def
-/FA {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arc tMatrix setmatrix stroke} bdef
-/PA {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arc closepath tMatrix setmatrix fill} bdef
-/FAn {newpath tMatrix currentmatrix pop translate scale
-  0 0 1 5 -2 roll arcn tMatrix setmatrix stroke} bdef
-/PAn {newpath tMatrix currentmatrix pop translate 0 0 moveto scale
-  0 0 1 5 -2 roll arcn closepath tMatrix setmatrix fill} bdef
-/vradius 0 def /hradius 0 def /lry 0 def
-/lrx 0 def /uly 0 def /ulx 0 def /rad 0 def
-/MRR {/vradius xdef /hradius xdef /lry xdef /lrx xdef /uly xdef
-  /ulx xdef newpath tMatrix currentmatrix pop ulx hradius add uly
-  vradius add translate hradius vradius scale 0 0 1 180 270 arc 
-  tMatrix setmatrix lrx hradius sub uly vradius add translate
-  hradius vradius scale 0 0 1 270 360 arc tMatrix setmatrix
-  lrx hradius sub lry vradius sub translate hradius vradius scale
-  0 0 1 0 90 arc tMatrix setmatrix ulx hradius add lry vradius sub
-  translate hradius vradius scale 0 0 1 90 180 arc tMatrix setmatrix
-  closepath} bdef
-/FRR {MRR stroke } bdef
-/PRR {MRR fill } bdef
-/MlrRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lry uly sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 90 270 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 270 90 arc tMatrix setmatrix
-  closepath} bdef
-/FlrRR {MlrRR stroke } bdef
-/PlrRR {MlrRR fill } bdef
-/MtbRR {/lry xdef /lrx xdef /uly xdef /ulx xdef /rad lrx ulx sub 2 div def
-  newpath tMatrix currentmatrix pop ulx rad add uly rad add translate
-  rad rad scale 0 0 1 180 360 arc tMatrix setmatrix lrx rad sub lry rad
-  sub translate rad rad scale 0 0 1 0 180 arc tMatrix setmatrix
-  closepath} bdef
-/FtbRR {MtbRR stroke } bdef
-/PtbRR {MtbRR fill } bdef
-/stri 6 array def /dtri 6 array def
-/smat 6 array def /dmat 6 array def
-/tmat1 6 array def /tmat2 6 array def /dif 3 array def
-/asub {/ind2 exch def /ind1 exch def dup dup
-  ind1 get exch ind2 get sub exch } bdef
-/tri_to_matrix {
-  2 0 asub 3 1 asub 4 0 asub 5 1 asub
-  dup 0 get exch 1 get 7 -1 roll astore } bdef
-/compute_transform {
-  dmat dtri tri_to_matrix tmat1 invertmatrix 
-  smat stri tri_to_matrix tmat2 concatmatrix } bdef
-/ds {stri astore pop} bdef
-/dt {dtri astore pop} bdef
-/db {2 copy /cols xdef /rows xdef mul dup 3 mul string
-  currentfile exch readhexstring pop
-  dup 0 3 index getinterval /rbmap xdef
-  dup 2 index dup getinterval /gbmap xdef
-  1 index dup 2 mul exch getinterval /bbmap xdef pop pop}bdef
-/it {gs np dtri aload pop moveto lineto lineto cp c
-  cols rows 8 compute_transform 
-  rbmap gbmap bbmap true 3 colorimage gr}bdef
-/il {newpath moveto lineto stroke}bdef
-currentdict end def
-%%EndProlog
-
-%%BeginSetup
-MathWorks begin
-
-0 cap
-
-end
-%%EndSetup
-
-%%Page: 1 1
-%%BeginPageSetup
-%%PageBoundingBox:    44   195   550   602
-MathWorks begin
-bpage
-%%EndPageSetup
-
-%%BeginObject: obj1
-bplot
-
-/dpi2point 8.33333 def
-portraitMode 0150 5100 csm
-
-  223    79  4218  3390 MR c np
-76 dict begin %Colortable dictionary
-/c0 { 0.000000 0.000000 0.000000 sr} bdef
-/c1 { 1.000000 1.000000 1.000000 sr} bdef
-/c2 { 0.900000 0.000000 0.000000 sr} bdef
-/c3 { 0.000000 0.820000 0.000000 sr} bdef
-/c4 { 0.000000 0.000000 0.800000 sr} bdef
-/c5 { 0.910000 0.820000 0.320000 sr} bdef
-/c6 { 1.000000 0.260000 0.820000 sr} bdef
-/c7 { 0.000000 0.820000 0.820000 sr} bdef
-c0
-1 j
-1 sg
-   0    0 4803 3605 PR
-4.16667 w
-0 2938 3722 0 0 -2938 624 3208 4 MP
-PP
--3722 0 0 2938 3722 0 0 -2938 624 3208 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
- 624 3208 mt 4346 3208 L
- 624 3208 mt  624  270 L
- 624 3208 mt  624 3189 L
- 624  270 mt  624  288 L
- 624 3208 mt  624 3170 L
- 624  270 mt  624  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 563 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 655 3283 mt 
-(0) s
- 904 3208 mt  904 3189 L
- 904  270 mt  904  288 L
-1067 3208 mt 1067 3189 L
-1067  270 mt 1067  288 L
-1184 3208 mt 1184 3189 L
-1184  270 mt 1184  288 L
-1274 3208 mt 1274 3189 L
-1274  270 mt 1274  288 L
-1348 3208 mt 1348 3189 L
-1348  270 mt 1348  288 L
-1410 3208 mt 1410 3189 L
-1410  270 mt 1410  288 L
-1464 3208 mt 1464 3189 L
-1464  270 mt 1464  288 L
-1511 3208 mt 1511 3189 L
-1511  270 mt 1511  288 L
-1554 3208 mt 1554 3189 L
-1554  270 mt 1554  288 L
-1554 3208 mt 1554 3170 L
-1554  270 mt 1554  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-1493 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-1585 3283 mt 
-(1) s
-1834 3208 mt 1834 3189 L
-1834  270 mt 1834  288 L
-1998 3208 mt 1998 3189 L
-1998  270 mt 1998  288 L
-2114 3208 mt 2114 3189 L
-2114  270 mt 2114  288 L
-2204 3208 mt 2204 3189 L
-2204  270 mt 2204  288 L
-2278 3208 mt 2278 3189 L
-2278  270 mt 2278  288 L
-2340 3208 mt 2340 3189 L
-2340  270 mt 2340  288 L
-2394 3208 mt 2394 3189 L
-2394  270 mt 2394  288 L
-2442 3208 mt 2442 3189 L
-2442  270 mt 2442  288 L
-2485 3208 mt 2485 3189 L
-2485  270 mt 2485  288 L
-2485 3208 mt 2485 3170 L
-2485  270 mt 2485  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2424 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-2516 3283 mt 
-(2) s
-2765 3208 mt 2765 3189 L
-2765  270 mt 2765  288 L
-2928 3208 mt 2928 3189 L
-2928  270 mt 2928  288 L
-3045 3208 mt 3045 3189 L
-3045  270 mt 3045  288 L
-3135 3208 mt 3135 3189 L
-3135  270 mt 3135  288 L
-3209 3208 mt 3209 3189 L
-3209  270 mt 3209  288 L
-3271 3208 mt 3271 3189 L
-3271  270 mt 3271  288 L
-3325 3208 mt 3325 3189 L
-3325  270 mt 3325  288 L
-3372 3208 mt 3372 3189 L
-3372  270 mt 3372  288 L
-3415 3208 mt 3415 3189 L
-3415  270 mt 3415  288 L
-3415 3208 mt 3415 3170 L
-3415  270 mt 3415  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-3354 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-3446 3283 mt 
-(3) s
-3695 3208 mt 3695 3189 L
-3695  270 mt 3695  288 L
-3859 3208 mt 3859 3189 L
-3859  270 mt 3859  288 L
-3975 3208 mt 3975 3189 L
-3975  270 mt 3975  288 L
-4065 3208 mt 4065 3189 L
-4065  270 mt 4065  288 L
-4139 3208 mt 4139 3189 L
-4139  270 mt 4139  288 L
-4201 3208 mt 4201 3189 L
-4201  270 mt 4201  288 L
-4255 3208 mt 4255 3189 L
-4255  270 mt 4255  288 L
-4303 3208 mt 4303 3189 L
-4303  270 mt 4303  288 L
-4346 3208 mt 4346 3189 L
-4346  270 mt 4346  288 L
-4346 3208 mt 4346 3170 L
-4346  270 mt 4346  307 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-4285 3334 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
-4377 3283 mt 
-(4) s
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 2328 mt  661 2328 L
-4346 2328 mt 4308 2328 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 400 2358 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 492 2307 mt 
-(0.65) s
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 1024 mt  661 1024 L
-4346 1024 mt 4308 1024 L
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 400 1054 mt 
-(10) s
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 55.5556 FMSR
-
- 492 1003 mt 
-(0.66) s
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt  624 3208 L
- 624 3208 mt 4346 3208 L
- 624  270 mt 4346  270 L
- 624 3208 mt  624  270 L
-4346 3208 mt 4346  270 L
-gs 624 270 3723 2939 MR c np
-/c8 { 0.000000 0.000000 1.000000 sr} bdef
-c8
-68 0 100 0 65 0 95 0 107 0 100 0 97 0 93 0 
-87 0 72 -1 55 0 84 -1 55 -1 64 -1 60 -3 60 -3 
-62 -5 45 -7 62 -10 48 -14 51 -19 71 -28 51 -38 90 -53 
-71 -82 104 -164 122 -184 190 -333 254 -480 560 -1229 1184 3208 31 MP stroke
-gr
-
-c8
-gs 1133 501 3146 2759 MR c np
-  25   25 1184 3208 FO
-  25   25 1744 1979 FO
-  25   25 1998 1499 FO
-  25   25 2188 1166 FO
-  25   25 2310  982 FO
-  25   25 2414  818 FO
-  25   25 2485  736 FO
-  25   25 2575  683 FO
-  25   25 2626  645 FO
-  25   25 2697  617 FO
-  25   25 2748  598 FO
-  25   25 2796  584 FO
-  25   25 2858  574 FO
-  25   25 2903  567 FO
-  25   25 2965  562 FO
-  25   25 3025  559 FO
-  25   25 3085  556 FO
-  25   25 3149  555 FO
-  25   25 3204  554 FO
-  25   25 3288  553 FO
-  25   25 3343  553 FO
-  25   25 3415  552 FO
-  25   25 3502  552 FO
-  25   25 3595  552 FO
-  25   25 3692  552 FO
-  25   25 3792  552 FO
-  25   25 3899  552 FO
-  25   25 3994  552 FO
-  25   25 4059  552 FO
-  25   25 4159  552 FO
-  25   25 4227  552 FO
-gr
-
-gs 624 270 3723 2939 MR c np
-DD
-/c9 { 1.000000 0.000000 0.000000 sr} bdef
-c9
-68 0 100 0 65 0 95 0 107 0 100 0 97 0 93 0 
-87 0 72 0 55 0 84 0 55 0 64 0 60 0 60 0 
-62 0 45 0 62 0 48 0 51 0 71 0 51 0 90 0 
-71 0 104 0 122 0 190 0 254 0 560 0 1184 547 31 MP stroke
-gr
-
-c9
-DD
-0 sg
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
-2226  206 mt 
-(Energie Norm) s
-2308 3434 mt 
-(Elemente) s
- 343 1883 mt  -90 rotate
-(eNorm) s
-90 rotate
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 66.6667 FMSR
-
- 302 1631 mt  -90 rotate
-(2) s
-90 rotate
-%%IncludeResource: font Helvetica
-/Helvetica /ISOLatin1Encoding 83.3333 FMSR
-
- 612 3239 mt 
-( ) s
-4335  300 mt 
-( ) s
-SO
-1 sg
-0 232 1000 0 0 -232 3305 3168 4 MP
-PP
--1000 0 0 232 1000 0 0 -232 3305 3168 5 MP stroke
-2.77778 w
-DO
-SO
-4.16667 w
-0 sg
-3305 3168 mt 4305 3168 L
-3305 2936 mt 4305 2936 L
-3305 3168 mt 3305 2936 L
-4305 3168 mt 4305 2936 L
-3305 3168 mt 4305 3168 L
-3305 3168 mt 3305 2936 L
-3305 3168 mt 4305 3168 L
-3305 2936 mt 4305 2936 L
-3305 3168 mt 3305 2936 L
-4305 3168 mt 4305 2936 L
-3628 3031 mt 
-(testAA Analytisch) s
-gs 3305 2936 1001 233 MR c np
-c8
-248 0 3354 3000 2 MP stroke
-gs 3427 2949 103 103 MR c np
-  25   25 3478 3000 FO
-gr
-
-gr
-
-c8
-0 sg
-3628 3134 mt 
-(extrapoliert) s
-gs 3305 2936 1001 233 MR c np
-DD
-c9
-248 0 3354 3103 2 MP stroke
-SO
-gr
-
-c9
-
-end %%Color Dict
-
-eplot
-%%EndObject
-
-epage
-end
-
-showpage
-
-%%Trailer
-%%EOF
diff --git a/src/mex_build_AU.cpp b/src/mex_build_AU.cpp

deleted file mode 100644 (file)

index 7bc237e..0000000
--- a/src/mex_build_AU.cpp
+++ /dev/null
@@ -1,287 +0,0 @@
-/*\r
- * Art der Berechnung durch Parameter bestimmt...\r
- * Analytisch/Semianalytisch/...\r
- *\r
- *\r
- */\r
-//#include <iostream>\r
-#include <cmath>\r
-//#include <cassert>\r
-#include "mex.h"\r
-/*/#include <stdlib.h>\r
-\r
-//#include "gauss.hpp"\r
-*/\r
-#define M_EPS 10e-8\r
-/*\r
-//#include "tbb/parallel_for.h"\r
-*/\r
-#include "slpRectangle.hpp"\r
-/*\r
-//using namespace std;\r
-//using namespace slpR;\r
-*/\r
-int dimOfVec(double* vec) {\r
-       if (vec[2] != 0)\r
-               return 2;\r
-       else if (vec[1] != 0)\r
-               return 1;\r
-       else if (vec[0] != 0)\r
-               return 0;\r
-       else {/*\r
-//             cerr << "dimOfVec : (" << vec[0] << " " << vec[1] << " " << vec[2]\r
-//                             << ") alle Werte 0 \n";*/\r
-               return -1;\r
-       }\r
-\r
-}\r
-\r
-inline int dimOfThird(int a, int b) {\r
-       return ((-(a + b) % 3) + 3) % 3;\r
-}\r
-\r
-void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {\r
-/*\r
-//     initeQuad();\r
-//     cout << Quad[1].nod[0];\r
-*/\r
-       int i, j, k; //Schleifenindizes\r
-       double tmp; //Zwischenspeicherung der Einzelnen Werte\r
-       int count;\r
-\r
-       //Sicherheitsabfragen zu Datengroessen\r
-       if ((nrhs != 4))\r
-               mexErrMsgTxt(\r
-                               "expected (coordinates(Nx3),elements(Mx4),mu(double),type(int))");\r
-       if (nlhs > 1)\r
-               mexErrMsgTxt("has only one output argument");\r
-\r
-       int cm = mxGetM(prhs[0]);\r
-       int cn = mxGetN(prhs[0]);\r
-       if (cn != 3)\r
-               mexErrMsgTxt("expected coordinates (Nx3)");\r
-//     cout << "  Koordinaten:" << cm << endl;\r
-\r
-       int em = mxGetM(prhs[1]);\r
-       int en = mxGetN(prhs[1]);\r
-       if (en != 4)\r
-               mexErrMsgTxt("expected elements (Mx4)");\r
-//     cout << "  Elemente:" << em << endl;\r
-\r
-       //Auslesen der Parameter\r
-\r
-       plhs[0] = mxCreateDoubleMatrix(em, em, mxREAL);\r
-       double * A = mxGetPr(plhs[0]);\r
-       double * C = mxGetPr(prhs[0]);\r
-       double * E = mxGetPr(prhs[1]);\r
-\r
-       int type = (int) *(mxGetPr(prhs[3]));\r
-       double * mu = mxGetPr(prhs[2]);\r
-\r
-       //Initialisieren der Hilfsvektoren\r
-\r
-       double * x0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * x1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * x2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * x3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * xa = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * xb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-\r
-       double * y0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * y1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * y2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * y3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * ya = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * yb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-\r
-       double * d = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-       double * dt = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
-\r
-       // Welche Funktion soll verwendet werden\r
-       double(*ctypeP)(double, double, double, double, double, double, double,\r
-                       double*);\r
-       double(*ctypeO)(double, double, double, double, double, double, double,\r
-                       double*);\r
-\r
-       //Art der Berechnung bestimmen\r
-//     cout << "  Typ:" << type << endl;\r
-       switch (type) {\r
-       default:\r
-               ctypeP = cParO1;\r
-               ctypeO = NULL;\r
-               break;\r
-       case 0:\r
-               ctypeP = cParO0;\r
-               ctypeO = NULL;\r
-               break;\r
-       case 2:\r
-               ctypeP = cParO2;\r
-               ctypeO = NULL;\r
-               break;\r
-       case 3:\r
-               ctypeP = cParO3;\r
-               ctypeO = NULL;\r
-               break;\r
-       }\r
-\r
-       //LageInformationen\r
-       int rx, rxa, rxb, ry, rya, ryb;\r
-\r
-//     cout << "  Progress: #";\r
-       count = 0;\r
-       //Ausrechnen\r
-       for (j = 0; j < em; ++j) {\r
-               x0[0] = C[(int) E[j] - 1];\r
-               x0[1] = C[cm + (int) E[j] - 1];\r
-               x0[2] = C[2 * cm + (int) E[j] - 1];\r
-\r
-               x1[0] = C[(int) E[em + j] - 1];\r
-               x1[1] = C[cm + (int) E[em + j] - 1];\r
-               x1[2] = C[2 * cm + (int) E[em + j] - 1];\r
-\r
-               x2[0] = C[(int) E[2 * em + j] - 1];\r
-               x2[1] = C[cm + (int) E[2 * em + j] - 1];\r
-               x2[2] = C[2 * cm + (int) E[2 * em + j] - 1];\r
-\r
-               x3[0] = C[(int) E[3 * em + j] - 1];\r
-               x3[1] = C[cm + (int) E[3 * em + j] - 1];\r
-               x3[2] = C[2 * cm + (int) E[3 * em + j] - 1];\r
-\r
-               for (i = 0; i < 3; ++i)\r
-                       xa[i] = x1[i] - x0[i];\r
-\r
-               for (i = 0; i < 3; ++i)\r
-                       xb[i] = x3[i] - x0[i];\r
-\r
-               // Lageeigenschaften des Flaechenstuecks\r
-\r
-               rxa = dimOfVec(xa);\r
-               rxb = dimOfVec(xb);\r
-               rx = dimOfThird(rxa, rxb);\r
-\r
-               //kleinste Ecke finden und fuer \delta verwenden\r
-\r
-               if (xa[rxa] > 0) {\r
-                       if (xb[rxb] > 0) {\r
-                               for (i = 0; i < 3; ++i)\r
-                                       dt[i] = -x0[i];\r
-                       } else {\r
-                               for (i = 0; i < 3; ++i)\r
-                                       dt[i] = -x3[i];\r
-                       }\r
-               } else {\r
-                       if (xb[rxb] > 0) {\r
-                               for (i = 0; i < 3; ++i)\r
-                                       dt[i] = -x1[i];\r
-                       } else {\r
-                               for (i = 0; i < 3; ++i)\r
-                                       dt[i] = -x2[i];\r
-                       }\r
-               }\r
-\r
-               for (k = 0; k <= j; ++k) {\r
-                       y0[0] = C[(int) E[k] - 1];\r
-                       y0[1] = C[cm + (int) E[k] - 1];\r
-                       y0[2] = C[2 * cm + (int) E[k] - 1];\r
-\r
-                       y1[0] = C[(int) E[em + k] - 1];\r
-                       y1[1] = C[cm + (int) E[em + k] - 1];\r
-                       y1[2] = C[2 * cm + (int) E[em + k] - 1];\r
-\r
-                       y2[0] = C[(int) E[2 * em + k] - 1];\r
-                       y2[1] = C[cm + (int) E[2 * em + k] - 1];\r
-                       y2[2] = C[2 * cm + (int) E[2 * em + k] - 1];\r
-\r
-                       y3[0] = C[(int) E[3 * em + k] - 1];\r
-                       y3[1] = C[cm + (int) E[3 * em + k] - 1];\r
-                       y3[2] = C[2 * cm + (int) E[3 * em + k] - 1];\r
-\r
-                       for (i = 0; i < 3; ++i)\r
-                               ya[i] = y1[i] - y0[i];\r
-\r
-                       for (i = 0; i < 3; ++i)\r
-                               yb[i] = y3[i] - y0[i];\r
-\r
-                       // Lageeigenschaften des Flaechenstuecks\r
-\r
-                       rya = dimOfVec(ya);\r
-                       ryb = dimOfVec(yb);\r
-                       ry = dimOfThird(rya, ryb);\r
-\r
-                       //kleinste Ecke finden und fuer \delta verwenden\r
-\r
-                       if (ya[rya] > 0) {\r
-                               if (yb[ryb] > 0) {\r
-                                       for (i = 0; i < 3; ++i)\r
-                                               d[i] = dt[i] + y0[i];\r
-                               } else {\r
-                                       for (i = 0; i < 3; ++i)\r
-                                               d[i] = dt[i] + y3[i];\r
-                               }\r
-                       } else {\r
-                               if (yb[ryb] > 0) {\r
-                                       for (i = 0; i < 3; ++i)\r
-                                               d[i] = dt[i] + y1[i];\r
-                               } else {\r
-                                       for (i = 0; i < 3; ++i)\r
-                                               d[i] = dt[i] + y2[i];\r
-                               }\r
-                       }\r
-\r
-                       if (rx == ry) {\r
-                               if (rxa == rya) {\r
-                                       tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rxa]),\r
-                                                       fabs(yb[rxb]), d[rxa], d[rxb], d[rx], mu);\r
-\r
-                               } else {\r
-                                       tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[rxa]),\r
-                                                       fabs(ya[rxb]), d[rxa], d[rxb], d[rx], mu);\r
-                               }\r
-\r
-                       } else {\r
-\r
-                               if (rxa == rya) {\r
-                                       tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(ya[rya]),\r
-                                                       fabs(yb[ryb]), d[rxb], d[rxa], d[rx], mu);\r
-                               } else if (rxa == ryb) {\r
-                                       tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(yb[ryb]),\r
-                                                       fabs(ya[rya]), d[rxb], d[rxa], d[rx], mu);\r
-                               } else if (rxb == rya) {\r
-                                       tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rya]),\r
-                                                       fabs(yb[ryb]), d[rxa], d[rxb], d[rx], mu);\r
-                               } else {\r
-                                       tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[ryb]),\r
-                                                       fabs(ya[rya]), d[rxa], d[rxb], d[rx], mu);\r
-                               }\r
-\r
-                       }\r
-                       A[(k * em) + j] = 1. / (4 * M_PI) * tmp;\r
-                       if (k != j)\r
-                               A[(j * em) + k] = 1. / (4 * M_PI) * tmp;\r
-/*                     if (count++ > ((em * (em + 1)) / 2) / 10) {\r
-                               count = 0;\r
-//                             mexPrintf("#");\r
-//                             cout << "#";\r
-//                             cout.flush();\r
-                       }*/\r
-\r
-               }\r
-\r
-       }\r
-//     cout << endl;\r
-       //Rueckgabe (eventuell zurueck schreiben)\r
-/*     mxFree(x0);\r
-       mxFree(x1);\r
-       mxFree(x3);\r
-       mxFree(xa);\r
-       mxFree(xb);\r
-       mxFree(y0);\r
-       mxFree(y1);\r
-       mxFree(y3);\r
-       mxFree(ya);\r
-       mxFree(yb);\r
-       mxFree(d);\r
-       mxFree(dt);\r
-*/\r
-       return;\r
-}\r
diff --git a/src/mex_build_V.cpp b/src/mex_build_V.cpp

new file mode 100644 (file)

index 0000000..6f11e60
--- /dev/null
+++ b/src/mex_build_V.cpp
@@ -0,0 +1,290 @@
+/*\r
+ * Art der Berechnung durch Parameter bestimmt...\r
+ * Analytisch/Semianalytisch/...\r
+ *\r
+ *\r
+ */\r
+//#include <iostream>\r
+#include <cmath>\r
+//#include <cassert>\r
+#include "mex.h"\r
+/*/#include <stdlib.h>\r
+\r
+//#include "gauss.hpp"\r
+*/\r
+//#define M_EPS 10e-8\r
+/*\r
+//#include "tbb/parallel_for.h"\r
+*/\r
+#include "slpRectangle.hpp"\r
+/*\r
+//using namespace std;\r
+//using namespace slpR;\r
+*/\r
+int dimOfVec(double* vec) {\r
+       if (vec[2] != 0)\r
+               return 2;\r
+       else if (vec[1] != 0)\r
+               return 1;\r
+       else if (vec[0] != 0)\r
+               return 0;\r
+       else {\r
+               mexErrMsgTxt("length of Site is zero");\r
+               /*\r
+//             cerr << "dimOfVec : (" << vec[0] << " " << vec[1] << " " << vec[2]\r
+//                             << ") alle Werte 0 \n";*/\r
+               return -1;\r
+       }\r
+\r
+}\r
+\r
+// Gibt von [0 1 2] die Fehlende Zahl zurueck\r
+inline int dimOfThird(int a, int b) {\r
+       return ((-(a + b) % 3) + 3) % 3;\r
+}\r
+\r
+void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) {\r
+/*\r
+//     initeQuad();\r
+//     cout << Quad[1].nod[0];\r
+*/\r
+       int i, j, k; //Schleifenindizes\r
+       double tmp; //Zwischenspeicherung der Einzelnen Werte\r
+       int count;\r
+\r
+       //Sicherheitsabfragen zu Datengroessen\r
+       if ((nrhs != 4))\r
+               mexErrMsgTxt(\r
+                               "expected (coordinates(Nx3),elements(Mx4),zeta(double),type(int))");\r
+       if (nlhs > 1)\r
+               mexErrMsgTxt("has only one output argument");\r
+\r
+       int cm = mxGetM(prhs[0]);\r
+       int cn = mxGetN(prhs[0]);\r
+       if (cn != 3)\r
+               mexErrMsgTxt("expected coordinates (Nx3)");\r
+//     cout << "  Koordinaten:" << cm << endl;\r
+\r
+       int em = mxGetM(prhs[1]);\r
+       int en = mxGetN(prhs[1]);\r
+       if (en != 4)\r
+               mexErrMsgTxt("expected elements (Mx4)");\r
+//     cout << "  Elemente:" << em << endl;\r
+\r
+       //Auslesen der Parameter\r
+\r
+       plhs[0] = mxCreateDoubleMatrix(em, em, mxREAL);\r
+       double * A = mxGetPr(plhs[0]);\r
+       double * C = mxGetPr(prhs[0]);\r
+       double * E = mxGetPr(prhs[1]);\r
+\r
+       int type = (int) *(mxGetPr(prhs[3]));\r
+       double * zeta = mxGetPr(prhs[2]);\r
+\r
+       //Initialisieren der Hilfsvektoren\r
+\r
+       double * x0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * x1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * x2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * x3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * xa = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * xb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+\r
+       double * y0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * y1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * y2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * y3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * ya = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * yb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+\r
+       double * d = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+       double * dt = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL));\r
+\r
+       // Welche Funktion soll verwendet werden\r
+       double(*ctypeP)(double, double, double, double, double, double, double,\r
+                       double*);\r
+       double(*ctypeO)(double, double, double, double, double, double, double,\r
+                       double*);\r
+\r
+       //Art der Berechnung bestimmen\r
+//     cout << "  Typ:" << type << endl;\r
+       switch (type) {\r
+       default:\r
+               ctypeP = cParO1;\r
+               ctypeO = NULL;\r
+               break;\r
+       case 0:\r
+               ctypeP = cParO0;\r
+               ctypeO = NULL;\r
+               break;\r
+       case 2:\r
+               ctypeP = cParO2;\r
+               ctypeO = NULL;\r
+               break;\r
+       case 3:\r
+               ctypeP = cParO3;\r
+               ctypeO = NULL;\r
+               break;\r
+       }\r
+\r
+       //LageInformationen\r
+       int rx, rxa, rxb, ry, rya, ryb;\r
+\r
+//     cout << "  Progress: #";\r
+       count = 0;\r
+       //Ausrechnen\r
+       for (j = 0; j < em; ++j) {\r
+               x0[0] = C[(int) E[j] - 1];\r
+               x0[1] = C[cm + (int) E[j] - 1];\r
+               x0[2] = C[2 * cm + (int) E[j] - 1];\r
+\r
+               x1[0] = C[(int) E[em + j] - 1];\r
+               x1[1] = C[cm + (int) E[em + j] - 1];\r
+               x1[2] = C[2 * cm + (int) E[em + j] - 1];\r
+\r
+               x2[0] = C[(int) E[2 * em + j] - 1];\r
+               x2[1] = C[cm + (int) E[2 * em + j] - 1];\r
+               x2[2] = C[2 * cm + (int) E[2 * em + j] - 1];\r
+\r
+               x3[0] = C[(int) E[3 * em + j] - 1];\r
+               x3[1] = C[cm + (int) E[3 * em + j] - 1];\r
+               x3[2] = C[2 * cm + (int) E[3 * em + j] - 1];\r
+\r
+               for (i = 0; i < 3; ++i)\r
+                       xa[i] = x1[i] - x0[i];\r
+\r
+               for (i = 0; i < 3; ++i)\r
+                       xb[i] = x3[i] - x0[i];\r
+\r
+               // Lageeigenschaften des Flaechenstuecks\r
+\r
+               rxa = dimOfVec(xa);\r
+               rxb = dimOfVec(xb);\r
+               rx = dimOfThird(rxa, rxb);\r
+\r
+               //kleinste Ecke finden und fuer \delta verwenden\r
+\r
+               if (xa[rxa] > 0) {\r
+                       if (xb[rxb] > 0) {\r
+                               for (i = 0; i < 3; ++i)\r
+                                       dt[i] = -x0[i];\r
+                       } else {\r
+                               for (i = 0; i < 3; ++i)\r
+                                       dt[i] = -x3[i];\r
+                       }\r
+               } else {\r
+                       if (xb[rxb] > 0) {\r
+                               for (i = 0; i < 3; ++i)\r
+                                       dt[i] = -x1[i];\r
+                       } else {\r
+                               for (i = 0; i < 3; ++i)\r
+                                       dt[i] = -x2[i];\r
+                       }\r
+               }\r
+\r
+               for (k = 0; k <= j; ++k) {\r
+                       y0[0] = C[(int) E[k] - 1];\r
+                       y0[1] = C[cm + (int) E[k] - 1];\r
+                       y0[2] = C[2 * cm + (int) E[k] - 1];\r
+\r
+                       y1[0] = C[(int) E[em + k] - 1];\r
+                       y1[1] = C[cm + (int) E[em + k] - 1];\r
+                       y1[2] = C[2 * cm + (int) E[em + k] - 1];\r
+\r
+                       y2[0] = C[(int) E[2 * em + k] - 1];\r
+                       y2[1] = C[cm + (int) E[2 * em + k] - 1];\r
+                       y2[2] = C[2 * cm + (int) E[2 * em + k] - 1];\r
+\r
+                       y3[0] = C[(int) E[3 * em + k] - 1];\r
+                       y3[1] = C[cm + (int) E[3 * em + k] - 1];\r
+                       y3[2] = C[2 * cm + (int) E[3 * em + k] - 1];\r
+\r
+                       for (i = 0; i < 3; ++i)\r
+                               ya[i] = y1[i] - y0[i];\r
+\r
+                       for (i = 0; i < 3; ++i)\r
+                               yb[i] = y3[i] - y0[i];\r
+\r
+                       // Lageeigenschaften des Flaechenstuecks\r
+\r
+                       rya = dimOfVec(ya);\r
+                       ryb = dimOfVec(yb);\r
+                       ry = dimOfThird(rya, ryb);\r
+\r
+                       //kleinste Ecke finden und fuer \delta verwenden\r
+\r
+                       if (ya[rya] > 0) {\r
+                               if (yb[ryb] > 0) {\r
+                                       for (i = 0; i < 3; ++i)\r
+                                               d[i] = dt[i] + y0[i];\r
+                               } else {\r
+                                       for (i = 0; i < 3; ++i)\r
+                                               d[i] = dt[i] + y3[i];\r
+                               }\r
+                       } else {\r
+                               if (yb[ryb] > 0) {\r
+                                       for (i = 0; i < 3; ++i)\r
+                                               d[i] = dt[i] + y1[i];\r
+                               } else {\r
+                                       for (i = 0; i < 3; ++i)\r
+                                               d[i] = dt[i] + y2[i];\r
+                               }\r
+                       }\r
+\r
+                       if (rx == ry) {\r
+                               if (rxa == rya) {\r
+                                       tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rxa]),\r
+                                                       fabs(yb[rxb]), d[rxa], d[rxb], d[rx], zeta);\r
+\r
+                               } else {\r
+                                       tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[rxa]),\r
+                                                       fabs(ya[rxb]), d[rxa], d[rxb], d[rx], zeta);\r
+                               }\r
+\r
+                       } else {\r
+\r
+                               if (rxa == rya) {\r
+                                       tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(ya[rya]),\r
+                                                       fabs(yb[ryb]), d[rxb], d[rxa], d[rx], zeta);\r
+                               } else if (rxa == ryb) {\r
+                                       tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(yb[ryb]),\r
+                                                       fabs(ya[rya]), d[rxb], d[rxa], d[rx], zeta);\r
+                               } else if (rxb == rya) {\r
+                                       tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rya]),\r
+                                                       fabs(yb[ryb]), d[rxa], d[rxb], d[rx], zeta);\r
+                               } else {\r
+                                       tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[ryb]),\r
+                                                       fabs(ya[rya]), d[rxa], d[rxb], d[rx], zeta);\r
+                               }\r
+\r
+                       }\r
+                       A[(k * em) + j] = 1. / (4 * M_PI) * tmp;\r
+                       if (k != j)\r
+                               A[(j * em) + k] = 1. / (4 * M_PI) * tmp;\r
+/*                     if (count++ > ((em * (em + 1)) / 2) / 10) {\r
+                               count = 0;\r
+//                             mexPrintf("#");\r
+//                             cout << "#";\r
+//                             cout.flush();\r
+                       }*/\r
+\r
+               }\r
+\r
+       }\r
+//     cout << endl;\r
+       //Rueckgabe (eventuell zurueck schreiben)\r
+/*     mxFree(x0);\r
+       mxFree(x1);\r
+       mxFree(x3);\r
+       mxFree(xa);\r
+       mxFree(xb);\r
+       mxFree(y0);\r
+       mxFree(y1);\r
+       mxFree(y3);\r
+       mxFree(ya);\r
+       mxFree(yb);\r
+       mxFree(d);\r
+       mxFree(dt);\r
+*/\r
+       return;\r
+}\r
diff --git a/src/test_sol.m b/src/test_sol.m

index 76bcb4c09d460ad6a1a0f7123605dc685031fd68..23a1fa7ad0ed165a6d2fd09bde923739bb59a852 100644 (file)
--- a/src/test_sol.m
+++ b/src/test_sol.m
@@ -2,7 +2,7 @@
  format longG
  
  % Matrix MEX Funktion neu Compilieren
-mex mex_build_AU.cpp slpRectangle.cpp
+mex mex_build_V.cpp slpRectangle.cpp
  
  % Test ausführen
author	Peter Schaefer <peter.schaefer@tuwien.ac.at>
	Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)
committer	Peter Schaefer <peter.schaefer@tuwien.ac.at>
	Wed, 28 Mar 2012 20:24:37 +0000 (22:24 +0200)
doc/doc.pdf		patch \| blob \| history
doc/doc.tex		patch \| blob \| history
src/A_plot.m	[deleted file]	patch \| blob \| history
src/A_plots.m		patch \| blob \| history
src/A_run.m		patch \| blob \| history
src/A_step.m		patch \| blob \| history
src/Aelements.m	[deleted file]	patch \| blob \| history
src/area.m	[deleted file]	patch \| blob \| history
src/areaQuad.m	[new file with mode: 0644]	patch \| blob
src/computeEstSlpMuTilde.m	[deleted file]	patch \| blob \| history
src/expmlAA_2DQuad_error.eps	[deleted file]	patch \| blob \| history
src/expmlAA_2DQuad_hminmax.eps	[deleted file]	patch \| blob \| history
src/expmlAA_2DQuad_norm.eps	[deleted file]	patch \| blob \| history
src/mex_build_AU.cpp	[deleted file]	patch \| blob \| history
src/mex_build_V.cpp	[new file with mode: 0644]	patch \| blob
src/test_sol.m		patch \| blob \| history