数学プログラムを Haskell で書くべき 6 の理由

数学プログラムを  Haskell で書くべき６の理由早稲田大学基幹理工学部  数学科四年石井大海

自己紹介石井大海（いしい・ひろみ）早稲田大学数学科四年（楫元研究室）専門：集合論、計算機代数四月から筑波の数学研究科に進学塩谷研・照井研

発表内容純粋関数型言語 Haskell の紹介今年一年間 Haskell で計算機代数ライブラリを書いてきた経験を踏まえて実装や言語の詳細は後程捕まえていて下さい

純粋関数型言語 Haskell

純粋関数型言語？  Haskell？？

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意

要するに？

Haskell は  記述力が凄い

何がどのように？

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意より数学らしいプログラム

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意より数学らしいプログラム簡潔な記述・見掛け上無限を扱える

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意より数学らしいプログラム簡潔な記述・見掛け上無限を扱える安全性の保証・数学的構造の抽象化

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意より数学らしいプログラム簡潔な記述・見掛け上無限を扱える安全性の保証・数学的構造の抽象化アルゴリズムに適した記述が可能

Haskell の特徴宣言的・純粋関数型代数的データ型遅延評価強力な型システム柔軟な副作用の制御並列・並行処理が得意より数学らしいプログラム簡潔な記述・見掛け上無限を扱える安全性の保証・数学的構造の抽象化アルゴリズムに適した記述が可能効率的なプログラムの記述

宣言的・純粋関数型

宣言的・純粋関数型 C, Ruby……命令的・手続型何をどうするか記述する（命令する）ループが基本 def fib(n) a = b = 1 for i in 1..n a, b = b, a + b end return a end

宣言的・純粋関数型 C, Ruby……命令的・手続型何をどうするか記述する（命令する）ループが基本 Haskell：宣言的、関数型何がどうであるか記述する（宣言する）再帰が基本 def fib(n) a = b = 1 for i in 1..n a, b = b, a + b end return a end fib :: Int → Int → Int fib 0 = 1 fib 1 = 1 fib n = fib (n-2) + fib (n-1)

宣言的・純粋関数型 C, Ruby……命令的・手続型何をどうするか記述する（命令する）ループが基本 Haskell：宣言的、関数型何がどうであるか記述する（宣言する）再帰が基本 def fib(n) a = b = 1 for i in 1..n a, b = b, a + b end return a end fib :: Int → Int → Int fib 0 = 1 fib 1 = 1 fib n = fib (n-2) + fib (n-1) 宣言型の方が数学の「定義」に近い

パターンマッチ帰納的に定義した型の値を、その形によって場合分けできる再帰的にデータを走査する際に便利例：構文木の評価 eval :: Expr -> Int eval (Num i) = i eval (Neg i) = negate (eval i) eval (Plus n m) = eval n + eval m eval (Mult n m) = eval n * eval m eval (If p t f) = if p == 0 then f else t

遅延評価 Haskell は必要になるまで値を評価しない C, Ruby の && や ¦¦ の振舞いに似ているこれらは第一引数の値で結果が決まる場合、第二引数は評価しない Haskell ではデフォルトで全てがこの評価戦略

遅延評価の例可算無限の対象を扱える（列挙できるなら）下：素数から成る無限リスト（注：これはとても非効率）全部読もうとすると当然止まらないが、一部分だけ使う分には問題ない sieve [] = [] sieve (n : ns) = n : sieve [ m | m ← ns, m `mod` n ≠ 0] ! primes = 2 : sieve [3, 5 ..] ! primesUpto100 = takeWhile (≦ 100) primes

遅延評価の応用例有限次元代数とは限らない多項式環の剰余環イデアルが零次元なら剰余環は有限次元、有限個の元の乗算表を用意すれば良い ➡ 乗算表の値を必要になるまで遅延しておけば、無限次元であっても扱うことが出来る！遅延評価は思わぬメモリ漏れを引き起こしもするので、注意が必要

強力な型システム（1）静的型付き：コンパイル時に型をチェック（これはCも同じ） Ruby など：実行時に型エラーで不正終了することがあるこうしたバグが防げるだけで、デバッグが格段に楽になる型推論：コンパイラが型を推定してくれるので、大抵型を自分で書かなくても大丈夫型安全型が付いたプログラムは（何かの意味で）絶対に刺さらない C, C++, ... : 型安全ではない：cast や printf で簡単にSEGV Haskell, OCaml：型安全である

型安全性と型の意味そもそも型とは？ ○ 値の種類を表すタグ int は自然数、char は文字、[int] はint値のリスト etc, etc... ◎ 値や関数の不変条件の表明実際に不変条件として振る舞うには、システムが型安全性を満たす必要がある Haskell の型は十分強力なので、配列の境界条件や変数条件などを型の制約として書ける安全性を静的に保証出来る！

