fortran66のブログ

fortran について書きます。

【乞食速報】ACM Fortran forum 無料! 6月末まで

ACM SIGPLAN Fortran Forum が無料化

そろそろ4月号が出るころかと見に行ったら、いつの間にか free access になっていました。なお 4月号はまだのようです。無料は武漢肺炎に鑑みて ACM すべてにわたり期間は6月末までのようです。ちゃんコロナ様のおかげです。

今までも推薦してきましたが、とりあえず Gorelik の文法規格差表を落としておくとよいのではないかと思います。

Statements, data types and intrinsic procedures in the Fortran standards (1966-2008)

https://dl.acm.org/doi/abs/10.1145/2701654.2701655

The Fortran 2003 Handbook: The Complete Syntax, Features and Procedures

The Fortran 2003 Handbook: The Complete Syntax, Features and Procedures

新しめの free access 本もあるようです。 https://dl.acm.org/acmbooks/archive

【メモ帳】SciPy 中の Fortran ソースの割合 2位 25%

SciPy 中の Fortran ソース割合

Python に次いで多く四分の一。より詳しいソースがあった。

Language choices. According to analysis using the linguist library (https://github.com/github/linguist), SciPy is approximately 50% Python, 25% Fortran, 20% C, 3% Cython and 2% C++, with a dash of TeX, Matlab, shell script and Make.

www.nature.com

武漢肺炎】NY市人種毎の死亡率

NY市長による人種毎の死亡率の差。 アジア系(主に中国人だと思うが)の死亡率が低い。 f:id:fortran66:20200409012626p:plain

また男の方が女より死亡率が高くなってもいます

jp.reuters.com

これから導き出される結論は、武漢肺炎は人種差別的で性差別的な反リベラル勢力だということのようです。

サンダース選挙運動止める

予備選候補としては残るらしいので、サンダースもバイデンも年寄り同士コロナが危険なので、むしろ引きこもって相手を陽性にして死なそうという、ユダヤらしいヨセフス戦法とみたw

www.bbc.com

bbc の偏向振りも笑える。

【ニュース】llvm の flang 続報 その他

LLVM f18(旧名)の近況

何か障害があって llvm の正式ツールチェーンに入れなかった flang (旧f18) ですが、その後の動向が記事になっています。あまり大した進展もないようですが、悪いニュースもないようです。

www.phoronix.com

【悲報】ジョンソン首相ICUで酸素吸入

7~10日目に、坂道を転げ落ちるように急速に悪化するというパターンでヤバいです~

www.bbc.com

武漢肺炎】韓国も無・軽症者への検査制限をしている?w

blog.livedoor.jp

縦軸が狂っていて、2万, 5万, 8万 は、1.2万, 1.5万, 1.8万 の間違いのようです。 f:id:fortran66:20200408005214p:plain

6. 韓国人(スレ主)
>>5
実際に、検査を受けられる基準がぐんと上がった

7. 韓国人
>>6
もともとの基準が何で、それがどう上がったのか?

8. 韓国人(スレ主)
>>7
以前は、医師の所見だけあれば検査ができた
つまり、軽症でも予防できた
その後は、CTレントゲンで肺炎の兆候が見られたら検査できるようになった
つまり、重症にならないと検査ができないように変わった

11. 韓国人(スレ主)
>>10
実際、検査ができている基準が、3月10日から、医師の所見からCT X線で肺炎の症状がなければ検査できないようになった

13. 韓国人(スレ主)
>>12
今、病院行って検査を頼んでみろ
100%断られる

この内容が本当だとすれば、検査数は毎週どんどん絞られていて最盛期の半分になっていて、検査方針は日本式になっていることになります。

韓国のデータを見ると死亡者数は指数関数的に増えているのに感染者数は線形にしか増えていなくて、何かおかしいです。

91-divoc.com

ドイツ

そういえば、ドイツの死亡率が低めに出ているのは既往症があると、コロナ感染での武漢肺炎での死亡でも、既往症の病名での死亡にするからという話もありました。本当かどうかわかりません。

いずれにせよ今までのインフルエンザの時も、各地域ごとの速報値は基準がバラバラであてにならず、正確な値が出るのは流行が終わってからデータ整理をしてからで、最終的値ではオーダーが変わるのも普通らしいです。

【武漢肺炎】合言葉は、やってはいけない 集(シュウ)近(キン)閉 (ペイ)

合言葉は、やってはいけない 集(シュウ)近(キン)閉 (ペイ)

3: 2020/04/04(土) 18:49:09.33 id:KuBgJxDN

「集」まらない。

「近」づかない。

「閉」鎖空間を避ける。

「シュウキンペイ」と覚えましょう。

クッソw 三密より出来がいいwww なんとなくコンパクト集合というか有界閉集合を思わせるところもいいです。

『密閉・密集・密接の「3つの密」』これは覚えられないです。

そろそろアンチ 集・近・閉で蟄居・閉門して暮らしますかw

1:2020/04/02(Thu) 08:50:10ID:LDNKXE/Pa.net
やってはいけない事は
習近平 シューキンペイと覚えよう

「集」まるのは危険 →シュー
「近」くで話すのは危険 →キン
「閉」鎖された空間は危険 →ペイ

“3密”なんかより覚えやすいからみんなよかったら広めてや

alfalfalfa.com

【悲報】ジョンソン首相入院

池田隼人首相の『前がん症状』での入院を思い出します。

www.jiji.com

新型コロナ感染のジョンソン英首相入院 官邸「予防措置」

2020年04月06日07時01分

 【ロンドン時事】英首相官邸は5日、新型コロナウイルスに感染したジョンソン首相が入院したと明らかにした。感染判明後10日を経ても症状があるため「予防措置」として検査を受ける。

 首相は3月27日、検査で陽性だったと公表し、官邸で隔離生活に入った。官邸報道官によると、現在も高熱などの症状が続いており、検査のため5日夜に入院した。報道官は「医師の勧めであり、緊急入院ではない」と強調した。引き続き職務は行い、他の閣僚らとも連絡を取っているという。

 自らの病状について首相は3日、ツイッターに投稿した動画で「良くなってきたがまだ熱がある」と報告していた。英政権内では首相のほか、ハンコック保健相も同時期に感染したが、既に回復している。首相の妊娠中の婚約者も感染が疑われる症状を示していたが、本人が4日、ツイッターに「快方に向かっている」と書き込んだ。

【悲報】オナー・ブラックマン女史死亡

コロナとは書いていませんが・・・悼 www.bbc.com

The Avengers: The Inside Story

The Avengers: The Inside Story

【朗報】チャールズ皇太子のお言葉 & エリザベス女王のお言葉

www.bbc.com

www.youtube.com

www.bbc.com

www.youtube.com

Fortran Con

このご時世に、スイスで FortranCon が初開催されるようです。

FortranCon 2020 is the first International Conference targeting the development for and with the Fortran programming language and will take place on July 2 - 4 in Zurich, Switzerland.

tcevents.chem.uzh.ch

www.swissinfo.ch

連邦内閣は3月20日、公園、広場、遊歩道など公共のスペースで6人以上の集まりを4月19日まで禁止する措置を講じた。違反した場合は1人に付き100フラン(約1万1千円)の罰金を科す。5人以下で集まる場合も、互いに2メートル以上の「ソーシャル・ディスタンシング(社会的距離)」を保つことが望ましいとしている。

COBOL ネタにちなんで

【武漢肺炎】コロナでコボラー緊急募集!

ニュージャージー知事がボランティア求む!

但しゴールドマンサックス崩れの知事さんは「コバルト言語」と言ったようですw よく分かりませんが、失業手当申請者が激増してコボルのプログラムに支障が出ているようです。

josephsteinberg.com


New Jersey Needs COBOL Programmers

武漢肺炎の死者飽和について

以下のサイトで百万人当たりの死者や感染確認者をプロットできるのですが、西洋諸国の死者数が百万人当たり500程度で飽和しそうに見えます。

国民総人口で規格化しているようですが、中国では約一千万人規模の武漢市周辺の結果とみなすと、100倍して考えればよかろうかと思います。イタリアも北部地方中心だし、イランは人口ピラミッドが若いなど、地域人口の比率・老齢人口の比率などもありましょうが、なんとなく百万人当たり 500~1,000 人の死者になりそうです。さらに死亡率を感染者数の 1%~10% で考えると、百万人当たりの感染者は 5,000~100,000 人となります。

すると死亡者が飽和する感染者は人口の 0.5%~10% となって、集団免疫が出来上がるという 60%~90% という割合と大きな差があります。

移動制限を外すと再燃するのか、国家単位で比較するのが不適切で局所的に集団免疫レベルに達したのか、よく分かりません。

【追記: R2.4.7 】 NHK-BS でやっていた フランス2のニュースで、武漢肺炎が猖獗を極めたイタリアの北部地方で抗体検査を試験的に小数名に行ったところ、抗体を持っている比率が予想より少ない10%程度だった、と言っていました。ネット上ではイタリアで試験的に抗体検査が始まったというニュースしか見つかりませんでしたので、ソース的にイマイチですが、やはり集団免疫の理論値にはほど遠い状態のようです。

91-divoc.com

f:id:fortran66:20200405195737p:plain

日本の感染者数

上記死亡者グラフは 100 名以上死亡の国からのデータでしたので、日本のグラフは描けませんでした。しかし死者数に比べるとさらに不定要因が増えますが、人口百万人当たりに規格化した感染者数の方は描けます。

f:id:fortran66:20200405202338p:plain

対数プロットでの感染者数でみると、日本も直線に乗っているので指数関数的に増えていることは他国と変わりがありません。しかし日本の傾きは他の国々よりはるかに緩やかで、4~5倍の時間をかけてゆっくりと事態を推移させえいることが分かります。(直近はハッスルしてますが)となると、夏頃にピークがやってきて毎日数百人の死者が出る日々が数か月にわたって続くのではないかと思われます。

【メモ帳】Dijkstra, Wirth n 個の素数の表

Fortran への移植改変

ソース・プログラム

    program nprime
        implicit none
        integer, parameter :: np = 100, nroot = int(sqrt(real(np)))
        integer :: ip(np), iv(nroot)
        integer :: k, i, ix, isqr, lim
        logical :: qprime
        ip(1) = 2
        isqr = 2**2
        ix  = 1
        lim = 1
        write(*, '(i0, 1x)', advance = 'no') ip(1)
        do i = 2, np
            do
                ix = ix + 2
                if (isqr <= ix) then 
                    iv(lim) = isqr
                    lim = lim + 1
                    isqr = ip(lim)**2 
                end if
                k = 2
                qprime = .true.
                do k = 1, lim - 1
                    if (iv(k) < ix) iv(k) = iv(k) + ip(k)
                    qprime = ix /= iv(k)
                    if (.not. qprime) exit   
                end do    
                if (qprime) exit
            end do     
            ip(i) = ix
            write(*, '(i0, 1x)', advance = 'no') ip(i)
        end do
        print *
    end program nprime

実行結果

最初の 100 個の素数

2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59 61 67 71 73 79 83 89 97 101 103 107 109 113 127 131 137 139 149 151 157 163 167 173 179 1
81 191 193 197 199 211 223 227 229 233 239 241 251 257 263 269 271 277 281 283 293 307 311 313 317 331 337 347 349 353 359 367 373 379 3
83 389 397 401 409 419 421 431 433 439 443 449 457 461 463 467 479 487 491 499 503 509 521 523 541

アトム誕生の巻

アトム誕生の巻

  • 発売日: 2018/08/31
  • メディア: Prime Video

NHK 漢詩の時間

NHK 漢詩の時間も新年度が始まりました。講師は引き続き佐藤正光です。今期のテーマは『信』のようです。巧言令色で虚飾で嘘を騙る中共政府と相容れない言葉です。

www4.nhk.or.jp

しかし佐藤先生、進駐軍の民生局が露骨に選挙に干渉して、当選しそうな候補を片っ端から公職追放するなか、どさくさまぎれに総理大臣になった芦田均を称揚するなど、狂気を垂れ流しているので油断できません。終戦直後の選挙は、今の世の中だったら、国際選挙監視団によって選挙無効が宣言されるような有様です。

第51回「蓮を愛す」で出てきた「江南」(前漢後漢 無名氏)の詩で、リフレインの所が東西南北の順で方角が出てきたのが気になります。支那では甲骨文の頃から時計回りに東南西北の順で方角が出てくるのに、ここでは日本式に東西南北です。脚韻でも踏む都合かと思いきやそうでもないようです。

清明

清明節と言えば、『人面桃花』という唐代伝奇小説。その昔、どんな怪奇・奇談かと思ってワクワクして読んだら、ただのご都合主義のボーイミーツガール式小説で桃の花のように美しい顔とかふざけんな!と思った記憶が。

そんな中、中共唯物論のくせにお盆相当の清明節を国家規模で行事にしたようです。

武漢肺炎の犠牲者を悼んだらしいのですが、看板の文字を見る限り烈士wに祭り上げた連中と同胞とのみを悼んでいるようで、海外の国々に人的・経済的被害を与えたことに関する責任は、これっぽっちも感じていないようです。支那人の本質を見るようで、漢詩の理解も進みます。

f:id:fortran66:20200405211001j:plain

news.tbs.co.jp

中国、“警鐘”の医師ら14人の「烈士」を哀悼

 中国は4日、祖先を供養する伝統的な祭日「清明節」を迎えています。

 中国政府は、新型コロナウイルスの対策にあたり犠牲となった医療関係者など14人に対し「烈士」の称号を与えていますが、清明節を迎えた4日、哀悼の意を表する活動が各地で行われました。

 国営メディアによりますと、烈士は「党と国家が授与する最高の栄誉ある称号」で、原因不明の肺炎にいち早く警鐘を鳴らしたにもかかわらず警察に処分され、その後、自らも感染し死亡した医師、李文亮さんにも与えられています。

 一方、中国の国家衛生当局によりますと、新型コロナウイルスによる中国本土での死者は3日、新たに湖北省で4人確認され、合わせて3326人となっています。

【メモ帳】おしっこは貴重な資源!

おしっこムーンベース 

gigazine.net

アフリカには牛糞で作った家とかありますが、月面基地がおしっこで出来るとすれば夢のワンダーランドですね!

f:id:fortran66:20200403013413j:plain

f:id:fortran66:20200403023507p:plain

構造化プログラミングの例題

 ダイクストラ、ホーア、ダール『構造化プログラミング』 

 N. Wirth『Systematic Programming』

ダイクストラ、ホーア、ダール共著の『構造化プログラミング』中のダイクストラ「構造化プログラミング論」の中にダイクストラ口述試験に出した問題として、

「0,1,2 でなる長さ1以上の列で、隣り合う部分列はどれも等しくないものを考え、そのような列を、長さ100のもの(つまり100個の数字の列)になるまで、次々とアルファベット順に創り出すプログラムを構成すること」

というのがあって、ダイクストラがウダウダ愚図っていて考える気もせず読み飛ばしていたのですが、この設問の考案者である N. Wirth の 『Systematic Programming  An Introduction』を眺めていたら最終章の例題の中に出てきて、stepwise program development の実例として考え方まで開陳してくれていたので Fortran で書いてみました。

Wirth の手にかかると簡単になって、かつ数字の対称性から数の置換は本質的に繰り返しになっていて本質的に違いのある配置数は 1/6 になるなど、啓蒙させられました。

なお、ダイクストラの出した別の問題

「2 つの数の n 乗和で、2通り以上の形で表わされるもので最小のものを求めよ。」

というのもあって、同じように解説してくれていました。3乗和の例題によって、ハーディの語るラマヌジャンの病床での逸話の数 1729 が求められています。

Fortran

ソース・プログラム

色々めんどくさいので、サブルーチン中でプログラム終了するお下品プログラムにしましたw なお Wirth は長さ3で数は 0,1,2 ではなく 1,2,3 でやっているのでこれに従いました。Wirth が表に与えている長さ20までの本質的に独立な列の数はチェックしました。

    program seq123
       implicit none
       integer, parameter :: n = 3, k = 3 ! n:length  k:1,2,3..k
       integer :: is(n), m
       
       is(1) = 1
       m = 1
       do 
           if (m == n) then
               if (is_ok()) print '(20i3)', is 
               call change()
           else 
               call extend()
           end if
           do 
               if (is_ok()) exit
               call change()
           end do
       end do    
    contains 

       subroutine extend()
           m = m + 1
           is(m) = 1
       end subroutine extend
       
       subroutine change()
           do 
               if (is(m) == k) then 
                   m = m - 1
               else
                   exit
               end if
               if (m <= 0) stop 'normal end'
           end do
           is(m) = is(m) + 1
       end subroutine change
       
       logical function is_ok()
           integer :: i
           is_ok = .true.
           do i = 1, m / 2
               is_ok = any(is(m-2*i+1:m-i) /= is(m-i+1:m)) 
               if (.not. is_ok) exit 
           end do    
       end function is_ok
       
    end program seq123

出力

  1  2  1
  1  2  3
  1  3  1
  1  3  2
  2  1  2
  2  1  3
  2  3  1
  2  3  2
  3  1  2
  3  1  3
  3  2  1
  3  2  3
normal end

長さ 100 の場合の条件を満たす最初の列

出力は 0,1,2 にしました。

  0  1  0  2  0  1  2  0  2  1  0  1  2  0  1  0  2  0  1  2
  0  2  1  0  2  0  1  0  2  1  0  1  2  0  1  0  2  0  1  2
  0  2  1  0  1  2  0  1  0  2  1  0  1  2  0  2  1  0  2  0
  1  0  2  1  0  1  2  0  1  0  2  0  1  2  0  2  1  0  1  2
  1  0  2  0  1  0  2  1  0  1  2  0  1  0  2  0  1  2  0  2

Systematic Programming: An Introduction

Systematic Programming: An Introduction

  • 作者:Wirth, Niklaus
  • 発売日: 1973/04/01
  • メディア: ハードカバー

プーチン氏 自己隔離や否や

www.asahi.com

プーチン氏が「自主隔離」 接触した医師がコロナ感染 新型肺炎コロナウイルス

モスクワ=石橋亮介 2020年4月1日 21時37分

 ロシア大統領府は1日、新型コロナウイルスの感染拡大防止のため、プーチン大統領が側近らとの接触を避け、会議などの業務を遠隔で行うと発表した。プーチン氏と面会した医師の感染が確認されたためで、ウイルスの検査も毎日行うという。

 感染が確認された医師はデニス・プロツェンコ氏。3月24日、モスクワ市内の病院に視察に来たプーチン氏らを案内していた。31日になってこの医師の感染が明らかになった。大統領府は同日には「大統領は定期的に検査を受けており、問題ない」とし、「自主隔離」には触れていなかった。

 タス通信によると、大統領府のペスコフ報道官は1日、「大統領は当面、(テレビ電話を使った)オンラインで執務する」と説明。「自分も含め、病院で大統領と一緒にいた人は毎日検査を行い、会議などを避けて可能な限り距離をとる」と強調した。(モスクワ=石橋亮介)

www.yomiuri.co.jp

視察し握手もした病院長感染、プーチン氏は隔離措置取らず 2020/04/01 22:42 新型コロナ

 【モスクワ=工藤武人】新型コロナウイルスの感染拡大を受け、ロシアのプーチン大統領が3月24日に視察した病院の院長が、ウイルスに感染していたことが31日、明らかになった。プーチン氏は、院長と素手で握手もしていたが、自主的な隔離措置は取らず、通常の公務を続けている。

 プーチン氏の視察先となったのは、新型コロナウイルス感染者の治療の最前線となっているモスクワ市内の病院だ。プーチン氏は防護服を着て治療現場を視察したが、院長との会談ではマスクを着けずに至近距離で話し込み、エレベーターにも同乗した。院長も白衣をはおっただけだった。視察には、ゴリコワ副首相やモスクワ市のソビャーニン市長も同行した。

 タス通信によると、ペスコフ大統領報道官は31日、「プーチン氏は定期的に検査を受けており、全て順調だ」と述べたが、4月1日にはプーチン氏が「握手を断っている」とし、感染予防に努めていることを強調した。プーチン氏は1日、政権幹部らとの会合を初めてテレビ会議に切り替えた。

【メモ帳】どうして PFLOTRAN は C++ より Modern Fortran を選んだのは何故?

PFLOTRAN

参考までに。

PFLOTRAN is written in object oriented, free formatted Fortran 2003. The choice of Fortran over C/C++ was based primarily on the need to enlist and preserve tight collaboration with experienced domain scientists, without which PFLOTRAN's sophisticated process models would not exist.

分野のおっさんたちは Fortran しかしゃべらないんだす。

www.pflotran.org

f:id:fortran66:20200402021156p:plain

p.20 The Origins Evolution and Applications of PFLOTRAN. (Conference) | OSTI.GOV https://www.osti.gov/biblio/1594643-origins-evolution-applications-pflotran

Modern Fortran: Building Efficient Parallel Applications

Modern Fortran: Building Efficient Parallel Applications

  • 作者:Curcic, Milan
  • 発売日: 2020/07/28
  • メディア: ペーパーバック

【寝言】 NHK 聴き逃し再放送に 音楽遊覧飛行 エキゾチック・クルーズ 

NHK-FM 放送の音楽遊覧飛行はサラーム海上のエキゾチック・クルーズだけが好きなんですが、録音してたまに聞いたりしています。エキゾチック・クルーズは、月の第一週だけなのでしょっちゅう忘れて聴き逃してしまいますが、まぁそのくらいがちょうどいい塩梅かも。わりと現代ものの民族音楽をやっていて、伝統民族音楽に囚われず、かといって若者ポップスばかりにもならずと結構面白い巻が多いです。

それが今月からネットの聴き逃し放送に対応するようで有難いです。今月のテーマはアフリカポップスでイマイチですw

www.nhk.or.jp

合唱の時間も聴き逃し放送に対応してくれるとうれしいです。

民族音楽などを紹介する人は、若いうちは聞きやすい/なじみやすい曲を流して呉れるのですが、年を取るにつれ煮詰まってきて、何が面白いのか分からないようなマイナー曲を延々と流し続けるようになるのが不思議です。最近では土曜の朝のピーター・バラカンさんもいい感じに煮詰まってくれて、全然面白くないのですが、結構お便りが送られているところを見ると需要はあるようです。

イスパニアFortran

f:id:fortran66:20200402024124j:plain

The Xavier Cugat Collection 1933-58 [Clean]

The Xavier Cugat Collection 1933-58 [Clean]

  • 発売日: 2013/07/08
  • メディア: MP3 ダウンロード


そよ風と私 (ザビア・クガート)


Patricia