# HG changeset patch
# User Jeff Hammel <k0scist@gmail.com>
# Date 1500666415 25200
# Node ID c2f545f320251beda8bb5f6f2ad9f5bca7421343
# Parent  411db53cb9fb9f705dc75b418788e5828ca1ec24
move conformity ensurance to separate function

diff -r 411db53cb9fb -r c2f545f32025 numerics/__init__.py
--- a/numerics/__init__.py	Fri Jul 21 09:37:28 2017 -0700
+++ b/numerics/__init__.py	Fri Jul 21 12:46:55 2017 -0700
@@ -1,4 +1,4 @@
 """
-personal experiments in plotting
+personal experiments in plotting and numerics
 """
 
diff -r 411db53cb9fb -r c2f545f32025 numerics/conformity.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/numerics/conformity.py	Fri Jul 21 12:46:55 2017 -0700
@@ -0,0 +1,19 @@
+"""
+ensures data is what we assert it to be
+"""
+
+
+class NonConformantRowLengths(Exception):
+    """nested arrays have different lengths"""
+
+
+def ensure_row_length(data):
+    """
+    ensures that all rows of array `data` are the same
+    If so, return that length.
+    If not, raise NonConformantArrayLengths
+    """
+    lengths = [len(i) for i in data]
+    if len(set(lengths)) != 1:
+        raise NonConformantRowLengths("Different lengths to array_mean: {}".format(' '.join(lengths)))
+    return lengths.pop()
diff -r 411db53cb9fb -r c2f545f32025 numerics/filters.py
--- a/numerics/filters.py	Fri Jul 21 09:37:28 2017 -0700
+++ b/numerics/filters.py	Fri Jul 21 12:46:55 2017 -0700
@@ -1,20 +1,24 @@
 """
-filter functions for stats
+filter functions for statistics
 """
 
+from .conformity import ensure_row_lengths
+
 __all__ = ['mean', 'array_mean', 'median']
 
 def mean(data):
+    """return arithemetic mean of a vector"""
+
     return sum(data)/float(len(data))
 
+
 def array_mean(data):
     if not data:
         return []
-    lengths = [len(i) for i in data]
-    if len(set(lengths)) != 1:
-        raise AssertionError("Different lengths to array_mean: {}".format(' '.join(lengths)))
+    ensure_row_lengths(data)
     return [mean(i) for i in zip(*data)]
 
+
 def median(data):
     length = len(data)
     index = length/2
diff -r 411db53cb9fb -r c2f545f32025 tests/test_conformity.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tests/test_conformity.py	Fri Jul 21 12:46:55 2017 -0700
@@ -0,0 +1,36 @@
+#!/usr/bin/env python
+
+"""
+test conformity
+"""
+
+import unittest
+from numerics import conformity
+
+class TestConformity(unittest.TestCase):
+    """tests for ensuring data conformity"""
+
+    def test_equal_lengths(self):
+
+        data = [[1,2,3],
+                [4,5,6]
+                [7,8,9]]
+
+        assert conformity.ensure_row_length(data) == 3
+
+    def test_nonequal_lengths(self):
+        data = [[1,2,3],
+                [4,5,6]
+                [7,8,9, 10]  # oops!
+                ]
+
+        e = None
+        try:
+            conformity.ensure_row_length(data)
+        except conformity.NonformantRowLengths as e:
+            pass
+        assert e is not None
+        assert isinstance(e, NonConformantRowLengths)
+
+if __name__ == '__main__':
+    unittest.main()