[Tau_0]

Это то самое набившее оскомину cтруктурное программирование, к созданию которого и приложил руку Эдсгер. »

Тогда мне совсем непонятен полёт Вашей высокой мысли...
Продолжим однако, оно того стОИт… Заодно и проверим фортранщика на зуб…
Вот Вам код подпрограммы всего на двух страничках…
LU_73
LU_74

Это не деЦЦкий код. Это из одного серьёзного пакета для FEM/МКЭ, написанного ещё на легендарном FORTRAN 66 для DEC/VAX. По крайней мере код профилирован, и поверьте мне на слово, более чем удовлетворяет требованиям структурного программирования… Как сверху вниз, так и снизу вверх, а также в любом направлении…

Речь идёт в целом море “квадратов и квадратиков” и подпрограмм в пакете полно… Но ни в одной компоненте не найдёте юниорских, так милых Вашему сердцу проверок…

Добавлено
=======
Нашёл pdf оригинала -- его проще обозреть...
ComputerSolution of Sparse Linear Systems