久しぶりに

書いたら、インターフェースが大分変わってる。 MacBook Air 11インチ欲しい! 最近は汎用な感じにUNIX設定ファイルを書きなおしたり、 MUA を mutt に変えたり、RSS Reader を newsbeuter に変えたり、 tmux を使ってみたり、端末大好きっ子です。 MacBook …

配列

d:id:konnyakmannan:20091212 をみて。 彼が FORTRAN の配列をどこで便利に感じているのかはわからないのだけれど、 C で任意次元配列をサブルーチンで扱うために一次元配列を利用するのを思いだした。 やり方は Row- and column-major order - Wikipedia や…

現実逃避として、

OpenBSD インストールちゅう。

めもりのせたけど、

起動したらECC えらーがでてる。 ECC えらーってそんなにでるもんなのかしら。 だとすると、Non-ECC ってこわい。さて、メモリをいっぱいに使った計算でどのくらいかかるか知りたいけれど、 そのまえにいろいろなおさないといけない。

昨日は

セミナー後つかれて大学で寝てしまった。

なんか結果おかしい

時間かかるんで後まわしにしてた、伏見表示してみっか。 というか projection した方が時間的に多少ましか。

intel compiler

だと、jansson が一発 make できなかった。 icc が $ CC=icc ./configure --prefix=$HOME $ make ... Making all in src /bin/sh ../libtool --tag=CC --mode=compile icc -DHAVE_CONFIG_H -I. -I.. -Wall -Wextra -Werror -g -O2 -MT hashtable.lo -MD -MP …

並列化かんりょ

メモリも無いので、とりあえず格子数 16777216、100000 ステップ流しとく。来週はセミナ発表準備につかいましょう。 結果は厳しいので、文献紹介かな。 文献はやっただけっぽいから難しいことはないはず。

xmalloc, xfree

みたいなのを用意するか、どうかな。

研究室の

Wiki に書いてたんだけども、 物理学会の予稿は http://ci.nii.ac.jp/vol_issue/nels/AA11439205_ja.html から見られる。 先輩や先生から見せてもらうのってもいいんですけどね。しかし、Web から一般公開されているわけなので、下手なこと書けないわけです…

いまだに、

プログラミング言語の機能なんてものがもてはやされるのは理解できない。 だいたい、原理的にできることは同じなんだから新機能なんかよりも犯しやすいミスを防ぐ、読みやすいコードを書くといったような制約を課す方がよっぽどまともに思える。そういう機能…

紆余、曲折を経て、

たぶん、これでほとんど仕上がったはず。久しぶりに、C 言語で書いたもんだからいろいろ忘れていたりした。 例えば、 void do_something(struct spam* s, int N){ for(int i=0; i<N; ++i) (s+i)->member = i; } みたいなのにひっかかった。なんという基本的なところだろう。</n;>…

つかまらない。

ので、復習してして考えてみる。やっぱり間違ってますね。 結局自分で考えたやつでやったやつでよかったというのは悲しい。 ので、御利益ありそうな方法をさぐる。ちゃんと御利益はありそうですね。計算過程が若干減らせるし、 前のアルゴリズムの微妙な点を…

プログラムを

つめている最中色々考えたが、相互作用表示のアルゴリズムが誤っているような。 初等的ミスかなぁ。 とりあえず、明日議論できればいいけど。

プログラムの

ロジックがすっきりするなぁ。 研究のプログラム、メインのソルバーはほぼ書けてしまった。C++ をやめる理由。 昨日が多すぎる。これであれが実現できるんじゃないか?とか思うと大抵泥沼で バッドノウハウを積み上がっていく。 結果として元々の実現したかっ…

version の

Mac

古い、Xcode の入手法。 https://connect.apple.com/ にアクセス。 ログイン、Downloads、右メニュー Developer Tools

もう、

C++ は使わないと決めたので、json を読みこむのに Jansson — C library for working with JSON data を使ってみました*1。 大分すっきりしたコードになりました。結局使う機能がしぼられている方が、集中しやすいし、ロジックを簡単なものにしやすい、 C++…

