t-nissieの日記: Fortran
日記 by
t-nissie
久しぶりにFortranでプログラミング。
call DGEMV('N', 5, 5, 1.0d0, O, 5, Dzz_prime(1,i), 1, 0.0d0, xi, 1)
と書くべきところを
call DGEMV('N', 5, 5, 1.0, O, 5, Dzz_prime(1,i), 1, 0.0, xi, 1)
と書いてしまい、ハマりました。
DGEMVはBLASのサブルーチンです。
g77には「浮動小数点数定数はぜんぶ倍精度ですよ」という
オプションはないのでしょうか。
BLASとLAPACKの主要サブルーチンのinterface文 (Cでいえば
関数の宣言が書いてあるヘッダファイルみたいなもの)、
どこかに落ちていないでしょうか。
というか、Fortranでプログラミングしたくないです。
call DGEMV('N', 5, 5, 1.0d0, O, 5, Dzz_prime(1,i), 1, 0.0d0, xi, 1)
と書くべきところを
call DGEMV('N', 5, 5, 1.0, O, 5, Dzz_prime(1,i), 1, 0.0, xi, 1)
と書いてしまい、ハマりました。
DGEMVはBLASのサブルーチンです。
g77には「浮動小数点数定数はぜんぶ倍精度ですよ」という
オプションはないのでしょうか。
BLASとLAPACKの主要サブルーチンのinterface文 (Cでいえば
関数の宣言が書いてあるヘッダファイルみたいなもの)、
どこかに落ちていないでしょうか。
というか、Fortranでプログラミングしたくないです。
Fortran More ログイン