C++で前進代入、後退代入、gauss消去

Question

//メインのコード
#include <stdio.h>
#include <iostream>
#include <math.h>
#include "LinearAlgebra.hpp"

int main(int argc,const char * argv[]){
    using namespace std;
    vector<vector<double> > A(3,vector<double>(3));//行列Aの設定
    A= {{3,-1,-1},{-1,3,-1},{-1,-1,3}};
    vector<double> b(3);
    b= {1,1,1};//ベクトルbの設定
    cout<<"A="<<endl;
    printMatrix(A);//Aを表示
    cout<<"b="<<endl;
    printVector(b);//bを表示
    vector<double> x(3);
    x=GaussianElimination(A,b);//gauss消去
    cout<<"x="<<endl;
    printVector(x);//Ax=bの解を出力```ここに言語を入力
    return 0;
}

//関数ファイルLinearAlgebra.cpp
#include "LinearAlgebra.hpp"
#include <vector>
#include <math.h>
#include <iostream>
using namespace::std;

//ベクトルaを出力
void printVector(vector<double> a){
    for(int i=0;i<a.size();i++){
        printf("%.2e\t",a[i]);
    }
    printf("\n");
}

//行列Aを出力
void printMatrix(vector<vector<double> > A){
    for(int i=0;i<A.size();i++){
        for(int j=0;j<A.size();j++){
            printf("%2e\t",A[i][j]);
        }
        printf("\n");
    }
    printf("\n");
}
//後退代入Ax=b
vector<double> Backward_sub(vector<vector<double> > A,vector<double> b){
    int N=int(b.size());
    vector<double> x(N);
    for(int k=N-1;k>=0;k--){
        float C=0;
        for(int j=k+1;j<N;j++){
            C=A[k][j]*x[j];
        }
        x[k]=(b[k]-C)/A[k][k];
    }
    return x;
}
//gause消去
vector<double> GaussianElimination(vector<vector<double> > a,vector<double> c){
    unsigned long N=c.size();
    vector<vector<double> > A(N,vector<double>(N));
    vector<double> b(N);
    A=a;
    b=c;
    //前進消去
    for(int k=0;k<N-1;k++){
        for (int i=k+1;i<=N;i++){
            int B=a[i][k]/a[k][k];
            for(int j=k+1;j<N;j++){
                A[i][j]=A[i][j]-A[i][k]*A[k][j];
            }
            b[i]=b[i]-A[i][k]*b[k];
        }
    }
    cout<<"前進消去\n";
    cout<<"A="<<endl;
    printMatrix(A);
    cout<<"b="<<endl;
    printVector(b);
    //後退代入
    cout<<"後退代入\n";
    vector<double> x(N);
    x=Backward_sub(A,b);
    return x;
}

//ヘッダーファイル　LinearAlgebra.hpp
#ifndef LinearAlgebra_hpp
#define LinearAlgebra_hpp

#include <stdio.h>
#include <vector>
#include<math.h>
using namespace::std;

void printVector(vector<double>);
void printMatrix(vector<vector<double> >);
vector <double> VectorSubtract(vector <double>,vector <double>);
vector <double> MatrixVector(vector<vector<double> >,vector <double>);
vector <double> ResidualError(vector<vector <double> >,vector <double>,vector <double>);
double VectorNorm1(vector <double>);
double VectorNorm2(vector <double>);
double VectorNormInfty(vector <double>);
double MatrixNorm1(vector<vector <double> >);
double MatrixNormIntfty(vector<vector <double> >);
vector<double> GaussianElimination(vector<vector<double> > ,vector<double> );
vector<double> Backward_sub(vector<vector<double> > ,vector<double> );

#endif /* LinearAlgebra_hpp */

コードを全て記載しておきます。行列AとベクトルbのAx＝bの解xを求めるコードを作っています。
出力するとエラーが出てしまい、そのエラーコードが
main.cpp:9:8: error: expected expression
A= {{3,-1,-1},{-1,3,-1},{-1,-1,3}};
^
main.cpp:11:8: error: expected expression
b= {1,1,1};
^
となってます。教授にも質問しましたが、おそらくLinearAlgebra.cppのどこかが間違っていると思われます。
どのように修正すればいいか教えてくれると幸いです.
使ってるのはXcodeのversion１１.５です

Comments

[episteme]

https://teratail.com/questions/273681
でベストアンサーつけたろ? 解決したんだろ?

Answer 1

コンパイルエラーは出ませんが、GaussianElimination()の以下の箇所でiが3のときにa[i],A[i],b[i]が範囲外アクセスしていますよ。

		for (int i = k + 1; i <= N; i++) {
			int B = a[i][k] / a[k][k];
			for (int j = k + 1; j < N; j++) {
				A[i][j] = A[i][j] - A[i][k] * A[k][j];
			}
			b[i] = b[i] - A[i][k] * b[k];
		}

Answer 2

問題を再現できません。
Visual Studio 2019 でコンパイル/実行を確認しました。

A=
3.000000e+00    -1.000000e+00   -1.000000e+00
-1.000000e+00   3.000000e+00    -1.000000e+00
-1.000000e+00   -1.000000e+00   3.000000e+00

b=
1.00e+00        1.00e+00        1.00e+00

[追記]

教授にも質問しましたが、おそらくLinearAlgebra.cppのどこかが間違っていると思われます。

それが教授の返答だとしたら能無し教授だ。