[cig-commits] [commit] master: Added Test_debye (b405b14)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/burnman

On branch  : master
Link       : https://github.com/geodynamics/burnman/compare/6f75d3766c3c8e3bf6599b894e7bbe2993f57419...e7f78ff83a8092bc07fea4398232e73406e7cbbb

>---------------------------------------------------------------

commit b405b1435bea8ec0c9818f00ddcada67c202d728
Author: I <kmanunterborn at gmail.com>
Date:   Fri Dec 12 11:07:35 2014 -0800

    Added Test_debye
    
    Added Test_debye


>---------------------------------------------------------------

b405b1435bea8ec0c9818f00ddcada67c202d728
 tests/test_debye.py | 47 +++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 47 insertions(+)

diff --git a/tests/test_debye.py b/tests/test_debye.py
new file mode 100644
index 0000000..b1e2556
--- /dev/null
+++ b/tests/test_debye.py
@@ -0,0 +1,47 @@
+import unittest
+import os, sys
+
+sys.path.insert(1, os.path.abspath('..'))
+import warnings
+
+import burnman
+from burnman import minerals
+
+from util import BurnManTest
+
+
+class mypericlase(burnman.Mineral):
+    """
+    Stixrude & Lithgow-Bertelloni 2005 and references therein 
+    """
+
+    def __init__(self):
+        self.params = {
+            'equation_of_state': 'slb3',
+            'V_0': 11.24e-6,
+            'K_0': 161.0e9,
+            'Kprime_0': 3.8,
+            'G_0': 131.0e9,
+            'Gprime_0': 2.1,
+            'molar_mass': .0403,
+            'n': 2,
+            'Debye_0': 773.,
+            'grueneisen_0': 1.5,
+            'q_0': 1.5,
+            'eta_s_0': 2.8}
+
+
+class Debye(BurnManTest):
+    def test_temperature(self):
+        #rock = mypericlase()
+        pressure = 0.
+        temperature = 300.
+        eoses = [burnman.slb.SLB2(), burnman.slb.SLB3(), burnman.birch_murnaghan.BM2(), burnman.birch_murnaghan.BM3()]
+
+        test_debye = burnman.debye.debye_fn(temperature)
+        test_debye_cheb = burnman.debye.debye_fn_cheb(temperature)
+        self.assertFloatEqual(test_debye, test_debye_cheb)
+
+
+if __name__ == '__main__':
+    unittest.main()