秋分 23 日 Equinox を迎える喜び。国民の休日になっているのももっとも。χ を崇める。 DevCloud intel の 無料試用 cloud server で fortran が利用できますが、python 用の kernel しか用意されていません。 昔、MS Azure 用に試みた手法で fortran 用の k…

白ロシアのルカシェンコ大統領のジョークセンス twitter.com euroweeklynews.com ウクライナ侵攻への協力を頼まれて「ロシア軍は強いから助太刀は不要」と断るセンスもいいし、それ以前も中東方面で経済移民目的の「難民」を飛行機で集めて EU 圏に送り込む…

Gustafson 語る Gustafson によると RISC-V 用に新浮動小数点数フォーマットが開発されているとのこと。ヨーロッパでも IEEE754 以外のフォーマットが作られているとか。 また昔のように、ベンダー毎に数値フォーマットが入り乱れて、互いに罵り合って欲しい…

FP8 短浮動小数点数もここまで来たかｗ 6 bit まで縮めて昔の 6 bit コンピュータ 3+3 の octal 時代へ帰ろう. We're pleased to support FP8 alongside @Intel and @NVIDIA to drive forward AI development. This allows developers to leverage existing …

BBC がソ連の原子力電池のレポート 以前、北方領土にある灯台ビーコンのエネルギー源として使われているソ連の原子力電池について触れましたが、英ＢＢＣが何故かそれがらみのビデオレポートを出してきました。 fortran66.hatenablog.com The nuclear lighth…

Fortran 遊び場！感度三千倍 play.fortran-lang.org 今年の Google Summer of Code の成果のようです。 fortran-lang.discourse.group Go 言語の Playground のように例題講座もついているようです。 中のコンパイラは gfortran 11.2 のようです。 Modern Fo…

Ternary computer ソ連で三進法のコンピュータが作られていたようです。名前は Setun 。 en.wikipedia.org en.wikipedia.org 三進法だと補数表現をとっても、負の絶対値最大値と正の最大値が二進法のようにずれることが無くなっていいかもしれません。 数値…

並列Fortranに関するシンポジウムの（第6回） 発表資料見られます。 sites.google.com 依頼講演動画は苦労を厭わず字幕を入れられたようです。 The State of Fortran 2021（ 日本語字幕入りビデオ上映 45分）Laurence Kedward, University of Bristol, UK LF…

Seminar: Generic Programming in Fortran 202Y www.youtube.com あまり例が整理されておらず、やや分かり難い点もありますが、まぁ勉強になります。 C++ 等で問題になった部分をうまく避けるようにデザインしているようです。 マクロ的なアプローチは避けて…

父親を狙ったテロで死亡か？ www.jiji.com www.bbc.com ロシアのウラジーミル・プーチン大統領の世界観に大きく影響したとされる国家主義思想家の娘が20日夜、モスクワ近郊で車が爆発したため死亡したという。ロシアの複数メディア報道を、ロシア連邦捜査委…

わたみんの論語 以前、和民社長の渡邊美樹氏の愛読論語が、絶版になっている木村英一編の講談社文庫版のものであること書きましたが、その後復刊がなったようです。現在電子版の半額セールが行われています。 fortran66.hatenablog.com お尻の方の解説を斜め…

Absoft 廃業 かつては、Cray Fortran とのフロントエンド共通を謳い、その後は Mac 方面に力を入れていた Absoft でしたが、昨今のコンパイラ無料化の流れを受け廃業することになったようです。 無料の gfortran に加え、intel fortran も無償化されたことで…

多倍長計算ルーチン GNU MPFR の Fortran binding GNU MPFR Library github.com Fortran 次期規格の Generics について The State of Fortran Genericseverythingfunctional.wordpress.com 考えているのは単なる型引数でしょうかね。 ぷにるはかわいい…

計算法 以前メビウス関数の記事を書きましたが、Fortran での計算法を別のやり方でやってみます。 fortran66.hatenablog.com ここでは素数の因数ごとに符号反転、素数の二乗以上のべきが混じるものは 0 にするという、愚直な方法でまず一覧表を作ってそれを…

(1/4)! Youtube のおすすめで (1/4)! を計算する動画が出てきたのですが、 (1/4)! は に等しいので見てみたところレムニスケート周率で表わされることが学べたので、メモっておくことにします。 www.youtube.com exp(-x4) の積分 で置換すると なので、 が示…

