import platform
platform.platform()
'Darwin-10.8.0-i386-64bit'
import sys
print sys.version
2.7.6 (default, Feb 24 2014, 09:53:59) [GCC 4.2.1 (Apple Inc. build 5666) (dot 3)]
import numpy
numpy.__version__
'1.8.0'
%run test_mc_compute_stationary.py
Elements in each row of P sum to one ... ok Elements of v sum to one ... ok All the elements of v are nonnegative ... FAIL v is a left eigen vector ... ok Elements in each row of P sum to one ... ok Elements of v sum to one ... ok All the elements of v are nonnegative ... FAIL v is a left eigen vector ...
==== Testing with prameter values set 0 N = 27 , epsilon = 0.01 , move = sequential v = [ -1.88568573e-02 -2.55846807e-03 -1.67136105e-04 -6.99899937e-06 -2.11025107e-07 -4.87796717e-09 -8.98786730e-11 -1.35483738e-12 -1.69042488e-14 -6.39194809e-17 -7.40955589e-17 -8.55722397e-17 -9.85973897e-17 -1.13501133e-16 -1.30026937e-16 -1.49343124e-16 -1.71629892e-16 -1.21574765e-16 8.31176299e-15 8.05182291e-13 6.41140988e-11 4.25292338e-09 2.30817772e-07 9.98537756e-06 3.31181689e-04 7.90861874e-03 1.21062702e-01 8.92276953e-01] TOL = 0.01 ==== Testing with prameter values set 1 N = 3 , epsilon = 1e-14 , move = sequential P = [[ 1.00000000e+00 5.00000000e-15 0.00000000e+00 0.00000000e+00] [ 3.33333333e-01 4.99600361e-15 6.66666667e-01 0.00000000e+00] [ 0.00000000e+00 3.33333333e-15 6.66666667e-01 3.33333333e-01] [ 0.00000000e+00 0.00000000e+00 5.00000000e-15 1.00000000e+00]] v = [ -1.14942529e-02 -8.51809045e-17 1.49066583e-14 1.01149425e+00] TOL = 0.01 ====
ok Elements in each row of P sum to one ... ok Elements of v sum to one ... ok All the elements of v are nonnegative ... FAIL v is a left eigen vector ...
An exception has occurred, use %tb to see the full traceback. SystemExit: True
Testing with prameter values set 2 N = 5 , epsilon = 1e-15 , move = simultaneous P = [[ 1.00000000e+00 2.50000000e-15 2.50000000e-30 1.25000000e-45 3.12500000e-61 3.12500000e-77] [ 1.00000000e+00 2.50000000e-15 2.50000000e-30 1.25000000e-45 3.12500000e-61 3.12500000e-77] [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30 2.77555756e-15 1.00000000e+00] [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30 2.77555756e-15 1.00000000e+00] [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30 2.77555756e-15 1.00000000e+00] [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30 2.77555756e-15 1.00000000e+00]] v = [ -1.00000000e+00 -2.60902411e-15 -5.55111512e-17 -5.55111512e-17 5.49560397e-15 2.00000000e+00] TOL = 0.01
ok ====================================================================== FAIL: All the elements of v are nonnegative ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/local/Cellar/python/2.7.6/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest self.test(*self.arg) File "/Volumes/HD1/Users/oyama/Dropbox/documents/Teaching/Zemi/2014/test_mc_compute_stationary/test_mc_compute_stationary.py", line 161, in __call__ eq_(np.prod(v >= 0-TOL), 1) AssertionError: 0 != 1 ====================================================================== FAIL: All the elements of v are nonnegative ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/local/Cellar/python/2.7.6/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest self.test(*self.arg) File "/Volumes/HD1/Users/oyama/Dropbox/documents/Teaching/Zemi/2014/test_mc_compute_stationary/test_mc_compute_stationary.py", line 161, in __call__ eq_(np.prod(v >= 0-TOL), 1) AssertionError: 0 != 1 ====================================================================== FAIL: All the elements of v are nonnegative ---------------------------------------------------------------------- Traceback (most recent call last): File "/usr/local/Cellar/python/2.7.6/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/nose/case.py", line 197, in runTest self.test(*self.arg) File "/Volumes/HD1/Users/oyama/Dropbox/documents/Teaching/Zemi/2014/test_mc_compute_stationary/test_mc_compute_stationary.py", line 161, in __call__ eq_(np.prod(v >= 0-TOL), 1) AssertionError: 0 != 1 ---------------------------------------------------------------------- Ran 12 tests in 0.068s FAILED (failures=3)