testing chunker_interior in demo_axisysym_neumann

fastalgorithms · Oct 28, 2023 · ddad202 · ddad202
1 parent 39529c2
commit ddad202
Showing 1 changed file with 16 additions and 0 deletions.
diff --git a/chunkie/demo/demo_axissym_neumann.m b/chunkie/demo/demo_axissym_neumann.m
@@ -22,6 +22,22 @@
 
 fprintf('Done building geometry\n');
 
+
+ngrid = 100;
+grid = linspace(-1,1,ngrid);
+[xt,yt] = meshgrid(grid, grid);
+targs = [xt(:).'; yt(:).'];
+
+opts_in = [];
+opts_in.axissym = 1;
+opts_in.fmm = true;
+[in] = chunkerinterior(chnkr, targs, opts_in);
+
+
+plot(chnkr,'kx'); hold on;
+plot(targs(1,in), targs(2,in), 'b.');
+
+
 % source strengths
 strengths = randn(ns, 1);