fortran66のブログ

fortran について書きます。

【メモ帳】R ソース中 Fortran

R における Fortran ソースの割合 その2 以前 2011 年の何某氏の blog をソースとして、R ソース中の Fortran の割合を < 25% と見積もりましたが、最近また別の何某氏の bolg にてさらに別の何某氏が github 上のソース割合表示から Fortran 割合を 約 23%…

【ニュース】NAG Fortran が Apple Silicon 対応 他

NAG Fortran アップルの秘密主義のせいか噂も出ていませんでしたが、NAG が Apple silicon 用 Fortran を提供するようです。意外でした。NAG は文法に厳しい正しい言葉遣いが特徴で、実行時はそんなに速くないイメージですがどうなるでしょうか。ともかくめ…

【寝言】LGBTQ Fortran

Loli, Gachi-muchi, BBA, TKTT, オバQ gfortran は LBGT 止まりだが、 Intel Fortran は一歩進んでQも可。 Q Returns the number of characters remaining in an input record. とあるので、以下の print 文では全く意味を持たない。 fortran66.hatenablog.c…

【ニュース】 R 界隈からの Apple silicon での Fortran 情報

gfortran 以前、R 言語中で Fortran ソースの割合がそこそこあるというメモを書きましたが、R 開発者界隈の方から Apple Silicon での fortran 最新情報が出てきました。llvm, gnu とも公式対応は以前と変わらず芳しくないようですが、gfortran に関しては個…

【メモ帳】Fortran から Julia 呼び出し その4 (SymPy)

SymPy で記号演算した結果の Fortran からの利用 Fortran から Julia を呼んで、そこで PyCall + SymPy して記号処理します。出てきた結果を Julia の関数として定義して、それを Fortran から呼び出します。 実行結果 まず julia で普通に f(x) = x2 と定義…

【メモ帳】Fortran から Julia 呼び出し その3

Fortran から julia 呼んでそこから python 呼んでそこから Fortran 呼び出し Python から呼び出す Fortran subroutine sub.f90 subroutine test() bind(c, name = 'test') print *, 'subroutine test: Fortran calls Julia calls Python calls Fortran subr…

【メモ帳】Fortran から Julia 呼び出し その2

Julia from Fortran 前回の続きで、関数呼び出し、配列受け渡しなどの例題をやってみました。 一応、配列構造体を定義してみました。しかしコンパイル時の定数で構造が変わるようなので、よく分かりません。 しかし構造体などをきちんと定義していないので、…

【メモ帳】Fortran から Julia 呼び出し

Embedding Julia in Fortran スタバのソファーでふんぞり返って、タラタラとネットサーフィンwしていたら C から Julia を呼び出す本家のページにたどり着きましたので眺めやると、それほど面倒くさそうでもなさそうなので Fortran から呼び出すことに挑戦…

【メモ帳】道教の成立

武内義雄「支那思想史」第18章 道教の成立 そこで私は、道教は最初張陵から張魯に至る符呪祈祷の俗信から起って、次に葛玄の神仙養生術を取り入れ、又次に魏晋の際一世を風靡した老荘の哲学を牢籠し、最後に仏教の教理で修飾して出来上がった宗教だと結論…

【メモ帳】昔のネタ再訪

BBC Radio 4 FORTRAN 5年前の放送ですが podcast となって聞けていたようです。今も聞けます。 fortran66.hatenablog.com Using FORTRAN was absolute bliss! Codes that Changed the WorldBBC Radio 4ScienceUSD 0podcasts.apple.com 途中、動画の音声が流…

【寝言】鼻歌検索

Google Assistant の鼻歌検索 すごく昔から気になっていた曲があって、どの鼻歌検索でも引っかからなかったのに、Google Assistant が最近始めたらしい鼻歌検索してみたらすぐ出てきました。Google Assistant はウザくて嫌っていましたが、少し評価したいで…

【メモ帳】IBM llvm Fortran 他

alias 問題 背景として最近の流れでは、Fortran などのコンパイラで生成した中間形式 (IR) を、llvm のバックエンドで最適化をかけて、加速装置など個別のハードウェア構成に適合したバイナリを吐くようにしたいのではないかと思われます。 副プログラムの引…

【ニュース】Intel Fortran update

Intel® Parallel Studio XE Composer Edition 2020 update 4 が来ました。バグ修正のみのようです。 software.intel.com software.intel.com 大室家 (4) 特装版 (4) (百合姫コミックス)作者:なもり発売日: 2020/12/23メディア: コミック ゆるゆり (19) 特装…

【メモ帳】分割数の漸化式 オイラーの五角数版

オイラーの五角数を用いた分割数の漸化式 但し とする。 Julia の BigInt の試用目的で昔の計算をやり直してみます。 fortran66.hatenablog.com fortran66.hatenablog.com fortran66.hatenablog.com 整数の分割作者:ジョージ・アンドリュース,キムモ・エリク…

【メモ帳】第五元先という奇妙な名前

鄭玄の先生 秋も深まってきて、またぞろ武内義雄の支那思想史を読んでおりました。面白いので何度読んでも飽きません。戦前の学術の水準の高さよ。でも何度も読んでいると細かいところが気になってきます。 第十二章 後漢の経学 に、鄭玄ははじめ第五元先に…

