as_interger_ratio()では正確に元の分数には戻せない a = 0.1 を区間型に代入する際にうまく包含できないので、色々実験してみた。 何が問題なのか? import pint as pn x=pn.interval(0.1) と、2進数小数で表現不可能な数字なので下限と上限は format(x.inf,…

twosumアルゴリズムの研究 from pint import roundfloat as rf from pint import roundmode as rdm a=0.1 b=0.2 誤差が発生するパターン rf.rdadd(a, b, rdm.nearest) 0.30000000000000004 最近点丸め rf.rdadd(a, b, rdm.up) 0.30000000000000004 上方向丸…

成果 github.com ソースコードが長くなりそうだったので分離した。 本当は一つのファイルで行きたかったけど、管理や読み込みなどを考えたら分離せざるを得なくなった。 github.com github.com list型を継承したinterval型も入る便利な配列を作る。 に該当す…

成果 PEP8に対応した。 Githubにコードを公開した。 github.com 反省 PEP8対応のため、手動で途中まで直したが、autopep8コマンドで十分だった。(勉強にはなったけど) 今後の方針 list型を継承したinterval型も入る便利な配列を作る。 標準のmathを継承し、i…

成果 numpyのndarray型に対して、numpyのndarrayの加減乗っぽく計算する、仮想丸め方向付き計算クラスroundarrayを実装 反省 ndarrayの除算っぽい計算が実装できなかった(数値計算にとって必要なのかもよくわからない……) double型の行列に対してtwosum,twopr…

なんか他に書くこと溜まってるような気がするけど、便利そうだったので書いておく。 codenvy.com 数年前にニュースサイトで目にしたのでお気に入りに放り込んでおいてたが、全く使っていなかった。 お気に入りの整理をしている最中に目についたので、アクセ…

成果 interval型の配列をplotする plotinfsup,plotlengthを作成した。 反省 しかし、実装がとても汚いので、公開には遠い 理想はplotinfsup(x,y)でplotできるようにすること とある場所で進捗を発表したら、フィードバックがあったので。 今後の方針 行列の…

interval型の四則演算は加減乗まで完成。 丸め方向問題は柏木雅英先生の無誤差変換を用いた方向付き丸めのエミュレートを利用することにした。クラス roundfloat型の中に丸め方向を指定した加算乗除の計算関数を実装。 str( )などの関数は未実装 肝心要のplo…

今期はMATLAB(Octave)で作業する時間が多かったので、すぐプロットしたりプレビューできる環境にハマってしまった。そこでpythonを用いて、似たようなものが作れないかとやってみる。 目指すべき特徴 python3.4.3をサポート RaspberryPi上で動けば完璧だが、…