CakePHPを学ぶ(1)

　本研究室ではもっぱらスクラッチからPHPとMySQLを使ってWebアプリを作成してきましたが，ボチボチMVC(Model-View-Control)思想に基づくPHPフレームワークにも取り組む必要が出てきた感じています。スクラッチビルドも悪くないですが，正直同じ教材を毎年使うのも飽きたので更新したいという事情も手伝っています。

　ということで昨年12月からCakePHPで遊んでみましたが，Tutorialは手っ取り早く勉強できてありがたいものの，作りたいWebアプリを自在に組み立てるのは思いのほか難儀で，Bakeが作った雛形から脱していない代物が多いのも頷けます。年寄なので，Visual Basicが流行った頃を思い出しますね。結局，作り込むためには地道なプログラミングをしないといかんと。

　ということで本日反省がてら作ったPowerPoint資料を上げておきます。

　ということで，

ユーザ認証機能のカスタマイズ（任意のテーブル名でユーザ管理できるように）
全体的なデザインの変更とカスタマイズ
Webサーバへのアップロード方法の模索（ユーザ単位のWebアプリとして外部公開したい）
第5章で制作する機能を全部盛り込んだサイトを構築

というところまで，今年度内に何とかしたいと考えています。ここに書いたメモを元に教材作成までできれば御の字ですね。

2018年11月15日2019年3月25日

遠江総合高等学校：オープンデータとWebアプリケーション

「オープンデータとWebアプリケーション」
静岡理工科大学　情報学部　コンピュータシステム学科
幸谷　智紀（こうや　とものり）

０．【宣伝】静岡理工科大学情報学部について
静岡理工科大学情報学部のWebページ

１．オープンデータとは？
「オープンデータ(Open Data)」とは，誰でも閲覧できる状態で広く公開されているデータのことです。今の時代はデータを扱う主体はコンピューターですので，オープンデータと言っても紙媒体ではなく，電子化されているもの，即ち，コンピューターが直接読み書きできるファイルの形式になっているものを指すことが殆どです。
電子化されたデータは，コンピューター同士が繋がっているインターネット(The Internet)の上でやり取りしやすい形式ですので，公開するにしても，インターネット上で誰しも読める状態になっている方が，読む方としても都合が良いわけです。インターネットでは様々なサービスが展開されていますが，現在ではWeb(World Wide Web)を通じたサービスが普通になっていますので，オープンデータもWeb上で，つまり，データを持っている会社，行政組織，個人，各種団体のWebページ(ホームページ）上に置き，誰でもブラウザを通じて読むことができるようにしています。こうしたオープンデータを活用したWeb上のサービス，いわゆるWebアプリケーションが様々な形で展開されています。

以下，オープンデータとは，Web上で誰でもブラウザから閲覧できるようになっている電子ファイル状態のデータである，と定義します。最近では，行政機関が積極的にデータを公開するようになっています。静岡県も積極的にオープンデータの活用を後押ししていますし，静岡県内の市，静岡市，浜松市では下記のところでオープンデータにアクセスできるようになっています。

オープンデータは，単なる公開情報という存在を超えるものです。コンピューターで利用しやすいファイルの形式になっているということは，そのままソフトウェアで加工できることになります。これをうまく活用することで，より便利なWebアプリケーションを作ることができます。

２．オープンデータをどのように活用するか？
ここでは，簡単なオープンデータの活用方法を具体例で示します。

皆さんは，「国民の祝日」，例えば正月，成人の日，建国記念日，春分の日，がどこでどのように決まっているか，ご存知ですか？　検索すればすぐに分かりますが，これは政府が法律に基づいて決めているもので，内閣府に解説ページがあります。ただ，これは情報を人間に見やすく表示しているだけで，余り使い勝手の良いものではありません。そこで，近年の国民の祝日を表計算ソフトで扱いやすい形式で公開しています。これが国民の祝日のオープンデータです。この形式であれば，表計算ソフトだけでなく，様々なソフトウェアに利用することができるようになります。例えば，下記の図は，FullcalendarというJavaScriptで構築されたスケジュール帳に，国民の祝日を上のURLからファイルを読み取り，自動的に挿入したものです(サンプル)。

但し，これはあまり使い勝手の良いものではありません。元々の国民の祝日ファイルは，単なるテキストファイルなので，データそのものを使いやすい形で形成し，流し込む処理が必要になります。例えば，上記のFullcalendarに取り込むためには，次のようなPHPスクリプト(プログラム)の手助けが必要です。

<?php
	// holiday.php: J国民の祝日を日本政府サイトから取得する
	// Copyright (c) 2017 Tomonori Kouya

	// 初期設定&データベースオープン
	require "database.php";

	// URL
	$holiday_url = 'http://www8.cao.go.jp/chosei/shukujitsu/syukujitsu.csv';
	$holiday_url_charset = 'SJIS-win';

	// validate_date
	// https://stackoverflow.com/questions/19271381/correctly-determine-if-date-string-is-a-valid-date-in-that-format
	function validate_date($date)
	{
		$d = DateTime::createFromFormat('Y-m-d', $date);
		return $d && $d->format('Y-m-d') === $date;
	}

	// http://php.net/manual/ja/function.str-getcsv.php
	$array = array_map('str_getcsv', file($holiday_url));

	// 休日配列化
	$holiday_array = array();
	foreach($array as $index => $tmp_array)
	{
		if(validate_date($tmp_array[0]))
		{
			//echo $tmp_array[0] . ' -> ' . mb_convert_encoding($tmp_array[1], 'UTF-8', $holiday_url_charset) . "\n";
			$title = mb_convert_encoding($tmp_array[1], 'UTF-8', $holiday_url_charset);
			$holiday = new Event([
				//'id'     => $index,
				'title'  => $title,
				'allDay' => true,
				'start'  => $tmp_array[0],
				'end'    => $tmp_array[0],
			], $default_timezone);
			$holiday_array[] = $holiday->toArray();
		}
			//echo $holiday_array[0] . ' -> ' . mb_convert_encoding($holiday_array[1], 'UTF-8', $holiday_url_charset) . "\n";
	}

	// JSON形式で送り込み
	echo json_encode($holiday_array);

?>

<?php

// holiday.php: J国民の祝日を日本政府サイトから取得する

// 初期設定&データベースオープン

require "database.php";

// URL

$holiday_url = 'http://www8.cao.go.jp/chosei/shukujitsu/syukujitsu.csv';

$holiday_url_charset = 'SJIS-win';

// validate_date

// https://stackoverflow.com/questions/19271381/correctly-determine-if-date-string-is-a-valid-date-in-that-format

function validate_date($date)

{

$d = DateTime::createFromFormat('Y-m-d', $date);

return $d && $d->format('Y-m-d') === $date;

}

// http://php.net/manual/ja/function.str-getcsv.php

$array = array_map('str_getcsv', file($holiday_url));

// 休日配列化

$holiday_array = array();

foreach($array as $index => $tmp_array)

{

if(validate_date($tmp_array[0]))

{

//echo $tmp_array[0] . ' -> ' . mb_convert_encoding($tmp_array[1], 'UTF-8', $holiday_url_charset) . "\n";

$title = mb_convert_encoding($tmp_array[1], 'UTF-8', $holiday_url_charset);

$holiday = new Event([

//'id' => $index,

'title' => $title,

'allDay' => true,

'start' => $tmp_array[0],

'end' => $tmp_array[0],

], $default_timezone);

$holiday_array[] = $holiday->toArray();

}

//echo $holiday_array[0] . ' -> ' . mb_convert_encoding($holiday_array[1], 'UTF-8', $holiday_url_charset) . "\n";

}

// JSON形式で送り込み

echo json_encode($holiday_array);

このように，直接データを取得して自動的に加工できるようになると，人間を介しての作業が不要になります。

こういうWeb技術を利用するプログラミングの詳細については，例えば本研究室の教材などを参照して下さい。

オープンデータの利用のためには，もう少し，利用者の立場に立って，使い勝手の良い形式で必要な分量のデータを受け取れるようにする必要があります。最悪な形式は，いわゆる「ネ申Excel（三重大学・奥村晴彦）」と呼ばれる，書類をそのまま表計算ソフトウェアで成形しただけのシロモノで，人間が手作業でチェックするには便利かもしれませんが，コンピューターにとっては最も加工の難しいものです。オープンデータは人間よりも，コンピューターにとって使い勝手の良いシンプルな形式であることが望まれます。

しかし，もっと使い勝手の良いオープンデータは，情報を加工したい側が欲しいデータだけをいつでも取得できるようになっているものです。そのためには，上記のFullcalendarに祝日データを渡すプログラムのように，オープンデータと利用者の間で仲介し，利用者の要求に応じた返答を返すような仕組みを作らなくてはなりません。これを実現しているのが，Web API (Application Programming Interface)と呼ばれるものです。

３．Web APIの実例
最も良く知られたweb APIは，GoogleやYahoo! Japan，Bingといった検索エンジンです。

どの検索エンジンも，URLの中に検索すべき単語を埋め込んであります。例えばGoogleの場合は

https://www.google.co.jp/#q=遠江総合高等学校

となっています。このように，ブラウザからアクセスする際に，利用者が欲しいデータの範囲（この場合は「遠江総合高等学校」を含むWebページの一覧）を指定することで，必要のないデータまで取り込むことがなくなります。

現在では様々なデータを，Web APIを通じて取り込むことができるようになっています。Google map（地図情報）やぐるなび（グルメ情報）なども，細かくデータを指定して取り込むことができるようになっています。

その他，下記のようなWeb APIもあります。

リクルート：Small Talk API
Yahoo! Japan: テキスト解析API

４．openBDを基盤とする書誌データの活用
オープンデータとWeb API活用の事例として，現在本研究室で取り組んでいるopenBDについて紹介します。
2017年初めに広く公開された本のオープンデータ化を進めるプロジェクトで，2017年6月現在では約90万件の本のデータが取得できます。今のところ，本のISBN番号一覧と，ISBN番号に紐づけされた書誌データ（表紙画像も含む）を取得できるWeb APIが提供されており，基本的にはユーザーがダウンロードして使用する全データを提供することを主目的としています。
本研究室でも，簡単なWebのインターフェースを提供しています。

openBD API

本一冊一冊ごとに異なるISBN番号が振られており，それが分かればその本のデータだけを入手することが可能です。しかし，普通本を買うときにISBN番号で指定することはありません。普通は，本のタイトル，著者名，出版社名，定価等で探すものです。

それを可能にするため，本研究室では全件検索ができるような仕組みを作り，公開しています。使い勝手はまだ改良の余地がありますが，これも改良されたopenBDのWeb API化の一つと言える仕事になっています。

５．おわりに
インターネットが当たり前になった今では，様々なデータを自由に取得でき，もっと便利に使うためのWeb APIの利用が盛んになっています。今後は，これらのオープンデータやWeb APIを通じてより高度なデータ解析を自動的に，例えばAIを使って行うことが普通になっていきます。
人間の側はもっと賢く，これらのオープンデータやWeb API，そしてその分析結果を使って高度な知を使いこなしていくことになります。プログラミングを行うだけでなく，それをより良い方向に生かす技術を，今後の大学生活で身につけて下さい。

2018年10月31日2018年10月31日

PHP 5.6 security update終了

　WordPress用のセキュリティプラグインを提供しているWordfenceから，PHP 5.6セキュリティアップデート終了のアナウンスがありました。

　記事によれば，PHPの最新バージョンは7.2で，7.0のアップデートもあと一月で終了とのこと。一つのバージョン(Version x.x相当)につき3年のスパンでサポートしており，最初の2年が過ぎるとセキュリティアップデートのみ1年行ってサポート完全終了，という流れだそうです。

　で，慌てて自分の管理サーバも見ましたが，このblogのあるサーバはPHP 5.4という更に古い奴を使っており，もう一つの方も7.0系統でした。どちらもLinux distributionサイトからのアップデートが出ていないので，すぐにセキュリティの危機になるわけじゃありませんが，OSごとアップグレードしてやろうかという気分になってきましたが，そうなると，次のような手順でやりたいものです。

バックアップ用のサーバを一台用意し，そちらを更新時の新しいOS環境にセッティングする
古いファイルを丸ごとバックアップ用サーバに移し，コンテンツがきちんと動くことを確認（特にWeb, MySQLサーバ)する。その後，現サーバマシンのDNS情報をバックアップ用サーバに書き換えておく
現サーバマシンに新しいOS環境を構築し，バックアップ用サーバからコンテンツを書き戻す
きちんと動作することを確認したら，DNS情報を現サーバマシンに書き換え。

