fortran66のブログ

fortran について書きます。

【ニュース】Julia ヨイショ記事と思いきや、珍しく Fortran 褒められる

意外に公正で草 どうせ julia 上げ記事かと思ってみてませんでしたが、覗いてみると意外に Fortran も上げられてました。 お金が絡むとシビアです。 towardsdatascience.com 要旨 During the last months, I have been working intensively on quantitative …

pointer function で連想配列 改

以前、Fortran2008 の pointer function を用いて、簡単な連想配列を作ってみましたが、ハッシュ法を N.Wirth の『アルゴリズム+データ構造=プログラム』 中の方法を用いてみました。キーは散らばっているでしょうか? 関数名を配列名のように使うので、こ…

【ニュース】金のトイレ盗まる!

トランプ大統領に貸与の申し出があった 18 金トイレ盗まる! 英美術館で使用可能状態で展示中に盗まれたようです。 www.kitco.com America According to media reports, the working toilet was reported stolen from the Blenheim Palace in Oxfordshire, E…

【ニュース】「人気の高いプログラミング言語2019」

Fortran さんじゅうはち位 数日前の記事ですが、メモ帳代わりに。 IEEE Spectrum、「The Top Programming Languages 2019」を発表 news.infoseek.co.jp 最後に紹介しているのは、38位にランクインしている「Fortran」。科学計算言語としての莫大な遺産により…

【メモ帳】暇つぶしに Julia で Camplbell-Baker-Hausdorff 計算してみる

Camplbell-Baker-Hausdorff の公式 Camplbell-Baker-Hausdorff の公式という割に、実際の展開式を導いたのはソ連の Dynkin らしいので、この名でいいのかとも思いますが、アカと間違えられて FBI に通報されても困るので、長い物には巻かれましょうw exp(A)…

【メモ帳】Fortran 2018 C言語の相互運用 MFE §21.6 その2

MFE §21.6.5 Allocatable objects C 側で Fortran の allocatable 変数を deallocate/allocate します。 これを使えば、ZMQ の zmq_msg_init_data も出来ると思います。 fortran66.hatenablog.com 実行結果 gfortran-9/gcc-9 hp8@HP8:~/f2018$ gfortran-9 mf…

【メモ帳】○○的の「的」について

杉山孚富 著「漢文故事熟語詳解」より 明治36年 (1903年)初版 明治38年(1905年)訂正6版 面倒なので当用漢字でw カタカナもひらがなに。その他適当に濁点づけ。 てノ部 的 「近来人好んで的の字を用ゆ何々的など云へり此は元来助字にして意味あるに…

【寝言】ペンス副大統領の反中演説ようやく!?

ペンス演説第二段 当初、天安門事件30周年の日、六月四日にする予定で、その後延び延びになっていた、マイク・ペンス米副大統領の反中演説が、十月一日の中共の建国記念日に向けて行われるようです。 www.youtube.com 去年十月にもハドソン研究所で画期的…

【メモ帳】Fortran 2018 C言語の相互運用 MFE §21.6

MFE §21.6 Accessing Fortran objects ~21.6.4 実行結果 gfortran-9 hp8@HP8:~/f2018$ gfortran-9 c_mfe216.c mfe216.f90 hp8@HP8:~/f2018$ ./a.out GCC version 9.1.0 -mtune=generic -march=x86-64 x contiguous associated ip => null() discontiguous di…

【メモ帳】Fortran 2018 C言語の相互運用 MFE §21.5

Modern Fortran Explained - 2018 §21.5 C descriptors 暑いので適当w interface に type(*) :: x(..) を使ったのが新奇。 attribute がうまく取れていないのが謎。 実行結果 gfortran-9/gcc-9 hp8@HP8:~/f2018$ gfortran-9 c_mfe215b.c mfe215b.f90 hp8@HP…

【寝言】覗き魔 google 通販履歴貯めてるw その他

Googleの「自分のデータをダウンロード」 お盆休みに、google の「自分のデータをダウンロード」を、ほんの一部だけ見てみたのですが、通販での買い物履歴が json データになってたまっていてムカつきました。amazon などは別メアドなのに、スパムフィルター…

【メモ帳】Fortran 2018 C言語の相互運用 MFE §21.4

Modern Fortran Explained - 2018 §21.4 assumed character length Fortran の文字列は文字の配列ではなく、文字列長などの付加的な情報を持っています。その情報は、C言語側で Fortran 2018 が用意したヘッダーファイルを include することで、構造体の要素…

【寝言】「支那」という語について

宋代の『翻訳名義集』に現れた支那 国会図書館デジタルライブラリ『翻訳名義集 : 支那撰述. 巻第3』17 コマ dl.ndl.go.jp 国立国会図書館オンライン | National Diet Library Online 「脂那」こっちの字面の方が、下水油で脂ぎった支那料理のようで、お似合…

【メモ帳】Fortran 2018 C言語の相互運用 MFE §21.2

Modern Fortran Explained - 2018 Fortran2018 で導入された Fortran と C の interoperability について、少しづつ試してゆきたいと思います。 MFE にはコード断片しかないので、残りを適当に埋めて、実行可能な形にして試していみることにします。 § 21.2 …

【メモ帳】メモリー階層レイテンシー

