FORTRAN

はじめまして !! 神奈川県内（川崎らへん）の大学に通う現在修士１年（４月から２年です！）のkenkenです！ 今は高い強度と高い延性を兼ね備えた金属材料、TWIP材料というものを研究対象に、変形シミュレーションを行っています！（PCカタカタしてます） Fortranというなんとも古き時代に取り残されたような言語（数値計算の分野ではまだまだ現役！）を用いて日々コードとにらめっこしています、、 今このブログを書いているPCはLinux mintというOSを入れたもので、せっかく入れたならなら何かやってみようとブログを書き始めることにしました！！ 研究でもLinux環境を用いて開発することが多いの…

５日前の10月10日に全銀システムに障害が発生し、一部の取引が混乱して銀行間送金が２日間にわたり停止した。今日は全銀システム障害について調べてみた。 まず、全銀システムは、銀行間取引を仲介し、銀行間送金を可能にしているシステムであるが、これを可能にすべく当システムは、各銀行のシステムとネットワークを介して接続する形となっている。今回の障害は、全銀システムの第７次システム更改（2019～）の中で起きたもので、各銀行との間のトランザクション（取引処理）を中継するＲＣ（リレーコンピュータ）の更改において、新規ソフトウェアにエラーが発生し障害に至ったと言われている。 上図は、今回の障害の概要を表してい…

学生時代に、コンピュータプログラムをFORTRAN言語で書いた。社会人になって、パソコンのプログラムをBASIC言語で書いた。その後は、ビジネスの自動化のためにLOTUS 1-2-3やExcel、Word、ファイルメーカーなどではスクリプトで処理をしていた。 スマホでアプリを作ってみたいという気持ちがあり、定年後にプログラム言語をもう一度ゆっくり勉強するつもりだった。なにしろC言語などの本格的なプログラム言語は筆者にとっては異次元の世界で，何度見てもあの記号だらけのプログラムが理解できない。定年年齢を過ぎたが，平の契約社員などでだらだらと仕事を続けているうちに、2022年末に登場したのがCha…

例のウイルスにかかる人がたくさんいます。基礎疾患を持っている人は特に気をつけなければいけません。 そして、基礎をしっかりと学ぶことで応用力を育むことができます。勉強の仕方は人それぞれいろいろあると思うのですが、やはり基礎がしっかりしていないと応用問題を解くときに化けの皮がはがれます。めちょんこ昔、なぜか１日だけ、とあるお子さまの家庭教師をすることになりました。その子は決まった形の計算問題はすらすらと解けるのですが、ちょっとひねった問題を出すと「わからない……」と言って泣きかけるのでどうにかこうにかなだめすかした記憶があるのですが、もしかしたらそんな事実はなかったのかもしれません。とにかく、問題…

試した環境 mac OS 12.6 ifort (IFORT) 2021.7.0 20220726 Python 3.9.12 :: Intel Corporationgfortranを使っているQiitaの記事を参考に試してみたところ、いくつか対策が必要になった。 ※Qiitaの記事 qiita.com具体的には、次の3点。 forpyのモジュールをリンクする際に、呼ぶPythonのライブラリのパスとして-L/opt/intel/oneapi/intelpython/latest/libを追加しないと、libpython3.9.dynalibがリンクできないために最終的な実行形式ファイルがで…

普通はf2pyでPythonからFortranを呼ぶことが多いのだろうが、その逆で、FortranからPythonを呼ぶためのFotranモジュールが公開されている。 公開サイトの例ではGNUのコンパイラが使われているがIntelのコンパイラ（ifort）でも動くらしい。 但し、自分では確認していない。 github.com 実際に試してみた人のリポート： qiita.com

