fortran66のブログ

fortran について書きます。

2017-01-01から1年間の記事一覧

スティーブン・バノン氏来日!

招いたのは東トルキスタンやチベットなどの後援組織でしょうか。www.sankei.com トランプ氏の“参謀”バノン氏「アジア歴訪は素晴らしかった」、「覇権広げる」と中国批判 東京で講演 トランプ米大統領の有力側近で首席戦略官兼上級顧問を8月まで務めたスティ…

メモ帳:最近の Fortran 記事 など

Happy 60th birthday, Fortran 珍しく、オープンソース系のサイトでの記事です。 大したことは書いてありませんw スヌーピー・カレンダーが出てきます。 opensource.com Hybrid Fortran Pulls Legacy Codes into Acceleration Era 日本の論文が元ネタの記事…

Gustafson の面白記事

SC17 での Gustafson Gustafson & Mullin の共著記事が結構面白かったです。 www.hpcwire.comPFD版 Gustafson & Mullin https://arxiv.org/ftp/arxiv/papers/1709/1709.09108.pdf話題は昔ばなしから始まり、近頃流行りの AI テンソル計算、結合則を満たす厳…

英 coarray ビデオ講義

イギリスのビデオ講義集ですが、視聴回数が2年で10回程度のものが並んでいて悲しいですw イギリス人はラジオか電信で。 俺も見てないけどそのうち見る。Anton Shterenlikht - Writing parallel programs with Fortran 2008 and 2015 coarraysAnton Shter…

Apery の定数 その2

加速 正攻法で和を取るとは収束が遅いので、加速します。より、という関係式を使うと、nの4乗に比例して項が小さくなるので、少し加速されます。 実行結果 1000項まで足してみました。項の大きさは n^4 の大きさになっていますが、和の収束は n^3 くらいに…

ARM Fortran 正式版出る!・・・ が、ページがバタ臭くてw

ARM が HPC 用コンパイラツール一式を正式出荷したようです。 www.hpcwire.comハイライト孫引き引用 Arm Fortran Compiler,released as a beta compiler in June 2017, is now a fully supported commercial Fortran compiler with Fortran 2003 and prior s…

ζ(3) アペリーの定数

中央二項係数でアペリーの定数 中央二項係数の逆数和(その6) より、 ζ(3)=Σ1/n^3=5/2Σ(-1)^(n-1)/n^3(2n,n) を用いて、ζ(3)アペリーの定数を計算してみます。 実行結果 Central Binomial Coefficients 2mCm 1 2. 2 6. 3 20. 4 70. 5 252. 6 924.…

最近の Fortran 本

第3版 有限要素法による流れのシミュレーション OpenMPに基づくFortranソースコード付作者: 日本計算工学会出版社/メーカー: 丸善出版発売日: 2017/07/29メディア: 単行本(ソフトカバー)この商品を含むブログを見る建築構造設計・解析入門 Fortran 解析プ…

【ニュース遅報】Forcheck 社買収されてた

新聞によりますと、Fortran ソースコード静的解析プログラムの Forcheck を開発販売していた 蘭 Forcheck社が今年の1月に買収されたそうです。そしてその買収会社がまた買収された模様です。 www.tomihisa.co.jp Forcheck は地味に Fortran2008 まで対応し…

中央二項係数 再び

中央二項係数 中央二項係数とは、二項係数の丁度真ん中の値になります。2nCn とも書き表せます。 以前も計算したのですが、32bit 整数の範囲だったのであまり多くは求められませんでした。 fortran66.hatenablog.com 今回 64bit 整数で計算してもいいのです…

Fortran 利用状況情報サイト

Companies using Fortran 商利用の情報のサイトのようです。 idatalabs.com利用者の多くは、科学技術研究機関・重厚長大産業の模様です。

Wilson の定理で素数を求める

Wilson の定理 計算量が多くて役に立たないとされる Wilson の定理をたわむれに計算して見ます。Wilson の定理とは、素数の時 (p-1)! ≡ -1 (mod p) が成り立ち、合成数では 0 になります。 実行結果 10^3 以下の素数リストと個数 2 3 5 7 11 13 17 19 23 29 …

NEC の 次期ス―パーコンピュータ SX-Aurora TSUBASA はキャプテン翼から

PCIe に差すボードの形で、完全に独立にジョブを投げる構成の模様です。最近ブームの AI も、線形代数の問題のようなので、なんとなく成功しそうな気がします。コンパイラは従来のまま Fortran、C、C++ で問題もなさそうです。 TSUBASAは、「欧州から和名で…

デバッグの第二法則??

第一法則はなんじゃい?w The second law of debugging - Loeser - 1976 - Software: Practice and Experience - Wiley Online LibraryThe second law of debugging Authors R. Loeser, E. M. Gaposchkin これによると、第二法則は以下のようです。 If you …

メモ帳:週末小ネタ集

"you cannot do all of these at the same time." パッと見、良さげな約束をしてくれますが、同時には成り立たないのは基本w The problem with Common LispNow I think that one of the main reasons for this is that while you can write scientific code…

Using OpenMP 新刊

