直列加算と並列加算でdoubleの足し算の結果が変わる話 - きしだのHatena
Javaに限った話ではないのだけど、Javaで並列加算が気軽にできるようになったので、気に留めておいたほうがいい話。 まず、次のようなコードを動かしてみます。 public static void main(String[] args){ double[] data = { 1.234E80, -1.234E80, 2, 3}; System.out.println(Arrays.stream(data).sum()); System.out.println(Arrays.stream(data).parallel().sum()); } 1.234×10^80と-1.234×10^80という、桁が大きくて符号の違う数を並べて、そのあとに2と3という1桁の数値を置いています。 これらを加算すると、1.234×10^80と-1.234×10^80は符号が違うだけなので、当然結果は0になります
そうなのか。 レガシーでよくあんな複雑なもん作れるなあ。 自動化どころか..
そうなのか。 レガシーでよくあんな複雑なもん作れるなあ。 自動化どころか、テストそのものすら絶望的に難しい気がするがどうしてるの？ 特定のアイテムの組み合わせを持った状態で特定の場所で特定の行動を取るとバグるとか、確かめようがない気がする…。 ちなみに数値計算系もテスト無理。ゲームだと流体計算とかシェーディングとかか。あれが正しい計算結果になってるかどうかとか、特定の境界条件で変な解が出ないかとか、無理だろ。 解析解が得られるケースで試すくらいしかないと思ってるけど、他に何かあるなら教えて欲しい。
Numerical integration trick
Suppose you want to evaluate the following integral: We’d like to do a change of variables to make the range of integration finite, and we’d like the transformed integral to be easy to evaluate numerically. The change of variables t = 1/x2 transforms the range of integration from (30, ∞) to (0, 1/900). I’ll explain where this choice came from and why it results in a nice function to integrate. The
zunda framework - neuron simulator
zunda frameworkは，コンダクタンスベースの神経細胞モデルを構築し4次のルンゲ・クッタ法またはオイラー法でシミュレーションできるようにするフレームワークです。このフレームワークはC#で書かれていて.NET Frameworkで実行されます。 モデルの構築にはこのフレームワークで定義されているクラスを継承することで新しい細胞を定義し，いくつかのイオンチャネルやシナプスを追加します。それぞれのチャネルの特性やシナプスの特性はプログラミング上で数式などで定義します。 ライセンス ライセンスは，zlib/libpngです。日本語訳。 ダウンロード zunda framework プロジェクト本体（サンプル付き） zundaframework.zip チュートリアル Hodgkin Huxley モデル リファレンス zunda framework API Referenceをご覧くださ
NArray and scientific computing with Ruby - RubyKaigi2010
NArray is an n-dimensional numerical array library developed for use in science fields, but it is usable for handling binary data efficiently. I talk about NArray for the first time after 10 years from the first release. I will also present topics including the next version of NArray, and distributed scientific computing using Ruby for multi-core era.
python - 機械学習の「朱鷺の杜Wiki」
python一般† python.org：公式サイト Wiki Package Index iOS用Python環境 Computable Python for iOS Python Math PyPad Pythonista 英語資料 Python Course：Python のいろいろなトピックについての講義を集めたサイト Python Quick Reference @ Richard Gruet (旧版 日本語訳) Google's Python Class OLamp.com -- Python Dev Center：O'Reillyのpython関連ニュース Python Miro Community：チュートリアル講演ビデオのポータル Wikipedia:Python_(programming_language) Python tools that everyone shou
NArrayの特徴
NArrayは、大規模な多次元数値配列の計算を、 簡単かつ高速に実現するRubyのクラスです。 同様な数値配列には、 Python/NumPy, Perl/PDL, Yorick などがあります。 NArrayはそのRuby版で、特徴としては同じなんですが、 実装は独自です。 数値およびオブジェクトの多次元配列 NArray は、多次元配列のクラスです。 配列の要素の型として、 8,16,32 bit 整数 単精度/倍精度の実数/複素数、 Rubyオブジェクト をサポートしています。 1要素の取出しは、a[ 12, 34 ] のように各次元のインデックスを指定します。 部分配列は、a[ 10..15, 5..-1 ] のように範囲指定で取り出せます。 演算の記述 +,-,*,/,%,** の演算は、要素ごとに適用します。 たとえば、2つの NArray クラスのインスタンス a,b が
Ruby/GSL
GFD Dennou Club / Dennou Ruby / Tutorials Title: Ruby/GSL 編集者: 塚原大輔, 中野満寿男 最終更新: 2006/03/08 (中野満寿男) 新規作成: 2005/02/12 (塚原大輔) Contents Ruby/GSL とは バイナリパッケージでインストール ソースコードからインストール plotutils のインストール GSL のインストール rb-gsl のインストール 使用例 EOF 解析 1 Ruby/GSL とは GSL(GNU Scientific Library) はその名の通り科学技術計算ライブラリで さまざまな数値計算法の関数がたくさん集められています．<URL:http://www.gnu.org/software/gsl/> ANSI C で記述されていて，C や C++ から呼び出せます． Ruby/
Numerical Ruby NArray
多次元 数値 配列クラスです。 1,2,4 byte 整数、単/倍 精度 実数/複素数、 およびRubyオブジェクトを要素に持つことができます。 これにより Ruby でも 大量の数値を扱う計算が、簡単かつ高速にできるようになります。 NArray 特徴 NArray メソッド一覧 NMatrix & NVector - LU分解により線形方程式を解く NArray デモ: マンデルブロ 画像のリサイズ 画像のスムージング(FFTW使用) ライフゲーム ベンチマーク Future rewrite plan ソース: narray-0.5.9p5.tar.gz (2008-06-10) mswin32バイナリへのリンク: narray-0.5.9-i386-mswin32-1.8.zip Known bugs NArray の拡張 FFTW3 FFTW (高速 DFT (Discrete
計算物理のためのC/C++言語入門
Last revised on May 25, 2000 UNIX & XWindow 環境における C/C++プログラミング言語の基礎とその計算物理への応用を解説します。 このホームページの古い版の偽ミラーサイトが存在するのでご注意ください。 本物のサイトは http://www-cms.phys.s.u-tokyo.ac.jp/~naoki/CIPINTRO/ です。 前書き 推奨図書一覧 C言語とC++言語の基礎文法の初心者向け解説 (プログラム初心者を対象としたC/C++言語の入門です) 第1章 C言語とC++言語の役割 第2章 簡単な計算 第3章 ループで繰り返し計算 第4章 配列で大量データ処理 第5章 関数で計算の分担 第6章 文字列の操作 第7章 ファイルの操作 第8章 グラフィック(XWindow System利用者向け) 高校物理の計算例 (高校生を対象として初歩的な物
ローレンツアトラクタ (The Lorenz Attractor)
— ローレンツアトラクタを描画する, あるいはそれに類似することをするときのために. — ソースコード l1.c /* ローレンツ方程式のカオスのストレンジアトラクタを描く. */ /* ** ローレンツ方程式の説明についてはlorenz.hのコメント文を参照せよ. ** アトラクタとは時間とともに体積が縮小する散逸力学系において ** 過渡状態経過後の定常状態のことである. ** カオスのアトラクタを特に, ストレンジアトラクタト呼ぶ. ** 常微分方程式の数値計算にはルンゲクッタ法を用いた. ** ルンゲクッタ法の説明についてはrk.hのコメント文を参照せよ. ** 描画にはgnuplotを用いた. ** このプログラムはlinuxで動作確認している. */ #include <stdio.h> #include <stdlib.h> #include "rk.h" #include
数値計算 - MacWiki
数値計算ライブラリ[編集] ATLAS[編集] Automatically Tuned Linear Algebra Software (ATLAS) http://math-atlas.sourceforge.net/ https://sourceforge.net/project/showfiles.php?group_id=23725 Absoft Fortran には、ATLAS の 64bit 版が付属しているそうです。 http://www.hulinks.co.jp/software/pf_mac/section01.html Bayes++[編集] Bayesian Filtering Classes for C++ Kalman filter や各種の非線型フィルタの数値計算ライブラリ http://bayesclasses.sourceforge.net/Bayes++
1
リリース、障害情報などのサービスのお知らせ
最新の人気エントリーの配信
j次のブックマーク
k前のブックマーク
lあとで読む
eコメント一覧を開く
oページを開く
[iPhone][iPad][無料] Linpack | iPad iPhone Wire
iPad で Linpack - 最適化問題に対する超高速＆安定計算
Ruby科学データ処理ツールの開発 NArrayとPwrake
Num0610b.dvi
Apple Developer
巫女 GNYO/Linux - Seesaa Wiki（ウィキ）