From: treecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Date: Tue, 8 Nov 2011 10:41:44 +0000 (+0000)
Subject: [src] update fürst test_calcInt1
X-Git-Url: https://git.leopard-lacewing.eu/?a=commitdiff_plain;h=e210a2d0b62fd3f597e198c6dd94c82dff44bf25;p=bacc.git

[src] update fürst test_calcInt1

git-svn-id: https://drops.fb12.tu-berlin.de/svn/bacc/trunk@65 26120e32-c555-405d-b3e1-1f783fb42516
---

diff --git a/src/meshSave/Iso320_1.mat b/src/meshSave/Iso320_1.mat
new file mode 100644
index 0000000..c6984d8
Binary files /dev/null and b/src/meshSave/Iso320_1.mat differ
diff --git a/src/meshSave/Iso320_2.mat b/src/meshSave/Iso320_2.mat
new file mode 100644
index 0000000..03bb063
Binary files /dev/null and b/src/meshSave/Iso320_2.mat differ
diff --git a/src/meshSave/Iso320_3.mat b/src/meshSave/Iso320_3.mat
new file mode 100644
index 0000000..53caa52
Binary files /dev/null and b/src/meshSave/Iso320_3.mat differ
diff --git a/src/meshSave/Iso320_4.mat b/src/meshSave/Iso320_4.mat
new file mode 100644
index 0000000..3ad8967
Binary files /dev/null and b/src/meshSave/Iso320_4.mat differ
diff --git a/src/meshSave/Iso320_5.mat b/src/meshSave/Iso320_5.mat
new file mode 100644
index 0000000..8d5e632
Binary files /dev/null and b/src/meshSave/Iso320_5.mat differ
diff --git a/src/meshSave/anIso320_1.mat b/src/meshSave/anIso320_1.mat
new file mode 100644
index 0000000..40119e5
Binary files /dev/null and b/src/meshSave/anIso320_1.mat differ
diff --git a/src/meshSave/anIso320_2.mat b/src/meshSave/anIso320_2.mat
new file mode 100644
index 0000000..0ede8b8
Binary files /dev/null and b/src/meshSave/anIso320_2.mat differ
diff --git a/src/meshSave/anIso320_3.mat b/src/meshSave/anIso320_3.mat
new file mode 100644
index 0000000..25eb052
Binary files /dev/null and b/src/meshSave/anIso320_3.mat differ
diff --git a/src/meshSave/anIso320_4.mat b/src/meshSave/anIso320_4.mat
new file mode 100644
index 0000000..beaa16c
Binary files /dev/null and b/src/meshSave/anIso320_4.mat differ
diff --git a/src/t.mat b/src/t.mat
new file mode 100644
index 0000000..9808609
Binary files /dev/null and b/src/t.mat differ
diff --git a/src/test_calcInt1.m b/src/test_calcInt1.m
index 89515a5..ad6da47 100644
--- a/src/test_calcInt1.m
+++ b/src/test_calcInt1.m
@@ -1 +1 @@
-

coo = @(h,diff)[0 0 0;1 0 0; 1 1 0;0 1  0; 0 0 0 ; 0+h 0 0; 0+h h 0 ; 0 h 0]+[zeros(4,3);repmat(diff,4,1)];
elements=[1 2 3 4;5 6 7 8];
neigh = zeros(2,8);
dat = [];

diff = [2 0 0];

