{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "import platform\n", "platform.platform()" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 1, "text": [ "'Darwin-13.3.0-x86_64-i386-64bit'" ] } ], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "import sys\n", "print sys.version" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "2.7.6 |Anaconda 1.9.2 (x86_64)| (default, Jan 10 2014, 11:23:15) \n", "[GCC 4.0.1 (Apple Inc. build 5493)]\n" ] } ], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "import numpy\n", "numpy.__version__" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 3, "text": [ "'1.8.0'" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "%run test_mc_compute_stationary.py" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "Elements in each row of P sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "Elements of v sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "All the elements of v are nonnegative ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "FAIL\n", "v is a left eigen vector ... " ] }, { "output_type": "stream", "stream": "stdout", "text": [ "====\n", "Testing with prameter values set 0\n", "\n", "N = 27 , epsilon = 0.01 , move = sequential \n", "\n", "v =\n", "[ -2.48255180e-02 -3.36828636e-03 -2.20038808e-04 -9.21435543e-06\n", " -2.77819761e-07 -6.42196414e-09 -1.18327463e-10 -1.78366133e-12\n", " -2.22487348e-14 -8.34484590e-17 -9.59425880e-17 -1.07363310e-16\n", " -1.20471395e-16 -1.35327791e-16 -1.52499877e-16 -1.72571646e-16\n", " -1.95756734e-16 -1.47184432e-16 8.34181730e-15 8.10541310e-13\n", " 6.45429448e-11 4.28137221e-09 2.32361769e-07 1.00521722e-05\n", " 3.33397045e-04 7.96152144e-03 1.21872520e-01 8.98245614e-01] \n", "\n", "TOL = 0.01 \n", "\n", "====" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "Elements in each row of P sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "Elements of v sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "All the elements of v are nonnegative ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "FAIL\n", "v is a left eigen vector ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "Elements in each row of P sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "Elements of v sum to one ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "All the elements of v are nonnegative ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "v is a left eigen vector ... " ] }, { "ename": "SystemExit", "evalue": "True", "output_type": "pyerr", "traceback": [ "An exception has occurred, use %tb to see the full traceback.\n", "\u001b[0;31mSystemExit\u001b[0m\u001b[0;31m:\u001b[0m True\n" ] }, { "output_type": "stream", "stream": "stdout", "text": [ "\n", "Testing with prameter values set 1\n", "\n", "N = 3 , epsilon = 1e-14 , move = sequential \n", "\n", "P =\n", "[[ 1.00000000e+00 5.00000000e-15 0.00000000e+00 0.00000000e+00]\n", " [ 3.33333333e-01 4.99600361e-15 6.66666667e-01 0.00000000e+00]\n", " [ 0.00000000e+00 3.33333333e-15 6.66666667e-01 3.33333333e-01]\n", " [ 0.00000000e+00 0.00000000e+00 5.00000000e-15 1.00000000e+00]] \n", "\n", "v =\n", "[ -1.14942529e-02 -8.51809045e-17 1.48233916e-14 1.01149425e+00] \n", "\n", "TOL = 0.01 \n", "\n", "====\n", "Testing with prameter values set 2\n", "\n", "N = 5 , epsilon = 1e-15 , move = simultaneous \n", "\n", "P =\n", "[[ 1.00000000e+00 2.50000000e-15 2.50000000e-30 1.25000000e-45\n", " 3.12500000e-61 3.12500000e-77]\n", " [ 1.00000000e+00 2.50000000e-15 2.50000000e-30 1.25000000e-45\n", " 3.12500000e-61 3.12500000e-77]\n", " [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30\n", " 2.77555756e-15 1.00000000e+00]\n", " [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30\n", " 2.77555756e-15 1.00000000e+00]\n", " [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30\n", " 2.77555756e-15 1.00000000e+00]\n", " [ 5.27109897e-77 4.74778387e-61 1.71056941e-45 3.08148791e-30\n", " 2.77555756e-15 1.00000000e+00]] \n", "\n", "v =\n", "[ 1.00000000e+00 2.60902411e-15 0.00000000e+00 0.00000000e+00\n", " 0.00000000e+00 0.00000000e+00] \n", "\n", "TOL = 0.01 \n", "\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "\n", "======================================================================\n", "FAIL: All the elements of v are nonnegative\n", "----------------------------------------------------------------------\n", "Traceback (most recent call last):\n", " File \"//anaconda/lib/python2.7/site-packages/nose/case.py\", line 197, in runTest\n", " self.test(*self.arg)\n", " File \"/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/test_mc_compute_stationary.py\", line 161, in __call__\n", " eq_(np.prod(v >= 0-TOL), 1)\n", "AssertionError: 0 != 1\n", "\n", "======================================================================\n", "FAIL: All the elements of v are nonnegative\n", "----------------------------------------------------------------------\n", "Traceback (most recent call last):\n", " File \"//anaconda/lib/python2.7/site-packages/nose/case.py\", line 197, in runTest\n", " self.test(*self.arg)\n", " File \"/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/test_mc_compute_stationary.py\", line 161, in __call__\n", " eq_(np.prod(v >= 0-TOL), 1)\n", "AssertionError: 0 != 1\n", "\n", "----------------------------------------------------------------------\n", "Ran 12 tests in 0.108s\n", "\n", "FAILED (failures=2)\n" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }