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) ...
-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,:)',:);
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;
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
--- /dev/null
+\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