fortran66のブログ

fortran について書きます。

【寝言】夏至

至点 今日は夏至でしたが、珍しく晴れ間が広がって太陽が輝いていました。 普段は梅雨時で太陽を見ないことが多いのでありがたいことです。 これも天皇陛下の神徳か! ksk 主義w 加速主義もイタリア未来派の再来を見るようで既視感がある。右派左派両方が沸…

【ニュース】Intel さん 今から LLVM 求人

Intel Fortran LLVM 版 intel は Fortran compiler のバックエンドを llvm のものに置き換える構想で、来年の頭には製品版を出すと抱負を述べていましたが、開発チーム募集の求人が出ていました。この調子だと気長に待つことになりそうです。 求人場所が東海…

【ニュース】TIOBE 月旦評6月号で Fortran 17位

亂世之英雄 Fortran なんか今年に入って集計のアルゴリズムが変わったのか、Fortran の上昇が止まりません。Swift 越えまであと一歩w amazon.co.jp も集計対象らしいので、古本屋の皆さんは図書館除籍のカビの生えた FORTRAN 本を大量に1円出品してくださ…

【メモ帳】M1 Mac で fpm

arm 版 gfortran での場合 公式のバイナリ(α版)は intel 版で、arm 版の gfortran しかないとランタイムの問題で動きません。 github.com github.com しかし、Fortran のソースがあるので、それをコンパイルすれば大丈夫です。 メイクもスクリプトも要らず、…

【寝言】レ・ジーン