%% Laage Übersicht
figure(1)
h = 1;
coordinates=coo(h,diff)
current = coordinates(elements(1,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'g');
hold on
current = coordinates(elements(2,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'b');

h = h/2;
coordinates=coo(h,diff);
current = coordinates(elements(2,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'y');
hold off
legend('Element1(h)','Element2(h)','Element2(h/2)');
title('Laage der Elemente');

%% Integrale bei kleiner werdenden Elementen

h = 2;
for I = 1:50

  h = h/2;
  coordinates=coo(h,diff);
  
  A0 = mex_build_AU(coordinates,elements,0,0);
  A2 = mex_build_AU(coordinates,elements,1,2);
  A1 = mex_build_AU(coordinates,elements,1,1);
  I
  dat(I,1:4) = [h A0(1,2) A2(1,2) A1(1,2)];
    dat(dat<0)=0;
end

figure(2)
loglog(dat(:,1),abs(dat(:,2)),dat(:,1),abs(dat(:,3)),dat(:,1),abs(dat(:,4)));

legend('Analytisch','Quad Element','Element vertauschen','location','southeast');
xlabel('Elementgroesse (kürzeste Seite)');
ylabel('Integral');
title('Integral bei kleiner werdenden Element');

%% Netzverfeinerung ab einem Wert

h =  (1/2)^25;
coordinates=coo(h,diff);

A_loadMesh(coordinates,elements,neigh);
datA=[];
for I = 1:5
  datA(I,:) = A_stepIso(1,[0 1 3 2],1,0);
  I
end
figure(3)
loglog(datA(2:end,1),datA(2:end,3),datA(2:end,1),datA(2:end,6),datA(2:end,1),datA(2:end,12),datA(2:end,1),datA(2:end,9))

legend('Analytisch','Element vertauschen','Quad Element','Quad Achse','location','northeast');
xlabel('Anzahl der Elemente');
ylabel('mu Schaetzer');
title('mu Schaetzer mit "schlechtem" Startnetz');

figure(4)
loglog(datA(:,1),datA(:,4),datA(:,1),datA(:,7),datA(:,1),datA(:,13),datA(:,1),datA(:,10))

legend('Analytisch','Element vertauschen','Quad Element','Quad Achse','location','southeast');
xlabel('Anzahl der Elemente');
ylabel('EnergieNorm ^2 ');
title('EnergieNorm ^2 mit "schlechtem" Startnetz');

datA
\ No newline at end of file
+

coo = @(h,diff)[0 0 0;1 0 0; 1 1 0;0 1  0; 0 0 0 ; 0+h 0 0; 0+h h 0 ; 0 h 0]+[zeros(4,3);repmat(diff,4,1)];
elements=[1 2 3 4;5 6 7 8];
neigh = zeros(2,8);
dat = [];

diff = [2 0 0];

%% Laage Übersicht
figure(1)
h = 1;
coordinates=coo(h,diff)
current = coordinates(elements(1,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'g');
hold on
current = coordinates(elements(2,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'b');

h = h/2;
coordinates=coo(h,diff);
current = coordinates(elements(2,[1:4,1])',:);
fill3(current(:,1),current(:,2),current(:,3),'y');
hold off
legend('Element1(h)','Element2(h)','Element2(h/2)');
title('Laage der Elemente');

%% Integrale bei kleiner werdenden Elementen

h = 2;
for I = 1:50

  h = h/2;
  coordinates=coo(h,diff);
  
  A0 = mex_build_AU(coordinates,elements,0,0);
  A2 = mex_build_AU(coordinates,elements,1,2);
  A1 = mex_build_AU(coordinates,elements,1,1);
  I
  dat(I,1:4) = [h A0(1,2) A2(1,2) A1(1,2)];
    dat(dat<0)=0;
end

figure(2)
loglog(dat(:,1),abs(dat(:,2)),dat(:,1),abs(dat(:,3)),dat(:,1),abs(dat(:,4)));

legend('Analytisch','Quad Element','Element vertauschen','location','southeast');
xlabel('Elementgroesse (kürzeste Seite)');
ylabel('Integral');
title('Integral bei kleiner werdenden Element');

%% Netzverfeinerung ab einem Wert

h =  (1/2)^25;
coordinates=coo(h,diff);

A_loadMesh(coordinates,elements,neigh);
datA=[];
for I = 1:5
  datA(I,:) = A_stepIso(5,[0 1 3 2]);
%   datA(I,:) = A_stepAniso(6,[0 1 3 2],1,0);
  I
end
figure(3)
loglog(datA(2:end,1),datA(2:end,3),datA(2:end,1),datA(2:end,6),datA(2:end,1),datA(2:end,12),datA(2:end,1),datA(2:end,9))

legend('Analytisch','Element vertauschen','Quad Element','Quad Achse','location','northeast');
xlabel('Anzahl der Elemente');
ylabel('mu Schaetzer');
title('mu Schaetzer mit "schlechtem" Startnetz');

figure(4)
loglog(datA(:,1),datA(:,4),datA(:,1),datA(:,7),datA(:,1),datA(:,13),datA(:,1),datA(:,10))

legend('Analytisch','Element vertauschen','Quad Element','Quad Achse','location','southeast');
xlabel('Anzahl der Elemente');
ylabel('EnergieNorm ^2 ');
title('EnergieNorm ^2 mit "schlechtem" Startnetz');

datA
\ No newline at end of file