2017年5月10日水曜日

開発環境

行列プログラマー(Philip N. Klein (著)、 松田 晃一 (翻訳)、 弓林 司 (翻訳)、 脇本 佑紀 (翻訳)、 中田 洋 (翻訳)、 齋藤 大吾 (翻訳)、オライリージャパン)の3章(ベクトル空間)、3.2(線形包)、3.2.4(線形結合の線形結合)のクイズ 3.2.12 を取り組んでみる。

クイズ 3.2.12

2 a 1 + a 2 +2 a 3 =3 2 a 3 =0 a 1 +2 a 2 +2 a 3 =0 2 a 1 + a 2 =3 a 1 +2 a 2 =0 3 a 2 =3 a 2 =1 a 1 =2 2 a 1 + a 2 +2 a 3 =0 2 a 3 =2 a 1 +2 a 2 +2 a 3 =0 a 3 =1 2 a 1 + a 2 =2 a 1 +2 a 2 =2 3 a 2 =2 a 2 = 2 3 a 1 =2+ 4 3 = 2 3 2 a 1 + a 2 +2 a 3 =0 2 a 3 =0 a 1 +2 a 2 +2 a 3 =1 a 3 =0 2 a 1 + a 2 =0 a 1 +2 a 2 =1 3 a 2 =2 a 2 = 2 3 a 1 =1 4 3 = 1 3

コード(Emacs)

#!/usr/bin/env python3
# -*- coding: utf-8 -*-

import unittest
from vec import Vec


def list2vec(l):
    return Vec(set(range(len(l))), {i: v for i, v in enumerate(l)})


def gf2mul(u, v):
    def mul(x, y):
        if x == 0 or y == 0:
            return 0
        return 1
    return [mul(x, y) for x, y in zip(u, v)]


def gf2add(u, v):
    def add(x, y):
        if x == y:
            return 0
        return 1
    return [add(x, y) for (x, y) in zip(u, v)]


class Test(unittest.TestCase):

    def setUp(self):
        self.v1 = list2vec([2, 0, 1])
        self.v2 = list2vec([1, 0, 2])
        self.v3 = list2vec([2, 2, 2])

    def tearDown(self):
        pass

    def test_1(self):
        v = list2vec([3, 0, 0])
        a1 = 2
        a2 = -1
        a3 = 0

        self.assertEqual(v, a1 * self.v1 +
                         a2 * self.v2 + a3 * self.v3)

    def test_2(self):
        v = list2vec([0, 2, 0])
        a1 = -2 / 3
        a2 = -2 / 3
        a3 = 1

        self.assertEqual(v, a1 * self.v1 +
                         a2 * self.v2 + a3 * self.v3)

    def test_3(self):
        v = list2vec([0, 0, 1])
        a1 = -1 / 3
        a2 = 2 / 3
        a3 = 0

        self.assertEqual(v, a1 * self.v1 +
                         a2 * self.v2 + a3 * self.v3)

if __name__ == '__main__':
    unittest.main()

入出力結果(Terminal, IPython)

$ ./sample3.py -v
test_1 (__main__.Test) ... ok
test_2 (__main__.Test) ... ok
test_3 (__main__.Test) ... ok

----------------------------------------------------------------------
Ran 3 tests in 0.001s

OK
$

0 コメント:

コメントを投稿

Comments on Google+: