csv — обратная матрица из текстового файла с использованием исключения Гаусса. Переполнение стека

Мне нужно написать программу, которая читает матрицу из файла значений, разделенных запятыми, а затем использует метод исключения Гаусса для вычисления обратного значения и записи этого обратного в новый файл.

Чтение этого в порядке, так же как и выписывание обратно. Я думаю, что понимаю, как работает метод исключения Гаусса, и могу сделать это с помощью массива.

#include <iostream>
#include <fstream>
#include <iomanip>
#include <vector>
#include <string>

using namespace std;

//  main program starts here
int main() {
//  create a vector of a vector to store original matrix
//  and matrices after each calculation
vector<vector<double>> data;
int n_com = 0;

//  try to read input file
ifstream readFile("test_data.txt");

if (readFile.is_open()){
while (!readFile.eof()){
int i;

//  declare temporary line vector and line string
vector<double> vline;
string aLine;

//  assign line of file to line string
getline(readFile, aLine);

//  count number of commas in line string
n_com = count(aLine.begin(), aLine.end(), ',');

//  define integer for start of each element
int start = 0;

//  loop over all but final element
for (i=0; i < n_com; i++) {

//  declare and find position of next comma
int comma_pos;
comma_pos = aLine.find(',', start);

//  declare string for element and assign substring to it
string elems;
elems = aLine.substr(start, comma_pos - start);

//  convert string to double
double elemd = atof(elems.c_str());

//  push back double to temporary vector
vline.push_back(elemd);

//  redefine start for next iteration
start = comma_pos + 1;
}
//  assign final element to string
string final_elems = aLine.substr(start, aLine.length() - start);

//  convert final element to double and push back to vector
double final_elemd = atof(final_elems.c_str());
vline.push_back(final_elemd);

//  push back line vector to data vector
data.push_back(vline);
}
}
else {
//  print error if unable to open file
printf("Error unable to open input file!\n");
//  exit program
exit(1);
}
//  close input file
readFile.close();

Так вот как я прочитал в исходной матрице.

Это код, который я имею для исключения Гаусса для массива

//  calculate width and length of original data (no. of rows and columns)
int length = data.size();
int width = n_com + 1;

//  create new file to write to
ofstream writeFile ("tranpose.txt");

//  check outfile is open
if (writeFile.is_open()){

//  declare indices
//  width (columns)
int i;
// length (rows)0
int j;
//  k
int k;

//  declare a float?
float data[10][10] = {0},d;

//  identity matrix
for (i=1; i <= length; i++){
for (j=1; j <= 2 * length; j++){
if (j == (i + length)){
data[i][j] = 1;
}
}
}

//  partial pivoting
for (i=length; i > 1; i--){
if (data[i-1][1] < data[i][1]){
for(j=1;j <= length * 2; j++){
d = data[i][j];
data[i][j] = data[i-1][j];
data[i-1][j] = d;
}
}
}
cout<<"Augmented Matrix: "<<endl;
for (i=1; i <= length; i++){
for (j=1;j <= length * 2; j++){
cout<<data[i][j]<<"  ";
}
cout<<endl;
}

//  reducing to diagonal matrix
for (i=1; i <= length; i++){
for (j=1; j <= length * 2; j++){
if (j != i){
d = data[j][i] / data [i][i];
for (k=1; k<= length * 2; k++){
data[j][k] = data[j][k] - (data[i][k] * d);
}
}
}
}

//  reducing to unit matrix
for (i=1; i <= length; i++){
d = data[i][i];
for (j=1; j <= length * 2; j++){
data[i][j] = data[i][j] / d;
}
}
//  print inverse matrix in console
cout<<"Inverse Matrix "<<endl;
for (i=1; i <= length; i++){
for (j = length + 1; j <= length * 2; j++){
cout<<data[i][j]<<"  ";
}
cout<<endl;
}
//  loop over all rows
for (i=1; i <= length; i++){
//  loop over all columns
for (j = length + 1; j <= length * 2; j++){
//  print data in transposed positions excluding last value
//  i.e. [j][i] instead of [i][j]
writeFile << setw(4) << fixed << setprecision(2) << data[i][j] << ",";
}
// print onto new line
int i = length - 1;
writeFile << setw(4) << fixed << setprecision(2) << data[i][j] << "\n";

}

//  close written file
writeFile.close();

Как я могу написать это так, чтобы он использовал данные, хранящиеся в моем векторе или матричном файле, вместо обычного массива чисел?