48 bit への回帰 昔の CDC の実数は 48 bit 長だったと思います。 www.hpcwire.com IBM も 360 シリーズに統一される前は、現在の 8 bit = 1byte ではなく、1 文字 6 bit で科学技術用は単精度浮動小数点数が 1 語 36 bit 2 進型で精度が今より良くなってい…

黙祷 これほどの人物の命を、大切に守ってゆこうとしない社会にショックです。残念でなりません。 妬み僻み嫉みで足の引っ張り合いしかしないサヨク・メディアが憎しみを煽ったことが巡り巡った結果だと思います。 疫病、戦争、暗殺、次は飢餓でしょうか？グ…

COBOL/Fortran が Latin 語扱いｗ www.nytimes.com ウクライナ侵攻後の Yandex に関する記事 Nginx はどうなのか知りたい。 www.nytimes.com BL エロ同人界隈 イギリス政界の中心地、ロンドン・ウェストミンスターにある会員制クラブ「カールトン」での出来…

Wendy Carlos demonstrates her Moog Synthesizer in 1970 昔、Walter Carlos がいつの間にか Wendy Carlos に性転換していてびっくりしましたが、1970 年の時点で十分カマっぽいです。 www.youtube.com 89 年には限界突破していたようです。 www.youtube.co…

ガンマ関数 Stirling の公式 をまた考えてみます。 で、また であるので、 を考えることにします。 この積分の被積分関数は が大きくなると増大する とそれを減衰させて抑え込む の積になっていて、 の増加につれてまず増えた後に尻すぼみに減る山のような形…

Programming in America in the 1950s – Some Personal Impressions 同名の記事について、以前記事にしていました。 fortran66.hatenablog.com This lecture’s transcript was included in the edited volume from the conference, viz. Backus, J., “Progra…

新しい月 六月三十日の夕方、夕焼け雲が綺麗なので写真を撮ったのですが、後から拡大して見てみると月齢 1 台の月が写っていました。ネットの月齢サイトで調べてみると月齢 1.3~1.4 程度と思われます。 こんなに細い月を見るのはしばらくぶりです。うれしい…

Moonglow / thema from picnic picnic の映画は見てないんですが曲が好きで。メドレーと書かれたりもしますが、今風に言うとマッシュアップのようになっていて、同時に二つの曲の旋律が流れるところが良かったりします。 この前、月が綺麗に輝いていたし。 …

導出メモ すぐ忘れるｗ これを の二次方程式と見て解くと 但しプラスの解だけを取った。ここで両辺の対数を取ると、 となる。 さて と置き換えることを考えると なので、 が求まる。 check program arcsine implicit none real :: x complex :: c, ci = (0.0…

至点 www.youtube.com Euler の γ 収束が遅い（～O(1/n) ） Euler の γ ですが、引くlog(n) を定義からずらして log(n+1/2) にすると 定義より早く O(1/24n2) で収束するのが草。 計算による比較 小さい項からの和 program euler_gamma implicit none intege…

Fortran 講座 メモリー確保位置 メモリーの割り付けについて詳しく説明されているようですが、ロシア語なのでよく分かりません。 よくみるとウクライナの方のようです。 www.youtube.com ファイヤーフォックス エクステンデッド・エディション(字幕版)クリン…

zeta 関数の非自明零点の計算 昨日求めた von Mangoldt 関数 の定義から、zeta 関数の対数微分が のように得られるようです。この対数微分を逆さまに見れば Newton 法での解の更新を連想させます。 は 0 が多いのでこのまま行けないかと思ったのですが、やは…

Riemann zeta function の非自明零点 先日、「ビジュアル リーマン予想入門」という本を眺めていたのですが、色々やさしく分かりやすく書いてあるので、その結果の一部を再現してみることにします。 ビジュアル リーマン予想入門 ~グラフで解き明かす素数と…

FORTRAN in 100 Seconds www.youtube.com 100秒まとめシリーズに Fortran が来るとか。 The State of Fortran が publish さる ieeexplore.ieee.org Brad Richardson on Twitter: "An article I helped author titled The State of #Fortran was finally rel…

朝日リナさんが三角形を描いています ビデオドライバをフルスクラッチで開発しているようです。 ３D おっさんのライブコーディングは余程興味が無いと見る気がしませんが、２D 美少女 V-tuber は流しておけます。 pic.twitter.com/dz68hvalLv— Asahi Lina /…