　現マシンのOS環境を毎書き換えるには時間がかかります（数時間程度）し，スッピンのOS環境を無法The Internet空間にさらし者にするのは避けるべきです。また，新しいOS環境構築の練習としてバックアップ用サーバを構築するのは役立ちます。使い終わったバックアップ用サーバの転用も考えていますので，初期投資として無駄になることはないでしょう。ちなみに，サーバ構築手順の解説は「ネコでもわかる！さくらのVPS講座」がお勧めです。

　ということで，来年3月ぐらいにはバックアップ用サーバを契約し，新年度に入る前にさっさと移行しようと計画しています。ま，見かけは変わらないでしょうが，本サーバの説明書きは大幅に変わりそうです。

2018年10月18日2018年10月18日

Update or perish

　先週の「数値解析2」で，実習環境で使っているLinuxサーバのSamba（ネットワーク共有フォルダを提供するサーバソフト）に接続できない（「\\サーバIPアドレスが見つかりません」というエラーが出る）という現象が一部の受講生に出ました。不思議なのは，私やもう一人の受講生は普通に接続できること。接続できるグループとできないグループに分かれたことで，その時間内にはトラブルが解決できませんでした。TCP/IP接続には問題がないことは，Ping応答やTeraTerm接続ができることで確認済みです。

　後であれこれ検索してみたところ，実習用Linuxサーバの設定が古く（6年前に設定したまま放置状態），Windows 10 1709で無効化されたSMB 1.0プロトコルしか使えない状態になっていることが判明しました。Microsoftに言わせれば

SMBv1 のみをサポートする製品でこの問題を回避するには、製品の製造元に連絡し、SMBv2.02 以降のバージョンをサポートするソフトウェアまたはファームウェアの更新プログラムを依頼してください。
https://support.microsoft.com/ja-jp/help/4034314/smbv1-is-not-installed-by-default-in-windows

だそうで，ぐうの音も出ない正論です。

　ということで，Sambaの設定(/etc/samba/smb.conf)に”max protocol = SMB2″を追記して，SMB2接続ができるように変更しました。これで受講生全員が接続できるようになって一安心です。実際，PowerShellの管理者モードでGet-SmbConnectionコマンドを使うとSMB 2.0.2で接続されていることが確認できます。

　古くからSambaやNASを使っていたり，Windows 8以下から10にアップデートしたマシンはSMB 1.0が最初から有効になっているので，問題なく接続できたのはそれが理由でしょう。あるいはWindows Updateをサボっているか，です。

　ICT技術はすさまじい勢いで進化しており，セキュリティ確保のための基盤ソフトウェアやネットワークの更新も頻繁に行われます。更新に伴う不都合が起きた時には，古い環境をそのまま維持するのではなく，Microsoftが主張するように，新しい環境に適応していく方向でトラブル対処を行うようにしましょう。”Update or perish”（アップデートせよ，できなきゃ滅びよ）はICTの世界で生きていく技術者が肝に銘じておくべき格言です。

2018年4月18日2018年4月18日

PCケースファンの交換

　本研究室のデスクトップマシンは大概古ぼけたPCケースに収まっています。既に３世代以上更新を重ねているものも少なくありません。見かけは古いですが，中身はマザーボードごと新品に入れ替えつつ，大事に使っています。

　私がメインに使っているマシンも，最初はIntel Core 2 Quadプロセッサの入ったシステムでしたが，今はCore i7になっています。10年以上使っていることになりますね。さすがにケースファンは軸のグリスが固まってしまっているせいか，全く動かなくなっていました。そのままでも使えますが，夏場の動作が心配ですので，付け替えることにしました。緑色に発行する120mmタイプ。

　LEDに電力使われるのも癪なので，発光しないものを探したのですが見つかりませんでしたので，仕方なくこれにしました。まぁ，起動していることが後ろから見ても分かるというのはメリットなんでしょうが，派手に見せて何かいいことがあるのかなと，年寄としては不満であります。

　光物といえば，マザーボードもGPUボードも，光らないものを見つけるのが面倒なほど，LED付きのものが当たり前になってますね。このマシンも既にそうなっています。

　透明アクリル板付きのケースにして，PCショップの店頭に並べておくのにはふさわしいのかもしれませんが，純然たる仕事用マシンとしては無駄以外の何物でもありません。光らない地味なマシンを追究していきたいところですが，現状難しいとなると，本研究室の卒研用マシンも全部光物になる日も近そうです。

2018年1月18日

日経新聞電子版の普及

　2018年1月18日現在，日経新聞電子版の会員総数は3,843,506だそうです。そのうち有料版の購読者数は558,900，紙版を購読していない純粋な電子有料版購読者数は345,743とのこと。私もこの仲間です。新聞紙不要派。

