]> git.leopard-lacewing.eu Git - bacc.git/commitdiff
[src] Geschwindigkeitstest zwischen Matlab und Mex
authortreecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Fri, 22 Apr 2011 11:48:06 +0000 (11:48 +0000)
committertreecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Fri, 22 Apr 2011 11:48:06 +0000 (11:48 +0000)
git-svn-id: https://drops.fb12.tu-berlin.de/svn/bacc/trunk@15 26120e32-c555-405d-b3e1-1f783fb42516

src/test_SolveSpeed1.fig [new file with mode: 0644]
src/test_solve.m
src/test_solveSpeed.m [new file with mode: 0644]
src/test_solveSpeed1.mat [new file with mode: 0644]

diff --git a/src/test_SolveSpeed1.fig b/src/test_SolveSpeed1.fig
new file mode 100644 (file)
index 0000000..24401f8
Binary files /dev/null and b/src/test_SolveSpeed1.fig differ
index 422b32360b3c8efd89fca3e73b1b88deef7284bf..1de19e0ee8293db082d9eb7aa47b8c870d450194 100644 (file)
@@ -5,6 +5,7 @@ A1 = zeros(N);
 
 x = zeros(N,2);
 xe = zeros(1,2);
+t = zeros(1,2);
 
 % untere schranke s t obere schranke k l
 intF = @(f,s1,s2,k1,k2,t1,t2,l1,l2) ...
@@ -12,7 +13,7 @@ intF = @(f,s1,s2,k1,k2,t1,t2,l1,l2) ...
     -f(k1,s2,l1,l2)+f(k1,s2,l1,t2)+f(k1,s2,t1,l2)-f(k1,s2,t1,t2)...
     -f(s1,k2,l1,l2)+f(s1,k2,l1,t2)+f(s1,k2,t1,l2)-f(s1,l2,t1,t2)...
     +f(s1,s2,l1,l2)-f(s1,0,l1,t2)-f(s1,s2,t1,l2)+f(s1,s2,t1,t2);
-
+tic
 for j = 1:N
     for k = 1:N
         ej = coordinates(elements(j,:)',:);
@@ -31,10 +32,11 @@ for j = 1:N
 
     end
 end
+t(1) = toc;
 disp ' '
-
+tic
 A2 = build_A(coordinates,elements);
-
+t(2) = toc;
 b = sqrt(sum(quad_norm(coordinates,elements,'w').^2,2));
 
 x(:,1) = A1\b;
@@ -43,5 +45,6 @@ xe(1) = x(:,1)'*A1*x(:,1);
 x(:,2) = A2\b;
 xe(2) = x(:,2)'*A2*x(:,2);
 
-x
-xe
\ No newline at end of file
+x;
+xe
+t
\ No newline at end of file
diff --git a/src/test_solveSpeed.m b/src/test_solveSpeed.m
new file mode 100644 (file)
index 0000000..3ccfa70
--- /dev/null
@@ -0,0 +1,23 @@
+\r
+rt_times = 5;\r
+\r
+time = zeros(rt_times,3);\r
+\r
+    test_solve\r
+    \r
+    time(1,:) = [size(x,1) t];\r
+\r
+for l = 2:rt_times\r
+    test_refine\r
+    test_solve\r
+    \r
+    time(l,:) = [size(x,1) t];\r
+end\r
+\r
+\r
+loglog(time(:,1),time(:,2),time(:,1),time(:,3));\r
+\r
+legend('Matlab','Mex')\r
+xlabel('Elemente')\r
+ylabel('Zeit')\r
+\r
diff --git a/src/test_solveSpeed1.mat b/src/test_solveSpeed1.mat
new file mode 100644 (file)
index 0000000..904568a
Binary files /dev/null and b/src/test_solveSpeed1.mat differ