Added testing framework with some tests.

Author: tobydriscoll

+sctool/runTests.m

@@ -0,0 +1,11 @@
+function result = runTests
+
+import matlab.unittest.TestSuite;
+
+dirname = which('lapsolve');
+[pathstr,fname,fext] = fileparts(dirname);
+dirname = fullfile(pathstr,'tests');
+suiteClass = TestSuite.fromFolder(dirname);
+result = run(suiteClass);
+
+end

tests/testDisk.m

@@ -0,0 +1,60 @@
+classdef testDisk < matlab.unittest.TestCase
+    
+    properties
+        map
+    end
+    
+    
+    %     methods(TestMethodSetup)
+    %             function createFigure(testCase)
+    %             end
+    %     end
+    %
+    %     methods(TestMethodTeardown)
+    %     end
+    
+    methods (TestClassSetup)
+        function createMap(testCase)
+            opt = sctool.scmapopt('trace',0,'tol',1e-12);
+            p = polygon([4 2i -2+4i -3 -3-1i 2-2i]);
+            testCase.map = diskmap(p,opt);
+        end
+    end
+    
+    methods (Test)
+        
+        function testForwardMap(testCase)
+            result = testCase.map([0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)]);
+            expected = [
+                -2.301479291389453 + 0.891455618349974i,...
+                -2.959017053517382 - 0.004724964608807i,...
+                -2.920229222824237 + 0.110570172682907i,...
+                -2.699042997340806 + 1.203828010636846i,...
+                ];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testInverseMap(testCase)
+            val = testCase.map([0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)]);
+            result = evalinv( testCase.map, val );
+            expected = [0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testCenter(testCase)
+            f = center(testCase.map,0);
+            result = f(0);
+            testCase.verifyEqual(result,0,'abstol',1e-11);
+        end
+        
+        function testPlot(testCase)
+            fig = figure;
+            plot(testCase.map,4,3)
+            close(fig)
+        end
+        
+    end
+    
+    
+    
+end

tests/testExterior.m

@@ -0,0 +1,54 @@
+classdef testExterior < matlab.unittest.TestCase
+    
+    properties
+        map
+    end
+    
+    
+    %     methods(TestMethodSetup)
+    %             function createFigure(testCase)
+    %             end
+    %     end
+    %
+    %     methods(TestMethodTeardown)
+    %     end
+    
+    methods (TestClassSetup)
+        function createMap(testCase)
+            opt = sctool.scmapopt('trace',0,'tol',1e-12);
+            p = polygon([4 2i -2+4i -3 -3-1i 2-2i]);
+            testCase.map = extermap(p,opt);
+        end
+    end
+    
+    methods (Test)
+        
+        function testForwardMap(testCase)
+            result = testCase.map([0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)]).';
+            expected = [
+                3.383320944105799 - 2.338017988574543i
+                -3.257095120413423 + 0.536117032603197i
+                -3.428512520416633 - 0.641446812228358i
+                2.571844815094428 - 1.428155184905573i
+                ];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testInverseMap(testCase)
+            val = testCase.map([0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)]);
+            result = evalinv( testCase.map, val );
+            expected = [0.5+0.5i -0.9 -0.8+0.3i (1+1i)/sqrt(2)];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testPlot(testCase)
+            fig = figure;
+            plot(testCase.map,4,3)
+            close(fig)
+        end
+        
+    end
+    
+    
+    
+end

tests/testHalfplane.m

@@ -0,0 +1,45 @@
+classdef testHalfplane < matlab.unittest.TestCase
+    
+    properties
+        map
+    end
+   
+    methods (TestClassSetup)
+        function createMap(testCase)
+            opt = sctool.scmapopt('trace',0,'tol',1e-12);
+            p = polygon([4 2i -2+4i -3 -3-1i 2-2i]);
+            testCase.map = hplmap(p,opt);
+        end
+    end
+    
+    methods (Test)
+        
+        function testForwardMap(testCase)
+            result = testCase.map([-1+0.01i 2i 4+0.5i Inf]);
+            expected = [
+                3.718839996085665 - 0.046084791699413i,...
+                1.734612962216089 - 0.777136490010106i,...
+                1.179285609480821 - 1.753573737693204i,...
+                2.000000000000000 - 2.000000000000000i,...
+                ];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testInverseMap(testCase)
+            val = testCase.map([-1 -1+0.01i 2i 4+0.5i]);
+            result = evalinv( testCase.map, val );
+            expected = [ -1 -1+0.01i 2i 4+0.5i ];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testPlot(testCase)
+            fig = figure;
+            plot(testCase.map,2,3)
+            close(fig)
+        end
+        
+    end
+    
+    
+    
+end

tests/testRectangle.m

@@ -0,0 +1,45 @@
+classdef testRectangle < matlab.unittest.TestCase
+    
+    properties
+        map
+    end
+    
+    methods (TestClassSetup)
+        function createMap(testCase)
+            opt = sctool.scmapopt('trace',0,'tol',1e-12);
+            p = polygon([4 2i -2+4i -3 -3-1i 2-2i]);
+            testCase.map = rectmap(p,1:4,opt);
+        end
+    end
+    
+    methods (Test)
+        
+        function testForwardMap(testCase)
+            result = testCase.map([1.5 1.4+3i -0.6+1i 1].');
+            expected = [
+                3.641550444027862 - 0.358449555972138i
+                -0.005336970451055 + 1.988135598448822i
+                -1.643459212104280 + 0.428597577267735i
+                1.646072527976422 - 1.929214505595285i
+                ];
+            testCase.verifyEqual(result,expected,'abstol',1e-7);
+        end
+        
+        function testInverseMap(testCase)
+            val = testCase.map([1.5 1.4+3i -0.6+1i 1]);
+            result = evalinv( testCase.map, val );
+            expected = [1.5 1.4+3i -0.6+1i 1];
+            testCase.verifyEqual(result,expected,'abstol',1e-7);
+        end
+                
+        function testPlot(testCase)
+            fig = figure;
+            plot(testCase.map,4,3)
+            close(fig)
+        end
+        
+    end
+    
+    
+    
+end

tests/testStrip.m

@@ -0,0 +1,45 @@
+classdef testStrip < matlab.unittest.TestCase
+    
+    properties
+        map
+    end
+    
+    methods (TestClassSetup)
+        function createMap(testCase)
+            opt = sctool.scmapopt('trace',0,'tol',1e-12);
+            p = polygon([4 2i -2+4i -3 -3-1i 2-2i]);
+            testCase.map = stripmap(p,[1 4],opt);
+        end
+    end
+    
+    methods (Test)
+        
+        function testForwardMap(testCase)
+            result = testCase.map([1+1i Inf -2+0.5i 0].');
+            expected = [
+                -3.000000000000000 - 0.312856946533931i
+                -3.000000000000000 + 0.000000000000000i
+                3.534971699300866 - 0.074755569347609i
+                0.000000000000000 + 2.000000000000000i
+                ];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testInverseMap(testCase)
+            val = testCase.map([4 -2+0.5i 1+0.5i]);
+            result = evalinv( testCase.map, val );
+            expected = [4 -2+0.5i 1+0.5i];
+            testCase.verifyEqual(result,expected,'abstol',1e-10);
+        end
+        
+        function testPlot(testCase)
+            fig = figure;
+            plot(testCase.map,4,3)
+            close(fig)
+        end
+        
+    end
+    
+    
+    
+end