From: Peter Schaefer Date: Wed, 28 Mar 2012 20:24:37 +0000 (+0200) Subject: [src] nutzlose Funktionen entfernt X-Git-Url: https://git.leopard-lacewing.eu/?a=commitdiff_plain;h=7f19cbe523b6408c38fe372b2cdfb11d62d92dd5;p=bacc.git [src] nutzlose Funktionen entfernt [src] area und mex_build_AU umbenannt [src] Variablen hoffentlich komplett einheitlich --- diff --git a/doc/doc.pdf b/doc/doc.pdf index 8597ce2..add9a02 100644 Binary files a/doc/doc.pdf and b/doc/doc.pdf differ diff --git a/doc/doc.tex b/doc/doc.tex index 5f20f08..3d663ab 100644 --- 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 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 5c694fe..0b603b7 100644 --- 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 be5943d..1476eef 100644 --- 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 27ef6d5..71657c1 100644 --- 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 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 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 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 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 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 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 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 index 7bc237e..0000000 --- a/src/mex_build_AU.cpp +++ /dev/null @@ -1,287 +0,0 @@ -/* - * Art der Berechnung durch Parameter bestimmt... - * Analytisch/Semianalytisch/... - * - * - */ -//#include -#include -//#include -#include "mex.h" -/*/#include - -//#include "gauss.hpp" -*/ -#define M_EPS 10e-8 -/* -//#include "tbb/parallel_for.h" -*/ -#include "slpRectangle.hpp" -/* -//using namespace std; -//using namespace slpR; -*/ -int dimOfVec(double* vec) { - if (vec[2] != 0) - return 2; - else if (vec[1] != 0) - return 1; - else if (vec[0] != 0) - return 0; - else {/* -// cerr << "dimOfVec : (" << vec[0] << " " << vec[1] << " " << vec[2] -// << ") alle Werte 0 \n";*/ - return -1; - } - -} - -inline int dimOfThird(int a, int b) { - return ((-(a + b) % 3) + 3) % 3; -} - -void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { -/* -// initeQuad(); -// cout << Quad[1].nod[0]; -*/ - int i, j, k; //Schleifenindizes - double tmp; //Zwischenspeicherung der Einzelnen Werte - int count; - - //Sicherheitsabfragen zu Datengroessen - if ((nrhs != 4)) - mexErrMsgTxt( - "expected (coordinates(Nx3),elements(Mx4),mu(double),type(int))"); - if (nlhs > 1) - mexErrMsgTxt("has only one output argument"); - - int cm = mxGetM(prhs[0]); - int cn = mxGetN(prhs[0]); - if (cn != 3) - mexErrMsgTxt("expected coordinates (Nx3)"); -// cout << " Koordinaten:" << cm << endl; - - int em = mxGetM(prhs[1]); - int en = mxGetN(prhs[1]); - if (en != 4) - mexErrMsgTxt("expected elements (Mx4)"); -// cout << " Elemente:" << em << endl; - - //Auslesen der Parameter - - plhs[0] = mxCreateDoubleMatrix(em, em, mxREAL); - double * A = mxGetPr(plhs[0]); - double * C = mxGetPr(prhs[0]); - double * E = mxGetPr(prhs[1]); - - int type = (int) *(mxGetPr(prhs[3])); - double * mu = mxGetPr(prhs[2]); - - //Initialisieren der Hilfsvektoren - - double * x0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * x1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * x2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * x3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * xa = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * xb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - - double * y0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * y1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * y2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * y3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * ya = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * yb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - - double * d = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - double * dt = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); - - // Welche Funktion soll verwendet werden - double(*ctypeP)(double, double, double, double, double, double, double, - double*); - double(*ctypeO)(double, double, double, double, double, double, double, - double*); - - //Art der Berechnung bestimmen -// cout << " Typ:" << type << endl; - switch (type) { - default: - ctypeP = cParO1; - ctypeO = NULL; - break; - case 0: - ctypeP = cParO0; - ctypeO = NULL; - break; - case 2: - ctypeP = cParO2; - ctypeO = NULL; - break; - case 3: - ctypeP = cParO3; - ctypeO = NULL; - break; - } - - //LageInformationen - int rx, rxa, rxb, ry, rya, ryb; - -// cout << " Progress: #"; - count = 0; - //Ausrechnen - for (j = 0; j < em; ++j) { - x0[0] = C[(int) E[j] - 1]; - x0[1] = C[cm + (int) E[j] - 1]; - x0[2] = C[2 * cm + (int) E[j] - 1]; - - x1[0] = C[(int) E[em + j] - 1]; - x1[1] = C[cm + (int) E[em + j] - 1]; - x1[2] = C[2 * cm + (int) E[em + j] - 1]; - - x2[0] = C[(int) E[2 * em + j] - 1]; - x2[1] = C[cm + (int) E[2 * em + j] - 1]; - x2[2] = C[2 * cm + (int) E[2 * em + j] - 1]; - - x3[0] = C[(int) E[3 * em + j] - 1]; - x3[1] = C[cm + (int) E[3 * em + j] - 1]; - x3[2] = C[2 * cm + (int) E[3 * em + j] - 1]; - - for (i = 0; i < 3; ++i) - xa[i] = x1[i] - x0[i]; - - for (i = 0; i < 3; ++i) - xb[i] = x3[i] - x0[i]; - - // Lageeigenschaften des Flaechenstuecks - - rxa = dimOfVec(xa); - rxb = dimOfVec(xb); - rx = dimOfThird(rxa, rxb); - - //kleinste Ecke finden und fuer \delta verwenden - - if (xa[rxa] > 0) { - if (xb[rxb] > 0) { - for (i = 0; i < 3; ++i) - dt[i] = -x0[i]; - } else { - for (i = 0; i < 3; ++i) - dt[i] = -x3[i]; - } - } else { - if (xb[rxb] > 0) { - for (i = 0; i < 3; ++i) - dt[i] = -x1[i]; - } else { - for (i = 0; i < 3; ++i) - dt[i] = -x2[i]; - } - } - - for (k = 0; k <= j; ++k) { - y0[0] = C[(int) E[k] - 1]; - y0[1] = C[cm + (int) E[k] - 1]; - y0[2] = C[2 * cm + (int) E[k] - 1]; - - y1[0] = C[(int) E[em + k] - 1]; - y1[1] = C[cm + (int) E[em + k] - 1]; - y1[2] = C[2 * cm + (int) E[em + k] - 1]; - - y2[0] = C[(int) E[2 * em + k] - 1]; - y2[1] = C[cm + (int) E[2 * em + k] - 1]; - y2[2] = C[2 * cm + (int) E[2 * em + k] - 1]; - - y3[0] = C[(int) E[3 * em + k] - 1]; - y3[1] = C[cm + (int) E[3 * em + k] - 1]; - y3[2] = C[2 * cm + (int) E[3 * em + k] - 1]; - - for (i = 0; i < 3; ++i) - ya[i] = y1[i] - y0[i]; - - for (i = 0; i < 3; ++i) - yb[i] = y3[i] - y0[i]; - - // Lageeigenschaften des Flaechenstuecks - - rya = dimOfVec(ya); - ryb = dimOfVec(yb); - ry = dimOfThird(rya, ryb); - - //kleinste Ecke finden und fuer \delta verwenden - - if (ya[rya] > 0) { - if (yb[ryb] > 0) { - for (i = 0; i < 3; ++i) - d[i] = dt[i] + y0[i]; - } else { - for (i = 0; i < 3; ++i) - d[i] = dt[i] + y3[i]; - } - } else { - if (yb[ryb] > 0) { - for (i = 0; i < 3; ++i) - d[i] = dt[i] + y1[i]; - } else { - for (i = 0; i < 3; ++i) - d[i] = dt[i] + y2[i]; - } - } - - if (rx == ry) { - if (rxa == rya) { - tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rxa]), - fabs(yb[rxb]), d[rxa], d[rxb], d[rx], mu); - - } else { - tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[rxa]), - fabs(ya[rxb]), d[rxa], d[rxb], d[rx], mu); - } - - } else { - - if (rxa == rya) { - tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(ya[rya]), - fabs(yb[ryb]), d[rxb], d[rxa], d[rx], mu); - } else if (rxa == ryb) { - tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(yb[ryb]), - fabs(ya[rya]), d[rxb], d[rxa], d[rx], mu); - } else if (rxb == rya) { - tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rya]), - fabs(yb[ryb]), d[rxa], d[rxb], d[rx], mu); - } else { - tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[ryb]), - fabs(ya[rya]), d[rxa], d[rxb], d[rx], mu); - } - - } - A[(k * em) + j] = 1. / (4 * M_PI) * tmp; - if (k != j) - A[(j * em) + k] = 1. / (4 * M_PI) * tmp; -/* if (count++ > ((em * (em + 1)) / 2) / 10) { - count = 0; -// mexPrintf("#"); -// cout << "#"; -// cout.flush(); - }*/ - - } - - } -// cout << endl; - //Rueckgabe (eventuell zurueck schreiben) -/* mxFree(x0); - mxFree(x1); - mxFree(x3); - mxFree(xa); - mxFree(xb); - mxFree(y0); - mxFree(y1); - mxFree(y3); - mxFree(ya); - mxFree(yb); - mxFree(d); - mxFree(dt); -*/ - return; -} diff --git a/src/mex_build_V.cpp b/src/mex_build_V.cpp new file mode 100644 index 0000000..6f11e60 --- /dev/null +++ b/src/mex_build_V.cpp @@ -0,0 +1,290 @@ +/* + * Art der Berechnung durch Parameter bestimmt... + * Analytisch/Semianalytisch/... + * + * + */ +//#include +#include +//#include +#include "mex.h" +/*/#include + +//#include "gauss.hpp" +*/ +//#define M_EPS 10e-8 +/* +//#include "tbb/parallel_for.h" +*/ +#include "slpRectangle.hpp" +/* +//using namespace std; +//using namespace slpR; +*/ +int dimOfVec(double* vec) { + if (vec[2] != 0) + return 2; + else if (vec[1] != 0) + return 1; + else if (vec[0] != 0) + return 0; + else { + mexErrMsgTxt("length of Site is zero"); + /* +// cerr << "dimOfVec : (" << vec[0] << " " << vec[1] << " " << vec[2] +// << ") alle Werte 0 \n";*/ + return -1; + } + +} + +// Gibt von [0 1 2] die Fehlende Zahl zurueck +inline int dimOfThird(int a, int b) { + return ((-(a + b) % 3) + 3) % 3; +} + +void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { +/* +// initeQuad(); +// cout << Quad[1].nod[0]; +*/ + int i, j, k; //Schleifenindizes + double tmp; //Zwischenspeicherung der Einzelnen Werte + int count; + + //Sicherheitsabfragen zu Datengroessen + if ((nrhs != 4)) + mexErrMsgTxt( + "expected (coordinates(Nx3),elements(Mx4),zeta(double),type(int))"); + if (nlhs > 1) + mexErrMsgTxt("has only one output argument"); + + int cm = mxGetM(prhs[0]); + int cn = mxGetN(prhs[0]); + if (cn != 3) + mexErrMsgTxt("expected coordinates (Nx3)"); +// cout << " Koordinaten:" << cm << endl; + + int em = mxGetM(prhs[1]); + int en = mxGetN(prhs[1]); + if (en != 4) + mexErrMsgTxt("expected elements (Mx4)"); +// cout << " Elemente:" << em << endl; + + //Auslesen der Parameter + + plhs[0] = mxCreateDoubleMatrix(em, em, mxREAL); + double * A = mxGetPr(plhs[0]); + double * C = mxGetPr(prhs[0]); + double * E = mxGetPr(prhs[1]); + + int type = (int) *(mxGetPr(prhs[3])); + double * zeta = mxGetPr(prhs[2]); + + //Initialisieren der Hilfsvektoren + + double * x0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * x1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * x2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * x3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * xa = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * xb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + + double * y0 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * y1 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * y2 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * y3 = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * ya = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * yb = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + + double * d = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + double * dt = mxGetPr(mxCreateDoubleMatrix(3, 1, mxREAL)); + + // Welche Funktion soll verwendet werden + double(*ctypeP)(double, double, double, double, double, double, double, + double*); + double(*ctypeO)(double, double, double, double, double, double, double, + double*); + + //Art der Berechnung bestimmen +// cout << " Typ:" << type << endl; + switch (type) { + default: + ctypeP = cParO1; + ctypeO = NULL; + break; + case 0: + ctypeP = cParO0; + ctypeO = NULL; + break; + case 2: + ctypeP = cParO2; + ctypeO = NULL; + break; + case 3: + ctypeP = cParO3; + ctypeO = NULL; + break; + } + + //LageInformationen + int rx, rxa, rxb, ry, rya, ryb; + +// cout << " Progress: #"; + count = 0; + //Ausrechnen + for (j = 0; j < em; ++j) { + x0[0] = C[(int) E[j] - 1]; + x0[1] = C[cm + (int) E[j] - 1]; + x0[2] = C[2 * cm + (int) E[j] - 1]; + + x1[0] = C[(int) E[em + j] - 1]; + x1[1] = C[cm + (int) E[em + j] - 1]; + x1[2] = C[2 * cm + (int) E[em + j] - 1]; + + x2[0] = C[(int) E[2 * em + j] - 1]; + x2[1] = C[cm + (int) E[2 * em + j] - 1]; + x2[2] = C[2 * cm + (int) E[2 * em + j] - 1]; + + x3[0] = C[(int) E[3 * em + j] - 1]; + x3[1] = C[cm + (int) E[3 * em + j] - 1]; + x3[2] = C[2 * cm + (int) E[3 * em + j] - 1]; + + for (i = 0; i < 3; ++i) + xa[i] = x1[i] - x0[i]; + + for (i = 0; i < 3; ++i) + xb[i] = x3[i] - x0[i]; + + // Lageeigenschaften des Flaechenstuecks + + rxa = dimOfVec(xa); + rxb = dimOfVec(xb); + rx = dimOfThird(rxa, rxb); + + //kleinste Ecke finden und fuer \delta verwenden + + if (xa[rxa] > 0) { + if (xb[rxb] > 0) { + for (i = 0; i < 3; ++i) + dt[i] = -x0[i]; + } else { + for (i = 0; i < 3; ++i) + dt[i] = -x3[i]; + } + } else { + if (xb[rxb] > 0) { + for (i = 0; i < 3; ++i) + dt[i] = -x1[i]; + } else { + for (i = 0; i < 3; ++i) + dt[i] = -x2[i]; + } + } + + for (k = 0; k <= j; ++k) { + y0[0] = C[(int) E[k] - 1]; + y0[1] = C[cm + (int) E[k] - 1]; + y0[2] = C[2 * cm + (int) E[k] - 1]; + + y1[0] = C[(int) E[em + k] - 1]; + y1[1] = C[cm + (int) E[em + k] - 1]; + y1[2] = C[2 * cm + (int) E[em + k] - 1]; + + y2[0] = C[(int) E[2 * em + k] - 1]; + y2[1] = C[cm + (int) E[2 * em + k] - 1]; + y2[2] = C[2 * cm + (int) E[2 * em + k] - 1]; + + y3[0] = C[(int) E[3 * em + k] - 1]; + y3[1] = C[cm + (int) E[3 * em + k] - 1]; + y3[2] = C[2 * cm + (int) E[3 * em + k] - 1]; + + for (i = 0; i < 3; ++i) + ya[i] = y1[i] - y0[i]; + + for (i = 0; i < 3; ++i) + yb[i] = y3[i] - y0[i]; + + // Lageeigenschaften des Flaechenstuecks + + rya = dimOfVec(ya); + ryb = dimOfVec(yb); + ry = dimOfThird(rya, ryb); + + //kleinste Ecke finden und fuer \delta verwenden + + if (ya[rya] > 0) { + if (yb[ryb] > 0) { + for (i = 0; i < 3; ++i) + d[i] = dt[i] + y0[i]; + } else { + for (i = 0; i < 3; ++i) + d[i] = dt[i] + y3[i]; + } + } else { + if (yb[ryb] > 0) { + for (i = 0; i < 3; ++i) + d[i] = dt[i] + y1[i]; + } else { + for (i = 0; i < 3; ++i) + d[i] = dt[i] + y2[i]; + } + } + + if (rx == ry) { + if (rxa == rya) { + tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rxa]), + fabs(yb[rxb]), d[rxa], d[rxb], d[rx], zeta); + + } else { + tmp = ctypeP(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[rxa]), + fabs(ya[rxb]), d[rxa], d[rxb], d[rx], zeta); + } + + } else { + + if (rxa == rya) { + tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(ya[rya]), + fabs(yb[ryb]), d[rxb], d[rxa], d[rx], zeta); + } else if (rxa == ryb) { + tmp = ctypeO(fabs(xb[rxb]), fabs(xa[rxa]), fabs(yb[ryb]), + fabs(ya[rya]), d[rxb], d[rxa], d[rx], zeta); + } else if (rxb == rya) { + tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(ya[rya]), + fabs(yb[ryb]), d[rxa], d[rxb], d[rx], zeta); + } else { + tmp = ctypeO(fabs(xa[rxa]), fabs(xb[rxb]), fabs(yb[ryb]), + fabs(ya[rya]), d[rxa], d[rxb], d[rx], zeta); + } + + } + A[(k * em) + j] = 1. / (4 * M_PI) * tmp; + if (k != j) + A[(j * em) + k] = 1. / (4 * M_PI) * tmp; +/* if (count++ > ((em * (em + 1)) / 2) / 10) { + count = 0; +// mexPrintf("#"); +// cout << "#"; +// cout.flush(); + }*/ + + } + + } +// cout << endl; + //Rueckgabe (eventuell zurueck schreiben) +/* mxFree(x0); + mxFree(x1); + mxFree(x3); + mxFree(xa); + mxFree(xb); + mxFree(y0); + mxFree(y1); + mxFree(y3); + mxFree(ya); + mxFree(yb); + mxFree(d); + mxFree(dt); +*/ + return; +} diff --git a/src/test_sol.m b/src/test_sol.m index 76bcb4c..23a1fa7 100644 --- 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