行列式のプログラミング

Question

行列式のプログラミングを書いたつもりなのですが、行の入れ替えを行う際につくはずの負の符号が出力するとへんな位置に出されます。下のコードの場合、det=-12のはずがdet=1-2 となります。
何がまずいのでしょうか？返答お待ちしております。

コード
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
int main(void){
 double a[4][4]={{0.0,0.0,0.0,1.0},
                    {0.0,-2.0,0.0,0.0},
                    {0.0,0.0,-3.0,0.0},
                    {2.0,0.0,0.0,0.0}};
 double sum, tmp, absmax;
 n=4;
 int imax;

 //部分ピボット選択付きの前進消去のプログラム
 int k;
 int i;
 int j;
 double d;

 for(k=0;k<=n-2;k++){
 absmax=fabs(a[k][k]);  //とりあえずa[k][k]が最大としておく
 imax=k;
 for(i=k+1;i<n;i++){
  if(fabs(a[i][k])>absmax){  //もっと大きいものがあれば
   absmax=fabs(a[i][k]);    //最大の候補を入れ替え
   imax=i;
   }
  }
   if(imax!=k){ // imaxがkでない場合だけ交換作業
    for(j=k;j<n;j++){
     tmp=(-1.0)*a[k][j];          //a[k][j]とa[imax][j]を交換
     a[k][j]=a[imax][j];   //tmpは交換のための
     a[imax][j]=tmp;       //作業用変数
     }
    }
   for(i=k+1;i<=n-1;i++){
    for(j=k+1;j<=n-1;j++){
     a[i][j]=a[i][j]-a[k][j]*a[i][k]/a[k][k];
     }
    }
  for(i=0;i<=n-2;i++){
   d=a[0][0];
   d*=a[i+1][i+1];
   }
 }
 cout << "det=" << showpoint << setprecision << d << "\n";
}

Answer 1

setprecision
->
setprecision(2)

cout << d << "\n";
としても2が出てきます

Comments

[TakumiMiyamoto]

ありがとうございます。一箇所凡ミスも見つかりました。