From: treecity <treecity@26120e32-c555-405d-b3e1-1f783fb42516>
Date: Fri, 22 Apr 2011 11:48:06 +0000 (+0000)
Subject: [src] Geschwindigkeitstest zwischen Matlab und Mex
X-Git-Url: https://git.leopard-lacewing.eu/?a=commitdiff_plain;h=0a25d808f48667ad3a99607407d00c4c3352c898;p=bacc.git

[src] Geschwindigkeitstest zwischen Matlab und Mex

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

diff --git a/src/test_SolveSpeed1.fig b/src/test_SolveSpeed1.fig
new file mode 100644
index 0000000..24401f8
Binary files /dev/null and b/src/test_SolveSpeed1.fig differ
diff --git a/src/test_solve.m b/src/test_solve.m
index 422b323..1de19e0 100644
--- a/src/test_solve.m
+++ b/src/test_solve.m
@@ -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
index 0000000..3ccfa70
--- /dev/null
+++ b/src/test_solveSpeed.m
@@ -0,0 +1,23 @@
+
+rt_times = 5;
+
+time = zeros(rt_times,3);
+
+    test_solve
+    
+    time(1,:) = [size(x,1) t];
+
+for l = 2:rt_times
+    test_refine
+    test_solve
+    
+    time(l,:) = [size(x,1) t];
+end
+
+
+loglog(time(:,1),time(:,2),time(:,1),time(:,3));
+
+legend('Matlab','Mex')
+xlabel('Elemente')
+ylabel('Zeit')
+
diff --git a/src/test_solveSpeed1.mat b/src/test_solveSpeed1.mat
new file mode 100644
index 0000000..904568a
Binary files /dev/null and b/src/test_solveSpeed1.mat differ