　ざっとグラフにしてみると，電子版購読者数のうち，大半は無料で見れられる範囲のみで満足しているようで，有料購読者のうち電子版だけ読んでいるのは1/10程度ですね。ただ，月々の購読価格は紙版と変わらないので，配送印刷紙代諸経費を考えると，大変利益率の高そうな購読者層ということになります。紙版を併用して読んでいる人より多いというのも最近のメディア状況を犯得していると言えそうですね。

　グラフにしてみると一目瞭然ですが，きちんとした情報を得るためにはそれなりの対価を払わないとダメという時代になったと言えそうです。

2017年9月5日2019年4月4日

[模擬講義] コンピュータの計算は正しくない！

オープンキャンパス模擬講義

「コンピュータの計算は正しくない！」

静岡理工科大学　情報学部
コンピュータシステム学科
幸谷　智紀（こうや　とものり）

１．「数」とは？

　現在のICT社会で使用される「データ(data)」は例外なく全て「数」を表現したものです。小学校では自然数，分数（有理数），小数を，中学校では負の数と整数を，高校では無理数を習いますが，それらの数を総称して「実数(real number)」と呼びます。実数を幾何学的に表現する時には直線＝数直線を使用します。全ての数は，この数直線上のどこかに「点(point)」として存在しています。

　自然数，整数，有理数は全て分数として表現できる数です。下記の図に示す通り，分数は無数に存在します。

　ちょっと脱線しますが，下記のような問題はご存知でしょうか？

　解答については検索して調べて頂ければ何らかの手かがりが得られるでしょうが，実際の分数の計算はこのようにしてはいけませんね。どのように計算できるのか，ちょっと試してみましょう。

分数の計算フォーム

　実数にはもう一つ，分数としては表現できない「無理数(irrational number)」があります。有名なものとしては円周率 $\pi = 3.1415...$ や平方根 $\sqrt{2} = 1.4142...$ があります。これも実数の一種なので，数直線上の点として実在しています。

　有理数(分数)を小数で表現すると，次のように必ず同じ数のパターンが繰り返し出てきます。

\[\begin{split}
\frac{1}{2} &= 0.50000\cdots = 0.5\dot{0} \\
\frac{1}{3} &= 0.33333\cdots = 0.\dot{3}
\end{split}\]

　無理数では同じパターンの数の羅列が永久に続くことはありません。ちょっと計算してみましょう。

Try! MPFR

　従って，実数を一言で表現すると「循環したりしなかったりする）無限桁の小数として表わすことのできる数」ということになります。何故無限桁が必要になるのかは，微分積分（数学II, III)に登場する「極限」というものと深い関係がありますが，それは大学に入ってからじっくり考えてみて下さい。

　現在の科学技術では実数を土台とする計算が不可欠です。そしてほとんどの計算はコンピューターの中で行います。では，次にコンピューターの中身について見ていくことにしましょう。

２．コンピュータの仕組みと「データ」

　今のコンピューターの筐体とぱかっと開けると，マザーボードと呼ばれるデジタル回路基板の上に様々な集積回路が載っていることが分かります。下記は本研究室で故障したマシンのマザーボードの写真です。

　マザーボードの上に乗っている集積回路のうち，重要なものがCPU(Central Processing Unit)とメモリ(RAM)です。コンピューターで行われる処理は，全てのデータが「数」として表現されてメモリに格納され，処理内容に応じてメモリ内のデータが読み出されて処理され，またメモリに書き戻されます。

　このCPUとメモリとのやり取りに要する時間，CPU内の処理（計算）に要する時間が短ければ短いほど，コンピューターの性能が高いと言えます。とても高速な処理が可能になっていますが，困ったことに，メモリの大きさは有限で，処理時間も短いとはいえゼロにはなりません。「無限の長さを持つ実数」をそのままコンピューターのメモリに格納することは不可能ですし，仮にできたとしても，無限の長さの処理時間を要することになります。「無限大(∞)」とは果てのないことを意味し，それ故に完ぺきな数学，特に微分積分のような解析学という体系が成り立つわけですが，これではコンピューターには実数の計算が出来ないことになってしまいます。

３．「誤差」を含まざるを得ない科学技術計算

　仕方がないので，コンピューターの中では無限桁の小数を適当な長さに切って，有限桁の小数（＝有理数の一種）として扱うことにします。この適当長さに切る操作のことを「丸め（Round-off)」と呼び，丸めに伴って起きる元の正しい値（真値）とのズレを「誤差(error)」と呼びます。

　この丸めに伴う誤差を，関数グラフアプリを使ってみていくことにしましょう。

MPFRgraph

コンピューターは計算が高速なので，いくらでも長い桁の計算ができそうですが，長くなればなるほど計算時間はかかるようになりますし，科学技術開発競争を世界中で行っている昨今，あまり誤差ばかり気にして長い桁の計算ばかりすることは難しいのが現状です。従って，最近のAIの中核技術であるDeep Learningではなるべく短い桁の小数で計算しようとしています。

　しかし，いくら計算が速くても，途中で誤差が拡大されて，不正確な値しか出てこないようでは困ります。一般には，丸める桁数を長くすれば誤差も小さくできるわけですが，とてもたちの悪い問題があることは良く知られていて，いくら桁を長く取っても足りなくなるということもあります。

５．誤差評価付き計算の例

　性質の良くない問題の例として，「複雑系(Chaos system)」というものがあります。一番シンプルなものとして，ロジスティック写像で定義される次のような数列 $x_0, x_1, ..., x_i, x_{i+1}, ...$ の生成問題を取り上げましょう。

　出発値 $x_0 = 0.7501$ として，次の式の右辺を計算して次の $x_1$ を導出します。

\[ x_{i+1} := 4x_i (1 – x_i)　\]

これで実数列 $\{x_i\}^{100}_{i=1}$ を，下記のCプログラム

#include <stdio.h>

main()
{
    int i;
    float x[102]; /* 単精度 */

    x[0] = 0.7501;

    for(i = 0; i <= 100; i++)
    {
        x[i+1] = 4 * x[i] * (1 - x[i]);
        if(i%10 == 0) 
        	printf("%5d %25.17e\n", i, x[i]);
    }
}

#include <stdio.h>

main()