強力な型システム（2）型クラス（アドホック多相）似た操作・演算を持つ型を抽象化出来るオブジェクト指向の mix-in やトレイトに相当応用例：有理数型や有限体型に対して個別に線型代数やグレブナー基底計算を実装するのではなく、一般の体について実装出来る class Field k where zero :: k one :: k (+) :: k -> k -> k (*) :: k -> k -> k recip :: k -> Maybe k

強力な型システム（3）多相型（パラメトリック多相）リストや木など、中身に関係のないデータ構造や関数を総称的に定義出来る以下は二分木とその深さの例 data Tree a = Leaf a | Branch (Tree a) (Tree a) ! depth :: Tree a -> Int depth (Leaf _) = 0 depth (Branch t t') = 1 + max (depth t) (depth t')

柔軟な副作用の制御副作用とは？ Haskell は純粋関数型言語 Haskell の関数は全て数学的な意味での関数純粋 = 同じ引数に対して同じ値を返す最適化や並列化に有利・定義順を気にしなくてよい副作用……関数の純粋性を破るような作用入出力、乱数の生成、etc...

純粋だと何が嬉しいかデータの部分構造を共有出来る関数の定義順を気にする必要がない同じ計算を二度しなくて済む（答えが変わらないから）

Q. 副作用なしでどうやってプログラムを書くの？

モナドによる副作用の制御圏論のモナドの概念を使って副作用を抽象化する副作用を持ち得る値はモナドに包んで扱う型 a を持つけど何らかの意味で副作用を含むものをモナドで包んで m a 型の値として区別する代表的なモナド入出力：IO、非決定性計算：リスト、失敗計算：Maybe、  破壊的代入：ST などなど IO Int : 入出力を伴った整数、  [String]：文字列のリストまたは文字列を使った非決定性計算

モナドを理解せずとも  Haskell は使える型レベルで副作用が区別されることが大事 IO モナドや ST モナドの場合、副作用を外に取り出すことは出来ない C とか Ruby のように、手続型っぽく副作用を書くことが出来る自前の「モナド」を定義でき、DSL が手軽に作れる使う副作用の種類を限定でき、型によって副作用を使う場所が区別出来る代入や破壊的変更によるバグを防げる

例：クィックソート一般的に Haskell の例として出される qsort 実はこれは偽物 pivot してない、in-place じゃない、遅い…… qsort :: [Int] → [Int] qsort [] = [] qsort (a : as) = qsort [x | x ← xs, x ≦ a] ++ [a] ++ qsort [y | y ← ys, y > a]

例：クィックソート破壊的変更による本物のクィックソート（抜粋）（出典：Quicksort in Haskell） qsort' l r = if1 (r > l) $ do i ← auto l j ← auto (r+1) let v = a[l] :: E m a iLTj = i < (j :: E m i) while iLTj $ do while ((i += 1) < mub && a[i] < v) skip while (a[j -= 1] > v) skip if1 iLTj $ a[i] =:= a[j] a[l] =:= a[j] qsort' l (j-1) qsort' i r

副作用：まとめ型により副作用と純粋な関数を分離できる記述するアルゴリズムの性質に応じて、数学的な関数として記述するか、破壊的変更や可変状態を使って手続的に記述するかが選べる手軽に DSL が作れる

時間切れorz

今後の課題数値計算は LAPACK などへの高レベルなバインディングがある記号計算系のライブラリが圧倒的に不足効率的な計算の為にはモジュラー計算や記号線型代数のライブラリの開発が必要効率的なプログラムを書くのにはコツが必要どの言語でもそうだが、ノウハウの蓄積が必要

入門書すごい Haskell 楽しく学ぼう！他言語での経験がある人向けの、肩肘の張らない解説書の邦訳関数プログラミング入門 Haskell で学ぶ原理と技法 Haskell を使って、関数型プログラミングの考え方を基礎から学ぶ良書。型システム入門型システムに関する定番教科書の邦訳。この三冊にはレビュアーとして参加してたりします（自慢） Purely Functional Data Structures 純粋関数型を生かしたデータ構造の解説書

まとめ Haskell は数学的プログラムを記述するのに適す宣言的にも手続的にも書ける遅延評価による計算量の工夫リッチな型システムが安全性や表現力を増加詳細はシンポジウムで！

数学プログラムを Haskell で書くべき 6 の理由

by Hiromi Ishii, Undergraduate Student at Dept. of Math, School of Fundamenta Science and Engineering, Waseda Univ.l on Mar 13, 2014

Statistics

Views

Actions

0 Embeds 0

Accessibility

Categories

Upload Details

Usage Rights

Report content

数学プログラムを Haskell で書くべき 6 の理由 Presentation Transcript