Matlab use test class

+igrf/TestUnit.m

@@ -0,0 +1,28 @@
+classdef TestUnit < matlab.unittest.TestCase
+
+methods(TestMethodSetup)
+
+function tset(tc)
+cwd = fileparts(mfilename('fullpath'));
+run(fullfile(cwd, '../setup.m'))
+end
+end
+
+methods (Test)
+
+function test_basic(tc)
+
+%% basic
+mag = igrf.igrf(datetime(2010,7,2), 65, 85, 0);
+
+tc.verifyEqual(mag.north, 9295.415460, "RelTol", 0.001, 'north error')
+tc.verifyEqual(mag.east, 2559.7889298, "RelTol", 0.001, 'east error')
+tc.verifyEqual(mag.down, 59670.379598, "RelTol", 0.001, 'down error')
+tc.verifyEqual(mag.total, 60444.284008, "RelTol", 0.001, 'total error')
+
+tc.verifyEqual(mag.incl, 80.821575, "RelTol", 0.001, 'inclination error')
+tc.verifyEqual(mag.decl, 15.396590, "RelTol", 0.001, 'declination error')
+
+end
+end
+end

+igrf/igrf.m

@@ -19,16 +19,3 @@
 B.decl = xarray2mat(dat{'decl'});
 
 end
-
-
-function M = xarray2mat(V)
-% M = double(py.numpy.asfortranarray(V));
-V = V.values;
-S = V.shape;
-V = cell2mat(cell(V.ravel('F').tolist()));
-if length(S) == 1
-  M = V(1);
-else
-  M = reshape(V,[int64(S{1}), int64(S{2})]);
-end
-end

+igrf/private/xarray2mat.m

@@ -0,0 +1,11 @@
+function M = xarray2mat(V)
+% M = double(py.numpy.asfortranarray(V));
+V = V.values;
+S = V.shape;
+V = cell2mat(cell(V.ravel('F').tolist()));
+if length(S) == 1
+  M = V(1);
+else
+  M = reshape(V,[int64(S{1}), int64(S{2})]);
+end
+end

setup.m

@@ -1,3 +1,5 @@
+% this puts Python DLLs on Matlab shell path.
+
 cwd = fileparts(mfilename('fullpath'));
 addpath(cwd)