From 3d73d560ed753cad450d423987de33c8e9c003a6 Mon Sep 17 00:00:00 2001 From: Travis Askham Date: Mon, 28 Apr 2025 14:35:09 -0400 Subject: [PATCH 1/2] removing addpaths_loc, clc, clearvars, clear all from tests --- devtools/test/absconvgaussTest.m | 2 - devtools/test/adapgausswtsTest.m | 3 - devtools/test/addpaths_loc.m | 6 - devtools/test/arclengthfunTest.m | 4 +- .../test/auxiliary_testing_routs/green_sym.m | 52 --- .../test/auxiliary_testing_routs/kern_sym.m | 324 ------------------ devtools/test/axissymkernTest.m | 2 - devtools/test/axissymkernel_sphereTest.m | 2 - devtools/test/chunkerclassunitTest.m | 1 - devtools/test/chunkerfitTest.m | 2 - devtools/test/chunkerfuncTest.m | 2 - devtools/test/chunkerfuncuniTest.m | 2 - devtools/test/chunkerintegralTest.m | 2 - .../test/chunkerkerneval_correctionsTest.m | 2 - devtools/test/chunkerkerneval_gaussidTest.m | 2 - devtools/test/chunkerkerneval_greenhelmTest.m | 2 - devtools/test/chunkerkerneval_greenlapTest.m | 2 - .../test/chunkerkernevalmat_greenlapTest.m | 2 - devtools/test/chunkermatTest.m | 2 - devtools/test/chunkermat_axissymhelm2dTest.m | 2 - .../test/chunkermat_axissymhelm_neumannTest.m | 2 - .../chunkermat_axissymhelm_transmissionTest.m | 2 - ...unkermat_axissymhelm_transmissionTest_re.m | 2 - devtools/test/chunkermat_helm2dTest.m | 2 - devtools/test/chunkermat_l2scaleTest.m | 1 - devtools/test/chunkermat_quadadapTest.m | 2 - .../chunkermat_quadadap_closetotouchingTest.m | 2 - devtools/test/chunkermat_stok2dTest.m | 2 - devtools/test/chunkermat_stok_tractiontest.m | 2 - devtools/test/chunkermat_truepolygonTest.m | 2 - devtools/test/chunkermatapplyTest.m | 2 - devtools/test/chunkerpolyTest.m | 2 - devtools/test/chunkgraph_basicTest.m | 2 - devtools/test/chunkgrphconstructTest.m | 3 - devtools/test/chunkgrphrcipTest.m | 3 - devtools/test/chunkgrphrcipTransmissionTest.m | 1 - devtools/test/chunkgrphrcip_ignoreTest.m | 2 - devtools/test/chunkgrphregionTest.m | 4 - devtools/test/chunkrgrphOpdimTest.m | 2 - devtools/test/elastickernelsTest.m | 3 - devtools/test/flagnearTest.m | 2 - devtools/test/flagrectTest.m | 1 - devtools/test/flagselfTest.m | 1 - devtools/test/flamproxybylevelTest.m | 2 - devtools/test/flamutilitiesTest.m | 2 - devtools/test/helm1d_greenTest.m | 1 - devtools/test/helm2d_greenTest.m | 5 +- devtools/test/kernel_interleaveTest.m | 2 - devtools/test/kernelclassTest.m | 2 - devtools/test/legeexpsunitTest.m | 1 - devtools/test/mixedbcTest.m | 1 - devtools/test/rcipTest.m | 2 - devtools/test/singularkernelTest.m | 2 - devtools/test/stokes_dtracTest.m | 1 - devtools/test/trappermatTest.m | 2 - 55 files changed, 2 insertions(+), 486 deletions(-) delete mode 100644 devtools/test/addpaths_loc.m delete mode 100644 devtools/test/auxiliary_testing_routs/green_sym.m delete mode 100644 devtools/test/auxiliary_testing_routs/kern_sym.m diff --git a/devtools/test/absconvgaussTest.m b/devtools/test/absconvgaussTest.m index 5de9e8dc..4d223e59 100644 --- a/devtools/test/absconvgaussTest.m +++ b/devtools/test/absconvgaussTest.m @@ -2,8 +2,6 @@ % % This file tests the absconvgauss routine -clearvars; close all; -addpaths_loc(); a = 0.75; b = -0.75; diff --git a/devtools/test/adapgausswtsTest.m b/devtools/test/adapgausswtsTest.m index 64ba1b0b..6f9f7cd9 100644 --- a/devtools/test/adapgausswtsTest.m +++ b/devtools/test/adapgausswtsTest.m @@ -3,14 +3,11 @@ % % define geometry and test adaptive integration routine -clearvars; close all; iseed = 8675309; rng(iseed); zk = randn() + 1i*randn(); -addpaths_loc(); - cparams = []; cparams.eps = 1.0e-10; cparams.nover = 2; diff --git a/devtools/test/addpaths_loc.m b/devtools/test/addpaths_loc.m deleted file mode 100644 index e3544a41..00000000 --- a/devtools/test/addpaths_loc.m +++ /dev/null @@ -1,6 +0,0 @@ -function addpaths_loc() - -addpath('../../chunkie'); -addpath('../../chunkie/FLAM') -addpath(genpath_ex('../../chunkie/FLAM')) -addpath('../../chunkie/fmm2d/matlab') \ No newline at end of file diff --git a/devtools/test/arclengthfunTest.m b/devtools/test/arclengthfunTest.m index 1d11d70a..7465925e 100644 --- a/devtools/test/arclengthfunTest.m +++ b/devtools/test/arclengthfunTest.m @@ -2,10 +2,8 @@ % along length of curve % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); % geometry parameters and construction @@ -33,4 +31,4 @@ [s, nchs, ~] = arclengthfun(chnkrtotal); assert(norm(s(:,1:nchs(1)) - ts) < 1e-12); -assert(norm(s(:,nchs(1)+1:end) - rfac*ts) < 1e-12); \ No newline at end of file +assert(norm(s(:,nchs(1)+1:end) - rfac*ts) < 1e-12); diff --git a/devtools/test/auxiliary_testing_routs/green_sym.m b/devtools/test/auxiliary_testing_routs/green_sym.m deleted file mode 100644 index 548ee623..00000000 --- a/devtools/test/auxiliary_testing_routs/green_sym.m +++ /dev/null @@ -1,52 +0,0 @@ -function [val,grad,hess] = green_sym(k,src,targ) -%CHNK.HELM2D.GREEN evaluate the helmholtz green's function -% for the given sources and targets - -[~,ns] = size(src); -[~,nt] = size(targ); - -xs = repmat(src(1,:),nt,1); -ys = repmat(src(2,:),nt,1); - -xt = repmat(targ(1,:).',1,ns); -yt = repmat(targ(2,:).',1,ns); - -rx = xt-xs; -ry = yt-ys; - -rx2 = rx.*rx; -ry2 = ry.*ry; - -r2 = rx2+ry2; - -r = sqrt(r2); - - -if nargout > 0 - h0 = besselh(0,1,k*r); - val = 0.25*1i*h0; -end - -[m,n] = size(xs); - -if nargout > 1 - h1 = besselh(1,1,k*r); - grad = zeros(m,n,2); - - grad(:,:,1) = -1i*k*0.25*h1.*rx./r; - grad(:,:,2) = -1i*k*0.25*h1.*ry./r; - -end - -if nargout > 2 - - hess = zeros(m,n,3); - - r3 = r.^3; - - h2 = 2*h1./(k*r)-h0; - - hess(:,:,1) = 0.25*1i*k*((rx-ry).*(rx+ry).*h1./r3 - k*rx2.*h0./r2); - hess(:,:,2) = 0.25*1i*k*k*rx.*ry.*h2./r2; - hess(:,:,3) = 0.25*1i*k*((ry-rx).*(rx+ry).*h1./r3 - k*ry2.*h0./r2); -end diff --git a/devtools/test/auxiliary_testing_routs/kern_sym.m b/devtools/test/auxiliary_testing_routs/kern_sym.m deleted file mode 100644 index bd19e973..00000000 --- a/devtools/test/auxiliary_testing_routs/kern_sym.m +++ /dev/null @@ -1,324 +0,0 @@ - -function submat= kern_sym(zk,srcinfo,targinfo,type,varargin) -%CHNK.HELM2D.KERN standard Helmholtz layer potential kernels in 2D -% -% Syntax: submat = chnk.heml2d.kern(zk,srcinfo,targingo,type,varargin) -% -% Let x be targets and y be sources for these formulas, with -% n_x and n_y the corresponding unit normals at those points -% (if defined). Note that the normal information is obtained -% by taking the perpendicular to the provided tangential deriviative -% info and normalizing -% -% Kernels based on G(x,y) = i/4 H_0^{(1)}(zk |x-y|) -% -% D(x,y) = \nabla_{n_y} G(x,y) -% S(x,y) = G(x,y) -% S'(x,y) = \nabla_{n_x} G(x,y) -% -% Input: -% zk - complex number, Helmholtz wave number -% srcinfo - description of sources in ptinfo struct format, i.e. -% ptinfo.r - positions (2,:) array -% ptinfo.d - first derivative in underlying -% parameterization (2,:) -% ptinfo.d2 - second derivative in underlying -% parameterization (2,:) -% targinfo - description of targets in ptinfo struct format, -% if info not relevant (d/d2) it doesn't need to -% be provided. sprime requires tangent info in -% targinfo.d -% type - string, determines kernel type -% type == 'd', double layer kernel D -% type == 's', single layer kernel S -% type == 'sprime', normal derivative of single -% layer S' -% type == 'c', combined layer kernel D + i eta S -% varargin{1} - eta in the combined layer formula, otherwise -% does nothing -% -% Output: -% submat - the evaluation of the selected kernel for the -% provided sources and targets. the number of -% rows equals the number of targets and the -% number of columns equals the number of sources -% -% see also CHNK.HELM2D.GREEN - -src = srcinfo.r; -targ = targinfo.r; - -[~,ns] = size(src); -[~,nt] = size(targ); - -if strcmpi(type,'d') - srcnorm = srcinfo.n; - [~,grad] = green_sym(zk,src,targ); - nx = repmat(srcnorm(1,:),nt,1); - ny = repmat(srcnorm(2,:),nt,1); - submat = -(grad(:,:,1).*nx + grad(:,:,2).*ny); - src(2,:) = - src(2,:); - [~,grad] = green_sym(zk,src,targ); - nx = repmat(srcnorm(1,:),nt,1); - ny = repmat(srcnorm(2,:),nt,1); - submat =submat -(grad(:,:,1).*nx - grad(:,:,2).*ny); -end - -if strcmpi(type,'sprime') - targnorm = targinfo.n; - [~,grad] = green_sym(zk,src,targ); - nx = repmat((targnorm(1,:)).',1,ns); - ny = repmat((targnorm(2,:)).',1,ns); - - submat = (grad(:,:,1).*nx + grad(:,:,2).*ny); -end - -if strcmpi(type,'sdtau') - targtan = targinfo.d; - [~,grad] = green_sym(zk,src,targ); - dx = repmat((targtan(1,:)).',1,ns); - dy = repmat((targtan(2,:)).',1,ns); - dn = sqrt(dx.*dx+dy.*dy); - submat = (grad(:,:,1).*dx./dn + grad(:,:,2).*dy)./dn; -end - -if strcmpi(type,'s') - submat = green_sym(zk,src,targ); -end - -if strcmpi(type,'dprime') - targnorm = targinfo.n; - srcnorm = srcinfo.n; - [~,~,hess] = green_sym(zk,src,targ); - nxsrc = repmat(srcnorm(1,:),nt,1); - nysrc = repmat(srcnorm(2,:),nt,1); - nxtarg = repmat((targnorm(1,:)).',1,ns); - nytarg = repmat((targnorm(2,:)).',1,ns); - submat = -(hess(:,:,1).*nxsrc.*nxtarg + hess(:,:,2).*(nysrc.*nxtarg+nxsrc.*nytarg)... - + hess(:,:,3).*nysrc.*nytarg); -end - -if strcmpi(type,'c') - srcnorm = srcinfo.n; - eta = varargin{1}; - [s,grad] = green_sym(zk,src,targ); - nx = repmat(srcnorm(1,:),nt,1); - ny = repmat(srcnorm(2,:),nt,1); - submat = -(grad(:,:,1).*nx + grad(:,:,2).*ny) + 1i*eta*s; -end - -if strcmpi(type,'all') - %targnorm = chnk.normal2d(targinfo); - %srcnorm = chnk.normal2d(srcinfo); - % added by Shidong Jiang to avoid O(N^2) calculation of normals - targnorm = targinfo.n; - srcnorm = srcinfo.n; - - submat = zeros(2*nt,2*ns); - % S - [submats,grad,hess] = green_sym(zk,src,targ); - %[submat(1:2:2*nt,2:2:2*ns),grad,hess] = green_sym(zk,src,targ); - - nxsrc = repmat(srcnorm(1,:),nt,1); - nysrc = repmat(srcnorm(2,:),nt,1); - nxtarg = repmat((targnorm(1,:)).',1,ns); - nytarg = repmat((targnorm(2,:)).',1,ns); - % D' - %submat(2:2:2*nt,1:2:2*ns) = -(hess(:,:,1).*nxsrc.*nxtarg ... - submatdp = -(hess(:,:,1).*nxsrc.*nxtarg ... - + hess(:,:,2).*(nysrc.*nxtarg+nxsrc.*nytarg)... - + hess(:,:,3).*nysrc.*nytarg); - % S' -% submat(2:2:2*nt,2:2:2*ns) = (grad(:,:,1).*nxtarg + grad(:,:,2).*nytarg); - submatsp = (grad(:,:,1).*nxtarg + grad(:,:,2).*nytarg); - % D -% submat(1:2:2*nt,1:2:2*ns) = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - submatd = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - - - submat(1:2:2*nt,1:2:2*ns) = submatd; - submat(1:2:2*nt,2:2:2*ns) = submats; - submat(2:2:2*nt,1:2:2*ns) = submatdp; - submat(2:2:2*nt,2:2:2*ns) = submatsp; -end - -if strcmpi(type,'eval') - coef = varargin{1}; - - srcnorm = srcinfo.n; - - submat = zeros(nt,2*ns); - % S - [submats,grad] = green_sym(zk,src,targ); - - nxsrc = repmat(srcnorm(1,:),nt,1); - nysrc = repmat(srcnorm(2,:),nt,1); - % D - submatd = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - - submat(:,1:2:2*ns) = coef*submatd; - submat(:,2:2:2*ns) = submats; -end - - -if strcmpi(type,'evalg') - coef = varargin{1}; - - srcnorm = srcinfo.n; - - submat = zeros(nt,ns,6); - % S - [submats,grad,hess] = green_sym(zk,src,targ); - - nxsrc = repmat(srcnorm(1,:),nt,1); - nysrc = repmat(srcnorm(2,:),nt,1); - % D - submatd = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - - submat(:,:,1) = coef*submatd; - submat(:,:,2) = submats; - submat(:,:,3) = -coef*(hess(:,:,1).*nxsrc + hess(:,:,2).*nysrc); - submat(:,:,4) = grad(:,:,1); - submat(:,:,5) = -coef*(hess(:,:,2).*nxsrc + hess(:,:,3).*nysrc); - submat(:,:,6) = grad(:,:,2); -end - - - -if strcmpi(type,'trans1') - %targnorm = chnk.normal2d(targinfo); - %srcnorm = chnk.normal2d(srcinfo); - % added by Shidong Jiang to avoid O(N^2) calculation of normals - targnorm = targinfo.n; - srcnorm = srcinfo.n; - - - nxsrc = repmat(srcnorm(1,:),nt,1); - nysrc = repmat(srcnorm(2,:),nt,1); - nxtarg = repmat((targnorm(1,:)).',1,ns); - nytarg = repmat((targnorm(2,:)).',1,ns); - - - submat = zeros(2*nt,2*ns); - - xs = repmat(src(1,:),nt,1); - ys = repmat(src(2,:),nt,1); - - [m,n] = size(xs); - xt = repmat(targ(1,:).',1,ns); - yt = repmat(targ(2,:).',1,ns); - - rx = xt-xs; - ry = yt-ys; - - rx2 = rx.*rx; - ry2 = ry.*ry; - - r2 = rx2+ry2; - - r = sqrt(r2); - k1 = repmat(targinfo.data(1,:).',1,ns); - k2 = repmat(targinfo.data(2,:).',1,ns); - c1 = repmat(targinfo.data(3,:).',1,ns); - c2 = repmat(targinfo.data(4,:).',1,ns); - - harg1 = k1.*r; - hargsz = size(harg1); - harg2 = k2.*r; - hargsz2= size(harg2); - harg = [harg1(:),harg2(:)]; - %szharg = size(harg) - [h00,h11] = chnk.helm2d.besselh01(harg); - %[h00,h11] = hankm103(harg); - %szh00 = size(h00) - h0 = reshape(h00(:,1),hargsz); - h1 = reshape(h11(:,1),hargsz); - %[h0,h1] = hankm103(k1.*r); - %h0 = besselh(0,1,k1.*r); - submats1 = 0.25*1i*h0; - - h1 = h1./r; - %h1 = besselh(1,1,k1.*r)./r; - - grad1 = zeros(m,n,2); - - ck = 0.25*1i*k1; - grad1(:,:,1) = -ck.*h1.*rx; - grad1(:,:,2) = -ck.*h1.*ry; - - hess1 = zeros(m,n,3); - - h2 = 2*h1./k1-h0; - tmp1 = (rx2-ry2).*h1./r2; - tmp2 = k1.*h0./r2; - - - hess1(:,:,1) = ck.*(tmp1 - rx2.*tmp2); - hess1(:,:,2) = ck.*k1.*rx.*ry.*h2./r2; - - hess1(:,:,3) = ck.*(-tmp1 - ry2.*tmp2); - - %[submat(1:2:2*nt,2:2:2*ns),grad,hess] = green_sym(zk,src,targ); - % D' - %submat(2:2:2*nt,1:2:2*ns) = -(hess(:,:,1).*nxsrc.*nxtarg ... - submatdp1 = -(hess1(:,:,1).*nxsrc.*nxtarg ... - + hess1(:,:,2).*(nysrc.*nxtarg+nxsrc.*nytarg)... - + hess1(:,:,3).*nysrc.*nytarg); - % S' -% submat(2:2:2*nt,2:2:2*ns) = (grad(:,:,1).*nxtarg + grad(:,:,2).*nytarg); - submatsp1 = (grad1(:,:,1).*nxtarg + grad1(:,:,2).*nytarg); - % D -% submat(1:2:2*nt,1:2:2*ns) = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - submatd1 = -(grad1(:,:,1).*nxsrc + grad1(:,:,2).*nysrc); - - h0 = reshape(h00(:,2),hargsz2); - h1 = reshape(h11(:,2),hargsz2); - %[h0,h1] = hankm103(k2.*r); - %h0 = besselh(0,1,k2.*r); - submats2 = 0.25*1i*h0; - - %h1 = besselh(1,1,k2.*r)./r; - h1 = h1./r; - - grad2 = zeros(m,n,2); - - ck = 0.25*1i*k2; - grad2(:,:,1) = -ck.*h1.*rx; - grad2(:,:,2) = -ck.*h1.*ry; - - hess2 = zeros(m,n,3); - - h2 = 2*h1./k2-h0; - tmp1 = (rx2-ry2).*h1./r2; - tmp2 = k2.*h0./r2; - - - hess2(:,:,1) = ck.*(tmp1 - rx2.*tmp2); - hess2(:,:,2) = ck.*k2.*rx.*ry.*h2./r2; - - hess2(:,:,3) = ck.*(-tmp1 - ry2.*tmp2); - - %[submat(1:2:2*nt,2:2:2*ns),grad,hess] = green_sym(zk,src,targ); - % D' - %submat(2:2:2*nt,1:2:2*ns) = -(hess(:,:,1).*nxsrc.*nxtarg ... - submatdp2 = -(hess2(:,:,1).*nxsrc.*nxtarg ... - + hess2(:,:,2).*(nysrc.*nxtarg+nxsrc.*nytarg)... - + hess2(:,:,3).*nysrc.*nytarg); - % S' -% submat(2:2:2*nt,2:2:2*ns) = (grad(:,:,1).*nxtarg + grad(:,:,2).*nytarg); - submatsp2 = (grad2(:,:,1).*nxtarg + grad2(:,:,2).*nytarg); - % D -% submat(1:2:2*nt,1:2:2*ns) = -(grad(:,:,1).*nxsrc + grad(:,:,2).*nysrc); - submatd2 = -(grad2(:,:,1).*nxsrc + grad2(:,:,2).*nysrc); - - - alpha1 = 2./(c1+c2); - alpha2 = 2./(1./c1+1./c2); - - - submat(1:2:2*nt,1:2:2*ns) = alpha1.*(c2.*submatd2-c1.*submatd1); - submat(1:2:2*nt,2:2:2*ns) = alpha1.*(submats2-submats1); - submat(2:2:2*nt,1:2:2*ns) = -alpha2.*(submatdp2-submatdp1); - submat(2:2:2*nt,2:2:2*ns) = -alpha2.*(1./c2.*submatsp2-1./c1.*submatsp1); -end - diff --git a/devtools/test/axissymkernTest.m b/devtools/test/axissymkernTest.m index 35e00419..13c09a77 100644 --- a/devtools/test/axissymkernTest.m +++ b/devtools/test/axissymkernTest.m @@ -1,6 +1,4 @@ -clearvars; close all; -addpaths_loc(); cparams = []; cparams.eps = 1.0e-4; pref = []; diff --git a/devtools/test/axissymkernel_sphereTest.m b/devtools/test/axissymkernel_sphereTest.m index 105d43c8..ea480084 100644 --- a/devtools/test/axissymkernel_sphereTest.m +++ b/devtools/test/axissymkernel_sphereTest.m @@ -1,8 +1,6 @@ -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); zk = 40.1; diff --git a/devtools/test/chunkerclassunitTest.m b/devtools/test/chunkerclassunitTest.m index 8c85bb5d..96a35a27 100644 --- a/devtools/test/chunkerclassunitTest.m +++ b/devtools/test/chunkerclassunitTest.m @@ -3,7 +3,6 @@ % setting up some tests for methods of the chunker class % -addpaths_loc(); % adversarial constructor tests isuccess = 0; try diff --git a/devtools/test/chunkerfitTest.m b/devtools/test/chunkerfitTest.m index a0cef29c..e57a9561 100644 --- a/devtools/test/chunkerfitTest.m +++ b/devtools/test/chunkerfitTest.m @@ -1,8 +1,6 @@ %CHUNKERFITTEST -clearvars; close all; -addpaths_loc(); % Sample a smooth curve at random points rng(0) diff --git a/devtools/test/chunkerfuncTest.m b/devtools/test/chunkerfuncTest.m index 4d1cac19..b9cee664 100644 --- a/devtools/test/chunkerfuncTest.m +++ b/devtools/test/chunkerfuncTest.m @@ -5,8 +5,6 @@ % -clearvars; close all; -addpaths_loc(); cparams = []; cparams.eps = 1.0e-4; pref = []; diff --git a/devtools/test/chunkerfuncuniTest.m b/devtools/test/chunkerfuncuniTest.m index 6f1d7c37..0ec18ea8 100644 --- a/devtools/test/chunkerfuncuniTest.m +++ b/devtools/test/chunkerfuncuniTest.m @@ -4,8 +4,6 @@ % This file tests the routine chunkerfunc on a couple of examples % as well as testing the plot, quiver, sort, and reverse utilities -clearvars; close all; -addpaths_loc(); nch = 10; diff --git a/devtools/test/chunkerintegralTest.m b/devtools/test/chunkerintegralTest.m index 43efdfcf..e53067f6 100644 --- a/devtools/test/chunkerintegralTest.m +++ b/devtools/test/chunkerintegralTest.m @@ -2,8 +2,6 @@ % % -clearvars; clear all; -addpaths_loc(); seed = 8675309; rng(seed); diff --git a/devtools/test/chunkerkerneval_correctionsTest.m b/devtools/test/chunkerkerneval_correctionsTest.m index a8eae558..9d7708d9 100644 --- a/devtools/test/chunkerkerneval_correctionsTest.m +++ b/devtools/test/chunkerkerneval_correctionsTest.m @@ -1,8 +1,6 @@ -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); % geometry parameters and construction diff --git a/devtools/test/chunkerkerneval_gaussidTest.m b/devtools/test/chunkerkerneval_gaussidTest.m index 7cb6d3e2..2d3455fd 100644 --- a/devtools/test/chunkerkerneval_gaussidTest.m +++ b/devtools/test/chunkerkerneval_gaussidTest.m @@ -3,10 +3,8 @@ % % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); doadap = true; diff --git a/devtools/test/chunkerkerneval_greenhelmTest.m b/devtools/test/chunkerkerneval_greenhelmTest.m index e4e2cef1..a4042432 100644 --- a/devtools/test/chunkerkerneval_greenhelmTest.m +++ b/devtools/test/chunkerkerneval_greenhelmTest.m @@ -3,10 +3,8 @@ % % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); % geometry parameters and construction diff --git a/devtools/test/chunkerkerneval_greenlapTest.m b/devtools/test/chunkerkerneval_greenlapTest.m index ecf5c926..8a08cddb 100644 --- a/devtools/test/chunkerkerneval_greenlapTest.m +++ b/devtools/test/chunkerkerneval_greenlapTest.m @@ -3,10 +3,8 @@ % % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); doadap = false; diff --git a/devtools/test/chunkerkernevalmat_greenlapTest.m b/devtools/test/chunkerkernevalmat_greenlapTest.m index 4ffa782b..c50210e1 100644 --- a/devtools/test/chunkerkernevalmat_greenlapTest.m +++ b/devtools/test/chunkerkernevalmat_greenlapTest.m @@ -3,10 +3,8 @@ % % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); doadap = false; diff --git a/devtools/test/chunkermatTest.m b/devtools/test/chunkermatTest.m index 7e36033f..d7386d7a 100644 --- a/devtools/test/chunkermatTest.m +++ b/devtools/test/chunkermatTest.m @@ -3,11 +3,9 @@ % % test the matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); cparams = []; cparams.eps = 1.0e-6; diff --git a/devtools/test/chunkermat_axissymhelm2dTest.m b/devtools/test/chunkermat_axissymhelm2dTest.m index 3b4985b2..bd25745d 100644 --- a/devtools/test/chunkermat_axissymhelm2dTest.m +++ b/devtools/test/chunkermat_axissymhelm2dTest.m @@ -3,11 +3,9 @@ % % test the matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); iftorus = 0; diff --git a/devtools/test/chunkermat_axissymhelm_neumannTest.m b/devtools/test/chunkermat_axissymhelm_neumannTest.m index ca630590..1b6482e6 100644 --- a/devtools/test/chunkermat_axissymhelm_neumannTest.m +++ b/devtools/test/chunkermat_axissymhelm_neumannTest.m @@ -1,8 +1,6 @@ -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); zk = 1.1; diff --git a/devtools/test/chunkermat_axissymhelm_transmissionTest.m b/devtools/test/chunkermat_axissymhelm_transmissionTest.m index a699359b..b53c1370 100644 --- a/devtools/test/chunkermat_axissymhelm_transmissionTest.m +++ b/devtools/test/chunkermat_axissymhelm_transmissionTest.m @@ -1,8 +1,6 @@ -clearvars; close all; iseed = 8675309; rng(iseed); -%addpaths_loc(); zk = 5.1; diff --git a/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m b/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m index 586c7ad3..c434917f 100644 --- a/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m +++ b/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m @@ -1,8 +1,6 @@ -clearvars; close all; iseed = 8675309; rng(iseed); -%addpaths_loc(); zk = 10.1; diff --git a/devtools/test/chunkermat_helm2dTest.m b/devtools/test/chunkermat_helm2dTest.m index 1c1aa8ca..1740e796 100644 --- a/devtools/test/chunkermat_helm2dTest.m +++ b/devtools/test/chunkermat_helm2dTest.m @@ -3,11 +3,9 @@ % % test the matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); cparams = []; cparams.eps = 1.0e-10; diff --git a/devtools/test/chunkermat_l2scaleTest.m b/devtools/test/chunkermat_l2scaleTest.m index 35d863c3..11ad1fdd 100644 --- a/devtools/test/chunkermat_l2scaleTest.m +++ b/devtools/test/chunkermat_l2scaleTest.m @@ -1,4 +1,3 @@ -clearvars; close all; zk0 = 2.0; % exterior zk1 = 3.0; % interior diff --git a/devtools/test/chunkermat_quadadapTest.m b/devtools/test/chunkermat_quadadapTest.m index db75f73a..09c8b601 100644 --- a/devtools/test/chunkermat_quadadapTest.m +++ b/devtools/test/chunkermat_quadadapTest.m @@ -3,11 +3,9 @@ % % define geometry and test adaptive matrix builder routine -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); zk = randn() + 1i*randn(); diff --git a/devtools/test/chunkermat_quadadap_closetotouchingTest.m b/devtools/test/chunkermat_quadadap_closetotouchingTest.m index 3f1bafaf..36f4aeed 100644 --- a/devtools/test/chunkermat_quadadap_closetotouchingTest.m +++ b/devtools/test/chunkermat_quadadap_closetotouchingTest.m @@ -6,8 +6,6 @@ % - compare accuracy using the adaptive routine to do nearby panel % evaluation vs smooth rule for everything but self and neighbors -clearvars; close all; -addpaths_loc(); iseed = 8675309; rng(iseed); diff --git a/devtools/test/chunkermat_stok2dTest.m b/devtools/test/chunkermat_stok2dTest.m index 0d1b2ede..4e843f64 100644 --- a/devtools/test/chunkermat_stok2dTest.m +++ b/devtools/test/chunkermat_stok2dTest.m @@ -3,11 +3,9 @@ % % test the matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); cparams = []; cparams.eps = 1.0e-10; diff --git a/devtools/test/chunkermat_stok_tractiontest.m b/devtools/test/chunkermat_stok_tractiontest.m index 7d820764..6651666b 100644 --- a/devtools/test/chunkermat_stok_tractiontest.m +++ b/devtools/test/chunkermat_stok_tractiontest.m @@ -3,11 +3,9 @@ % % test the matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); cparams = []; cparams.eps = 1.0e-10; diff --git a/devtools/test/chunkermat_truepolygonTest.m b/devtools/test/chunkermat_truepolygonTest.m index e12be30e..8bfdbc69 100644 --- a/devtools/test/chunkermat_truepolygonTest.m +++ b/devtools/test/chunkermat_truepolygonTest.m @@ -4,11 +4,9 @@ % test the matrix builder and do a basic solve for a true polygonal % domain. This is the Neumann problem... -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); % pre-defined vertices for a barbell shape diff --git a/devtools/test/chunkermatapplyTest.m b/devtools/test/chunkermatapplyTest.m index e32f0dba..8ed15120 100644 --- a/devtools/test/chunkermatapplyTest.m +++ b/devtools/test/chunkermatapplyTest.m @@ -3,8 +3,6 @@ % % -clearvars; clear all; -addpaths_loc(); seed = 8675309; rng(seed); diff --git a/devtools/test/chunkerpolyTest.m b/devtools/test/chunkerpolyTest.m index 473e1859..957e983c 100644 --- a/devtools/test/chunkerpolyTest.m +++ b/devtools/test/chunkerpolyTest.m @@ -1,7 +1,5 @@ %CHUNKERPOLYTEST -clearvars; close all; -addpaths_loc(); % pre-defined vertices for a barbell shape diff --git a/devtools/test/chunkgraph_basicTest.m b/devtools/test/chunkgraph_basicTest.m index 3fa0f585..204083e5 100644 --- a/devtools/test/chunkgraph_basicTest.m +++ b/devtools/test/chunkgraph_basicTest.m @@ -3,8 +3,6 @@ % test of basic chunkgraph class constructors, methods, etc. % -clearvars -addpaths_loc(); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % diff --git a/devtools/test/chunkgrphconstructTest.m b/devtools/test/chunkgrphconstructTest.m index 5c041118..d15b59c3 100644 --- a/devtools/test/chunkgrphconstructTest.m +++ b/devtools/test/chunkgrphconstructTest.m @@ -5,9 +5,7 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all -addpaths_loc(); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % @@ -16,7 +14,6 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a = -1.0; diff --git a/devtools/test/chunkgrphrcipTest.m b/devtools/test/chunkgrphrcipTest.m index c8a7601d..93b5fb06 100644 --- a/devtools/test/chunkgrphrcipTest.m +++ b/devtools/test/chunkgrphrcipTest.m @@ -5,9 +5,7 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all -addpaths_loc(); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % @@ -16,7 +14,6 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a = -1.0; diff --git a/devtools/test/chunkgrphrcipTransmissionTest.m b/devtools/test/chunkgrphrcipTransmissionTest.m index d05113a8..c17a96ea 100644 --- a/devtools/test/chunkgrphrcipTransmissionTest.m +++ b/devtools/test/chunkgrphrcipTransmissionTest.m @@ -5,7 +5,6 @@ % and tests the Helmholtz transmission problem % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -clearvars a = -1.0; b = 1.0; diff --git a/devtools/test/chunkgrphrcip_ignoreTest.m b/devtools/test/chunkgrphrcip_ignoreTest.m index c2c06eea..40f1d818 100644 --- a/devtools/test/chunkgrphrcip_ignoreTest.m +++ b/devtools/test/chunkgrphrcip_ignoreTest.m @@ -2,8 +2,6 @@ % artificial vertices are ignored -clearvars; close all; -addpaths_loc(); % set wave number zk = 1.1; diff --git a/devtools/test/chunkgrphregionTest.m b/devtools/test/chunkgrphregionTest.m index eba1771b..c2abad48 100644 --- a/devtools/test/chunkgrphregionTest.m +++ b/devtools/test/chunkgrphregionTest.m @@ -5,10 +5,7 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all -clearvars -addpaths_loc(); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % @@ -17,7 +14,6 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -%clear all %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a = -1.0; diff --git a/devtools/test/chunkrgrphOpdimTest.m b/devtools/test/chunkrgrphOpdimTest.m index 95957d63..928c418d 100644 --- a/devtools/test/chunkrgrphOpdimTest.m +++ b/devtools/test/chunkrgrphOpdimTest.m @@ -1,5 +1,3 @@ -clear all -addpaths_loc(); verts = [-3 1; 3 1; 3 5; -3 5; -3 -5; 3 -5; 3 -1; -3 -1].'; diff --git a/devtools/test/elastickernelsTest.m b/devtools/test/elastickernelsTest.m index 145df601..11725c7b 100644 --- a/devtools/test/elastickernelsTest.m +++ b/devtools/test/elastickernelsTest.m @@ -1,11 +1,8 @@ %ELASTIC_TEST % test the kernel definitions for elasticity -clear all iseed = 1234; rng(iseed,'twister'); -clc -%addpaths_loc(); % pde parameters diff --git a/devtools/test/flagnearTest.m b/devtools/test/flagnearTest.m index 71653eb2..7a9878d2 100644 --- a/devtools/test/flagnearTest.m +++ b/devtools/test/flagnearTest.m @@ -5,8 +5,6 @@ % -clearvars; close all; -addpaths_loc(); cparams = []; cparams.eps = 1.0e-6; pref = []; diff --git a/devtools/test/flagrectTest.m b/devtools/test/flagrectTest.m index 9a708253..76484109 100644 --- a/devtools/test/flagrectTest.m +++ b/devtools/test/flagrectTest.m @@ -3,7 +3,6 @@ % tests the flagging routine based on rectangles % -addpaths_loc(); ngrid = 100; chnkr = chunkerfunc(@(t) starfish(t)); chnkr = refine(chnkr); diff --git a/devtools/test/flagselfTest.m b/devtools/test/flagselfTest.m index d3981b8d..a3508d71 100644 --- a/devtools/test/flagselfTest.m +++ b/devtools/test/flagselfTest.m @@ -1,7 +1,6 @@ % FLAGSELFTEST test the routine for flagging when sources and targets % overlap -clear all; close all; nsrc = 40000; xtarg = linspace(0,1,floor(sqrt(nsrc))); diff --git a/devtools/test/flamproxybylevelTest.m b/devtools/test/flamproxybylevelTest.m index 88668b1c..109347d5 100644 --- a/devtools/test/flamproxybylevelTest.m +++ b/devtools/test/flamproxybylevelTest.m @@ -3,7 +3,6 @@ % % test the FLAM matrix builder with level-dependent proxy points -clearvars; close all; iseed = 8675309; rng(iseed); @@ -11,7 +10,6 @@ tol = 1e-14; proxybylevel = true; -% addpaths_loc(); cparams = []; cparams.maxchunklen = 4.0/abs(zk); diff --git a/devtools/test/flamutilitiesTest.m b/devtools/test/flamutilitiesTest.m index 5aa60f14..2cacf051 100644 --- a/devtools/test/flamutilitiesTest.m +++ b/devtools/test/flamutilitiesTest.m @@ -3,11 +3,9 @@ % % test the FLAM matrix builder and do a basic solve -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); nverts = 3; verts = exp(1i*2*pi*(0:(nverts-1))/nverts); diff --git a/devtools/test/helm1d_greenTest.m b/devtools/test/helm1d_greenTest.m index 64b11ee8..e7b9c665 100644 --- a/devtools/test/helm1d_greenTest.m +++ b/devtools/test/helm1d_greenTest.m @@ -1,4 +1,3 @@ -clear all close('all') ifsvdplot = false; diff --git a/devtools/test/helm2d_greenTest.m b/devtools/test/helm2d_greenTest.m index d6577bc2..617709ec 100644 --- a/devtools/test/helm2d_greenTest.m +++ b/devtools/test/helm2d_greenTest.m @@ -2,9 +2,6 @@ % % test that gradients are set up right in chnk.helm2d.green -addpaths_loc(); - -clearvars; close all; iseed = 8675309; rng(iseed); @@ -46,4 +43,4 @@ f = g1(2); g = h1(2:3); -end \ No newline at end of file +end diff --git a/devtools/test/kernel_interleaveTest.m b/devtools/test/kernel_interleaveTest.m index 9646dfe9..2622492b 100644 --- a/devtools/test/kernel_interleaveTest.m +++ b/devtools/test/kernel_interleaveTest.m @@ -1,8 +1,6 @@ -clearvars; close all; iseed = 8675309; rng(iseed); -addpaths_loc(); % dumb test diff --git a/devtools/test/kernelclassTest.m b/devtools/test/kernelclassTest.m index 8b436e78..7e507304 100644 --- a/devtools/test/kernelclassTest.m +++ b/devtools/test/kernelclassTest.m @@ -3,10 +3,8 @@ % % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); doadap = false; diff --git a/devtools/test/legeexpsunitTest.m b/devtools/test/legeexpsunitTest.m index 9b8d5aa6..4521dc43 100644 --- a/devtools/test/legeexpsunitTest.m +++ b/devtools/test/legeexpsunitTest.m @@ -1,6 +1,5 @@ % -addpaths_loc(); k = 19; diff --git a/devtools/test/mixedbcTest.m b/devtools/test/mixedbcTest.m index 84f90475..f72cf3e2 100644 --- a/devtools/test/mixedbcTest.m +++ b/devtools/test/mixedbcTest.m @@ -2,7 +2,6 @@ % mixed Dirichlet and Neumann BCs then test with Dirichlet and transmission % conditions, which has variable opdims -clearvars; close all; nverts = 4; vertsout = exp(1i*2*pi*(0:(nverts-1))/nverts); diff --git a/devtools/test/rcipTest.m b/devtools/test/rcipTest.m index 3de9a286..2b24862b 100644 --- a/devtools/test/rcipTest.m +++ b/devtools/test/rcipTest.m @@ -4,8 +4,6 @@ % This file tests the rcip routines for solving the exterior dirichlet % problem on a domain defined by two arcs of circles meeting at two vertices -clearvars; close all; -addpaths_loc(); ncurve = 2; chnkr(1,ncurve) = chunker(); diff --git a/devtools/test/singularkernelTest.m b/devtools/test/singularkernelTest.m index 7c7499e4..f07b01e7 100644 --- a/devtools/test/singularkernelTest.m +++ b/devtools/test/singularkernelTest.m @@ -2,10 +2,8 @@ % quadratures are working. % -clearvars; close all; seed = 8675309; rng(seed); -addpaths_loc(); doadap = false; diff --git a/devtools/test/stokes_dtracTest.m b/devtools/test/stokes_dtracTest.m index 7da991bf..71664924 100644 --- a/devtools/test/stokes_dtracTest.m +++ b/devtools/test/stokes_dtracTest.m @@ -1,4 +1,3 @@ -clear srcinfo = []; srcinfo.r = rand(2,1); srcinfo.n = rand(2,1); diff --git a/devtools/test/trappermatTest.m b/devtools/test/trappermatTest.m index 42f7463e..a595af94 100644 --- a/devtools/test/trappermatTest.m +++ b/devtools/test/trappermatTest.m @@ -3,13 +3,11 @@ % % test the matrix builder and do a basic solve -clearvars; close all; zk = 1.0; quadorder = 16; iseed = 8675309; rng(iseed); -addpaths_loc(); cparams = []; cparams.npt = 200; From 6c4891610c94dda02fa9543cd9ac82fb2eb4ee7b Mon Sep 17 00:00:00 2001 From: Travis Askham Date: Mon, 28 Apr 2025 15:44:11 -0400 Subject: [PATCH 2/2] encapsulates each test in a function --- devtools/test/absconvgaussTest.m | 9 +++++++++ devtools/test/adapgausswtsTest.m | 9 +++++++++ devtools/test/arclengthfunTest.m | 9 +++++++++ devtools/test/axissymkernTest.m | 9 +++++++++ devtools/test/axissymkernel_sphereTest.m | 9 +++++++++ devtools/test/chunker_nearestTest.m | 10 +++++++++- devtools/test/chunkerclassunitTest.m | 9 +++++++++ devtools/test/chunkerfitTest.m | 9 +++++++++ devtools/test/chunkerfuncTest.m | 9 +++++++++ devtools/test/chunkerfuncuniTest.m | 9 +++++++++ devtools/test/chunkerintegralTest.m | 10 +++++++++- devtools/test/chunkerinteriorTest.m | 9 +++++++++ devtools/test/chunkerkerneval_correctionsTest.m | 9 +++++++++ devtools/test/chunkerkerneval_gaussidTest.m | 9 +++++++++ devtools/test/chunkerkerneval_greenhelmTest.m | 9 +++++++++ devtools/test/chunkerkerneval_greenlapTest.m | 9 +++++++++ devtools/test/chunkerkernevalmat_greenlapTest.m | 9 +++++++++ devtools/test/chunkermatTest.m | 9 +++++++++ devtools/test/chunkermat_axissymhelm2dTest.m | 9 +++++++++ devtools/test/chunkermat_axissymhelm_neumannTest.m | 9 +++++++++ .../test/chunkermat_axissymhelm_transmissionTest.m | 9 +++++++++ .../test/chunkermat_axissymhelm_transmissionTest_re.m | 9 +++++++++ devtools/test/chunkermat_helm2dTest.m | 9 +++++++++ devtools/test/chunkermat_l2scaleTest.m | 9 +++++++++ devtools/test/chunkermat_quadadapTest.m | 9 +++++++++ .../test/chunkermat_quadadap_closetotouchingTest.m | 9 +++++++++ devtools/test/chunkermat_stok2dTest.m | 9 +++++++++ devtools/test/chunkermat_stok_tractiontest.m | 9 +++++++++ devtools/test/chunkermat_truepolygonTest.m | 9 +++++++++ devtools/test/chunkermatapplyTest.m | 9 +++++++++ devtools/test/chunkerpolyTest.m | 9 +++++++++ devtools/test/chunkgraph_basicTest.m | 9 +++++++++ devtools/test/chunkgrphconstructTest.m | 9 +++++++++ devtools/test/chunkgrphrcipTest.m | 9 +++++++++ devtools/test/chunkgrphrcipTransmissionTest.m | 9 +++++++++ devtools/test/chunkgrphrcip_ignoreTest.m | 9 +++++++++ devtools/test/chunkgrphregionTest.m | 9 +++++++++ devtools/test/chunkrgrphOpdimTest.m | 10 +++++++++- devtools/test/complexificationTest.m | 9 +++++++++ devtools/test/datafieldTest.m | 9 +++++++++ devtools/test/elastickernelsTest.m | 9 +++++++++ devtools/test/flagnearTest.m | 9 +++++++++ devtools/test/flagrectTest.m | 9 +++++++++ devtools/test/flagselfTest.m | 10 +++++++++- devtools/test/flamopdimsTest.m | 9 +++++++++ devtools/test/flamproxybylevelTest.m | 9 +++++++++ devtools/test/flamutilitiesTest.m | 9 +++++++++ devtools/test/gradient_check.m | 2 +- devtools/test/helm1d_greenTest.m | 9 +++++++++ devtools/test/helm2d_greenTest.m | 9 +++++++++ devtools/test/kernel_interleaveTest.m | 9 +++++++++ devtools/test/kernelclassTest.m | 9 +++++++++ devtools/test/legeexpsunitTest.m | 9 +++++++++ devtools/test/mixedbcTest.m | 9 +++++++++ devtools/test/pquadTest.m | 9 +++++++++ devtools/test/rcipTest.m | 9 +++++++++ devtools/test/singularkernelTest.m | 9 +++++++++ devtools/test/stokes_dtracTest.m | 9 +++++++++ devtools/test/trappermatTest.m | 9 +++++++++ 59 files changed, 523 insertions(+), 5 deletions(-) diff --git a/devtools/test/absconvgaussTest.m b/devtools/test/absconvgaussTest.m index 4d223e59..b5781f0b 100644 --- a/devtools/test/absconvgaussTest.m +++ b/devtools/test/absconvgaussTest.m @@ -1,3 +1,7 @@ +absconvgaussTest0(); + + +function absconvgaussTest0() %ABSCONVGAUSSTEST % % This file tests the absconvgauss routine @@ -44,6 +48,11 @@ assert(min(errsg) < 1e-6); end + + +end + + function [d,d2] = acgder(x,a,b,h) [~,d,d2] = chnk.spcl.absconvgauss(x,a,b,h); end diff --git a/devtools/test/adapgausswtsTest.m b/devtools/test/adapgausswtsTest.m index 6f9f7cd9..41318a48 100644 --- a/devtools/test/adapgausswtsTest.m +++ b/devtools/test/adapgausswtsTest.m @@ -1,3 +1,7 @@ +adapgausswtsTest0(); + + +function adapgausswtsTest0() %ADAPGAUSSWTSTEST % @@ -113,3 +117,8 @@ matcomp = mat1(istart:iend,jstart:jend); assert(norm(matcomp-mat,'inf') < 1e-11); + + +end + + diff --git a/devtools/test/arclengthfunTest.m b/devtools/test/arclengthfunTest.m index 7465925e..15b50951 100644 --- a/devtools/test/arclengthfunTest.m +++ b/devtools/test/arclengthfunTest.m @@ -1,3 +1,7 @@ +arclengthfunTest0(); + + +function arclengthfunTest0() %ARCLENGTHFUNTEST tests the arclengthfun routine for computing distance % along length of curve % @@ -32,3 +36,8 @@ assert(norm(s(:,1:nchs(1)) - ts) < 1e-12); assert(norm(s(:,nchs(1)+1:end) - rfac*ts) < 1e-12); + + +end + + diff --git a/devtools/test/axissymkernTest.m b/devtools/test/axissymkernTest.m index 13c09a77..2b963a51 100644 --- a/devtools/test/axissymkernTest.m +++ b/devtools/test/axissymkernTest.m @@ -1,3 +1,7 @@ +axissymkernTest0(); + + +function axissymkernTest0() cparams = []; cparams.eps = 1.0e-4; @@ -96,6 +100,11 @@ + + +end + + function [v] = get_exact_kernels(zk, srcinfo, targinfo, type) src = srcinfo.r; diff --git a/devtools/test/axissymkernel_sphereTest.m b/devtools/test/axissymkernel_sphereTest.m index ea480084..c51cb5e2 100644 --- a/devtools/test/axissymkernel_sphereTest.m +++ b/devtools/test/axissymkernel_sphereTest.m @@ -1,3 +1,7 @@ +axissymkernel_sphereTest0(); + + +function axissymkernel_sphereTest0() iseed = 8675309; rng(iseed); @@ -65,6 +69,11 @@ + + +end + + function [chnkobj, sources, targets] = get_geometry(type, pref, ns, nt, maxchunklen) if nargin == 0 diff --git a/devtools/test/chunker_nearestTest.m b/devtools/test/chunker_nearestTest.m index 54ab6c59..580a59b1 100644 --- a/devtools/test/chunker_nearestTest.m +++ b/devtools/test/chunker_nearestTest.m @@ -1,3 +1,7 @@ +chunker_nearestTest0(); + + +function chunker_nearestTest0() %chunker.nearest test % @@ -21,4 +25,8 @@ err = min(abs(th1-th2),abs(th1-th2+2*pi)); err = min(err,abs(th1-th2-2*pi)); assert(err < 1e-12); -end \ No newline at end of file +end + +end + + diff --git a/devtools/test/chunkerclassunitTest.m b/devtools/test/chunkerclassunitTest.m index 96a35a27..0635a20d 100644 --- a/devtools/test/chunkerclassunitTest.m +++ b/devtools/test/chunkerclassunitTest.m @@ -1,3 +1,7 @@ +chunkerclassunitTest0(); + + +function chunkerclassunitTest0() %chunkerclassunitTest % % setting up some tests for methods of the chunker class @@ -48,3 +52,8 @@ assert(chnkr.adj(1,i2) == j) end end + + +end + + diff --git a/devtools/test/chunkerfitTest.m b/devtools/test/chunkerfitTest.m index e57a9561..557f3f71 100644 --- a/devtools/test/chunkerfitTest.m +++ b/devtools/test/chunkerfitTest.m @@ -1,3 +1,7 @@ +chunkerfitTest0(); + + +function chunkerfitTest0() %CHUNKERFITTEST @@ -20,3 +24,8 @@ opts.ifclosed = false; chnkr = chunkerfit(r(:,1:10), opts); assert(checkadjinfo(chnkr) == 0); + + +end + + diff --git a/devtools/test/chunkerfuncTest.m b/devtools/test/chunkerfuncTest.m index b9cee664..d7e9e1fa 100644 --- a/devtools/test/chunkerfuncTest.m +++ b/devtools/test/chunkerfuncTest.m @@ -1,3 +1,7 @@ +chunkerfuncTest0(); + + +function chunkerfuncTest0() %CHUNKERFUNCTEST % @@ -110,3 +114,8 @@ % hold on % quiver(chnkr) % axis equal + + +end + + diff --git a/devtools/test/chunkerfuncuniTest.m b/devtools/test/chunkerfuncuniTest.m index 0ec18ea8..0987153e 100644 --- a/devtools/test/chunkerfuncuniTest.m +++ b/devtools/test/chunkerfuncuniTest.m @@ -1,3 +1,7 @@ +chunkerfuncuniTest0(); + + +function chunkerfuncuniTest0() %CHNKFUNCUNITEST % @@ -62,3 +66,8 @@ assert(abs(a - pi*r^2) < 1e-12,'area wrong for circle domain') + + +end + + diff --git a/devtools/test/chunkerintegralTest.m b/devtools/test/chunkerintegralTest.m index e53067f6..d2c26a3c 100644 --- a/devtools/test/chunkerintegralTest.m +++ b/devtools/test/chunkerintegralTest.m @@ -1,3 +1,7 @@ +chunkerintegralTest0(); + + +function chunkerintegralTest0() %CHUNKERINTEGRALTEST test the routines for integrating over chunks % % @@ -35,4 +39,8 @@ assert(abs(fscal_int1-fscal_int2)/abs(fscal_int2) < 1e-9); assert(abs(fscal_int3-fscal_int2)/abs(fscal_int2) < 1e-9); -assert(abs(fscal_int4-fscal_int2)/abs(fscal_int2) < 1e-9); \ No newline at end of file +assert(abs(fscal_int4-fscal_int2)/abs(fscal_int2) < 1e-9); + +end + + diff --git a/devtools/test/chunkerinteriorTest.m b/devtools/test/chunkerinteriorTest.m index 25b4b9bf..e52c0530 100644 --- a/devtools/test/chunkerinteriorTest.m +++ b/devtools/test/chunkerinteriorTest.m @@ -1,3 +1,7 @@ +chunkerinteriorTest0(); + + +function chunkerinteriorTest0() %CHUNKERINTERIORTEST tests the routines for checking whether a point is % inside a domain or not % @@ -130,3 +134,8 @@ assert(all(in0 == in)); + + +end + + diff --git a/devtools/test/chunkerkerneval_correctionsTest.m b/devtools/test/chunkerkerneval_correctionsTest.m index 9d7708d9..0c7e78e5 100644 --- a/devtools/test/chunkerkerneval_correctionsTest.m +++ b/devtools/test/chunkerkerneval_correctionsTest.m @@ -1,3 +1,7 @@ +chunkerkerneval_correctionsTest0(); + + +function chunkerkerneval_correctionsTest0() seed = 8675309; rng(seed); @@ -67,3 +71,8 @@ assert(norm(utrue-u_eval,inf)>1e-10) + + +end + + diff --git a/devtools/test/chunkerkerneval_gaussidTest.m b/devtools/test/chunkerkerneval_gaussidTest.m index 2d3455fd..6ce377fd 100644 --- a/devtools/test/chunkerkerneval_gaussidTest.m +++ b/devtools/test/chunkerkerneval_gaussidTest.m @@ -1,3 +1,7 @@ +chunkerkerneval_gaussidTest0(); + + +function chunkerkerneval_gaussidTest0() %CHUNKERKERNEVAL_GAUSSIDTEST test the routines for integrating over % chunks against Gauss' identity for the double layer potential % @@ -86,3 +90,8 @@ colorbar end + + +end + + diff --git a/devtools/test/chunkerkerneval_greenhelmTest.m b/devtools/test/chunkerkerneval_greenhelmTest.m index a4042432..10591b48 100644 --- a/devtools/test/chunkerkerneval_greenhelmTest.m +++ b/devtools/test/chunkerkerneval_greenhelmTest.m @@ -1,3 +1,7 @@ +chunkerkerneval_greenhelmTest0(); + + +function chunkerkerneval_greenhelmTest0() %CHUNKERKERNEVAL_GREENHELMTEST test the routines for integrating over % chunks against the Green's ID for Helmholtz % @@ -89,3 +93,8 @@ assert(relerr < 1e-11); + + +end + + diff --git a/devtools/test/chunkerkerneval_greenlapTest.m b/devtools/test/chunkerkerneval_greenlapTest.m index 8a08cddb..d2830a26 100644 --- a/devtools/test/chunkerkerneval_greenlapTest.m +++ b/devtools/test/chunkerkerneval_greenlapTest.m @@ -1,3 +1,7 @@ +chunkerkerneval_greenlapTest0(); + + +function chunkerkerneval_greenlapTest0() %CHUNKERKERNEVAL_GREENLAPTEST test the routines for integrating over % chunks against the Green's ID for Laplace % @@ -126,3 +130,8 @@ assert(relerr < 1e-11); + + +end + + diff --git a/devtools/test/chunkerkernevalmat_greenlapTest.m b/devtools/test/chunkerkernevalmat_greenlapTest.m index c50210e1..9e57d49e 100644 --- a/devtools/test/chunkerkernevalmat_greenlapTest.m +++ b/devtools/test/chunkerkernevalmat_greenlapTest.m @@ -1,3 +1,7 @@ +chunkerkernevalmat_greenlapTest0(); + + +function chunkerkernevalmat_greenlapTest0() %CHUNKERKERNEVAL_GREENLAPTEST test the routines for integrating over % chunks against the Green's ID for Laplace % @@ -102,3 +106,8 @@ assert(relerr < 1e-11); + + +end + + diff --git a/devtools/test/chunkermatTest.m b/devtools/test/chunkermatTest.m index d7386d7a..c2894e94 100644 --- a/devtools/test/chunkermatTest.m +++ b/devtools/test/chunkermatTest.m @@ -1,3 +1,7 @@ +chunkermatTest0(); + + +function chunkermatTest0() %CHUNKERMATTEST % @@ -108,3 +112,8 @@ assert(relerr < 1e-10,'low precision in chunkmat test at target'); + + +end + + diff --git a/devtools/test/chunkermat_axissymhelm2dTest.m b/devtools/test/chunkermat_axissymhelm2dTest.m index bd25745d..5582d2ce 100644 --- a/devtools/test/chunkermat_axissymhelm2dTest.m +++ b/devtools/test/chunkermat_axissymhelm2dTest.m @@ -1,3 +1,7 @@ +chunkermat_axissymhelm2dTest0(); + + +function chunkermat_axissymhelm2dTest0() %CHUNKERMAT_AXISSYMHELM2DTEST % @@ -139,3 +143,8 @@ fprintf('relative l_inf/l_1 error %5.2e\n',relerr2); assert(relerr < 1e-7); + + +end + + diff --git a/devtools/test/chunkermat_axissymhelm_neumannTest.m b/devtools/test/chunkermat_axissymhelm_neumannTest.m index 1b6482e6..4326bcbe 100644 --- a/devtools/test/chunkermat_axissymhelm_neumannTest.m +++ b/devtools/test/chunkermat_axissymhelm_neumannTest.m @@ -1,3 +1,7 @@ +chunkermat_axissymhelm_neumannTest0(); + + +function chunkermat_axissymhelm_neumannTest0() iseed = 8675309; rng(iseed); @@ -212,6 +216,11 @@ fprintf('difference between fast-direct and iterative %5.2e\n',err) + + +end + + function [chnkobj, sources, targets] = get_geometry(type, pref, ns, nt, maxchunklen) if nargin == 0 diff --git a/devtools/test/chunkermat_axissymhelm_transmissionTest.m b/devtools/test/chunkermat_axissymhelm_transmissionTest.m index b53c1370..a4bab772 100644 --- a/devtools/test/chunkermat_axissymhelm_transmissionTest.m +++ b/devtools/test/chunkermat_axissymhelm_transmissionTest.m @@ -1,3 +1,7 @@ +chunkermat_axissymhelm_transmissionTest0(); + + +function chunkermat_axissymhelm_transmissionTest0() iseed = 8675309; rng(iseed); @@ -222,6 +226,11 @@ fprintf('difference between fast-direct and iterative %5.2e\n',err) + + +end + + function [chnkobj, sources, targets] = get_geometry(type, pref, ns, nt, maxchunklen) if nargin == 0 diff --git a/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m b/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m index c434917f..786d526b 100644 --- a/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m +++ b/devtools/test/chunkermat_axissymhelm_transmissionTest_re.m @@ -1,3 +1,7 @@ +chunkermat_axissymhelm_transmissionTest_re0(); + + +function chunkermat_axissymhelm_transmissionTest_re0() iseed = 8675309; rng(iseed); @@ -225,6 +229,11 @@ fprintf('difference between fast-direct and iterative %5.2e\n',err) + + +end + + function [chnkobj, sources, targets] = get_geometry(type, pref, ns, nt, maxchunklen) if nargin == 0 diff --git a/devtools/test/chunkermat_helm2dTest.m b/devtools/test/chunkermat_helm2dTest.m index 1740e796..ae04a98e 100644 --- a/devtools/test/chunkermat_helm2dTest.m +++ b/devtools/test/chunkermat_helm2dTest.m @@ -1,3 +1,7 @@ +chunkermat_helm2dTest0(); + + +function chunkermat_helm2dTest0() %CHUNKERMAT_HELM2DTEST % @@ -114,3 +118,8 @@ fprintf('relative l_inf/l_1 error %5.2e\n',relerr2); assert(relerr < 1e-10); + + +end + + diff --git a/devtools/test/chunkermat_l2scaleTest.m b/devtools/test/chunkermat_l2scaleTest.m index 11ad1fdd..6a59edeb 100644 --- a/devtools/test/chunkermat_l2scaleTest.m +++ b/devtools/test/chunkermat_l2scaleTest.m @@ -1,3 +1,7 @@ +chunkermat_l2scaleTest0(); + + +function chunkermat_l2scaleTest0() zk0 = 2.0; % exterior zk1 = 3.0; % interior @@ -124,3 +128,8 @@ fprintf('error in density=%d\n',err1); + + +end + + diff --git a/devtools/test/chunkermat_quadadapTest.m b/devtools/test/chunkermat_quadadapTest.m index 09c8b601..ca390942 100644 --- a/devtools/test/chunkermat_quadadapTest.m +++ b/devtools/test/chunkermat_quadadapTest.m @@ -1,3 +1,7 @@ +chunkermat_quadadapTest0(); + + +function chunkermat_quadadapTest0() %CHUNKERMAT_QUADADAPTEST % @@ -82,3 +86,8 @@ + + +end + + diff --git a/devtools/test/chunkermat_quadadap_closetotouchingTest.m b/devtools/test/chunkermat_quadadap_closetotouchingTest.m index 36f4aeed..8e5e4170 100644 --- a/devtools/test/chunkermat_quadadap_closetotouchingTest.m +++ b/devtools/test/chunkermat_quadadap_closetotouchingTest.m @@ -1,3 +1,7 @@ +chunkermat_quadadap_closetotouchingTest0(); + + +function chunkermat_quadadap_closetotouchingTest0() %CHUNKERMAT_QUADADAP_CLOSETOTOUCHINGTEST % % - set up two discs which are nearly touching @@ -154,3 +158,8 @@ % figure(2) % plot3(chnkr,1) % + + +end + + diff --git a/devtools/test/chunkermat_stok2dTest.m b/devtools/test/chunkermat_stok2dTest.m index 4e843f64..9f0c8de2 100644 --- a/devtools/test/chunkermat_stok2dTest.m +++ b/devtools/test/chunkermat_stok2dTest.m @@ -1,3 +1,7 @@ +chunkermat_stok2dTest0(); + + +function chunkermat_stok2dTest0() %CHUNKERMAT_HELM2DTEST % @@ -162,3 +166,8 @@ assert(relerr < 1e-10) + + +end + + diff --git a/devtools/test/chunkermat_stok_tractiontest.m b/devtools/test/chunkermat_stok_tractiontest.m index 6651666b..d8ca7623 100644 --- a/devtools/test/chunkermat_stok_tractiontest.m +++ b/devtools/test/chunkermat_stok_tractiontest.m @@ -1,3 +1,7 @@ +chunkermat_stok_tractiontest0(); + + +function chunkermat_stok_tractiontest0() %CHUNKERMAT_HELM2DTEST % @@ -108,3 +112,8 @@ fprintf('relative frobenius error %5.2e\n',relerr); assert(relerr < 1e-10); + + +end + + diff --git a/devtools/test/chunkermat_truepolygonTest.m b/devtools/test/chunkermat_truepolygonTest.m index 8bfdbc69..7d835891 100644 --- a/devtools/test/chunkermat_truepolygonTest.m +++ b/devtools/test/chunkermat_truepolygonTest.m @@ -1,3 +1,7 @@ +chunkermat_truepolygonTest0(); + + +function chunkermat_truepolygonTest0() %CHUNKERMAT_TRUEPOLYGONTEST % @@ -136,3 +140,8 @@ assert(relerr < 1e-10); + + +end + + diff --git a/devtools/test/chunkermatapplyTest.m b/devtools/test/chunkermatapplyTest.m index 8ed15120..053ba3e1 100644 --- a/devtools/test/chunkermatapplyTest.m +++ b/devtools/test/chunkermatapplyTest.m @@ -1,3 +1,7 @@ +chunkermatapplyTest0(); + + +function chunkermatapplyTest0() %CHUNKERMATAPPLYTEST test the routines for a matrix free apply of the system %matrix % @@ -223,6 +227,11 @@ % assert(relerr < 1e-10) + + +end + + function [r,d,d2] = sinearc(t,amp,frq) xs = t; ys = amp*sin(frq*t); diff --git a/devtools/test/chunkerpolyTest.m b/devtools/test/chunkerpolyTest.m index 957e983c..9da20004 100644 --- a/devtools/test/chunkerpolyTest.m +++ b/devtools/test/chunkerpolyTest.m @@ -1,3 +1,7 @@ +chunkerpolyTest0(); + + +function chunkerpolyTest0() %CHUNKERPOLYTEST @@ -105,3 +109,8 @@ hold on quiver(chnkr_ref3) axis equal + + +end + + diff --git a/devtools/test/chunkgraph_basicTest.m b/devtools/test/chunkgraph_basicTest.m index 204083e5..37f16d2e 100644 --- a/devtools/test/chunkgraph_basicTest.m +++ b/devtools/test/chunkgraph_basicTest.m @@ -1,3 +1,7 @@ +chunkgraph_basicTest0(); + + +function chunkgraph_basicTest0() %chunkgraph_basicTest % % test of basic chunkgraph class constructors, methods, etc. @@ -143,6 +147,11 @@ % simpler calling sequence if all edges same function cgrph = chunkgraph(verts,endverts,fchnks); + + +end + + function idstrue = polygonids(cg,xx,yy) verts = cg.verts; diff --git a/devtools/test/chunkgrphconstructTest.m b/devtools/test/chunkgrphconstructTest.m index d15b59c3..5b216249 100644 --- a/devtools/test/chunkgrphconstructTest.m +++ b/devtools/test/chunkgrphconstructTest.m @@ -1,3 +1,7 @@ +chunkgrphconstructTest0(); + + +function chunkgrphconstructTest0() %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % . . . builds a simple pentagonal chunkergraph @@ -71,6 +75,11 @@ assert(all(cgrph1.verts(:) == cgrph2.verts(:))) assert(all(all(cgrph1.v2emat == cgrph2.v2emat))) + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/chunkgrphrcipTest.m b/devtools/test/chunkgrphrcipTest.m index 93b5fb06..a4dc280d 100644 --- a/devtools/test/chunkgrphrcipTest.m +++ b/devtools/test/chunkgrphrcipTest.m @@ -1,3 +1,7 @@ +chunkgrphrcipTest0(); + + +function chunkgrphrcipTest0() %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % . . . builds a simple pentagonal chunkergraph @@ -177,6 +181,11 @@ + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/chunkgrphrcipTransmissionTest.m b/devtools/test/chunkgrphrcipTransmissionTest.m index c17a96ea..04da2442 100644 --- a/devtools/test/chunkgrphrcipTransmissionTest.m +++ b/devtools/test/chunkgrphrcipTransmissionTest.m @@ -1,3 +1,7 @@ +chunkgrphrcipTransmissionTest0(); + + +function chunkgrphrcipTransmissionTest0() %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % @@ -135,6 +139,11 @@ + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/chunkgrphrcip_ignoreTest.m b/devtools/test/chunkgrphrcip_ignoreTest.m index 40f1d818..6e335ab3 100644 --- a/devtools/test/chunkgrphrcip_ignoreTest.m +++ b/devtools/test/chunkgrphrcip_ignoreTest.m @@ -1,3 +1,7 @@ +chunkgrphrcip_ignoreTest0(); + + +function chunkgrphrcip_ignoreTest0() % chunkgrphrcip_ignoreTest verify that solution is accurate when % artificial vertices are ignored @@ -160,6 +164,11 @@ + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/chunkgrphregionTest.m b/devtools/test/chunkgrphregionTest.m index c2abad48..4219f616 100644 --- a/devtools/test/chunkgrphregionTest.m +++ b/devtools/test/chunkgrphregionTest.m @@ -1,3 +1,7 @@ +chunkgrphregionTest0(); + + +function chunkgrphregionTest0() %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % . . . builds a simple pentagonal chunkergraph @@ -78,6 +82,11 @@ end end + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/chunkrgrphOpdimTest.m b/devtools/test/chunkrgrphOpdimTest.m index 928c418d..67c48103 100644 --- a/devtools/test/chunkrgrphOpdimTest.m +++ b/devtools/test/chunkrgrphOpdimTest.m @@ -1,3 +1,7 @@ +chunkrgrphOpdimTest0(); + + +function chunkrgrphOpdimTest0() verts = [-3 1; 3 1; 3 5; -3 5; -3 -5; 3 -5; 3 -1; -3 -1].'; @@ -48,4 +52,8 @@ end opts = struct('nonsmoothonly',false, 'rcip',true); -[sysmat] = chunkermat(cgrph, fkerns, opts); \ No newline at end of file +[sysmat] = chunkermat(cgrph, fkerns, opts); + +end + + diff --git a/devtools/test/complexificationTest.m b/devtools/test/complexificationTest.m index 0daf2849..276fb0d3 100644 --- a/devtools/test/complexificationTest.m +++ b/devtools/test/complexificationTest.m @@ -1,3 +1,7 @@ +complexificationTest0(); + + +function complexificationTest0() % test solve_sommerfeld_dens and eval_sommerfeld_dens % and compare to complexification solution @@ -74,6 +78,11 @@ + + +end + + function [f,fd,fdd] = flat_interface(t, a, b, t0, t1) phi = @(t,u,v,z) u*(t-z).*erfc(u*(t-z))*v - exp(-u^2*(t-z).^2)/sqrt(pi)*v; diff --git a/devtools/test/datafieldTest.m b/devtools/test/datafieldTest.m index f5b5e323..d705fb09 100644 --- a/devtools/test/datafieldTest.m +++ b/devtools/test/datafieldTest.m @@ -1,3 +1,7 @@ +datafieldTest0(); + + +function datafieldTest0() % A simple test to show data fields being used with a kernel. %% @@ -71,6 +75,11 @@ %% + + +end + + function kern = directional_der_S_kern() kern = kernel(); diff --git a/devtools/test/elastickernelsTest.m b/devtools/test/elastickernelsTest.m index 11725c7b..4a2dc688 100644 --- a/devtools/test/elastickernelsTest.m +++ b/devtools/test/elastickernelsTest.m @@ -1,3 +1,7 @@ +elastickernelsTest0(); + + +function elastickernelsTest0() %ELASTIC_TEST % test the kernel definitions for elasticity @@ -114,6 +118,11 @@ assert(err_neu < 1e-10); + + +end + + function [pde_errs,pdedalt_errs,div_errs,trac_errs,gid_err,daltgrad_errs] = ... test_kernels(chnkr,s,targ,lam,mu,f,niter) %TEST_KERNELS diff --git a/devtools/test/flagnearTest.m b/devtools/test/flagnearTest.m index 7a9878d2..d4fa187c 100644 --- a/devtools/test/flagnearTest.m +++ b/devtools/test/flagnearTest.m @@ -1,3 +1,7 @@ +flagnearTest0(); + + +function flagnearTest0() %FLAGNEARTEST % @@ -52,3 +56,8 @@ assert(nnz(flag2 ~= flag) == 0) + + +end + + diff --git a/devtools/test/flagrectTest.m b/devtools/test/flagrectTest.m index 76484109..75f0faae 100644 --- a/devtools/test/flagrectTest.m +++ b/devtools/test/flagrectTest.m @@ -1,3 +1,7 @@ +flagrectTest0(); + + +function flagrectTest0() %flagrectTest % % tests the flagging routine based on rectangles @@ -35,6 +39,11 @@ + + +end + + function rects = bounding_rects(convreg) % find minimal area bounding rectangle for each region diff --git a/devtools/test/flagselfTest.m b/devtools/test/flagselfTest.m index a3508d71..7a643c92 100644 --- a/devtools/test/flagselfTest.m +++ b/devtools/test/flagselfTest.m @@ -1,3 +1,7 @@ +flagselfTest0(); + + +function flagselfTest0() % FLAGSELFTEST test the routine for flagging when sources and targets % overlap @@ -25,4 +29,8 @@ errs = sum(vecnorm(srcs(:,flagslf(1,:)) - targs(:,flagslf(2,:)))>1e-10); fprintf('number of incorrectly flagged points: %d\n', errs) -assert(norm(srcs(:,flagslf(1,:)) - srcs(:,P(flagslf(2,:))),1)<1e-6) \ No newline at end of file +assert(norm(srcs(:,flagslf(1,:)) - srcs(:,P(flagslf(2,:))),1)<1e-6) + +end + + diff --git a/devtools/test/flamopdimsTest.m b/devtools/test/flamopdimsTest.m index fcb30e50..47d7063f 100644 --- a/devtools/test/flamopdimsTest.m +++ b/devtools/test/flamopdimsTest.m @@ -1,3 +1,7 @@ +flamopdimsTest0(); + + +function flamopdimsTest0() zk1 = 10; % exterior wave number zk2 = 15; % coating wave number @@ -46,6 +50,11 @@ + + +end + + function err1 = test_analytic_soln(zks, chnkr1, chnkr2, src1 ,src2, A, ifflam) zk1 = zks(1); zk2 = zks(2); diff --git a/devtools/test/flamproxybylevelTest.m b/devtools/test/flamproxybylevelTest.m index 109347d5..87eace99 100644 --- a/devtools/test/flamproxybylevelTest.m +++ b/devtools/test/flamproxybylevelTest.m @@ -1,3 +1,7 @@ +flamproxybylevelTest0(); + + +function flamproxybylevelTest0() %FLAMPROXYBYLEVELTEST % @@ -183,3 +187,8 @@ relerr2 = norm(utarg-sol_ifmm,'inf')/dot(abs(sig_rskelf(:)),wts(:)); fprintf('relative frobenius error %5.2e\n',relerr); fprintf('relative l_inf/l_1 error %5.2e\n',relerr2); + + +end + + diff --git a/devtools/test/flamutilitiesTest.m b/devtools/test/flamutilitiesTest.m index 2cacf051..89e9145e 100644 --- a/devtools/test/flamutilitiesTest.m +++ b/devtools/test/flamutilitiesTest.m @@ -1,3 +1,7 @@ +flamutilitiesTest0(); + + +function flamutilitiesTest0() %FLAMUTILITIESTEST % @@ -211,6 +215,11 @@ + + +end + + function [r,d,d2] = sinearc(t,amp,frq) xs = t; ys = amp*sin(frq*t); diff --git a/devtools/test/gradient_check.m b/devtools/test/gradient_check.m index feb25fa0..84b4af42 100644 --- a/devtools/test/gradient_check.m +++ b/devtools/test/gradient_check.m @@ -44,4 +44,4 @@ function cfprintf(ifprint,varargin) if (ifprint); fprintf(varargin{:}); end -end \ No newline at end of file +end diff --git a/devtools/test/helm1d_greenTest.m b/devtools/test/helm1d_greenTest.m index e7b9c665..13ed5f39 100644 --- a/devtools/test/helm1d_greenTest.m +++ b/devtools/test/helm1d_greenTest.m @@ -1,3 +1,7 @@ +helm1d_greenTest0(); + + +function helm1d_greenTest0() close('all') ifsvdplot = false; @@ -138,6 +142,11 @@ + + +end + + function [sol_gmres,vdens] = fast_solve_wrap(uu,m,E,chnkr,istart,iend,ab) figs(1) = 0; diff --git a/devtools/test/helm2d_greenTest.m b/devtools/test/helm2d_greenTest.m index 617709ec..a8c9f99f 100644 --- a/devtools/test/helm2d_greenTest.m +++ b/devtools/test/helm2d_greenTest.m @@ -1,3 +1,7 @@ +helm2d_greenTest0(); + + +function helm2d_greenTest0() %HELM2D_GREENTEST % % test that gradients are set up right in chnk.helm2d.green @@ -23,6 +27,11 @@ assert(min(errsfx) < 1e-8); assert(min(errsfy) < 1e-8); + + +end + + function [f,g] = fcn1(zk,src,trg) [f,g] = chnk.helm2d.green(zk,src,trg); diff --git a/devtools/test/kernel_interleaveTest.m b/devtools/test/kernel_interleaveTest.m index 2622492b..6dc20c22 100644 --- a/devtools/test/kernel_interleaveTest.m +++ b/devtools/test/kernel_interleaveTest.m @@ -1,3 +1,7 @@ +kernel_interleaveTest0(); + + +function kernel_interleaveTest0() iseed = 8675309; rng(iseed); @@ -197,3 +201,8 @@ fprintf('difference between fast-direct and iterative %5.2e\n',err) + + +end + + diff --git a/devtools/test/kernelclassTest.m b/devtools/test/kernelclassTest.m index 7e507304..7b0a26e7 100644 --- a/devtools/test/kernelclassTest.m +++ b/devtools/test/kernelclassTest.m @@ -1,3 +1,7 @@ +kernelclassTest0(); + + +function kernelclassTest0() %KERNELCLASSTEST check the use of kernel class for sending FMM % and singularity info to various routines % @@ -96,3 +100,8 @@ kerntmp = nan*kerns; assert(kerntmp.isnan); + + +end + + diff --git a/devtools/test/legeexpsunitTest.m b/devtools/test/legeexpsunitTest.m index 4521dc43..681f2545 100644 --- a/devtools/test/legeexpsunitTest.m +++ b/devtools/test/legeexpsunitTest.m @@ -1,3 +1,7 @@ +legeexpsunitTest0(); + + +function legeexpsunitTest0() % @@ -29,3 +33,8 @@ assert(abs(itrue-icoefs)< 1e-14); end + + +end + + diff --git a/devtools/test/mixedbcTest.m b/devtools/test/mixedbcTest.m index f72cf3e2..036bcc75 100644 --- a/devtools/test/mixedbcTest.m +++ b/devtools/test/mixedbcTest.m @@ -1,3 +1,7 @@ +mixedbcTest0(); + + +function mixedbcTest0() % MIXEDBCTEST test the code with mixed boundary conditions, first test with % mixed Dirichlet and Neumann BCs then test with Dirichlet and transmission % conditions, which has variable opdims @@ -286,3 +290,8 @@ relerr = max(relerr1, relerr2); fprintf('relative field error %5.2e\n',relerr); assert(relerr < 1e-4) + + +end + + diff --git a/devtools/test/pquadTest.m b/devtools/test/pquadTest.m index 7328a1cd..b18311b3 100644 --- a/devtools/test/pquadTest.m +++ b/devtools/test/pquadTest.m @@ -1,3 +1,7 @@ +pquadTest0(); + + +function pquadTest0() % testing product quadrature rule % planewave vec @@ -76,3 +80,8 @@ assert(rel_error < 1e-10) % profile viewer + + +end + + diff --git a/devtools/test/rcipTest.m b/devtools/test/rcipTest.m index 2b24862b..4592c21d 100644 --- a/devtools/test/rcipTest.m +++ b/devtools/test/rcipTest.m @@ -1,3 +1,7 @@ +rcipTest0(); + + +function rcipTest0() %RCIPTEST % @@ -249,6 +253,11 @@ % Auxiliary routines for generating boundary % + + +end + + function [r,d,d2] = circulararc(t,cpars) %%circulararc % return position, first and second derivatives of a circular arc diff --git a/devtools/test/singularkernelTest.m b/devtools/test/singularkernelTest.m index f07b01e7..aa272da3 100644 --- a/devtools/test/singularkernelTest.m +++ b/devtools/test/singularkernelTest.m @@ -1,3 +1,7 @@ +singularkernelTest0(); + + +function singularkernelTest0() %SINGULARKERNELTEST check that principal value and hypersingular type % quadratures are working. % @@ -72,3 +76,8 @@ relerr = norm(un-unbdry(:))/norm(unbdry); fprintf('%5.2e : error in un test (testing hs integration)\n',relerr); + + +end + + diff --git a/devtools/test/stokes_dtracTest.m b/devtools/test/stokes_dtracTest.m index 71664924..bd5e1107 100644 --- a/devtools/test/stokes_dtracTest.m +++ b/devtools/test/stokes_dtracTest.m @@ -1,3 +1,7 @@ +stokes_dtracTest0(); + + +function stokes_dtracTest0() srcinfo = []; srcinfo.r = rand(2,1); srcinfo.n = rand(2,1); @@ -33,3 +37,8 @@ f(2:2:end) = -p.*nty + (euxy.*ntx + euyy.*nty)*mu; norm(f - Kt) + + +end + + diff --git a/devtools/test/trappermatTest.m b/devtools/test/trappermatTest.m index a595af94..af9a84d4 100644 --- a/devtools/test/trappermatTest.m +++ b/devtools/test/trappermatTest.m @@ -1,3 +1,7 @@ +trappermatTest0(); + + +function trappermatTest0() %TRAPPERMATTEST % @@ -114,3 +118,8 @@ fprintf('relative frobenius error %5.2e\n',relerr); fprintf('relative l_inf/l_1 error %5.2e\n',relerr2); + + +end + +