fortran66 https://blog.hatena.ne.jp/fortran66/ fortran66のブログ https://fortran66.hatenablog.com/ Fortran2003 昔Householder法とQR法のサブルーチンを外人からもらいましたが、EISPACKの tred と tql で F66 で書かれていました。のちに Numerical Recipes の F77 版に置き換えたけれども、基本的に LAPACK のものだと思います。微妙に変えてあって、tql は旧版で収束しない時でも計算できるようになったけれど収束判定が厳しくなりすぎて遅くなった記憶が。tred は昔のプログラムなので、記憶領域を使わないように色々芸当を使っていますが、策士策に溺れるような感じがしなくもないような・・ワーク配列を取らないで、元の行列の空きスペースを利用するようにしています。… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Ffortran66.hatenablog.com%2Fentry%2F20111204%2F1322934815" title="Householder法による三重対角化 - fortran66のブログ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> http://cdn-ak.f.st-hatena.com/images/fotolife/f/fortran66/20111204/20111204025305.gif Hatena Blog https://hatena.blog 2011-12-04 02:53:35 rich https://fortran66.hatenablog.com/entry/20111204/1322934815 1.0 100%