1億年ぶりにfortranを使ったらコンパイルオプションの意味がわからずに整理していなかったことを後悔したので改めて整理しました。 前提として、 ソースはfortran2003をベースに一部fortran2008の機能を含んでいる ifortはLinux版。windows版ではハイフンがスラッシュになったり、区切り文字が違ったり...etc です。 後日CMakeについても同じようにまとめようと思います。 gfortran-debug gfortran-release ifort-debug ifort-release 参照サイト:2022-08-11時点 公式系 その他 gfortran-d…

Fortranで出力したバイナリをPython側で読み込むプログラムを書いていて、挙動が異なる事象に遭遇したのでメモしておきます。Fortran側のコードはこちら。 program main implicit none integer(4), parameter :: rp = kind(1.0d0) integer(4), parameter :: data1 = 10 real(rp), parameter :: data2 = 2.0_rp complex(rp), parameter :: data3 = cmplx(3.0_rp,4.0_rp,kind=rp) write(6,*) …

今回は僕の好きな動画をご紹介させていただきます。1965年〜2022年の間のプログラミング言語の人気の変遷が視覚的にわかり、面白いです。 www.youtube.com ２倍速でさくさく見ると3分半で栄枯盛衰を体験できます。 ざっくり言えば、 60~70年代：FORTRANの時代 〜コンピュータ黎明期〜80~90年代：Cの時代 〜コンピュータの普及〜00~10年代：Java, JSの時代 〜インターネットの普及〜20年代~：Pythonの時代 〜AIの普及〜 みたいな感じでしょうか。70年代にはCOBOLやBASIC、80年代にはPascal, Ada, Lisp、90年代にはC++、00年代…

アマダが量販できたAP40(Fortran版)の後継機として開発したAP40Ver.4(C言語版)が大失敗に終わり、当時アマダのC&Cの事業部長だった丹羽嘉男さんがこの失敗の開発費用を量販版AP40(Fortran版)の開発費用に紛れ込ますという誤魔化しをして以降、失敗が相次ぎ、ボロボロにできたAP60(C言語版)までに費やした外注費用が200億円を遥かに超えてしまった。 ボロボロにできたAP60(C言語版)は売るには売れたが販売台数が僅か数百台に止まり、とてもそれまでの開発費用を賄える額では無かった。しかし、アマダのC&C事業部の近藤課長は数百台売れたことに対して「元を取った」と豪語していた…

最近は自作 32-bit RISC-V CPU で Linux を動かしていた。その前段階として 32-bit RISC-V Linux をビルドして Spike で動かしたのでそのときの備忘録。 Spike でただ Linux を動かすだけなら Buildroot の Spike 用の設定ファイルを使えば簡単なのだが (ただし、ビルド時間が恐ろしく長い...)、今回は以下のような諸事情があって一からビルドした。 OpenSBI に手を加えて自作 CPU 向けの firmware を追加したい 命令セットとして RV32IMA_Zicsr_Zifencei_Zicntr を使いたい デフォルト…

Fortran の環境構築 ( VS Code ) Fortranの環境構築に必要なもの GFortran VSCode １．GFortranのインストール １．１ TDM-GCCのダウンロード jmeubank.github.io 上記のページにアクセスし、赤丸のリンク からTDM-GCCのダウンロードページにアクセスし、tdm-gcc-webdl.exe をダウンロード TDM-GCCのインストール １．２．TDM-GCCのインストール手順 インストーラーの実行: ダウンロードした tdm-gcc-webdl.exe をダブルクリックして実行 インストールウィザードの開始: インストールウ…

さぼっていた間にCUDAが何か全部忘れたので思い出す（せる？） こんにちわ、こんばんわ。かえるのクーの助手の「井戸中 聖」（いとなか セイ）でございます。 実験機のマザーボードetcを更新してWin10からWin11に完全新規の再インストールしたら、失敗を何度かしてライセンスがどっか行っちゃいました。当面、このままいきます。 さて、Win11にしたら環境も新規にインストールします。 とりあえずCUDAやVisualStudio Communityを入れて動かしてみます。 CUDAはサンプルプログラムを動かしたいので、現時点最新のCUDA 12.6ではなく、12.5をいれました。 いつものサンプ…

