Там человек всю жизнь писал на фортране 77. Почему он два абзаца документации должен писать на более новой версии фортрана? 😊 Это же отдаёт новой школой, а рулит старая. 😊

Если правильно писать на 90-м, то ошибки находятся гораздо раньше :-)
Да и если что-то чужое цеплять, то в 77 это можно выяснить графом вызовов, а в 90 графом модулей, который проще и короче

Какой граф вызовов? Это документация была в стиле

VAR0 SHOULD BE 0 WHEN N ITER > 100
KX SOULD BE < 0 FOR ....

Ну хотя бы такая есть)

Ну я просто тогда оценил, что "на любом языке можно писать, как на Фортране".

А к хаскелю нельзя прицепить в качестве выходного файла фортрановский код, который потом можно будет сктмпилировать?

GHC такого вроде не умеет. А зачем это вам? В Хаскеле значительно хуже с массивами, чем в Фортране. Даже в Клине хуже, хотя там они изменяемые и синтаксического сахару насыпали.

QuickJS только

Ему уже года три. Он не очень production ready был три года назад, и вряд ли оно двигается сильно.

Можно через плагины прицепиться к компилятору и на каком-то уровне выдавать фортран. Так делали, например, в clash, чтобы выдавать VHDL. Но задача эта не самая тривиальная и я не уверен, что найдётся команда, которая это хочет и может. Гораздо проще написать edsl, который будет генерировать фортран на выходе. С этим вполне может справиться средняя команда, и обычно так делают с Си, самое крупное такое решение это ivory и cryptol, но мелкие специализированные решения были во многих проектах (у меня минимум в 2 и ещё про несколько знаю)

Ещё есть FFI в фортран, если я правильно помню, что тоже может быть вариантом решения для задач. Про top-10, в принципе в Haskell есть режим сборки via-c, выдающий на выходе си код, который можно собирать любым сишным компилятором, можно и этим путём идти

“Do not trust this software because I honestly have no idea why it works.”

Моя новая любимая фраза

В принципе С тож пойдет. Главное, чтобы использование памяти оптимизировалось :)

Его уже юзают в проде🤔

Оптимизировалось под архитектуру в плане доступа с учетом топологии машины, или оптимизировалось - максимально уменьшалось?

Ну например чтобы лишний раз матрицы не сохраняло в отдельных переменных при их обновлении.
Например, чтобы не было так:
H_updated = 0
H_updated = H + H * dEdR * R
H = H_updated
Ну и лишний раз большие матрицы не транспонировало: это очень дорого по памяти

Ого, в фортране последнее на уровне компилятора?

*на уровне стандарта языка

Но капец дорогое по памяти

В имплементации GCC, по крайней мере