【メモ帳】分割数の漸化式

約数和を使った漸化式 Julia 言語は型を指定しないで書くと、色々な型に自動で対応してくれるというので、多倍長を試してみようと分割数を計算してみることにします。 まず Fortran で以前のプログラムを少し改良しておきます。Fortran のプログラムを移植し…

【寝言】小惑星にハンコ!

米宇宙探査機が小惑星ベンヌに押印! ハンコは古代シュメール文明から脈々と続く伝統なので大事にしたいものです。 米国の宇宙探査機が小惑星ベンヌに押印する映像を公開しましたが、ハンコを大事にしなさいという宇宙からのメッセージを読み取りました! yo…

【メモ帳】OpenMP メモ帳

ISUS セミナー動画とスライドあり。 www.isus.jp intel dev cloud OpenMP offload basics ログイン後 root 下の /data/oneapi_workshop に Jupyter Notebook 形式の Fortran/C 例題資料あり。 /data/oneapi_workshop/OMP_Offload.tar.gz software.intel.com …

【寝言】今期秋アニメ

今期の秋アニメ 今期の秋アニメで気に入ったのは「無能なナナ」です。ゆるゆりの腹黒ちなつちゃんが成長した姿かなと思いながら見ています。 予備知識なしで見たので楽しめてます。 新作ごちうさは、夏に放送予定だったのでしょうか?暑い暑い言っていて、違…

【寝言】米人 10・18→10^18 exascale day

謎の秋祭り 米国で10月18日を エクサ 1018 にちなんで exascale day として催しをやっております。 たぶん邪推するに Intel の Aurora が exaflops に達する予定だったのでそれに合わせたものだったのでしょう。出前は出ましたが、まだ到着しておりませ…

【メモ帳】Chapel 1.23 出る

Chapel 1.23 先日、Chapel の予感がして覗きに行きましたが、なんと同日アメリカ時間に 1.23 が上げられていたようです。早すぎたんだ・・・腐ってやがる。 fortran66.hatenablog.com コンパイラのコンパイル 単一機で multi locale のプログラムを試すには、qu…

【寝言】田中美智太郎も批判していた学術会議のありさま

暗黒世界 学術会議は共産党が私物化していることが問題なのに、誰も共産党の名すら出さない恐ろしさ。まさに封建的村社会ww アメリカ禁酒法時代のギャングのアル・カポネが役人・裁判官・マスコミを全部買収していて、復讐が恐ろしくて誰も証人として証言…

【メモ帳】Cray Chapel が Hewlett Packard Chapel に?

Cray Chapel が Hewlett Packard Chapel に? だんだん日暮れも早くなってきて、そろそろチャペルの更新時期かとのぞきに行ったところ、HP による Cray 買収の影響で、Brad Chamberlain のスライドが Hewlett Packard 印に変わっていました。買収の影響で今…

【メモ帳】flang 11 を試すつもりが 12 に   その他

flang f18 llvm 11 を試すつもりで llvm 11 一式のバイナリを落としてきましたが、flang が入っておらず、仕方が無いのでソースから build しようと思ったら llvm がらみのあれこれが足りぬと叱られるので、llvm 一式ソースから作るかと指示通りに git clone…

【ニュース】flang 11.0 出るが、まだフロントエンドのみ

ララベル速報 www.phoronix.com www.reddit.com j0mpz 1 day ago Is this still only the frontend (parser), or is LLVM code generation already available? zolamusic 1 day ago Only the frontend in this time. They call the gfortran by default. 【…

【メモ帳】LFortran 0.9  C++ の XEUS で Jupyter

lfortran LFortran が pre-alpha ながら Ver 0.9 くらいになったようです。C++ による書き換えが進んでいて、jupyter 対応も c++ 独自カーネルの XEUS にかわっているようです。機能としては、まだ組み込み関数も配列も使えないので、なかなか意味のあるプロ…

【寝言】スチーブン・バノン氏

無意識の解放! ブライトバートニュースのコアコンピタンス(他社にない自社の強み)は「ニュース」ではなく、「コメント欄」なのだとバノンは言う。「コメント欄」に人が集まり、コミュニティーを作る。そこから何かが起き、もう一人の自分たちがコミュニテ…

【メモ帳】英計算機協会 Fortran 分科会合 その他

BCS 2020 毎年、面白い内容が多いですが、今年はどうでしょうか。コロナでオンライン会合だった模様です。 fortran.bcs.org スライドあり。 Fortran Joint meeting with Computational Physics Group of Institute of Physics Introduction to Programming w…

【寝言】安倍前首相に関する foreignaffairs 記事

安倍前首相 www.foreignaffairs.com www.foreignaffairsj.co.jp 日本語版の出鱈目な表題と要旨はどうにかならんのかw 朝日の出版で中国人が編集ではしょうがないか。 それにつけても後任の菅首相も有能すぐる。六本木は乃木坂のアカの伏魔殿を燻し始めてワロタ…

【メモ帳】エラトステネスの篩

Sieve of Eratosthenes julia と fortran でおおむね同じなプログラムで比較。10**9 以下の素数をエラトステネスの篩で求める。 Julia 1.5.1 julia さん速くて Fortran -O3 と揺らぎの範囲で大差ない。 @time pr = prim(10^9) 33.018580 seconds (5 allocati…