合唱団 Les Djinns 昔、六本木 WAVE とか CD 屋で探しても見つからなかったのですが、今は検索ですぐ出てきますね。 ISIS の女神の加護により CD になっていたようです。 www.youtube.com Cole Porter www.youtube.com Du Music-Hall aux Années 60 (1950-19…

【メモ帳】WSLgで図を出す

gfortran 7 + gcc 7 の組み合わせでは block に問題あり WSL Ububtu 18.04 で我が xplot をコンパイルしようとしたところ、 block..end block 構文が複数回出てくると、c コンパイラが文句を垂れてコンパイルできないことがわかりました。 新版のコンパイラ…

【メモ帳】fpm で intel fortran

fpm で起動されるコンパイラを intel fortran に fpm のプロトタイプは gfortran 専用でしたが、いつの間にかコンパイラを指定できるようになっていました。コマンドラインのヘルプを出したら書いてありました。 fpm build --compiler ifort fpm run --compi…

【メモ帳】ブラウザからも阿部さん

canvas 形式での出力 postscript や HTML 形式での出力にはファイル名が必要になるので、インターフェースの改良が必要だがシンプルさが減る。昔作ったインターフェースでは考えてあった。 fortran66.hatenablog.com Modern Fortran: Style and Usage作者:Cl…

【メモ帳】abstract type の用い方について

abstract type に concrete な routine 以前、下記の様な記事を書きましたが、この手法はけっこう適用範囲が広いのではないかと思えてきましたので、メモっておきます。 fortran66.hatenablog.com はじめは abstract な interface に concrete な routine を…

【メモ帳】fpm に X11 用の c ファイル入れる

fpm に c プログラム入れて X window に図 fpm を利用して unicode plot と X-window で同じ図を出すことにします。 fortran66.hatenablog.com そのためには、fpm に c プログラムを認識させ、X11 のライブラリとリンクさせる必要があります。fpm 公式 gihub…

【メモ帳】点字プロットで作図

Logistic map & Lorenz attractor 昔のインターフェースを流用して、作図してみます。 fortran66.hatenablog.com fortran66.hatenablog.com 横に縞が入るのが、昔の電送写真のようでいいですね。 プログラム device.f90 module device implicit none type, a…

【メモ帳】点字 plot で pi

Monte Carlo で素朴に円周率を計算してみます 1000 点での見積もりなので二桁ほどの精度もありません。 でも点数を増やすと図がつぶれてしまいますw 実行結果 Monte Carlo: estimated pi = 3.19600010 ⡟⡛⠛⡛⠛⠛⠛⢿⠭⠭⣩⣉⣙⢉⠉⠙⠉⡩⠉⢉⠋⠉⠍⠉⠍⠝⠉⠉⠉⠉⠉⠉⠉⣙⠝⠉⠉⠍⠉⠉⠉⢙⠉⠩⠉⡉⠉⠉⠉⢹⠀…

【寝言】金星 水星 月食

月食 金星 水星? 男色家阿部さんの次は月食か。 天気が悪くて見えないかと思っていましたが、夕刻にかけて少し晴れ間が出ました。ふと西の空を見ると、金星がまさに沈もうとするところで、少し上の方に水星らしき星もかすかに見えていました。金星と水星が…

【メモ帳】unicode 点字で阿部さん

unicode 表示 Apple の terminal は unicode 文字を表示してくれるので、unicode の点字文字を使ったグラフィックスに挑戦してみることにします。ここで点字文字は 4 行 2 列の点で出来ています。 Windows の DOS 窓や WSL コンソールでは unicode 文字がう…

【メモ帳】intel fortran 次世代構想 roadmap

intel compiler 講演 intel の Fortran/C++ compiler は backend を llvm ベースに変えつつありますが、その移行期についての解説がなされています。 旧来の intel backend 版は intel fortran classic (ifort)となり、次世代 llvm backend 版が intel fortr…

【メモ帳】分割数計算 改訂

分割数 Hardy-Ramanujan-Rademacher の公式の exp の複素数和の所を、どうせ虚部は x 軸に対して対称の位置に出て打ち消しあうので、cos に置き換えて実数のみで計算します。 Fortran と Julia の以前のプログラムを改訂します。 Hardy-Ramanujan-Rademacher…

【メモ帳】補遺

ZDNet 記事の翻訳 先日紹介した記事の日本語訳が出て、少し話題になっておりました。 fortran66.hatenablog.com japan.zdnet.com Pyodide 以前、SciPy など Fortran 部分が WebAssembly 化できなくて困っていたのが一応解決したようです。 CBLAS/CLAPACK な…

【寝言】夕空に月水金!

三日月 水星 金星 夕空に三日月が出ているのを眺めていると、地平線近くに金星も見えていました。水星の方が高度は高かったですが、金星の方がはるかに見やすかったです。光度は水星が 0.1 等、金星が -3.9 等程のようです。 May, West earthsky.org Mae Wes…

【メモ帳】MATLAB の Moler の面白講演

Star Trek 映画に登場! LINPACK, EISPACK や MATLAB で有名な Moler の講演が結構面白いです。 www.youtube.com (以前 EISPACK の FORTRAN への移植についての Moler の講演を紹介しました。) fortran66.hatenablog.com MATLAB は元々は Moler によって F…

【寝言】水星が東方最大離角近し

夕空に水星 19時も大分過ぎると空も暗くなってきて肉眼でもよく見えました。 www.nao.ac.jp スマホのデジタルズームなのでガバガバです。 美少女戦士セーラームーン コスプレ水野亜美風 男性用 Lサイズメディア: 男性用 L サイズ・・・ MR. バッド・ガイ(ス…

【メモ】津波

自宅で津波に遭われた方々の動画 いずれも、足腰の悪い家族の方やペットなどと自宅に残られていて津波に遭遇しながらも、なんとか無事に助かった方々です。助かって本当に良かったです。 youtu.be youtu.be 現実とは思えない、家が船みたいになってる夢とか…

【メモ帳】1972年刊本に見る当時の事情

よそ様のページリンク 1972年刊行の Fortran techniques with special reference to non-numerical applications に当時の事情などを見て、色々感想が書かれています。 scruss.com internet archive に上げていますが、あそこも著作権無視の無法サイトだとい…

【ニュース】Fortran ブーム続くw

Tiobe 今月も 20 位 www.tiobe.com ZDnet 記事 www.zdnet.com まず最初に目につくのは、三枚モニターに同じ画面を出しているおじさんの写真が使われている事です。この三枚同画面おじさんは有名なフォトストック写真で、以前どこかで物笑いの種にされていた…

【メモ帳】アンドレイ・スルノフ

講演動画 www.youtube.com 前半は面白くない。後半は面白い。 夜の街をスケッチするのに、iPad だと明かりが要らないので便利というのが、デジタルデバイス固有の利点という事らしくて興味深いです。 インタビュー動画 www.youtube.com www.youtube.com チャ…

【メモ帳】Fortran から Swift を呼ぶ その4

swift JCL 説 swift JCL 説に則り pipeline は固定のまま、入力バッファも固定のまま、引数指定のみを繰り返して、2回平方根をとります。 import Metal @_cdecl("mymetal") public func mymetaln(n: Int32, px: UnsafeMutablePointer<Float32>) { let device = MTLCr</float32>…

【ネタ】わたみんの論語は木村英一

ワタミ社長の論語 写真を拡大してみると、和民社長の渡邉美樹氏の愛読論語は、講談社文庫版の木村英一編注の論語のようでした。今は絶版で私は見たことないのですが、少し標準的でない解釈の論語で一部で評価が高いようです。木村英一のハードカバー「孔子と…

【乞食速報】Humble Bundle に、やおいマンガ登場!

Humble Bundle アジア人を標的とした犯罪防止に寄付されるようですが、斯様な不健全なマンガはむしろ助長しそうなw 今の所、2千セットで4万ドルばかり売り上げています。 最近は、BL 一辺倒でやおいとかジュネ系とか言わなくなりましたね。 www.humblebun…

【メモ帳】Fortran から Swift を呼ぶ その3

Fortran の配列を Swift 経由で GPU に送り計算する Swift の文法もよくわからないまま、GPU へ突撃。GPU 用の機械語を吐き出すコンパイラの言語は C++14 になっています。 プログラム 以下では GPU 用のプログラム(拡張子 .metal)、 Swift 用のプログラム…

【メモ帳】Fortran から Swift を呼ぶ その2

関数・配列 関数の戻り値はスタックに積まれて普通に C 言語風にとれるようです。 Fortran の配列が Swift に渡せれば、いずれ Swift で GPU 計算させて Fortran に戻せるだろうという目論見です。FORTRAN 77 式の整合配列を渡すノリで C 風の生ポインタとし…

【メモ帳】Fortran から swift を呼ぶ

gfortran から簡単な swift routine 呼び出し 参考:gist.github.com Swift 側を dll にして呼べば良い様です、 プログラム swift @_cdecl("say") public func say() { print("Hey!") } fortran 短いのでインタフェースはメインプログラムに書いてもよかった…