revert pquadtest

Author: askhamwhat

devtools/test/pquadTest.m

@@ -7,15 +7,19 @@
 
 zk = norm(kvec);
 
-cparams = []; cparams.maxchunklen = 2.0/real(zk);
+% define geometry and boundary conditions
+% (vertices defined by another function)
+
+cparams = []; cparams.eps = 1e-9;
 chnkr = chunkerfunc(@(t) starfish(t),cparams);
+chnkr = refine(chnkr,struct('nover',1));
 
 % solve and visualize the solution
 
 % build laplace dirichlet matrix
 
 % fkern = kernel('laplace','d');
-fkern = kernel('helmholtz','c',zk,[1,-2i]);
+fkern = kernel('helmholtz','c',zk,[1.5,-2i]);
 opts = [];
 start = tic; C = chunkermat(chnkr,fkern,opts);
 t1 = toc(start);
@@ -27,7 +31,7 @@
 sys = -0.5*eye(chnkr.npt) + C;
 
 % rhs = chnkr.data(1,:); rhs = rhs(:);
-rhs = besselh(0,zk*abs((2+2.25*1i)-(chnkr.r(1,:)+1i*chnkr.r(2,:)))); rhs = rhs(:);
+rhs = besselh(0,zk*abs((1+1.25*1i)-(chnkr.r(1,:)+1i*chnkr.r(2,:)))); rhs = rhs(:);
 start = tic; sol = gmres(sys,rhs,[],1e-14,100); t1 = toc(start);
 
 fprintf('%5.2e s : time for dense gmres\n',t1)
@@ -61,16 +65,14 @@
 t1 = toc(start);
 fprintf('%5.2e s : time for kerneval (Helsing-Ojala for near)\n',t1);
 
-soltrue = besselh(0,zk*abs((2+2.25*1i)-(targets(1,in)+1i*targets(2,in))));
-
 start = tic;
 Csol = chunkerkerneval(chnkr,fkern,sol,targets(:,in)); t1 = toc(start);
 fprintf('%5.2e s : time for kerneval (adaptive for near)\n',t1);
 
 % Compare with reference solution Dsol
-err = max(abs(Csol-Csolpquad))/max(abs(Csol));
-err2 = max(abs(soltrue(:)-Csolpquad(:)))/max(abs(soltrue(:)));
-fprintf('%5.2e : Relative max error\n',err);
+rel_error = max(abs(Csol-Csolpquad))/max(abs(Csol));
+fprintf('%5.2e : Relative max error\n',rel_error);
 % 
 
-assert(err < 1e-10)
+assert(rel_error < 1e-10)
+% profile viewer