{

int i;

float x[102]; /* 単精度 */

x[0] = 0.7501;

for(i = 0; i <= 100; i++)

{

x[i+1] = 4 * x[i] * (1 - x[i]);

if(i%10 == 0)

printf("%5d %25.17e\n", i, x[i]);

}

を用いてIEEE754単精度(約7桁)及び倍精度(約15桁)で計算した結果は次のようになります。

    i      単精度          　倍精度
    0   7.501000e-01 7.50099999999999989e-01
   10   8.445168e-01 8.44495953602201199e-01
   20   1.229039e-01 1.42939724528399537e-01
   30   4.977781e-01 8.54296020314155413e-01
   40   5.816439e-01 7.74995884542777125e-01
   50   5.580120e-01 7.95132827423636751e-02
   60   2.287482e-02 2.73872762849587226e-01
   70   9.985481e-01 9.97021556611396687e-01
   80   9.421800e-01 3.52785425069070790e-01
   90   2.129389e-01 6.35558111134618908e-01
  100   7.568640e-01 5.15390006286616020e-01

i 単精度　倍精度

0 7.501000e-01 7.50099999999999989e-01

10 8.445168e-01 8.44495953602201199e-01

20 1.229039e-01 1.42939724528399537e-01

30 4.977781e-01 8.54296020314155413e-01

40 5.816439e-01 7.74995884542777125e-01

50 5.580120e-01 7.95132827423636751e-02

60 2.287482e-02 2.73872762849587226e-01

70 9.985481e-01 9.97021556611396687e-01

80 9.421800e-01 3.52785425069070790e-01

90 2.129389e-01 6.35558111134618908e-01

100 7.568640e-01 5.15390006286616020e-01

同じ値になる名図ですが， $x_{20}$ あたりでもう単精度の値は最初の一桁目しか合っていません。 $x_{30}$ 以降はもう全く違う数です。

　このような場合はもっと桁数を増やして計算する必要があります。50桁計算すると下記のようになります。

  0, 7.501e-01
 10, 8.4449595360221744753714870256154137267401952291229e-01
 20, 1.4293972451230765528428175723130628307376969200214e-01
 30, 8.5429600370442189166613118425888744374797589008230e-01
 40, 7.7497575311820124128022346126421168481511313675586e-01
 50, 9.3375332197703029055189668015168234762823013462054e-02
 60, 4.0822016829087813187102766181952686730359229911213e-01
 70, 7.1511999705058574897099346417593218350311240577887e-02
 80, 4.6325330290077571055993467458320747053550671068192e-01
 90, 1.3344050120868840464692012150021107621136608489291e-03
100, 7.8817989371509906806704770440086762650658767283966e-02

0, 7.501e-01

10, 8.4449595360221744753714870256154137267401952291229e-01

20, 1.4293972451230765528428175723130628307376969200214e-01

30, 8.5429600370442189166613118425888744374797589008230e-01

40, 7.7497575311820124128022346126421168481511313675586e-01

50, 9.3375332197703029055189668015168234762823013462054e-02

60, 4.0822016829087813187102766181952686730359229911213e-01

70, 7.1511999705058574897099346417593218350311240577887e-02

80, 4.6325330290077571055993467458320747053550671068192e-01

90, 1.3344050120868840464692012150021107621136608489291e-03

100, 7.8817989371509906806704770440086762650658767283966e-02

100桁計算すると，

  0, 7.501e-01
 10, 8.44495953602217447537148702561541372674019522912291280007848388118986290823419270975121157e-01
 20, 1.42939724512307655284281757231306283073769692002141121387637872878254442652405447154653323e-01
 30, 8.54296003704421891666131184258887443747975890082277146848716525830873157517999293072280730e-01
 40, 7.74975753118201241280223461264211684815113136725464353115716739919423777262939033931884043e-01
 50, 9.33753321977030290551896680151682347628230351487453020066411681729685131997586455460349178e-02
 60, 4.08220168290878131871027661819526867303629812881299069587191914174238088850476295365965331e-01
 70, 7.15119997050585748970993464175932183301720988330772170512558519903814407400108286900737962e-02
 80, 4.63253302900775710559934674583207430627755245547675072149154489741355095166291517037937750e-01
 90, 1.33440501208688404646920121499911906832767847107910033479094701619779865654038757784875352e-03
100, 7.88179893715099068067047704626992647240094450346038441154434366642728380004689090274960709e-02

0, 7.501e-01

10, 8.44495953602217447537148702561541372674019522912291280007848388118986290823419270975121157e-01

20, 1.42939724512307655284281757231306283073769692002141121387637872878254442652405447154653323e-01

30, 8.54296003704421891666131184258887443747975890082277146848716525830873157517999293072280730e-01

40, 7.74975753118201241280223461264211684815113136725464353115716739919423777262939033931884043e-01

50, 9.33753321977030290551896680151682347628230351487453020066411681729685131997586455460349178e-02

60, 4.08220168290878131871027661819526867303629812881299069587191914174238088850476295365965331e-01

70, 7.15119997050585748970993464175932183301720988330772170512558519903814407400108286900737962e-02

80, 4.63253302900775710559934674583207430627755245547675072149154489741355095166291517037937750e-01

90, 1.33440501208688404646920121499911906832767847107910033479094701619779865654038757784875352e-03

100, 7.88179893715099068067047704626992647240094450346038441154434366642728380004689090274960709e-02

となり，大分同じ数に近づいてきたことが分かります。

　このように，丸めによる誤差の影響が大きく出る「悪条件問題」というものが，科学技術計算にはたまに出てきます。計算途中で誤差がどの程度は行っているのかを自動的に調べる「区間演算(interval arithmetic)」という技術もありますが，このような悪条件問題の場合は大した役には立たず，結局桁数を多くしないとまともな値を得ることはできません。下記の例は，MPFIと呼ばれる多倍長計算を区間演算に利用できるソフトウェアを使ったロジスティック写像の掲載例です。一応，プログラムもつけておきましょう。

$ ./logistic_mpfi
prec(bits) = 128
    0, [7.5009999999999999e-1,7.5010000000000001e-1] 3.918e-39
   10, [8.4449595360221744e-1,8.4449595360221745e-1] 4.865e-33
   20, [1.4293972451230765e-1,1.4293972451230766e-1] 3.014e-26
   30, [8.5429600370442189e-1,8.5429600370442190e-1] 5.288e-21
   40, [7.7497575311819887e-1,7.7497575311820361e-1] 6.112e-15
   50, [9.3375329714176199e-2,9.3375334681229861e-2] 5.319e-08
   60, [4.0561739523590902e-1,4.1082572849550305e-1] 1.276e-02
   70, [-2.9373100713009889e43,2.1377489377594625e43] 5.075e+43
   80, [-1.2245050211453045e45127,8.9118419393669329e45126]2.116e+45127
   90, [-9.5278547330878072e46210753,6.9342905040200581e46210753]1.646e+46210754
  100, [-@Inf@,@Inf@]       inf

$ ./logistic_mpfi
prec(bits) = 256
    0, [7.5009999999999999e-1,7.5010000000000001e-1] 1.151e-77
   10, [8.4449595360221744e-1,8.4449595360221745e-1] 1.430e-71
   20, [1.4293972451230765e-1,1.4293972451230766e-1] 8.857e-65
   30, [8.5429600370442189e-1,8.5429600370442190e-1] 1.554e-59
   40, [7.7497575311820124e-1,7.7497575311820125e-1] 1.796e-53
   50, [9.3375332197703029e-2,9.3375332197703030e-2] 1.563e-46
   60, [4.0822016829087813e-1,4.0822016829087814e-1] 3.749e-41
   70, [7.1511999705058574e-2,7.1511999705058575e-2] 2.244e-34
   80, [4.6325330290077571e-1,4.6325330290077572e-1] 3.633e-29
   90, [1.3344050120868840e-3,1.3344050120868841e-3] 1.322e-20
  100, [7.8817989371509897e-2,7.8817989371509917e-2] 2.348e-16

$ ./logistic_mpfi

prec(bits) = 128

0, [7.5009999999999999e-1,7.5010000000000001e-1] 3.918e-39

10, [8.4449595360221744e-1,8.4449595360221745e-1] 4.865e-33

20, [1.4293972451230765e-1,1.4293972451230766e-1] 3.014e-26

30, [8.5429600370442189e-1,8.5429600370442190e-1] 5.288e-21

40, [7.7497575311819887e-1,7.7497575311820361e-1] 6.112e-15

50, [9.3375329714176199e-2,9.3375334681229861e-2] 5.319e-08

60, [4.0561739523590902e-1,4.1082572849550305e-1] 1.276e-02

70, [-2.9373100713009889e43,2.1377489377594625e43] 5.075e+43

80, [-1.2245050211453045e45127,8.9118419393669329e45126]2.116e+45127

90, [-9.5278547330878072e46210753,6.9342905040200581e46210753]1.646e+46210754

100, [-@Inf@,@Inf@] inf

$ ./logistic_mpfi

prec(bits) = 256

0, [7.5009999999999999e-1,7.5010000000000001e-1] 1.151e-77

10, [8.4449595360221744e-1,8.4449595360221745e-1] 1.430e-71

20, [1.4293972451230765e-1,1.4293972451230766e-1] 8.857e-65

30, [8.5429600370442189e-1,8.5429600370442190e-1] 1.554e-59

40, [7.7497575311820124e-1,7.7497575311820125e-1] 1.796e-53

50, [9.3375332197703029e-2,9.3375332197703030e-2] 1.563e-46

60, [4.0822016829087813e-1,4.0822016829087814e-1] 3.749e-41

70, [7.1511999705058574e-2,7.1511999705058575e-2] 2.244e-34

80, [4.6325330290077571e-1,4.6325330290077572e-1] 3.633e-29

90, [1.3344050120868840e-3,1.3344050120868841e-3] 1.322e-20

100, [7.8817989371509897e-2,7.8817989371509917e-2] 2.348e-16

// logistic 写像
// MPFR & MPFI版
#include <stdio.h>
#include "mpfr.h"
#include "mpfi.h"

int main()
{
	int i;
	unsigned long prec;
	mpfi_t x[102];
	mpfr_t relerr;

	printf("prec(bits) = "); scanf("%ld", &prec);
	mpfr_set_default_prec(prec);

	// 初期化
	mpfr_init(relerr);
	for(i = 0; i < 102; i++)
		mpfi_init(x[i]);

	// 初期値
	//x[0] = 0.7501;
	mpfi_set_str(x[0], "0.7501", 10);

	for(i = 0; i <= 100; i++)
	{
		if((i % 10) == 0)
		{
			printf("%5d, ", i);
			mpfi_out_str(stdout, 10, 17, x[i]);
			mpfi_diam(relerr, x[i]);
			mpfr_printf("%10.3RNe\n", relerr);
		}

		//x[i + 1] = 4 * x[i] * (1 - x[i]);
		mpfi_ui_sub(x[i + 1], 1UL, x[i]);
		mpfi_mul(x[i + 1], x[i + 1], x[i]);
		mpfi_mul_ui(x[i + 1], x[i + 1], 4UL);

	}

	// 消去
	mpfr_clear(relerr);
	for(i = 0; i < 102; i++)
		mpfi_clear(x[i]);

	return 0;
}

// logistic 写像

// MPFR & MPFI版

#include <stdio.h>

#include "mpfr.h"

#include "mpfi.h"

int main()

{

int i;

unsigned long prec;

mpfi_t x[102];

mpfr_t relerr;

printf("prec(bits) = "); scanf("%ld", &prec);

mpfr_set_default_prec(prec);

// 初期化

mpfr_init(relerr);

for(i = 0; i < 102; i++)

mpfi_init(x[i]);

// 初期値

//x[0] = 0.7501;

mpfi_set_str(x[0], "0.7501", 10);

for(i = 0; i <= 100; i++)

{

if((i % 10) == 0)

{

printf("%5d, ", i);

mpfi_out_str(stdout, 10, 17, x[i]);

mpfi_diam(relerr, x[i]);

mpfr_printf("%10.3RNe\n", relerr);

}

//x[i + 1] = 4 * x[i] * (1 - x[i]);

mpfi_ui_sub(x[i + 1], 1UL, x[i]);

mpfi_mul(x[i + 1], x[i + 1], x[i]);

mpfi_mul_ui(x[i + 1], x[i + 1], 4UL);

}

// 消去

mpfr_clear(relerr);

for(i = 0; i < 102; i++)

mpfi_clear(x[i]);

return 0;

}

６．まとめ

　以上，まとめますと

実数は無限桁の小数として表現される数である。
コンピューターのメモリは有限なので，無限桁の実数を格納することはできず，丸めて短い有限桁に収めなければならない。
丸めに伴う誤差が甚大な影響を及ぼす「悪条件問題」というものが存在する。
悪条件問題を解決するには桁を増やして対処するしかない。しかしこれは計算時間を要する解決策である。

となります。「高性能計算研究室」ではこのように長い桁の計算，多倍長計算を用いて様々な悪条件問題の解決に取り組んでいます。

2017年3月22日2017年3月22日

大学見学会「高性能計算研究室」

高性能計算とWebアプリケーション

静岡理工科大学　コンピュータシステム学科　幸谷こうや　智紀とものり

１．概要

　データの入出力をWeb経由で行うためのソフトウェアの開発，すなわち「Webプログラミング」の事例をいくつか紹介します。関数の計算式を入力することでグラフや数表を自動的に出力するアプリケーションや特別プログラムの製作物を体験して頂く予定です。

２．「高性能計算」研究室とは？

　現在の科学技術計算に欠かすことのできないのが，コンピュータ上で現象を再現するコンピュータシミュレーション技術です。このコンピュータシミュレーションの信頼性を高め，かつ，効率的に実行するためには，計算のための理論的枠組み，計算を効率よく行うためのアルゴリズム，そして堅牢かつ高性能なコンピュータが欠かせません。このようなコンピュータシミュレーションを下支えする学問分野を「高性能計算（High Performance Computing, HPC)」と呼びます。詳細はこちらの記事をご覧ください。

