「FORTRAN」を含むブログ

NAG の Fortran 検定 fortran66fortran66 2012-05-25

...tran/exam/index.html)をやってみました。文字列がらみの問題を軒並み間違ってしまいました。色々勘違いしている。ヤヴァイwww 文字変数は読み込み時は左詰め、書き出し時は右詰めw FORTRAN77 規格には倍精度複素数はないが、Fortran90 以降にはある。 <iframe src="http://rcm

sine and cosine curves fortran66fortran66 2012-05-24

昔よく見た文字グラフです。特に意味なしw Fortran2008 では、万能 G 型の Format がさらに拡張されて g0 によってあらゆる型を、よきに計らってくれます。また繰り返し指定子に * を使うことで無限反復となります。実行結果 ソース・コードprogram test implicit none real, parameter :: pi = 4 * atan(1.0) integer...

古書目録(06-02-06) kosyo1kosyo1 2012-05-23

古書目録(06-02-06) 詳細は「Bsearch書誌検索」で確認して下さい。 http://k-pj.com/ob/ http://k-pj.com/ob/ddata/search.cgi?mode=day&word=06-02-06 http://kosyo.net a3032 n32683 4-535-jp55008149-a 日本経済刑法概論 定塚道雄 日本評論社 - a3022...

Introduction to Modern Fortran fortran66fortran66 2012-05-23

ここの講義ノートが重箱の隅をつついていて結構面白いです。 http://people.pwf.cam.ac.uk/nmm1/Fortran/ 文字配列のメモリー上での置かれ方なんて考えたことがありませんでしたが、文字長が余分な配列の RANK として最内に加わった感じになっているそうです。これはもっともな気がします。実行結果 ソース・プログラムprogram Console6 implicit n...

エラー行を視覚的に教えてくれるプラグインsyntasticを導入してみた guyonguyon 2012-05-22

http://gyazo.com/6daa987520906ac7556413539dcdb911.png https://github.com/scrooloose/syntastic これを設定するとQuickFixで指定されたエラー行を視覚的に表示してくれます。 実際にコードを書いている時はエラーが表示されるようなケースは少ないないですが、lintチェックもできるみたいなので、そっち方面で...

Newton法 fortran66fortran66 2012-05-22