intel parallel universe magazine vol.37 p.7 より software.intel.com latency 1nsを1秒に換算すると、おおよそ 階層 換算時間 レジスタ 0.1秒 キャッシュ 1秒 DRAM 1分 不揮発メモリ 1時間 SSD 1日 HDD 100日 テープ 3年 分散メモリのネットワーク通信…

【寝言】【メモ帳】樗蒲の発明は老子

樗蒲(ちょほ) 博物志に『樗蒲は天竺より出づ、老子此に本づきて今の雙陸を造る』とか。 ja.wikipedia.org 『芸文類聚』巻74に引く馬融「樗蒲賦」では老子が西域に行ったときに暇潰しに樗蒲を遊んだとしている[1]。『太平御覧』巻754に引く西晋の張華『博物…

【寝言】ハドソン研究所創始者ハーマン・カーンと電子計算機

Herman Kahn 昨年のマイク・ペンス米副大統領の反中宣言会場に選ばれ、日高義樹でおなじみの米ワシントンのハドソン研究所ですが、最近ではトランプ大統領と不仲で辞職したマクマスター補佐官が、なぜか新設の日本部長になったとニュースに出ていました。ま…

【メモ帳】薔薇という名のコンパイラ族

ROSE compiler infrastructure rosecompiler.org ROSE is an open source compiler infrastructure to build source-to-source program transformation and analysis tools for large-scale C (C89 and C98), C++ (C++98 and C++11), UPC, Fortran (77/95/20…

【メモ帳】精度毎のルーチンで楽したい その2

submodule を使った場合(案) まだ何がいいかよく分かりません。 Modern Fortran Explained: Incorporating Fortran 2018 (Numerical Mathematics and Scientific Computation)作者: Michael Metcalf,John Reid,Malcolm Cohen出版社/メーカー: Oxford Univ …

【メモ帳】精度毎のルーチンで楽したい

generic の代用 薬局でもジェネリックを勧められても拒否するし、Fortran 文法でもジェネリクス導入を拒否しますw そもそも単精度では精度を保つため、積の結果だけ倍精度使うとか、完全に同じソースではすまないことがあります。 とはいえ、ジェネリクスの…

【メモ帳】gfortran-9 gif

binder gfortran-9 animation gif gfortran-9 は forall(integer::i ...) のローカル変数には対応していないが、findloc には対応していてやや進歩。 定数に関して、コンパイラオプションが必要 %fcflags:-fno-range-check f66blog.github.io ディレクトリ e…

【メモ帳】Backus の質疑応答から

1978年 "The History of Fortran I, II, and III" 講演の質疑応答より この質疑応答が中々いい所を突いていて面白いのです。 J.A.N. Lee. Transcript of question and answer session: "The History of Fortran I, II, and III". www.softwarepreservation.o…

【メモ帳】雑記

Backus 記憶違い www.semanticscholar.org Backus が 1950年代のコンピュータ事情を回顧していますが、手紙等を調べた結果、色々記憶違いがあって自分がそれまで公言していたことに誤りがあったと正直に述べています。 ワシントン (子どもの伝記全集 8)作者:…

【メモ帳】Julia 言語で BMP write

Julia 言語よく分からんw 型変換規則がよく分からない。構造体をまとめて書きだす方法も分からない。構造体に初期値を与えられるのかも分からない。ASCII 文字型も分からない。 頭部に定義部が無いと軽い感じになるのだなと思いました。xxx...end の block …

【メモ帳】mybinder 用 gfortran-9 Dockerfile

binder で gfortran-9 gfortran-9 は、まだ add-apt-repository ppa:ubuntu-toolchain-r/test しなければならないので、binder 備え付けの apt.txt では対応できず、Dockerfile を書く必要があります。しかし mybinder サイトの説明を読んでも幼女には難しい…

【メモ帳】分割数 p(721) チェック

分割数 p(721) ハーディーの「ラマヌジャン その生涯と業績に想起された主題による十二の講義 」に 721 の分割数計算の級数和が出ていたので(講義VIII, p.193)、検算のため比較してみました。 ラマヌジャン その生涯と業績に想起された主題による十二の講…

【寝言】支那人大暴れw

海外で中共派対香港派で争いw 最近ネットは朝鮮ニュースばかりで面白くないですが、一段落したのか、ようやく支那ネタも復活してきました。 阿片窟から這い出てきたような支那人がオースラリア・ニュージーランドで大暴れのようです。 youtu.be 傑作なのは…

【メモ帳】jupyter lab で IJulia 等の syntax coloring がされない件

jupyter lab での syntax coloring notebook で開いて色がつくのであれば、これを save すれば、lab の方でも色が残る。 IJulia, Wolfram language 等。

【メモ帳】Julia で分割数

分割数の直接計算 整数の分割数を Hardy-Ramanujan-Rademacher の公式で計算して見ます。 Fortran, Wolfram-Language 等でやったものの Julia 版です。 fortran66.hatenablog.com fortran66.hatenablog.com fortran66.hatenablog.com fortran66.hatenablog.c…

【ニュース】他言語の風の便り

Go 言語 例外処理の Try を拒否 www.theregister.co.uk N.Wirth の Oberon 言語の末裔として筋を通したのでしょう。 Fortran でも、例外構文とテンプレートを入れろという圧力が高まっていますが、バグるうえにもっさりしそうなので拒否してゆくべき?w Jul…