ようこそゲストさん

小さなそして広がっていく世界

2010/06/17(木) はじめての電子状態計算

2010/06/17 22:34 未分類

先輩,プログラムが動きません!

専門(より)の話.後輩にプログラムが動かないので助けて欲しいといわれました.
で,見てみたんですけど,どうやら,『はじめての電子状態計算』っていう本に付属しているプログラムらしい.
動作環境を見てみたら,「Windows 3.1,Windows95,Windows NT3.5 Windows NT4.0,MS-DOSも可」とある.
あれ,これだけ?
XPやVistaはどこでしょうか……?

さらに,記述言語はFORTRAN 77らしいです.うーむ古いプログラムですね.
これでは,コンパイル済みの実行ファイルでは動かないかも.
幸運なことに,ソースコードはついてました.
なので,「とりあえず,自分も試してみるから明日報告するよ」といって,後輩が帰ってからコンパイルを試みることに.

まず,バイナリをVistaで実行

本に書いてある通りに作業をして,一応それらしい動作しましたが,以下のようになりました…….
c:\DVXA>SETDVXA.BAT

c:\DVXA>cd CALC\CO

c:\DVXA\CALC\CO>makef05
DOS/4GW Protected Mode Run-time  Version 1.97
Copyright (c) Rational Systems, Inc. 1990-1994
DOS/4GW error (2001): exception 0Dh (general protection fault) at 21F:00006F12
TSF32: prev_tsf32 5284
SS       22F DS       22F ES       27F FS         0 GS         0
EAX        0 EBX      22F ECX        0 EDX        0
ESI      120 EDI     588C EBP     4F52 ESP     4448
CS:IP  21F:00006F12 ID 0D COD     7E40 FLG    13246
CS=  21F, USE16, byte granular, limit     914F, base  310F2D0, acc FB
SS=  22F, USE16, byte granular, limit     489F, base  3118420, acc F3
DS=  22F, USE16, byte granular, limit     489F, base  3118420, acc F3
ES=  27F, USE16, byte granular, limit     7CFF, base  31414C0, acc F3
FS=    0, USE16, byte granular, limit        0, base  31414C0, acc  0
GS=    0, USE16, byte granular, limit        0, base  31414C0, acc  0
CR0: unavailable

c:\DVXA\CALC\CO>dvscat n \dvxa\calc\co
DOS/4GW Protected Mode Run-time  Version 1.97
Copyright (c) Rational Systems, Inc. 1990-1994
DOS/4GW error (2001): exception 0Dh (general protection fault) at 21F:00006F12
TSF32: prev_tsf32 5288
SS       22F DS       22F ES       27F FS         0 GS         0
EAX        0 EBX      22F ECX        0 EDX        0
ESI      120 EDI     588E EBP     4F52 ESP     4448
CS:IP  21F:00006F12 ID 0D COD     7E40 FLG    13246
CS=  21F, USE16, byte granular, limit     914F, base  310F2D0, acc FB
SS=  22F, USE16, byte granular, limit     489F, base  3118420, acc F3
DS=  22F, USE16, byte granular, limit     489F, base  3118420, acc F3
ES=  27F, USE16, byte granular, limit     7CFF, base  31414C0, acc F3
FS=    0, USE16, byte granular, limit        0, base  31414C0, acc  0
GS=    0, USE16, byte granular, limit        0, base  31414C0, acc  0
CR0: unavailable
Could Not Find c:\DVXA\CALC\CO\f15

c:\DVXA\CALC\CO>
F05とかっていうファイルが出来るはずなのに出来ないし,EXにあったF05ファイルをコピーして実行してもエラー出たし…….
やっぱり完全には動かないってことですか?それとも設定ミスだろうか(書いてあることは全部やったんですが).
で,コンパイルしてみるかってことに.

これがFORTRANか!

ソースコードはGreekだった

ソースコードを見てみたんですが,わけわかんなかったです,はい…….
というか読む気が起きなかった……….
このコードにはコメント文のコの字もなかったので,理解するにはちゃんと勉強しないと駄目だということがわかりました.

gfortranでコンパイル

コンパイルの仕方については,本にはわずかに書いてありました.
なになに,バッチファイルがあるのでそれを実行しろ,とあります.
コンパイラは……Watcom FORTRAN?聞いたことがない名前です.
手持ちのFORTRANコンパイラはgfortranだったのでとりあえずそれで実行しました.

結果は……

出来たものと出来ないものがありました.
例えば,上のF05というファイルを作成するプログラムは成功し,実行して中身のあるファイルをはきました.
c:\DVXA\CALC\CO>MAKEF05
 f25 for nonsym calculation was generated

c:\DVXA\CALC\CO>
というわけで,実行ファイルがVistaやXPに対応していなかったらしいです.
ただ,全体的に失敗したものが多かったので,コンパイラを変える必要がありそうです.
g77(FORTRAN 77のGNUコンパイラ)をどこからともなく探してきて*1,コンパイルしてみるとか.
けど,今日は時間切れ.後輩に方針は示せそうなので後は頑張ってもらいましょう.

それにしても読みにくい本だ 加えて問題点

何となく読みづらいです.操作の手順も本文にコマンドラインの解説があるから見づらい見づらい.
また,アルゴリズムとかにはほとんど触れてないのも,理解しづらいですね.
それに版が古くて,最近のプログラムが公開されていないところを見ると,あまり開発とかは活発ではないようです.
プログラムも,あんまり大規模なものでもないし,高度なことをしようとすると壁に当たってしまうのではないかと思います.

何よりも問題は,ソースコードの改変とかが禁止されているということ.
これによって,ユーザ自身による改良という道がたたれています.
GPLとかにしておけばもう少し広まったかもしれないですね.

というわけでフリーのソフトを少し探してみる

結構いっぱいありますね.日本語ドキュメントのものでは,
文部科学省系の第一原理擬ポテンシャルバンド計算ソフトウエア「PAHSE」
Debianにパッケージングされているものでは,ABINITというものが見つかりました.
両者とも,ドキュメントが豊富で扱いやすそうです.
特に後者はGPLライセンスで再利用は(研究などに使う分には)容易です.


動画のエンコードなんかに使うより,こういうのに計算機を使った方が有用ですよね…….

Windows XPでは動作しました(2010年8月3日追記)

意外とこのページを見ている方がいるようでしたので,追記をしてみます.
『はじめての電子状態計算』っていう本に付属しているプログラムですが,
どうやら,XPであれば動作するようです.先輩と後輩,私が個別に実行して確認しました.
ただ,もう一度Vistaでやってみましたが,やはり動かず.NT6.0以上では駄目っぽそうです.

1: 青城光 2010年06月18日(金) 午後10時15分

いやあ、プログラムっていうのはその人のやりたい事を実現するひとつの道具に過ぎません。道具は好きなように使えばいいのですよ、誤った使い方をしなければ。

2: Hirayama Hirosugu 2010年06月18日(金) 午後10時22分

なるほど,正論ですね.ただ,エンコードが誤った使い方のような気もしてしまうのです…….


名前:  非公開コメント   

  • TB-URL(確認後に公開)  https://ac206223.ppp.asahi-net.or.jp/adiary/memo/adiary.cgi/tdiary2/037/tb/