少し改良。実行結果 ソース・プログラムmodule m_sub implicit none abstract interface pure elemental real function t_f(a, x) real, intent(in) :: a, x end function t_f end interface contains recursive function xiter(f, a, ...

Newton 法 fortran66fortran66 2012-05-21

最近関数プログラミングというのが流行っているわけですが、FORTRAN の父 John Backus は、関数プログラミングに対しても先駆者となっています。 昔 von Neumann モデルの限界がどうのこうの巷で騒がれたころに、Bakus のチューリング賞受賞講演 Can Programming Be Liberated from the von Neumann Style? http://w...

FizzBuss Fortran fortran66fortran66 2012-05-20

文字型変数と他の型の変数の変換には内部ファイルへの READ/WRITE を用います。内部ファイルへの I/O があった場合、pure 性が保たれるのか気になったので試してみました。結論からいうと pure でした。 内部ファイルの I/O は型変換の代入に近いので、もっともな仕様だと思います。66 時代には ENCODE/DECODE 関数で型変換してわけだし。 pure で elemental...

Minimal sample for reading a fortran unformatted data from Mathematica arakik10arakik10 2012-05-19

Reference : Google search: fortran unformatted file structure a sample fortran source code INTEGER,PARAMETER::NDAT=256; REAL(8),DIMENSION(0:NDAT-1)::dat=1D0 OPEN(1,FILE='dirname/filename',FOR...

eight queens fortran66fortran66 2012-05-19

昔、PASCALの再帰を使った例題でよく出ていた eight queens 問題を解いてみます。 前回の整数並べ替えの列挙を変形して利用します。配列の i 番目の要素が i 列目での queen の存在する行を指していると解釈します。 すでに飛車道を消した形になっているので、角道を消すことを考えます。これは対角線方向に新しい駒を置かないようにすれば実現されます。 なお回転と鏡像も重複して数えること...

順番付並べ替えの列挙 fortran66fortran66 2012-05-18

1からnまでの整数の順番付並べ替えを、再帰を使って列挙します。実行結果 ソース・プログラム , pack(list2, list2 /= list2(i)) この行では、list2 から要素を一個取り出して消して、list1 に加えています。配列をリストとして用いています。 は整数型の大きさ 0 の配列です。空リストの代用として使っています。 関数呼び出しのスタックを活用しているので、再帰呼び出し...

vi/vim 設定 コマンド 私的まとめ nightmare_timnightmare_tim 2012-05-17

vimの設定コマンドの私的まとめです。 基本は、各コマンドの入力はvimのスタンダードモード時に入力します。項目によりますが設定を永続化させたい場合は~/.vimrcに記述してください。また、論理値オプションと書かれた項目は、以下のような形で、オフにできます。:set no○○○○ ||-」に変更する。|:set listchars=tab:>-,trail:-| |コメント行の高さを3行に...

ソースコードのはてな記法 pospomepospome 2012-05-16

よく忘れるので、メモ。 a2ps a65 aap abap abaqus abc abel acedb ada aflex ahdl alsaconf amiga aml ampl ant antlr apache apachestyle arch art asm asm68k asmh8300 asn aspperl aspvbs asterisk asteriskvm atlas automa...

Napier's constant または Euler's Number fortran66fortran66 2012-05-16

Fortran2008 では、iso_fortran_env という intrinsic module が導入されています。その module には、単精度や倍精度実数の kind が定義されています。実行結果 ソース・プログラムprogram Napier use, intrinsic :: iso_fortran_env implicit none integer, parameter :: ...

関西Emacs勉強会 kenbeesekenbeese 2012-05-13

今更ながら関西Emacs勉強会に参加して、気になったelispを導入してみたので、その感想でも書こうかな。 quickrun 正直このさん作のelispは知っていたんだけど、Fortranがないし、まぁいっかと思ってた。 しかし、話を聞いてみると拡張が簡単らしい。 で、Fortranに対応させて、他の言語でも使ってみたら便利だった! 特にちょっとした試し書きにはめっちゃ便利!! 以下Fortran...

move_alloc を用いた整数の分割 fortran66fortran66 2012-05-13

昨日のバージョンでは、整数の分割は画面出力されるだけで、データとして値を返していませんでした。これを配列の配列として返すことを考えます。整数配列を要素として持つ構造体を返すことにします。 返り値となる割り付け配列の大きさを増加させていくために、Fortran2003 で導入された組み込みサブルーチン move_alloc を用いました。この方法では、サイズを増やす毎に、要素の全コピーによる移動が行...

整数の分割 fortran66fortran66 2012-05-12

以前整数の分割数 p(n) を漸化式によって求めましたが、今回は分割そのものを求めることにします。 再帰を使った導出法は、富田氏による 『Fortran90プログラミング』例題8.10を参考としました。 以下のプログラムでは、配列生成子 を用いて、サブルーチンの引数をスタックのように利用しています。とても短く記述できるので大変便利です。メインプログラムでは、大きさ0のダミー配列を与えればいいので...

Fortran のポインタがちょっと不便 ktakatamaktakatama 2012-05-10

Fortran にも 90 以降はポインタと名のつく機能があり,ある程度 C っぽいプログラムを書くことができます.幸いなことにポインタ演算機能はありませんが,構造体の連結リストを作成したり,下記のように部分配列を切り出すのに便利です.切り出しに際してデータのコピーを必要としないので,おそらく効率が良いのだろうと思います(未確認).〜前略〜 integer, parameter :: isize ...

10 May 2012 曇り一時雨 pokoritpokorit 2012-05-10

今日は一日中ジョブのセットアップ。やっと自分の研究の順番が回ってきた。午前中Fortran、午後C++。体感速度としてFortranの方がかなり速い。C++はまだ性能の出る「公式」が完全に解かっていない感じ。だいぶ掴めて来たような気もするし、まだまだなのかも知れない。今日のところはFortranとC++の結果が一致したことで良しとしよう。

プログラミングフォント・Rictyをインストールしてみた otome_pwotome_pw 2012-05-10

...で生産性がそんなに変わるわけがない……。と今でも思っている。 Ricty (リクティ) は Linux 環境での研究・開発に適したフリーのプログラミング用 TrueType フォントです。C、C++、FORTRAN、Python、Perl、Ruby、AWK、シェルスクリプト、Makefile、LaTeX など、UNIX 系のコーディングにおける使用を想定しています。以下の 2 つの等幅フォン...