fortran66のブログ

fortran について書きます。

メモ帳:週末晩報

インテル Parallel Studio XE 2018 Composer Edition update 1 が出た!

文書が日本語化されました。Help が使いにくくなったのが何とかなるようなことも書いてあった気が。

Fortran 言語的には、C 言語との相互運用部分が引き続き実装されたようです。

Fortran 2008 および Fortran 2015 機能の概要
インテル® Fortran コンパイラーは、Fortran 2008 標準のすべての機能をサポートします。また、 Fortran 2015 標準の草案のいくつかの機能もサポートします。その他の機能は将来のリリースで サポートされる予定です。このバージョンでサポートされる Fortran 2015 標準の草案の機能は、 次のとおりです。
• 「Technical Specification 29113 Further Interoperability with C」のすべての機能。 次の機能を含みます。 o 型引き継ぎ (TYPE(*)) o ランク引き継ぎ (DIMENSION(..)) o 互換性のある仮引数の制約の緩和 o Fortran で使用される C コード操作「C 記述子」を定義する C インクルード・ファ イル ISO_Fortran_binding.H

AMD の次世代 GPU に対応した ARM 用? GFortran

AMD 死んだかと思いきや、最近復活してきているようです。NVidia CUDA 憎しか IntelGPUAMD と組らしいです(http://www.geocities.jp/andosprocinfo/wadai17/20171111.htm)。日米欧 ARM スパコンAMD/ATI GPU がくるのかもしれません。
www.electronicsmedia.info

よく分かりませんが航空宇宙・軍事産業もやってるようなので Fortran が必要なのかも。
航空宇宙および軍事分野向けソリューション - Mentor Graphics

2016年の記事 スパコン故障原因

ASCI Q は1時間でクラッシュ、MacG5 クラスタはブートする前にクラッシュ、 IBM Blue Gene/L は1次キャッシュを切って利用とか、しばらく前に話題になった内容ですが、ただで読めるそうです。
spectrum.ieee.org

放電管も宇宙線による電離を最初の種火にしてるんだと言われてました。

2006年の記事 2,500万行の Fortran コードを使っているブルームバーグ

25 million lines of Fortran | Coding the markets

One of my colleagues is ex Bloomberg. He tells me there's only one thing I need to know about them: "25 million lines of Fortran !" Apparently they have 1500 coders firefighting, just keeping that legacy monster alive.

Sourcery Institute の Fortran ready な Virtual Box Image が更新 (2017.11.09)

と、思ったが、Guest session が存在せずログインできないので無事死亡
www.sourceryinstitute.org

比較的最近の雑誌における Fortran 記事

ここで比較的最近とは約四半世紀前までを指しますw

トラ技コンピュータ(CQ出版) 1994年7月、8月、10月号

http://www.cqpub.co.jp/try/shoseki/BN94.htm

解説 FORTRAN90  小西浩/大道学
第1回 派生型・ポインタ・動的割り当て
第2回 配列演算機能
第3回 数値計算の応用例
全三回

基本的な Fortran90 の文法を解説しています。第3回の応用例では77と90の両方のコードを出して比較解説しています。

Win3.1時代ゆえパソコン用 Fortran90 コンパイラが存在せず、MS-FORTRAN5.1 の 似たような機能をコラムで多く紹介してます。

月刊 C MAGAZINE (SBクリエイティブ) 1991年12月号

日本語版MS-FORTRAN Ver.5.1  中島信行
honto.jp

MS-FORTRAN Ver.5.1 の紹介記事です。コンパイラオプション、統合環境、簡易グラフィックス、デバッガ等が紹介されています。

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

招いたのは東トルキスタンチベットなどの後援組織でしょうか。

www.sankei.com

http://www.sankei.com/images/news/171115/wor1711150029-p1.jpg

トランプ氏の“参謀”バノン氏「アジア歴訪は素晴らしかった」、「覇権広げる」と中国批判 東京で講演

 トランプ米大統領の有力側近で首席戦略官兼上級顧問を8月まで務めたスティーブン・バノン氏が来日し、15日に都内で講演した。14日に終了したトランプ氏初のアジア歴訪を「素晴らしかった」と絶賛する一方、経済力などを背景に対外拡張を図る中国を批判した。

 バノン氏は、トランプ氏を「世界で最も賢い交渉人の一人だ」と持ち上げ、今回のアジア歴訪中のトランプ氏の言動を絶賛。各国指導者と「強い人間関係を築いた」と強調した。緊張が増している北朝鮮問題については「バランスを取った対応をした」との見方を示した。

 バノン氏は中国について「人民元を米ドルに変わる基軸通貨にしようとしている」と指摘し、金融面で覇権を広げることが目的だと分析した。「(そうなれば)自由民主主義や自由経済などが敗北する」と警鐘を鳴らし、米国が「中国の台頭にどう対応するか」が重要な課題と指摘した。

 辞任後もトランプ氏と頻繁に連絡を取っているとされるバノン氏は、9月には訪中し、習近平国家主席の腹心である王岐山氏との「秘密会談」を英メディアに報じられるなど、“トランプ氏の密使”としての動向も注目されている。

 バノン氏は中国の人権問題などについて研究する「諸民族青年リーダー研修会」で講演した。(岡田美月)

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

Happy 60th birthday, Fortran

珍しく、オープンソース系のサイトでの記事です。
大したことは書いてありませんw スヌーピー・カレンダーが出てきます。
opensource.com

Hybrid Fortran Pulls Legacy Codes into Acceleration Era

日本の論文が元ネタの記事です。Nvidia 利用
www.nextplatform.com

Can Vector Supercomputing Be Revived?

NEC SX Aurora 紹介記事です。SX シリーズ回顧も入ってます。
www.nextplatform.com

Flipping the Flops and Reading the Top500 Tea Leaves

Top 500 リストが中国のむちゃくちゃ登録のせいで歪んでいることを書いています。
一つの構成で測定がなされると、同じ構成の計算機は全部同じスペックで登録できるらしいです。

中国は政府の命令で、日本でいう楽天みたいな民間商業サイトまで動員かけられて linpack 測ってるという話が数年前ありましたが、まだやっているのかな?毛沢東の大躍進に倣ってどんどんやってくださいwww
www.hpcwire.com
 

Gustafson の面白記事

SC17 での Gustafson

Gustafson & Mullin の共著記事が結構面白かったです。
www.hpcwire.com

PFD版 Gustafson & Mullin
https://arxiv.org/ftp/arxiv/papers/1709/1709.09108.pdf

話題は昔ばなしから始まり、近頃流行りの AI テンソル計算、結合則を満たす厳密計算など話題を二転三転させながら、最後は IEEE754 の悪口で終わるワイドスクリーンバロック調で、読んで面白い記事でした。

八だよ、七だよ、六だよ、五
四だよ、三だよ、ニだよ、一
≪もっと引っぱる、≫いわくテンソル
≪もっと引っぱる、≫いわくテンソル
緊張、懸念、不和が来た

分解された男 (創元SF文庫)

分解された男 (創元SF文庫)


特に興味を引いたのは、「The Rise of the Installer Program」の段落で、従来のコンパイラ=リンカ=ローダー的な枠を崩して、コンパイラ(ないし「プリプロセッサ」)が中間言語を生成し、その後ハードウェアの構成に合わせたデータ構造を与える「インストーラ」が実行バイナリを構成するというアイデアです。

JIT 的な発想ですが、加速器ハードウェアが多様化しつつある今、特定のハードウェア構成やネットワーク・トポロジーに最適化された個別のソース・プログラムを書かざるおえない現状への面白アイデア提起かと思いました。

The End of Error: Unum Computing (Chapman & Hall/CRC Computational Science)

The End of Error: Unum Computing (Chapman & Hall/CRC Computational Science)

この他に気に入った部分として、データの順序を考えないで無構造な集合としてデータを与えるプログラムモデルもあります。言語側がデータをハードウェア上に最適に配置するのとセットになっているアイデアです。

最近イテレータにとらわれ過ぎて、本来順序のない(並列処理が自明な)構造にまでわざわざ順序を導入し全順序集合にしたあとで、逐次処理を結合則を利用して並列化出来ると言ったりしていますが、どこか狂ってると思います。

*本文で引用されていた記事

John Gustafson The Vector Gravy Train
ベクトル機はぼったくり 
The Vector Gravy Train

英 coarray ビデオ講義

イギリスのビデオ講義集ですが、視聴回数が2年で10回程度のものが並んでいて悲しいですw イギリス人はラジオか電信で。 俺も見てないけどそのうち見る。

Anton Shterenlikht - Writing parallel programs with Fortran 2008 and 2015 coarrays


Anton Shterenlikht - Writing parallel programs with Fortran 2008 and 2015 coarrays

関連DOC
http://coarrays.sourceforge.net/doc.pdf

The Fortran coarrays course


ARCHER Webinar: eCSE 0505 Open source exascale multi-scale framework

Apery の定数 その2

加速

正攻法で和を取ると\zeta(3)=\sum_{n=1}^\infty {1\over n^3}は収束が遅いので、加速します。

\sum{1\over n(n+1)(n+2)}={1\over4}より、\zeta(3)-{1\over4}=\sum{3n+2\over n^3(n+1)(n+2)}という関係式を使うと、nの4乗に比例して項が小さくなるので、少し加速されます。

実行結果

1000項まで足してみました。項の大きさは n^4 の大きさになっていますが、和の収束は n^3 くらいになっているようです。

中央二項係数による方法はとても収束が早かったですが、それに比べるとゆっくりしてます。

 Apery's constant =   1.20205690216284
続行するには何かキーを押してください . . .

ζ(3) = 1.202056903159594285399738161511449990764986292

ソース・プログラム

    program Apery
      implicit none
      integer(selected_int_kind(10)) :: i, n = 1000
      real(kind(1.0d0)) :: s, t, u
      s = 0.25d0
      do i = 1, n
        s = s + (3 * i + 2) / real(i**3 * (i + 1) * (i + 2), kind(1.0d0)) 
      end do    
      print *, 'Apery''s constant =', s
    end program Apery