PythonのPolarsでRのgroup_nestを実現する Polarsでgroup_nestがしたい R or Python 僕がRで分析するときはいつもこちらの記事を参考にしています。 https://qiita.com/kilometer/items/7184904765fbf0f33f04 tidyverseでデータを読み込んだ後は、なんらかの変数でグループ化して、ネストして、各データセットに対してmapを通じて処理をかけていくと。Rはスクリプト言語でありベクトル化された処理に最適化されているので、大容量のデータを扱う場合はこちらの方が速いですし、何より読みやすいです。 Rはユーザ…

OpenMP 6.0 承認 Task parallel 機能の充実と、GPU を見据えた計算デバイス利用の整備をうたっているようです。 Full support for C23, including C attribute syntax, Fortran 2023 and C++23; www.openmp.org 時々流れてくるベンチマークを見る限り、OpenMP はイマイチ性能が出てない感がありますが、デバイスが抽象化されているので長いスパンではベンダー固有のプログラミングより報われるかも知れません。まぁ AI がプログラミングの風景を激変させるだろうから、長期は意味ないのかも知れないけど…

AMD Hip Fortran オーストラリアのスパコンセンターが HIP Fortan の入門講座動画を出している様です。 HIP Fortran は AMD の GPU 用の CUDA Fortran の様なものです。確か外付け GPU と Linux でないと動かなかった気がします。Nvidia の GPU が高いので代替品としてこれから伸びるかも知れません。伸びないかも知れません。 まだ見てませんが、カンガルーやコアラに混じって勉強したいです。 引用しているもの以外にも沢山ありました。 www.youtube.com AMD AOCC Fortran AMD の CPU 向けのコンパイ…

祝 Fortran 8位 近年の最高位！ 復活！（多分）今世紀最高順位！ あの日 揚った Z旗を 父が仰いだ 波の上 今日はその子が その孫が 強く雄々しい 血を継いで 八重の潮路を 越えるのだ www.youtube.com プログラム言語の人気ランキングは Tiobe に限らずべき分布をしていることが知られていて、また時間が経過しても順位は入れ替わりはすれどもべき分布を保つことも観察されています。特定言語が時間軸方向にどのような人気曲線を描くのか一般的には断定はできないものの、プログラム言語人気の場合べき分布的な特徴を持つ可能性が高く、Fortran のようにある時点で高い人気を保ったもの…

FCのスーパーマリオブラザーズとか古いおもちゃなんかを見ると、あんな貧弱なハードウェアと低機能な言語であんなことできたのはすごいとしか思えない。 当時のパーソナルコンピュータでウインドウシステムを実現したXeroxもすさまじいと思う。 正直、FORTRANがこの世に現れてからプログラマーのしていることは何も変わっていない。 ただどんどん簡単になっていくので、誰でもなれるようになっただけだ。 私にスーパーマリオブラザーズをつくれといわれれば、そりゃつくれるが、１０００倍以上のリソースを食うのは間違いないし、組込や、当時の貧弱な仕様でやれといわれたら、無理だろ。現行PGの９９％は無理だろう。 だか…

アマダが量販できたAP40(Fortran版)の後継機として開発したAP40Ver.4(C言語版)が大失敗に終わり、当時アマダのC&Cの事業部長だった丹羽嘉男さんがこの失敗の開発費用を量販版AP40(Fortran版)の開発費用に紛れ込ますという誤魔化しをして以降、失敗が相次ぎ、ボロボロにできたAP60(C言語版)までに費やした外注費用が200億円を遥かに超えてしまった。 ボロボロにできたAP60(C言語版)は売るには売れたが販売台数が僅か数百台に止まり、とてもそれまでの開発費用を賄える額では無かった。しかし、アマダのC&C事業部の近藤課長は数百台売れたことに対して「元を取った」と豪語していた…