Using OpenMP の新刊が出るようです。 元々はベンダー毎にばらばらだった並列化 FORTRAN 用のコンパイラ指示行を統一する規格でしたが、20年の時を経て謎の進化を遂げています。Using OpenMP -- The Next Step: Affinity, Accelerators, Tasking, and SIMD…

1982 年の FORTRAN 誕生 25周年記念の集まりの記事

Annals of the History of Computing Issue 1 • Jan.-March 1984FORTRAN 誕生 25周年記念特集 IEEE Xplore: Annals of the History of Computing - ( Volume 6 Issue 1 )このビデオが作られた会のようです。The Beginnings of FORTRAN (Complete)全部見たわ…

メモ帳:Modern Fortran の構造について

Fortran の構造 Fortran の分類 Fortran は、約 60 年間に渡って変化・発展してきたプログラム言語で、言語内容からみて FORTRAN I から FORTRAN77 までの 古典 FORTRAN (Classical FORTRAN) と、Fortran90 以降の現代 Fortran (Modern Fortran) に大きく分…

メモ帳: N. Wirth の Modula-2 の function

N. Wirth の Programming in Modula-2 1988年版はフォントが小さいので、function の章が 2 ページに収まって、Springer の preview で 全部読めるwProgramming in Modula-2 | Niklaus Wirth | Springer 関数の副作用について まとめ 3.副作用について 4.…

function の side-effect について :追記あり

Fortran では、function の副作用は嫌われていて、文法的には許されていますが、道徳的に禁止になっています。Fortran95 以降では pure 指定子が出来たので、常に pure というのが行儀のよさだと思います。副作用のある場合は、サブルーチンを使います。この…

行列積を Cannon の algorithm で

昔、コネクションマシンの本を読んだときに、各並列プロセッサは 1bit 演算で、浮動小数点数も 1bit 演算を繰り返して計算すると知って、これは駄目ですわと思ったのですが、グスタフソンの動的可変長浮動小数点数フォーマットのことを思うと、任意長に出来…

メモ帳:パタヘネ付録 第3版  歴史展望より

Fortran90 では、多種多様な浮動小数点数フォーマットの乱立に対して、浮動小数点数フォーマットに関する情報を得る様々な関数を導入したのですが、ちょうどその頃 IEEE754 フォーマットにほぼ天下が一統されて、やや役立たずに終わりました。Linpack/Eispac…

RATFOR の Quick Reference Card

archive.orginternet archive に RATFOR の三つ折り Quick Reference Card がうpされておりました。昔はよくこの三つ折りの厚紙の Quick Reference Card がありましたが、あまり役に立ちませんでしたwRATFOR は古典 FORTRAN 時代に、構造化プログラミング…

ナウな言語とのベンチマーク

天文系の計算で Fortran と最近のナウい言語とのベンチマークがなされています。Fortran も 2008 水準のモダンなものを用いているようです。 Compiled languages: Fortran: GNU Fortran 6.3.1 C++: g++ 6.3.1 Java: JDK 1.8.0-74 Interpreted languages: Mat…

coding ground の gfortran が ver.7.1.1 になるもあまり動かず?

ブラウザ上で動くオンライン・コンパイラの coding ground が微妙にインターフェースを変えるとともに、コンパイラのバージョンをアップしてくれたのですが、パソコンからだとなかなか動いて呉れません。たまに動きます。またスマホ・タブレットだとやや動き…

メモ帳:FORTRAN 77 の制御構造が BLOCK IF 文だけになった経緯

昨日のエントリで日本での事情を調べましたが、米国での事情についてもメモしておきます。fortran66.hatenablog.com Martin N. Greenfield, History of FORTRAN standardization (1982) History of FORTRAN standardization時系列的には、FORTRAN 66 規格の…

FORTRAN の言語としての特性

FORTRAN の言語としての特性 菅忠義、プログラミング言語の最新の動向:3. FORTRAN 情報処理 22 (1981) 452 id.nii.ac.jpより引用。 FORTRAN II 1) IBM: The FORTRAN Automatic Coding System for the IBM 704 EDPM. 数値計算に関する任意の問題を表現可能。…

メモ帳:Knuth の An empirical study of FORTRAN programs

Donald E. Knuth が 1970年に "An empirical study of FORTRAN programs" という論文の中で、FORTRAN プログラムの解析をしています。年代からここでの FORTRAN は FORTRAN IV/66 水準のものだと分かります。すなわち制御が構造化される前です。結論からする…

メモ帳:情報学会誌『情報処理』の Fortran 規格関係記事 

■FORTRAN IV 言語の概略 1962年 Vol.3 No.3 西村恕彦 情報学広場:情報処理学会電子図書館■ALGOLとFORTRANのJIS原案 1967年 Vol.8 No.1 井上謙蔵 情報学広場:情報処理学会電子図書館■FORTRAN 77 の特徴:JIS FORTRAN との非互換性の観点から 1979年 Vol.20 …

bit 1972年8月号 グラビア COM による原子力発電所の地震シミュレーション解析(福島第一原発)

グラビアページにはコンピュータ・グラフィックスによる振動モード図が、記事には細かな解説がありました。 地震には勝てても、学生運動の反米左翼テロ活動で爆発