kashiの日記

• ライセンスを明記し、MITライセンスとした。
• interval.hppのproper_subsetを修正し、proper_subset([0,1],[0,2])が真になるようにした。
• beta.hppを追加。
• eig.hppを追加(高安 亮紀氏の提供による)。

(2015年10月6日 追記)

• kvライブラリの区間演算の部分を切り出して、簡単な区間演算ライブラリを作ってみました。
• テンプレートを排除しました。従って内部型はdoubleのみですが、プログラムが読みやすくなっています。
• テンプレートを排除したおかげでboostに依存しません。単体で使えます。
• 必要なファイルはinterval.hpp一つです。
• 内部型のカスタマイズが出来ないこと以外は、kvライブラリに含まれるinterval.hppと機能は同等です。数学関数も多数あります。
• 丸めの変更は単純にfesetroundを使っています。
• 文字列との方向丸め付き相互変換機能もあります。従って、文字列からその文字列の表す数値を含むように区間を初期化できるし、表示も真の値を含むように外側に丸めるようになっています。
• volatileを使った最適化対策はきちんとやってます。

• 原稿
• スライド

function testmm(n)
	disp('testing multiplication...')
	flag = 1;
	for k=1:n
		a=zeros(n);
		a(:,k)=ones(n,1)*1/3;
		b=zeros(n);
		b(k,:)=ones(1,n)*1/3;
		c = (rad(intval(a) * b) ~= 0) + 0.0;
		if prod(prod(c)) == 0
			disp('multiplication error!')
			flag = 0;
			break
		end
	end
	if flag == 1
		disp('multiplication OK!')
	end

	disp('testing addition...')
	flag = 1;
	for k=2:n
		a=zeros(n);
		a(:,1)=ones(n,1);
		a(:,k)=ones(n,1)*2^(-27);
		b=zeros(n);
		b(1,:)=ones(1,n);
		b(k,:)=ones(1,n)*2^(-27);
		c = (rad(intval(a) * b) ~= 0) + 0.0;
		if prod(prod(c)) == 0
			disp('addition error!')
			flag = 0;
			break
		end
	end
	if flag == 1
		disp('addition OK!')
	end
end

# CONSISTENT_FPCSR = 1

• sudo apt-get install devscripts
• 適当な作業directoryを作ってそこに移動
• sudo apt-get source openblas でソースコードを取得。
• sudo apt-get build-dep openblas でmakeするのに必要なものをインストール。
• cd openblas-0.2.8/ でソースdirectoryに入る。
• Makefile.ruleを編集し、コメントアウトされている
  

  CONSISTENT_FPCSR = 1
  

  を有効化する。
• dch -i で、version番号を例えば0.2.8-6ubuntu2のように進め、changelogを適当に書く。
• dpkg-buildpackage -b -uc -us でmake。
• これで必要なdebファイルが生成されたので、cd .. で一つ上に上がり、
• sudo dpkg -i libopenblas-base_0.2.8-6ubuntu2_amd64.deb
• sudo dpkg -i libopenblas-dev_0.2.8-6ubuntu2_amd64.deb

• libopenblas-ubuntu-rounding.zip