In [1]:

    

import platform
platform.platform()


    

    
Out[1]:





'Darwin-13.3.0-x86_64-i386-64bit'

In [2]:

    

import sys
print sys.version


    

    




2.7.6 |Anaconda 1.9.2 (x86_64)| (default, Jan 10 2014, 11:23:15) 
[GCC 4.0.1 (Apple Inc. build 5493)]

In [3]:

    

import numpy
numpy.__version__


    

    
Out[3]:





'1.8.0'

In [4]:

    

# Default values: N=27, epsilon=1e-2, move='sequential', tolerance=1e-2
%run unittest_mc_compute_stationary.py


    

    




test_left_eigen_vec (__main__.TestComputeStationary) ... ok
test_markov_matrix (__main__.TestComputeStationary) ... ok
test_nonnegative (__main__.TestComputeStationary) ... FAIL
test_sum_one (__main__.TestComputeStationary) ... 





    




N = 27 , epsilon = 0.01 

v =
[ -2.48255180e-02  -3.36828636e-03  -2.20038808e-04  -9.21435543e-06
  -2.77819761e-07  -6.42196414e-09  -1.18327463e-10  -1.78366133e-12
  -2.22487348e-14  -8.34484590e-17  -9.59425880e-17  -1.07363310e-16
  -1.20471395e-16  -1.35327791e-16  -1.52499877e-16  -1.72571646e-16
  -1.95756734e-16  -1.47184432e-16   8.34181730e-15   8.10541310e-13
   6.45429448e-11   4.28137221e-09   2.32361769e-07   1.00521722e-05
   3.33397045e-04   7.96152144e-03   1.21872520e-01   8.98245614e-01] 

TOL = 0.01 







    




ok

======================================================================
FAIL: test_nonnegative (__main__.TestComputeStationary)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/unittest_mc_compute_stationary.py", line 129, in test_nonnegative
    self.assertEqual(np.prod(self.v >= 0-TOL), 1)
AssertionError: 0 != 1

----------------------------------------------------------------------
Ran 4 tests in 0.006s

FAILED (failures=1)

In [5]:

    

%run unittest_mc_compute_stationary.py --N=3 --epsilon=1e-14


    

    




test_left_eigen_vec (__main__.TestComputeStationary) ... ok
test_markov_matrix (__main__.TestComputeStationary) ... ok
test_nonnegative (__main__.TestComputeStationary) ... FAIL
test_sum_one (__main__.TestComputeStationary) ... 





    




N = 3 , epsilon = 1e-14 

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.48233916e-14   1.01149425e+00] 

TOL = 0.01 







    




ok

======================================================================
FAIL: test_nonnegative (__main__.TestComputeStationary)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/unittest_mc_compute_stationary.py", line 129, in test_nonnegative
    self.assertEqual(np.prod(self.v >= 0-TOL), 1)
AssertionError: 0 != 1

----------------------------------------------------------------------
Ran 4 tests in 0.017s

FAILED (failures=1)

In [6]:

    

%run unittest_mc_compute_stationary.py --move='simultaneous' --N=5 --epsilon=1e-15


    

    




test_left_eigen_vec (__main__.TestComputeStationary) ... ok
test_markov_matrix (__main__.TestComputeStationary) ... ok
test_nonnegative (__main__.TestComputeStationary) ... ok
test_sum_one (__main__.TestComputeStationary) ... 





    




N = 5 , epsilon = 1e-15 

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   0.00000000e+00   0.00000000e+00
   0.00000000e+00   0.00000000e+00] 

TOL = 0.01 







    




ok

----------------------------------------------------------------------
Ran 4 tests in 0.003s

OK

In [6]: