diff --git a/CaseStudies/glass/src/Python/NewImplementation/test_All.py b/CaseStudies/glass/src/Python/NewImplementation/test_All.py
new file mode 100644
index 00000000..8dbb7188
--- /dev/null
+++ b/CaseStudies/glass/src/Python/NewImplementation/test_All.py
@@ -0,0 +1,110 @@
+import pytest
+from CurveADT import *
+from Data import *
+from SeqServices import *
+
+class TestAll:
+
+    @classmethod
+    def setup_class(cls):
+        """ setup any state specific to the execution of the given class (which
+        usually contains tests).
+        """
+
+    @classmethod
+    def teardown_class(cls):
+        """ teardown any state that was previously setup with a call to
+        setup_class.
+        """
+    
+    def setup_method(self):
+        X = [1, 2, 10]
+        Y = [1, 2, 10]
+        o = 1
+        self.c1 = CurveT(X, Y, o)
+
+        X = [0, 5, 9, 13]
+        Y = [0, 3, 30, 27]
+        o = 2
+        self.c2 = CurveT(X, Y, o)
+
+    def teardown_method(self, method):
+        self.c1 = None
+        self.c2 = None
+
+    # test SeqServices.py
+
+    @staticmethod
+    def test_isAscending():
+        assert (isAscending([1, 2, 3, 4, 5.5]))
+        assert (not (isAscending([1, 2, 3, 4, 0])))
+        assert (not (isAscending([1, 2, 0, 4, 5])))
+
+    @staticmethod
+    def test_isInBounds():
+        A = [1, 2, 3, 4, 5.5]
+        assert (isInBounds(A, 3.4))
+        assert (not(isInBounds(A, 100)))
+        assert (not(isInBounds(A, 0)))
+    
+    @staticmethod
+    def test_interpLin():
+        assert (interpLin(3 , 6 , 4.5 , 5.67 , 4) == 5.78)
+        
+    @staticmethod
+    def test_interpQuad():
+        assert (interpQuad(0, 0, 3 , 6 , 4.5 , 5.67 , 4) == (5.8533333333333335))
+
+    @staticmethod
+    def test_index():
+        A = [1, 2, 3, 4, 5.5]
+        assert (index(A, 1) == 0) # x is present in list, checks "=" req
+        assert (index(A, 5.5) == 3) # ensures that x < A[i+1]
+        assert (index([1, 2, 3, 4, 5.5], 3.12345678) == 2) # i can be found easily
+
+    # test CurveADT.py
+
+    @staticmethod
+    def test_CurveT_ssmErr():
+        with pytest.raises(SeqSizeMismatch):
+            c = CurveT([0,1,2], [0], 1)
+            c = CurveT([0], [0,1,2], 2)
+
+    @staticmethod
+    def test_CurveT_notAscendingErr():
+        with pytest.raises(IndepVarNotAscending):
+            c = CurveT([0,0,0], [0, 1, 3], 1) # x values are equal
+            c = CurveT([5,4,6], [0, 1, 3], 2)
+
+    @staticmethod
+    def test_CurveT_InvalidInterpErr():
+        with pytest.raises(InvalidInterpOrder):
+            c = CurveT([1, 2, 3], [0, 1, 3], 0)
+            c = CurveT([1, 2, 3], [0, 1, 3], 3)
+
+    def test_CurveT_oodErr(self):
+        with pytest.raises(OutOfDomain):
+            e = self.c1.eval(100)
+            e = self.c1.dfdx(0)
+            e = self.c2.dfdx(100.607)
+    
+    def test_minD(self):
+        assert (self.c1.minD() == 1)
+
+    def test_maxD(self):
+        assert (self.c1.maxD() == 10)
+
+    def test_order(self):
+        assert (self.c1.order() == 1)
+
+    def test_eval(self):
+        assert (self.c1.eval(3) == 3) # order 1 tested
+        assert (self.c2.eval(6) == 7.083333333333334) # order 2 tested
+
+    def test_dfdx(self):
+        assert (self.c1.dfdx(3) == 0.9999999999998899) # order 1 tested
+        assert (self.c2.dfdx(10) == 1.124062500004186) # order 2 tested
+
+    def test_d2fdx2(self):
+        assert (self.c1.d2fdx2(3) == 0) # order 1 tested
+        assert (self.c2.d2fdx2(7) == 1.4999999926601504) # order 2 tested
\ No newline at end of file