おー、

Mac mini に Snow Leopard Server 版が! Mac OS X Server 利用する選択肢としては、 Mac mini \104900〜 Mac Pro \278800+\53800=\332600〜 Xserve \338800〜 なので、価格設定としては非常にお得なモデルですね。 しかし、Snow Leopard Sever が unlimited…

つかれた、

夕食後、指導教官となぜか超流動のフォノンーロトン・スペクトルの話。 なかなか勉強になります。 岩波講座 現代物理学の基礎 物性 と Abrikosov Gor'kov Dzyaloshinskii は読んでおくべしとのこと。 こういう意外な引き出しもあって、近似法も何もかも空で…

それと Mac OS X で BLAS の関数、例えば、dgemm を呼び出すのにちょっとてこずったので、やり方を書いておく。 まず、ライブラリ含まれているかどうかを調べる。/usr/lib/libatlas.dylib に含まれているはず。 $ otool -v -t /usr/lib/libatlas.dylib|grep …

とりゃあっと

BLAS, LAPACK の subroutine がややこしいので wrapper を書いた。 しかし、BLAS や LAPACK は破壊的な subroutine なので、 operator だと内部でコピーするので余計な生成がおきてしまう。 *= += -= とかは、メンバに一時的な配列を加えておいて、copy して…

emacs を

やめてしまったのは設定を追いかける若さがなくなった。ということかもしれない。 バージョン、Carbon、Cocoa、X いろんなものがありすぎるし。Python 使っているのもその傾向かもしれない。

遊びで

実験装置作成。予想以上に間に合わせでできてしまって、工作としてはいまいち。 ただ、いくつか改良しないとまずい点があるが、どうするか。。。 とりあえず、手伝ってくれた卒研生の人は満足したようでよかった。 工作的にはスイッチ等のインターフェースや…

数値計算のパラメーターを

json にして読ませるようにした。例えば、 { "outputpath": "./data/", "eps": {"vector": {"expr": "1.0/sqrt(n)", "n": [1,10]}}, } みたいに書いておいて、python で eval するスクリプトを書いて、 { "outputpath": "./data/", "eps": [1.0, 0.707106781…

epwing 辞書の変換を

試みるもうまくいかず。以下ログ。ebd2html 実験バージョン を試してみるも、 途中で終了しちゃう。エラーメッセージくらいだしてほしい。 問題は、前処理の ebdump かもしれない。 どうしようもないので他の手段を考える、DDWin というアプリケーションでテ…

めも

Mac OS X の screen 問題関連。 ****** Returning nil _server ********** とか Terminal ででる問題。 2008-07-19T10:11:46+09:00 2008-07-25T23:19:42+09:00 screenでpbpasteできないのをどうにか - はこべにっき ♨ とりあえず patch をあてずに、最新版を…

なんとなく、

体調不良。Mac OS X Snow Leopard Up-to-Dateプログラム macOS Mojave - Apple(日本) そろそろ送ろう。 980円で update は安いですね。 行列の条件数とか、QR 分解、特異値分解とか。

数値計算の

パラメーター、計算条件、出力ディレクトリくらいをファイルにだしておこうと思いたって、 プログラム書き。 配列まであって、パーサーをよそからとってこれて、可読性がある ということで、フォーマットは json がいいのかな。

openmpで

Mac Pro をこきつかう。 1 ジョブで CPU 800%!!! #pragma omp parallel for おそるべし。だって時間かかるんだもの。 人がいるときはやってはいけない。しかし、これまでは new で呼び出されるコンストラクタが引数無しでないといけなかったので、global 変…

今週は

アルゴリズムを考えて、今日はプログラムのメイン部分を書いていた。 アルゴリズム練りをきっちりやったのでコーディングは速く進む。 最適化まで考えながらやってしまうけど、 並列化と一緒に後回しでやるべきだな。しかし、プログラムのテストが面倒だ。プ…