　本学の「高性能計算研究室」はその中でも特に，ソフトウェアを使ってシミュレーションの信頼性を効率的に上げるための研究を行っています。具体的には「多倍長数値計算(Multiple Precision Numerical Computation）」の研究が中心となります。現在のコンピュータが直接処理できるのは10進数で約16桁程度の有限桁小数ですが，どんな数も

$$ \frac{1}{3} \Longrightarrow 0.333333333333333 $$

と有限桁表現する必要があり，これを「丸め(rounding)」と呼びますが，この丸めによって生じる誤差が拡大することによって正しく計算されない「悪条件問題」というものがあり，シミュレーションの信頼性を損ねる一因となっています。それを防ぐ方法の一つとして，ソフトウェア的に桁数を増やして計算の精度を上げる，即ち，誤差を減らす手法が提案されており，それを「多倍長計算」と呼んでいます。

　多倍長計算の例として，単純なWebアプリケーションですが実行してみましょう。

Try MPFR!

　例えば，円周率$\pi = 3.1415\cdots $を10進約1000桁＝2進約3321bitsで計算するには

と入力すると

という結果が得られます。答えが出力されているテキストボックス左下のところをドラッグすると拡大することができます。

　本研究室では，この多倍長計算を効率的に行うために，たくさんのパソコンをネットワークで繋いだPCクラスタを作ってきました。その際培われたネットワークのノウハウを生かし，Webアプリケーションを開発するようになり，現在は

多倍長数値計算とその応用　＝高性能計算(HPC)
Webアプリケーションとの連携
数学ソフトウェアの教育利用　→　中間的

にも勤しんでいます。ここで紹介する「MPFRgraph」はHPCとしての多倍長数値計算を広報する目的で構築されたものです。

３．Webアプリの仕組み

　ここで，Webアプリケーション(Web application)，略して「Webアプリ」の構造を見ていくことにしましょう。先ほど紹介した「Try MPFR!」も一種のWebアプリです。これはもともと多倍長数値計算の土台を提供しているMPFRというライブラリ（沢山の機能を詰め込んだソフトウェアの素）を使って作った「mpfr_gexpr」という「ネイティブアプリ(native application)」をWeb経由で扱えるようにしたものです。

　ネイティブアプリは，パソコン（ハードウェア）を扱うための基本ソフトウェア（オペレーティングシステム, OS)の機能を利用して構築されたスタンドアローンタイプ(コンピュータ上で完結）のソフトウェアです。ハードウェア，OS，ネイティブアプリの関係性は，下記のようなレイヤー（層）で表現されます。

　これに対し，Webアプリは「Web （うぇぶ）」というインターネット上のサービスを土台としたアプリケーションで，「ブラウザ」と「サーバ」という，それぞれネットワーク上の独立したマシン上で協調して動作するものになっています。これらをレイヤー図で示すと下記のようになります。

　今のコンピュータは複数のプログラムを同時並行的に使用することができるので，ブラウザもサーバも同一マシンで動かすことは可能ですが，基本的にはインターネット上の独立したマシンでそれぞれ動くように作られています。ブラウザは今この記事をご覧になっている皆さんが直接目にしているソフトウェアです。代表的なブラウザとしては

があります。ちなみに，タブレットやスマートフォンではAndroid(Google)がChrome，iOSがSafariの利用が多いようです。Webサーバとしては

Apache
IIS (Microsoft)
NGiNX

の利用が多いようです（参照: NetCraft Web server servey）。

　Webアプリの開発は，このようなブラウザとサーバの連携を前提にしてソフトウェア構築をする必要があります。そのためには次のようなコンピュータ言語（ソフトウェアを記述するための開発言語）が使用されます。

HTML (HyperText Markup Language) ・・・静的Webページの構造を決める言語です。タグというもので括って文書構造を決定していきます。　→　サンプル（注意：ビデオと音声が流れます）
CSS (Cascading Style Sheet) ・・・静的Webページのデザイン（視覚効果）を決める言語です。
JavaScript ・・・主としてブラウザ上で動作することの多いインタプリタ言語（一行単位で逐次解釈・実行される）です。
PHP・・・主としてサーバ上で動作することの多いインタプリタ言語です。サーバ上のデータベースとの連携のためにも頻繁に利用されます。

　上記のうち，HTML, CSS, JavaScriptは一つの「HTMLファイル」（もしくは同等のもの）にまとめてサーバ側に保管され，２４時間いつでもインターネット上のクライアントマシンのブラウザからアクセスできるようになっています。サーバ側はHTMLファイルの内容を送り出しているだけです。それがブラウザに届いて初めてレンダリング（画面構築）が行われ，皆さんが現在目にしているような画面として表示されているわけです（下記参照）。

　従って，ブラウザからユーザが入力したデータを受け取るには，サーバ側にもその仕組みが必要になります。ブラウザにデータ入力のためのインターフェースを作るためにはHTMLでフォーム(form)を記述することで可能になります。例としてアンケートフォームを示します。

アンケートフォーム（HTMLのみ）

　しかし，フォームを作っただけではサーバにデータは届きません。そのための受け皿として，PHPで記述されたプログラム，PHPスクリプトを作り，フォームと連携させる必要があります。上記のアンケートフォームにPHPスクリプトを連携させたものを下記に示します。

　先に述べた，「Try MPFR!」もこれと同じ仕組みを使ってユーザからの情報「桁数(bits数）」と「数式(formula）」をサーバ側に伝えています。サーバ側ではそれをネイティブアプリである式パーサ「mpfr_gexpr」に与えて実行し，その結果をブラウザに送り込んでいます。

　この仕組みを拡張して関数の数表とグラフを同時に記述できるようにしたものが「MPFRgraph」です。

４．グラフ作成機能付き関数表作成アプリ「MPFRgraph」

　多倍長計算のためにはMPFRのような特殊なライブラリが必要ですので，Webアプリとしてネイティブアプリ多倍長精度式パーサmpfr_gexprを使って外部とのやりとりと連携させていました。しかしこの式パーサは数式に則って値の計算を行うことしかできません。多倍長計算のデモンストレーションとしては，様々な関数が使えることと，実際に関数 $y = f(x) $ の右辺式を与えて，独立変数$x$を区間 $[\alpha, \beta]$内で動かし，どのような値を取り，どのようなグラフになるのかを示すことが望ましいことになります。

　従って，区間$ [\alpha, \beta]$と分割数$n$，そして$f(x)$が$x$を含む数式として与えられていれば，分割された小区間幅を$ h = (\beta – \alpha) / n $として与え，分点$x_0, x_1, …, x_n$を

\[ x_0 = \alpha,\ x_1 = x_0 + h,\ …,\ x_{n-1} = x_0 + (n-1) h,\ x_n = \beta \]

として与えることで，関数の値$f(x_i)$ $ (i = 0, 1, …, n)$とセットで

\[ (x_0, f(x_0)), (x_1, f(x_1)), …, (x_n, f(x_n)) \]

をmpfr_gexprの機能を使うことだけで計算することができるようになります。こうして下記のように，グラフと数表を同時に出力できるようになりました。

　グラフについては，サーバ側でグラフィックスファイルととして表示させる方法と，ブラウザ側でHTML5のCanvas機能を使って表示させる方法があります。それぞれ作成してありますのでためしに様々な関数（多項式，三角関数，対数，指数等）を使って関数のグラフと数表を表示させてみてください。

数式入力＋PHPグラフ出力　（PHPlot使用）
数式入力＋Canvasグラフ出力　（jqPlot使用）

　使い方はどれも同じです。例えば

\[ y = \sin x^2 + \exp(\cos x),\ x \in [0, 10],\ \mbox{分割数: }\ n = 10 \]

の場合は下記のように入力します。

　関数式を入力する際は四則演算，関数は下記のように表記します。これらを組み合わせ，変数名は必ず「x」として正しい数式を入力して下さい。

加算：$ x + 3 $ → x + 3
減算：$ x – 3 $ → x – 3
乗算：$ 3x $ → 3 * x
除算：$ \frac{x}{3} $ → x / 3
平方根：$ \sqrt{x} $ → sqrt(x)
三角関数： $ \sin x,\ \cos x,\ \tan x $ → sin(x), cos(x), tan(x)
指数関数： $ e^x = \exp(x) $ → exp(x)
自然対数： $ \log(x) $ → log(x)

　これらのアプリケーションの一連の処理の流れを示したものが下記の図になります。

　詳細は煩雑になりますので省略しますが

関数のグラフ描画するための関数値の計算は，CanvasバージョンではJavaScriptが担っている
多倍長精度の計算はサーバ側でしか実行できないので，関数表は全てサーバ側のmpfr_gexprが計算したものを使用している

という役割分担がなされています。

５．最後に

　高性能計算研究室ではほかにもいくつかのWebアプリや，Webプログラミング開発のための教材も公開しています。

幸谷研究室＠静岡理工科大学

　また，「Webデザイン特別プログラム」では，過去の受講生の作品を公開しています。

Webデザイン特別プログラム

　24時間いつでも閲覧できますので，興味がありましたら色々見ていって下さい。特に優秀な作品は

Webデザイン特別プログラム優秀作品集

にまとめてあります。

参考URL

2016年9月9日2016年9月10日

模擬講義「コンピュータはどうやって計算しているのか？」（改題「数学と計算　～理論と実践～」）

概要

　数学（＝理論）と計算（＝実践）は全く性質の異なるもので，計算の効率性や正確性を追及すると数学とは別の方法論が必要になります。この講義では２次方程式や平均の計算を例に，「高性能な計算」をコンピューター上で行う手法を考えていきます。

１．数学＝理論，計算＝実践

　私の専門は「高性能計算(High Performance Computing)」です。読んで字のごとく，「高性能」な「計算法」を追求する学問の一分野ということになります。まだコンピューターが開発される前，20世紀までは人間があらゆる計算を行っていましたが，21世紀の現在，計算のうち単純なものはコンピューターが担うようになっています。従って，「高性能計算」も「コンピューターにとって」好都合なように計算方法を工夫するものになっています。この模擬講義では現在のコンピューターがどのような工夫のもとに「高性能計算」を実現しているかを具体的な例を挙げて解説します。

　さて，皆さんが思い浮かべる「計算」とはどのようなものでしょう？　「計算高い奴」と言うと，ずるがしこく立ち振る舞うという意味になりますが，ここで取り扱う計算とは「ごく単純な数の操作」とします。一番簡単な「数の操作」は，小学校で習う四則演算（加減乗除）ですね。ここで$a$, $b$は実数とします。

たし算（加算）・・・$a + b$
ひき算（減算）・・・$a – b$
かけ算（乗算）・・・$a \times b = ab$
わり算（除算）・・・$a / b$

　四則演算だけで実行できる計算といえば，平均値の計算があります。例えば$n$個の実数$c_1, c_2, …, c_n$の平均値$\mu$は，$n-1$の加算と1回の除算を組み合わせて
\[ \mu = \frac{1}{n}\sum^n_{i=1} c_i = \frac{c_1 + c_2 + \cdots + c_n}{n} \]
と計算することができます。

　もう少し複雑な計算になると，平方根（$\sqrt{a}$），三角関数（$\sin a, \cos a, \tan a $），指数関数（$\exp(x) = e^x$），対数関数（$\log a$）なども使う必要が出てきます。平方根を使う計算では，例えば2次方程式
\[ ax^2 + bx + c = 0 \]
の二つの解$x_1$, $x_2$は
\[ x_1 = \frac{-b + \sqrt{b^2 – 4ac}}{2a}, x_2 = \frac{-b – \sqrt{b^2 – 4ac}}{2a} \]
として計算することができます。
　こういう計算は今では普通に電卓やパソコンソフト(電卓ソフトや表計算)で簡単に実行することができます。

　ところで，こういう計算はどうしてこういう式として表わすことができるのでしょうか？　「平均値」の計算はどうして全ての数を足してその個数で割るのでしょうか？　2次方程式の解はどうしてこのような形になるのでしょうか？

　そのような疑問に答えるための学問体系が「数学(Mathematics)」です。

　例えば平均値。解釈の仕方はいろいろ考えられますが，例えばこう考えてはどうでしょうか？　いま，背丈の違う二人の人がいるとします。そのうち背の高い方から低い方へ背丈を分けて同じ身長にならしたとすれば，それがこの二人の「身長の平均値」を意味することになります。これが3人でも4人でも5人でも・・・$n$人でも，大きい人から小さい人へ背丈と分けて同じ身長に均した時の背丈が「$n$人の身長の平均値」であると言うことができます。

　2次方程式の解はどうでしょう？　教科書にはどうしてこういう解の公式が出てくるのか，ちゃんと解説が書いてあったはずですが皆さん覚えていますか？

　2次方程式の両辺を整理して，左側に$x$を含む1次式の2乗が来るようにすると
\[ \left(x + \frac{b}{2a}\right)^2 = \frac{b^2 – 4ac}{4a^2} \]
となりますので，
\[ x + \frac{b}{2a} = \pm \frac{\sqrt{b^2 – 4ac}}{2a} \]
から
\[ x = \frac{-b \pm \sqrt{b^2 – 4ac}}{2a} \]
となる訳です。

　こういう，「平均値の解釈」や「2次方程式の解の導出法」も「数学」の一部です。もちろん「平均値の計算」や「2次方程式の解の計算」も「数学」の一部ですが，これらの計算は「解釈」は「導出」の結果得られる「手続き」に過ぎません。計算の仕方が分かったとしても，どういう理由でそのような計算式になったのか，どういう理屈で計算式が導かれたのかは別に考える必要がある訳です。

　現在，コンピューターが担っている「計算」は数学的に導かれた解釈や導出法から得られる単純処理になっています。しかしこの先にもう一つ難しい問題が残っているのです。

　コンピューターはあらかじめ決められた手続きを教えてもらわなければ何も実行しません。手続きを記述した文書はプログラム（ソフトウェア）として人間が書いてやらなければなりません。平均値の計算や2次方程式の解の計算も，式に基づいてどういう手順で計算するのかを，四則演算や関数の組み合わせ方法まで指示しないと，効率的な計算ができないことがあります。このような手続きを「アルゴリズム(algorithm)」と呼びます。

　私が取り組んでいる「高性能計算」は，式を計算するためのアルゴリズムをどのように決定すれば効率が良くなるか，ということを考える学問分野なのです。

２．コンピューターの構造

　コンピュータは，全ての情報を電気信号として保持し，その信号を伝える回路（バス）と信号のやり取りを通じて何らかの処理を行う機器（デバイス）によって構築された金物（ハードウェア）と，ハードウェア内を伝わる電気信号の膨大な積み重ねによって構築された処理内容の蓄積，即ち軟体（ソフトウェア）によって構成されています。ハードウェアは堅い金属やプラスティックのカタマリですから，目で見ることができたり（可視物），触ってその存在を確認することができますが，ソフトウェアは単なる信号ですから，我々人間が感知するためにはデバイスを介して可視化する必要があります。処理を支持する信号を人間や外部の危機とやり取りするためのデバイスを入出力装置（I/O = Input & Output）と呼びます。こうした入出力装置を介しつつ，コンピューターに指令を与えるソフトウェアを作って初めてコンピューターを自在に操ることができるようになります。

　この原則は，先にあげた２次方程式の解の計算でも，平均値の計算でも変わりません。例えばC++(シープラスプラス)というコンピューター言語でこの計算を書くと次のようになります。これは計算部分のみの抜粋です（プログラム全体→quadratic_eq.cc)。

// print the given equation
cout << a << " * x^2 + (" << b << ") * x + (" << c << ") = 0 " << endl;

// Calculate solutions with analytic formulas
d = b * b - 4.0 * a * c;
sqrt_d = sqrt((complex<double>)d);
a2 = 2.0 * a;

x[0] = (- b + sqrt_d) / a2;
x[1] = (- b - sqrt_d) / a2;

// Print solutions
cout << "x[0] = " << setprecision(15) << x[0] << endl;
cout << "x[1] = " << setprecision(15) << x[1] << endl;

// print the given equation

cout << a << " * x^2 + (" << b << ") * x + (" << c << ") = 0 " << endl;

// Calculate solutions with analytic formulas

d = b * b - 4.0 * a * c;

sqrt_d = sqrt((complex<double>)d);

a2 = 2.0 * a;

x[0] = (- b + sqrt_d) / a2;

x[1] = (- b - sqrt_d) / a2;

// Print solutions

cout << "x[0] = " << setprecision(15) << x[0] << endl;

cout << "x[1] = " << setprecision(15) << x[1] << endl;

　解はx[0], x[1]という名前のメモリ上の領域に計算され格納されます。このプログラムはコンピュータが直接解釈できる形式に変換され（コンパイル），実行時に2次方程式の係数$a$, $b$, $c$を入力することで解$x_1$, $x_2$のをディスプレイに表示（出力）します。例えば
\[ x^2 + 3x + 2 = 0 \]
の解$-1, -2$を計算する場合は次のように実行します。
$ ./quadratic_eq 1 3 2 ← プログラム名，係数の入力 1 * x^2 + (3) * x + (2) = 0 ← 方程式の確認 x[0] = (-1,0)　←解1 x[1] = (-2,0)　←解2

　また，平均を計算するC++プログラム”average.cc“は次のようになります。

// print the given datum
for(i = 0; i < num_c; i++)
	cout << c[i] << " ";
cout << endl;

// Calculate average
average_mu = 0.0;
for(i = 0; i < num_c; i++)
	average_mu += c[i];

average_mu /= num_c;

// Print solutions
cout << "average = " << setprecision(15) << average_mu << endl;

// print the given datum

for(i = 0; i < num_c; i++)

cout << c[i] << " ";

cout << endl;

// Calculate average

average_mu = 0.0;

for(i = 0; i < num_c; i++)

average_mu += c[i];

average_mu /= num_c;

// Print solutions

cout << "average = " << setprecision(15) << average_mu << endl;

　このプログラムをコンパイルして，$(1 + 2 + \cdots + 10) / 10 = 5.5 $を計算させると，次のような結果が得られます。
$ ./average 10 1 2 3 4 5 6 7 8 9 10 ←プログラム名　データ数データ1 ・・・データ10 12 2 3 4 5 6 7 8 9 10 11 1 2 3 4 5 6 7 8 9 10 ←データの確認 average = 5.5　←平均値

４．浮動小数点数と桁落ち
　さて，例として挙げた2次方程式の解計算と平均値の計算，実はどちらも正確なものではありません。コンピュータ内部での実数は，有限桁の小数でしか表現できないため，どうしても真の値とのズレ（誤差）が発生してしまいます。

　例えば，円周率$\pi = 3.1415926535…$は，同じ数のパターンの繰り返しが現れない無理数として良く知られたものですが，正確に表現するためには無限の桁数が必要になります。しかしコンピュータのメモリは有限なので，どこかで打ち切って丸めなければ格納できません。例えば15桁で表現するためには
\[ \pi \approx 3.14159265358979 \]
としなければならず，
\[\pi – 3.14159265358979 = 0.0000000000000032384626433832795… \]
という誤差が生じます。これを丸め誤差と呼びます。例として挙げた二つのプログラムではどちらも約15桁程度の正確さを持っておりますので，16桁目以降は不正確ということになります。

　実用的には15桁も正確であればさほど問題発生しませんが，稀に「桁落ち」という現象が起き，結果が著しく不正確になることがあります。例えば2次方程式の例では
\[ x^2 – 12345x + 0.00001 = 0\]
を計算してみると
$ ./quadratic_eq 1 12345 0.00001 1 * x^2 + (12345) * x + (1e-05) = 0 x[0] = (-8.10359779279679e-10,0) x[1] = (-12344.9999999992,0)
という結果が得られますが，この場合$x_1 = -8.100445524504379… \times 10^{-10}$なので（正確な計算はこちらで実行できます），最初の3桁しか合っていません。誤差を含む，ほとんど同じ大きさの有限小数の加減算を行うとこのような「桁落ち」が起きます。これを防ぐには小数の桁数を増やすか，計算方法を変えるしかありません。

　ちなみに，桁数を約32桁に増やして計算してみると（プログラムはこちら）

$ ./quadratic_eq_dd 1 12345 0.00001 1.000000000000000000000000000000e+00 * x^2 + (1.234500000000000000000000000000e+04) * x + (1.000000000000000000000000000000e-05) = 0 x[0] = (-8.100445524504379240908904588516e-10,0.000000000000000000000000000000e+00) x[1] = (-1.234499999999918995544754956208e+04,0.000000000000000000000000000000e+00)

となり，桁落ち分をある程度フォローできていることが分かります。とはいえ，桁数を増やすより，もう少し賢いやり方はないのでしょうか？

　2次方程式の場合，次のように桁落ちしない正確な解$x_1’$を使って次の解$x_2’$を求めるという手順を取ることで，桁落ちを回避できることが知られています。
\begin{eqnarray*}
x_1′ &=& \frac{-b – \mbox{sign}(b)\sqrt{b^2 – 4ac}}{2a} \\
x_2′ &=& \frac{c}{2x_1′}
\end{eqnarray*}

　このような改良を加えたプログラムを”quadratic_eq2.cc“として実行してみます。
$ ./quadratic_eq2 1 12345 0.00001 1 * x^2 + (12345) * x + (1e-05) = 0 x[0] = (-12344.9999999992,0) x[1] = (-8.10044552450438e-10,-0)
　このように，正確な解を得ていることが分かります。

　平均値の計算でも桁落ちが発生します。絶対値の小さい数から足しこむことである程度，桁落ちを防ぐことはできますが，2次方程式のほどの抜本的効果はありませんので，小数の桁数を増やすことが一番の解決策となります。

５．マルチコアと並列化

　コンピュータに求められる能力は増えることはあれ，減ることはありません。人間の利便性を高めるために高性能化は常に求められています。
　コンピュータの性能を高める方法は，CPUの処理能力を上げることです。その一つの方法として，動作周波数(clock)を高めるという方法があります。実際に，20世紀最後の10年ほどで飛躍的に動作周波数は増加しました。

　しかし，21世紀に入ると打ち止め状態になっています。動作周波数を高めると電力を食うようになります。一般家庭でも使用できるようなパソコンが何千ワットも食うようでは使い物になりません。また，電力消費に伴う発熱も深刻で，安定的に動作させるための空冷ファンは大きくなるばかり。余計に電力消費が増えていきます。
　性能は上げなくてはならない，しかし，電力消費量の増加は極力抑えたい・・・という二律背反な要求に対し，CPUメーカーは電子回路の改良と共に，「マルチコア」という方向にかじを切りました。簡単に言うと，CPUの処理の中核を担う回路（コア）を複数（マルチ）持たせて，あたかも複数のパソコンを同時に使えるようにする技術です。もし6つのコアを持つCPUがあれば，6つの処理を同時にこなすことができ，短時間で終了することで多少電力消費は増えてもトータルでは少なくすることができるはずです。

　このように，複数の処理を同時並行にこなすことを「並列処理」と呼びます。対して，一つの処理だけを順次行う処理を「逐次処理」と呼びます。現在のパソコンやスマートフォン（スマホ）は2～4コアCPUが普通ですし，GPUのようにもっと大量のコアを備えてグラフィックスや計算を高速にこなす「メニー(many, 多数)コア」ハードウェアも登場しています。

　では，どうやったら複数のコアを当時に使うことができるのでしょうか？　並列処理を行うプログラム技法はいくつかありますが，一つのプログラムの中の処理を「スレッド」という単位に分割して並列に実行可能にさせる，というマルチスレッド化が良く用いられます。例えば，OpenMPという技法を使って2次方程式の解を同時に計算させるプログラムは次のようになります。

#pragma omp parallel sections
{
	// thread 0
	#pragma omp section
	{
		x[0] = (- b + sqrt_d) / a2; 
	}
	// thread 1
	#pragma omp section
	{
		x[1] = (- b - sqrt_d) / a2;
	}
}

#pragma omp parallel sections

{

// thread 0

#pragma omp section

{

x[0] = (- b + sqrt_d) / a2;

}

// thread 1

#pragma omp section

{

x[1] = (- b - sqrt_d) / a2;

}

$x_1, x_2$の計算は相互に依存していませんので，並行して計算させることができるわけです。逆に，桁落ちが起きないように改良してしまうと，片方の解が出るまでは次の解の計算に進めず，このような並列処理はできなくなります。

　しかし，2次方程式の計算は大変軽いものなので，並列化しても大したメリットはありません。桁落ちを防ぐ処理を行う方が実用的には有用と言えるでしょう。

６．並列処理の効果のほどは？

　では平均の計算はどうでしょう？　大量のデータ（数万，数億，・・・）の平均値の計算は普通に行われるものなので，並列処理するメリットはかなりありそうです。ではどのように行えばいいのでしょうか？　前述した

\[ \frac{1 + 2 + \cdots + 9 + 10}{10} = 5.5 \]

の計算を例に考えてみましょう。

　2コア（デュアルコア）のCPUでは，同時に２スレッドの処理が可能です。その際には，半分ずつ和の計算を行い，次にその部分和を足し合わせて1スレッドで合計値を求め，最後に10で割ります。

　4コア（クアッドコア）が使える場合はどうでしょう？　最初の部分和は4スレッドで計算し，次の部分和の計算は2スレッドで分担します。何故なら，最初の部分和の計算よりも計算量が増えると，ここで時間が余計にかかってしまうからです。

　最後は2個の和の計算を1スレッドで行って10で割ります。

　この場合，和の計算は最大10回で済みますから，5スレッド以上使える環境があっても意味がありません。

　こうして2スレッド，4スレッドそれぞれ並列実行できたとすると，計算時間は4スレッドの方が短くなることが期待されます。

　では更に大量のデータの平均の計算を実行しようとするとどうでしょうか？　データ数$n$に対して，同時に実行可能なスレッド数を$\alpha$ ($n > \alpha$)とすると最初の部分和の計算は$\lceil n / \alpha\rceil $回の和で済みます。ここで$\lceil x \rceil$は $x$を超えない最大の整数を意味します。

　次の段階で足すべきデータは$ n/2 $個になっています。これを繰り返すことで，最後は2個の和で済ますことができます。最終的に2個になったデータの和と$n$の除算は1スレッドで終了しますので，トータルでは和の計算の部分は最大
\[ \lceil \log_2 n \rceil \]
回の繰り返しで平均の計算が完了することになります。先にあげた10個のデータの場合は$\lceil\log_2 10\rceil$ $=\lceil 3.321928…\rceil = 4$となり，和の計算は4段階必要なことになります。

　このような集約を行う計算をリダクションと呼びます。これをOpenMPで記述すると下記のようになります。(プログラム全体: average_omp.cc)

// Calculate average
num_threads = 4;
if(num_threads > num_c)
	num_threads = num_c / 2;

omp_set_num_threads(num_threads);

cout << "#threads: " << omp_get_max_threads() << endl;
cout << "num_c   : " << num_c << endl;

#pragma omp parallel for reduction(+:average_mu)
for(i = 0; i < num_c; i++)
	average_mu += c[i];

average_mu /= num_c;

// Calculate average

num_threads = 4;

if(num_threads > num_c)

num_threads = num_c / 2;

omp_set_num_threads(num_threads);

cout << "#threads: " << omp_get_max_threads() << endl;

cout << "num_c : " << num_c << endl;

#pragma omp parallel for reduction(+:average_mu)

for(i = 0; i < num_c; i++)

average_mu += c[i];

average_mu /= num_c;

　使用するスレッド数はomp_set_num_threads関数で設定しています。データの数が少ない時はスレッド数を制限するようにしています。

　では実際に並列処理によってどの程度計算時間が軽減されるかを見ていきましょう。パソコンの環境は次の通りです。

CPU: Intel Core i7-3820 3.6GHz 4コア
メモリ: 32GB
OS: CentOS 6.5 x86_64
コンパイラ: Intel C++ 14.0.2
データ量: 16GB ($= 8\mbox{Bytes} \times 2 \times 1024^3$)

4コアありますので，最大4スレッドを並列に実行できることになります。そこで，スレッド数を1, 2, 4として実行してみました。”Time(s)”の所が処理時間(秒)になります。

$ export OMP_NUM_THREADS=4;./average_omp 1 #threads: 1 num_c : 2147483648 average = 1073738907.87731 Time(s) = 0.998027801513672 $ export OMP_NUM_THREADS=4;./average_omp 2 #threads: 2 num_c : 2147483648 average = 1073738907.87725 Time(s) = 0.568727016448975 $ export OMP_NUM_THREADS=4;./average_omp 4 #threads: 4 num_c : 2147483648 average = 1073738907.87729 Time(s) = 0.453949928283691

1スレッドでは約1秒でしたが，2スレッドで約0.57秒，4スレッドで0.45秒で終了しています。2スレッドで大体半減しますが，4スレッドにすると更に半分に，とはなりません。スレッド数が多いとその準備の手間で，理論的な見積りよりも処理時間が余計にかかります。

　ちなみに，太線部分のところがスレッド数によって異なっていることにお気づきでしょうか？　並列化してしまうと，特段指定しない限り，計算の順番は全く保障されないので，丸め誤差の差異が発生しているのです。これは桁落ちしていないのでせいぜい末尾の2桁程度で済んでます。

７．まとめ

　以上，駆け足でしたが，数学と計算の違い，そして計算そのものもきちんと考えることは意外と面倒なことが理解できたことでしょう。今まで述べてきたことのポイントを下記にまとめておきます。

理論式をそのまま逐次処理的に計算するにしても，桁落ちが起きないような工夫が可能であれば行う必要があります。
マルチコアCPUでは並列処理が可能なので，大量のデータの処理を行うときには威力を発揮できます。しかしそれは並列処理可能な部分に限られます。
並列処理を行って性能向上を図るためには，処理の流れ，即ちアルゴリズムを工夫する必要があります。

　現在の科学技術計算は，シミュレーションの理論化，計算処理，それぞれ分担して専門家が担うようになっています。ここで述べたのはごく簡単な例に基づいた計算処理の部分だけですが，数式をそのまま計算するというだけでは済まない面倒さがあることを覚えて頂ければ幸いです。

カテゴリー: エッセイ

情報学概論　第1回　レポート

CakePHPを学ぶ(1)

遠江総合高等学校：オープンデータとWebアプリケーション

PHP 5.6 security update終了

Update or perish

PCケースファンの交換

日経新聞電子版の普及

[模擬講義] コンピュータの計算は正しくない！

「コンピュータの計算は正しくない！」

静岡理工科大学　情報学部
コンピュータシステム学科
幸谷　智紀（こうや　とものり）

大学見学会「高性能計算研究室」

高性能計算とWebアプリケーション

模擬講義「コンピュータはどうやって計算しているのか？」（改題「数学と計算　～理論と実践～」）

2019年4月
日	月	火	水	木	金	土
« 3月
	1	2	3	4	5	6
7	8	9	10	11	12	13
14	15	16	17	18	19	20
21	22	23	24	25	26	27
28	29	30

共有:

共有:

共有:

共有:

共有:

共有:

共有:

「コンピュータの計算は正しくない！」

静岡理工科大学 情報学部 コンピュータシステム学科 幸谷 智紀（こうや とものり）

共有:

高性能計算とWebアプリケーション

共有:

共有:

静岡理工科大学　情報学部
コンピュータシステム学科
幸谷　智紀（こうや　とものり）