{ "metadata": { "name": "", "signature": "sha256:908f6eae2fbf5cf2bc27f67b5d151cd9d6c1a4376a38f8d1757db9dad5d861d5" }, "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.8 (default, Jul 2 2014, 10:14:46) \n", "[GCC 4.2.1 Compatible Apple LLVM 5.1 (clang-503.0.40)]\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.1'" ] } ], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "# Default values: N=27, epsilon=1e-2, move='sequential', tolerance=1e-2\n", "%run unittest_mc_compute_stationary.py" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "test_left_eigen_vec (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_markov_matrix (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_nonnegative (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "FAIL\n", "test_sum_one (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stdout", "text": [ "N = 27 , epsilon = 0.01 \n", "\n", "v =\n", "[ -1.80778809e-02 -2.45277781e-03 -1.60231716e-04 -6.70987087e-06\n", " -2.02307664e-07 -4.67645846e-09 -8.61657871e-11 -1.29886908e-12\n", " -1.62059348e-14 -6.12789737e-17 -7.10381353e-17 -8.20296809e-17\n", " -9.45029769e-17 -1.08779604e-16 -1.25281548e-16 -1.44569891e-16\n", " -1.66831270e-16 -1.17721634e-16 8.30861752e-15 8.04483364e-13\n", " 6.40581297e-11 4.24921050e-09 2.30616264e-07 9.97666012e-06\n", " 3.30892561e-04 7.90171435e-03 1.20957012e-01 8.91497977e-01] \n", "\n", "TOL = 0.01 \n", "\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "\n", "======================================================================\n", "FAIL: test_nonnegative (__main__.TestComputeStationary)\n", "----------------------------------------------------------------------\n", "Traceback (most recent call last):\n", " File \"/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/unittest_mc_compute_stationary.py\", line 129, in test_nonnegative\n", " self.assertEqual(np.prod(self.v >= 0-TOL), 1)\n", "AssertionError: 0 != 1\n", "\n", "----------------------------------------------------------------------\n", "Ran 4 tests in 0.002s\n", "\n", "FAILED (failures=1)\n" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "%run unittest_mc_compute_stationary.py --N=3 --epsilon=1e-14" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "test_left_eigen_vec (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_markov_matrix (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_nonnegative (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "FAIL\n", "test_sum_one (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stdout", "text": [ "N = 3 , epsilon = 1e-14 \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.49066583e-14 1.01149425e+00] \n", "\n", "TOL = 0.01 \n", "\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "\n", "======================================================================\n", "FAIL: test_nonnegative (__main__.TestComputeStationary)\n", "----------------------------------------------------------------------\n", "Traceback (most recent call last):\n", " File \"/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/unittest_mc_compute_stationary.py\", line 129, in test_nonnegative\n", " self.assertEqual(np.prod(self.v >= 0-TOL), 1)\n", "AssertionError: 0 != 1\n", "\n", "----------------------------------------------------------------------\n", "Ran 4 tests in 0.001s\n", "\n", "FAILED (failures=1)\n" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "%run unittest_mc_compute_stationary.py --move='simultaneous' --N=5 --epsilon=1e-15" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "test_left_eigen_vec (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_markov_matrix (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "test_nonnegative (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stderr", "text": [ "FAIL\n", "test_sum_one (__main__.TestComputeStationary) ... " ] }, { "output_type": "stream", "stream": "stdout", "text": [ "N = 5 , epsilon = 1e-15 \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 -5.55111512e-17 -5.55111512e-17\n", " 5.49560397e-15 2.00000000e+00] \n", "\n", "TOL = 0.01 \n", "\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "ok\n", "\n", "======================================================================\n", "FAIL: test_nonnegative (__main__.TestComputeStationary)\n", "----------------------------------------------------------------------\n", "Traceback (most recent call last):\n", " File \"/Users/oyama/Dropbox/Documents/Teaching/Zemi/2014/test_mc_compute_stationary/unittest_mc_compute_stationary.py\", line 129, in test_nonnegative\n", " self.assertEqual(np.prod(self.v >= 0-TOL), 1)\n", "AssertionError: 0 != 1\n", "\n", "----------------------------------------------------------------------\n", "Ran 4 tests in 0.001s\n", "\n", "FAILED (failures=1)\n" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 } ], "metadata": {} } ] }