至点 今日は夏至でしたが、珍しく晴れ間が広がって太陽が輝いていました。 普段は梅雨時で太陽を見ないことが多いのでありがたいことです。 これも天皇陛下の神徳か！ ksk 主義ｗ 加速主義もイタリア未来派の再来を見るようで既視感がある。右派左派両方が沸…
Intel Fortran LLVM 版 intel は Fortran compiler のバックエンドを llvm のものに置き換える構想で、来年の頭には製品版を出すと抱負を述べていましたが、開発チーム募集の求人が出ていました。この調子だと気長に待つことになりそうです。 求人場所が東海…
亂世之英雄 Fortran なんか今年に入って集計のアルゴリズムが変わったのか、Fortran の上昇が止まりません。Swift 越えまであと一歩ｗ amazon.co.jp も集計対象らしいので、古本屋の皆さんは図書館除籍のカビの生えた FORTRAN 本を大量に１円出品してくださ…
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…
gfortran 7 + gcc 7 の組み合わせでは block に問題あり WSL Ububtu 18.04 で我が xplot をコンパイルしようとしたところ、 block..end block 構文が複数回出てくると、c コンパイラが文句を垂れてコンパイルできないことがわかりました。 新版のコンパイラ…
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 に concrete な routine 以前、下記の様な記事を書きましたが、この手法はけっこう適用範囲が広いのではないかと思えてきましたので、メモっておきます。 fortran66.hatenablog.com はじめは abstract な interface に concrete な routine を…
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…
Monte Carlo で素朴に円周率を計算してみます 1000 点での見積もりなので二桁ほどの精度もありません。 でも点数を増やすと図がつぶれてしまいますｗ 実行結果 Monte Carlo: estimated pi = 3.19600010 ⡟⡛⠛⡛⠛⠛⠛⢿⠭⠭⣩⣉⣙⢉⠉⠙⠉⡩⠉⢉⠋⠉⠍⠉⠍⠝⠉⠉⠉⠉⠉⠉⠉⣙⠝⠉⠉⠍⠉⠉⠉⢙⠉⠩⠉⡉⠉⠉⠉⢹⠀…
月食 金星 水星？ 男色家阿部さんの次は月食か。 天気が悪くて見えないかと思っていましたが、夕刻にかけて少し晴れ間が出ました。ふと西の空を見ると、金星がまさに沈もうとするところで、少し上の方に水星らしき星もかすかに見えていました。金星と水星が…
unicode 表示 Apple の terminal は unicode 文字を表示してくれるので、unicode の点字文字を使ったグラフィックスに挑戦してみることにします。ここで点字文字は 4 行 2 列の点で出来ています。 Windows の DOS 窓や WSL コンソールでは unicode 文字がう…
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…
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年刊行の Fortran techniques with special reference to non-numerical applications に当時の事情などを見て、色々感想が書かれています。 scruss.com internet archive に上げていますが、あそこも著作権無視の無法サイトだとい…
Tiobe 今月も 20 位 www.tiobe.com ZDnet 記事 www.zdnet.com まず最初に目につくのは、三枚モニターに同じ画面を出しているおじさんの写真が使われている事です。この三枚同画面おじさんは有名なフォトストック写真で、以前どこかで物笑いの種にされていた…
講演動画 www.youtube.com 前半は面白くない。後半は面白い。 夜の街をスケッチするのに、iPad だと明かりが要らないので便利というのが、デジタルデバイス固有の利点という事らしくて興味深いです。 インタビュー動画 www.youtube.com www.youtube.com チャ…
swift JCL 説 swift JCL 説に則り pipeline は固定のまま、入力バッファも固定のまま、引数指定のみを繰り返して、２回平方根をとります。 import Metal @_cdecl("mymetal") public func mymetaln(n: Int32, px: UnsafeMutablePointer<Float32>) { let device = MTLCr</float32>…
ワタミ社長の論語 写真を拡大してみると、和民社長の渡邉美樹氏の愛読論語は、講談社文庫版の木村英一編注の論語のようでした。今は絶版で私は見たことないのですが、少し標準的でない解釈の論語で一部で評価が高いようです。木村英一のハードカバー「孔子と…
Humble Bundle アジア人を標的とした犯罪防止に寄付されるようですが、斯様な不健全なマンガはむしろ助長しそうなｗ 今の所、２千セットで４万ドルばかり売り上げています。 最近は、BL 一辺倒でやおいとかジュネ系とか言わなくなりましたね。 www.humblebun…
Fortran の配列を Swift 経由で GPU に送り計算する Swift の文法もよくわからないまま、GPU へ突撃。GPU 用の機械語を吐き出すコンパイラの言語は C++14 になっています。 プログラム 以下では GPU 用のプログラム（拡張子 .metal）、 Swift 用のプログラム…
関数・配列 関数の戻り値はスタックに積まれて普通に C 言語風にとれるようです。 Fortran の配列が Swift に渡せれば、いずれ Swift で GPU 計算させて Fortran に戻せるだろうという目論見です。FORTRAN 77 式の整合配列を渡すノリで C 風の生ポインタとし…
gfortran から簡単な swift routine 呼び出し 参考：gist.github.com Swift 側を dll にして呼べば良い様です、 プログラム swift @_cdecl("say") public func say() { print("Hey!") } fortran 短いのでインタフェースはメインプログラムに書いてもよかった…
