Add warnings when parameters are inconsistent with length of Y0

Author: Steven-Roberts

src/+otp/+allencahn/AllenCahnProblem.m

@@ -8,6 +8,16 @@
     end
 
     methods (Access = protected)
+        function validateNewState(obj, newTimeSpan, newY0, newParameters)
+            y0Len = length(newY0);
+            gridPts = newParameters.Size^2;
+            
+            if y0Len ~= gridPts
+                warning('Y0 has size %d, but there are %d grid points', ...
+                    y0Len, gridPts);
+            end
+        end
+        
         function onSettingsChanged(obj)
             n = obj.Parameters.Size;
             alpha = obj.Parameters.Alpha;

src/+otp/+cusp/CUSPProblem.m

@@ -19,6 +19,16 @@
     end
 
     methods (Access = protected)
+        function validateNewState(obj, newTimeSpan, newY0, newParameters)
+            y0Len = length(newY0);
+            gridPts = 3 * newParameters.Size;
+            
+            if y0Len ~= gridPts
+                warning('Y0 has size %d, but there are %d grid points', ...
+                    y0Len, gridPts);
+            end
+        end
+        
         function onSettingsChanged(obj)
             n = obj.Parameters.Size;
             epsilon = obj.Parameters.Epsilon;

src/+otp/+nbody/NBodyProblem.m

@@ -6,6 +6,18 @@
     end
 
     methods (Access = protected)
+        function validateNewState(obj, newTimeSpan, newY0, newParameters)
+            numMasses = length(newParameters.Masses);
+            expectedLen = numMasses * 2^newParameters.SpatialDim;
+            actualLen = length(newY0);
+            
+            if expectedLen ~= actualLen
+                warning( ...
+                    'With %d masses, Y0 should have length %d but has %d', ...
+                    numMasses, expectedLen, actualLen);
+            end
+        end
+        
         function onSettingsChanged(obj)
             spatialDim = obj.Parameters.SpatialDim;
             masses = obj.Parameters.Masses;

src/+otp/+pendulum/PendulumProblem.m

@@ -24,10 +24,25 @@
     end
 
     methods (Access = protected)
+        function validateNewState(obj, newTimeSpan, newY0, newParameters)
+            y0Len = length(newY0);
+            numMasses = length(newParameters.Masses);
+            numLens = length(newParameters.Lengths);
+            
+            if y0Len ~= 2 * numMasses
+                warning( ...
+                    'With %d masses, Y0 should have length %d but has %d', ...
+                    numMasses, 2 * numMasses, y0Len);
+            elseif numMasses ~= numLens
+                warning('Masses has length %d, but Lengths has length %d', ...
+                    numMasses, numLens);
+            end
+        end
+        
         function onSettingsChanged(obj)
             g = obj.Parameters.Gravity;
-            lengths = obj.Parameters.Lengths(:);
-            masses  = obj.Parameters.Masses(:);
+            lengths = obj.Parameters.Lengths;
+            masses  = obj.Parameters.Masses;
 
             numBobs = min(numel(lengths), numel(masses));
             lengths = lengths(1:numBobs);

src/+otp/+qg/QuasiGeostrophicProblem.m

@@ -100,6 +100,15 @@
     end
 
     methods (Access = protected)
+        function validateNewState(obj, newTimeSpan, newY0, newParameters)
+            y0Len = length(newY0);
+            gridPts = newParameters.Nx * newParameters.Ny;
+            
+            if y0Len ~= gridPts
+                warning('Y0 has size %d, but there are %d grid points', ...
+                    y0Len, gridPts);
+            end
+        end
 
         function onSettingsChanged(obj)