add lint test

Author: scivision

+matmap3d/MLint.txt

@@ -0,0 +1 @@
+

+matmap3d/TestALint.m

@@ -0,0 +1,50 @@
+classdef TestALint < matlab.unittest.TestCase
+
+properties
+TestData
+end
+
+properties (TestParameter)
+file = get_files()
+end
+
+methods (TestClassSetup)
+
+function get_files(tc)
+import matlab.unittest.constraints.IsFile
+import matlab.unittest.constraints.IsFolder
+
+cwd = fileparts(mfilename('fullpath'));
+cfg_file = fullfile(cwd, "MLint.txt");
+
+tc.assumeThat(cfg_file, IsFile)
+
+tc.TestData.cfg_file = cfg_file;
+
+end
+
+end
+
+
+methods (Test)
+
+function test_lint(tc, file)
+res = checkcode(file, "-config=" + tc.TestData.cfg_file, "-fullpath");
+
+for j = 1:length(res)
+  tc.verifyFail(append(file, ":", int2str(res(j).line), " ", res(j).message))
+end
+end
+
+end
+
+end
+
+
+function filenames = get_files()
+cwd = fileparts(mfilename('fullpath'));
+flist = dir(fullfile(cwd, '/**/*.m'));
+for i = 1:length(flist)
+ filenames{i} = fullfile(flist(i).folder, flist(i).name); %#ok<AGROW>
+end
+end

+matmap3d/referenceEllipsoid.m

@@ -22,9 +22,11 @@
 function e = referenceEllipsoid(name, lengthUnit)
 arguments
   name (1,1) string = "wgs84"
-  lengthUnit (1,1) string {mustBeMember(lengthUnit, ["m", "meters"])} = "m"
+  lengthUnit (1,1) string = "m"
 end
 
+mustBeMember(lengthUnit, ["m", "meters"])
+
 switch name
 case 'wgs84'
 %% WGS-84 ellipsoid parameters.

+matmap3d/unitTest.m

@@ -272,10 +272,6 @@ function test_ecef2eci(tc)
 utc = datetime(2019, 1, 4, 12,0,0);
 [x,y,z] = matmap3d.ecef2eci(utc, ecef(1), ecef(2), ecef(3));
 tc.verifyEqual([x,y,z], [-2.9818e6, 5.2070e6, 3.1616e6], 'RelTol', 0.01)
-
-[x,y,z] = matmap3d.geodetic2ecef([], 0, 0, 0);
-t = datetime(2000, 1, 1, 12, 0, 0, 'TimeZone', 'UTCLeapSeconds');
-[eci_x, eci_y, eci_z] = matmap3d.ecef2eci(t, x, y, z);
 end
 
 function test_ecef2eci_multiple(tc)

+matmap3d/vdist.m

@@ -125,7 +125,7 @@
   lambdaold(notdone) = lambda(notdone);
   sinsigma(notdone) = sqrt((cos(U2(notdone)).*sin(lambda(notdone)))...
       .^2+(cos(U1(notdone)).*sin(U2(notdone))-sin(U1(notdone)).*...
-      cos(U2(notdone)).*cos(lambda(notdone))).^2);
+      cos(U2(notdone)).*cos(lambda(notdone))).^2); %#ok<*AGROW>
   cossigma(notdone) = sin(U1(notdone)).*sin(U2(notdone))+...
       cos(U1(notdone)).*cos(U2(notdone)).*cos(lambda(notdone));
   % eliminate rare imaginary portions at limit of numerical precision:
@@ -177,7 +177,6 @@
     varargout{2} = reshape(rad2deg(a12),keepsize); % to degrees
 end
 if nargout > 2
-    a21=NaN*lat1;
     % From point #2 to point #1
     % correct sign of lambda for azimuth calcs:
     lambda = abs(lambda);