【SQL】NULLの上手な扱い方! IS NULL演算子とISNULL関数について
Fortran プログラミングの基礎知識 Fortran プログラミングの基礎知識 初版 1999. 31 最終改訂 2006. 20 今では Fortran はメジャーなプログラミング言語とは 言えませんが、数値計算の分野では、まだまだ現役の言語です。 複素数と整合配列を手軽に扱える言語は fortran しかありません。 Fortran90 も普及しつつあるようですが、過去の資産の活用、 FreeBSD や Linux で使えるコンパイラとなると、 fortran77 がまだまだ使われています。 今では書店に行ってもフォートランの本はほとんど置いておらず、 「本には載っていないけれども、実は基礎的なことがら」を知ることは 容易ではありません。 そこで、このページでは Fortran77 で プログラミングするにあたっての基礎的で実戦的な知識を 紹介致します。 以下では Fortran と表記した場合、Fortran77 を意味します。 暗黙の型宣言を使わない方法 Fortran はデフォルトで変数の型宣言が不要です。 i,j,k,l,m,n で 始まる変数は整数型、それ以外は実数型という暗黙の了解があります。 暗黙の型宣言を使いたくない場合は、サブルーチンの先頭で、 implicit none と宣言すると、変数の宣言が必ず必要になります。 暗黙の型宣言を使用すると、変数の宣言を逐一行う 煩わしさがなくなりますが、以下のようなデメリットが あります。 変数の名前をタイプミスすると、新しい変数と みなされてしまいます。 ftnchek などのツールを使うと この類のミスは発見してくれますが、コンパイラが発見して くれる方が楽です。 また、後述する include 文を使う場合、 ftnchek の表示は非常に見づらくなります。 変数名として自由な名前を付けることが できません。 たとえば最大値を max という変数で表したいとき、 rmax や dmax などのような変数名を付ける必要があります。 変数のタイプミスを防ぐというメリットは非常に大きいので、 implicit none は使うべきであると思います。 ちょっとした使い捨てプログラムを 書くとき以外は暗黙の型宣言は使ってはいけません。 この場合、約 15 桁の有効桁数が得られます。 0 4 byte 実数の精度で acos を求める として下さい。 ただし、ファイルやキーボードから読み込むときは、 単精度の表記をしても、倍精度の実数として認識されるようです。 1d0 と 書くのと同じ結果になります。 ただし、私の研究室で使っているグラフ描画ライブラリの 実数型引数は 4 byte real なので当研究室の人は注意して下さい。 暗黙の型宣言を使う場合、implicit 文の次のような使い方は 絶対に避けた方が良いでしょう。 変数名とサブルーチン名 Fortran77 の規格では変数名やサブルーチン名は 6 文字 以内となっていますが、 現在のほとんどの処理系では 6 文字以上の変数名やサブルーチン名 を許しています。 f2c, g77 をはじめとして 私が使ったことのある全てのコンパイラでこの機能はサポート されています。 定数は別ファイルで定義し、それをインクルードする これを行わない場合、あらゆるファイルに定数が散りばめられる ことになり、プログラムの保守が大変になります。 定数はインクルード ファイルで定義して、そのファイルをインクルードすべきです。 これを実現する方法として、次の 2 つの方法があります。 cpp C のプリプロセサ を使い、 define , include を使う Sun の純正 fortran や昔の Sony の News 付属の fortran 、g77 は 拡張子を. F とすると cpp を通した後、コンパイルします。 従って、 cpp の命令である define や include や ifdef が使えます。 FreeBSD 3. x の f77 コマンドなど. F のファイルを 扱えないコンパイラもあります。 これを解決するために、 私が作った という perl スクリプトがあります。 cpp を通した後、f77 に渡します。 また、多くの Linux のディストリビューションに 含まれている fort77 という perl スクリプトは拡張子が. F の 場合は cpp を通してから f2c を呼び出して、gcc でコンパイルするという perl スクリプトです。 parameter 文と include 文を使う 定数はインクルードファイル中に parameter 文で記述し、 include 文を使ってインクルードファイルを取り込みます。 プログラムは小文字で書く 小文字で書いた方が見やすいので小文字を使いましょう。 また、 大文字はマクロ define で定義する や parameter 文で定義する 定数にだけ使うと混乱が少なくて済みます。 do 〜 end do は使った方が良い do 〜 end do 構文は fortran77 規格には含まれていませんが、 FreeBSD, Linux 付属の f77, g77 や Sun 純正 fortran など、 ほぼ全てのコンパイラでサポートされているので 使った方が良いでしょう。 do while 条件 〜 end do 文も使った方が良いかもしれない do while 文は頻繁には使いませんが、FreeBSD , Linux の f77, g77 と Sun の純正 fortran では使えます。 無限ループ を組むには論理条件の部分を. true. として下さい。 do while 論理条件....... if 条件 exit! ループを抜ける : C の break 文に相当 if 条件 cycle! ループ先頭へ : C の continue 文に相当...... end do• 比較をする演算子について Fortran 77 では. などが if 文の中で 使用されますが、FreeBSD, Linux の f77, g77 は. しかし、 Sun の純正 fortran ではダメなようなので、これはお勧めできません。 使えるコンパイラもある、ということを知っておくにとどめましょう。 整合配列 Fortran ではサブルーチンに配列を渡し、サブルーチン側で 配列のサイズを定義することができます。 次の例を 見て下さい。 fortran の do ループはプログラムの文面に現れない iteration count によって制御されます。 そして、do ループ の range を実行するごとに、 1 do 変数に m3 を加える 2 iteration count から 1 を 引く 3 iteration count を調べ、0 より大きければ do ループを繰り返す という処理を行います。 1 の手順の後に、do ループを抜けるか どうかの 3 の判断するという点が重要です。 end do の do ループが終わった直後、i の値は 11 になります。
次の
