Файловое чтение и запись данных для решения СЛАУ

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <ctime>
#include <omp.h>
using namespace std;
 
int main() 
 
{
    //srand(time(0));
    //unsigned int start_time =  clock(); // начальное время
    double start_time, end_time, tick;
    start_time = omp_get_wtime();
    
 
 
    setlocale(LC_ALL,"Russian");
    int n=0;
    printf("Введите количество неизвестных\n");
    cin>>n;
    printf("Задайте матрицу коэффициентов\n");
//vstavka
//r(int i=0; i<n; i++)
//r(int j=0; j<m; j++)
   //mas[i][j] = 0+rand()%10;
 
//vstavka
 
 
 
 
    double ** A=new double*[n];
    //double ** L=new double*[n];
    for(int i=0;i<n;i++){
        A[i]=new double[n];
        //L[i]=new double[n];
        for(int j=0;j<n;j++){
        A[i][j] = 0+rand()%1000;
 
          //cin>>A[i][j];
        }
    }
    printf("Введите вектор свободных членов\n");
    double * B=new double[n];
    for(int i=0;i<n;i++){
      //cin>>B[i];
    B[i] = 0+rand()%1000;
    }
    //метод Гаусса
    //получаем верхнетреугольную матрицу
 
    for(int i=0;i<n-1;i++){
        omp_set_nested(1);
#pragma omp parallel num_threads(1) firstprivate(i)
        {
            int tid=omp_get_thread_num();
            int j=tid+i+1;
            double fix=A[j][i]/A[i][i];
            B[j]=B[j]-B[i]*fix;
#pragma omp parallel num_threads(1)
            {
                int k=omp_get_thread_num();
                k=k+i;
                A[j][k]=A[j][k]-A[i][k]*fix;
            }
        }
 
    }
 
    //вывод преобразованной расширенной матрицы
    printf("Матрица U\n");
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            printf("%f ",A[i][j]);
        }
        printf("  | %f \n",B[i]);
 
    }
 
    //определяем и выводим вектор решений
    double * X=new double[n];
//#pragma omp parallel num_threads(1) firstprivate(1) shared(count,X,A,B)
    {
#pragma omp for schedule(static,1)
        for(int i=n-1;i>=0;i--){
            X[i]=B[i];
 
            for(int j=n-1;j>i;j--){
 
                X[i]=X[i]-A[i][j]*X[j];
            }
#pragma omp critical
            X[i]=X[i]/A[i][i];
            printf("%f\n",X[i]);
        }
 
    }
    //system("pause");
    //unsigned int end_time = clock(); // конечное время
    //unsigned int search_time = end_time - start_time; // искомое время
 
    //cout << "vremya=" << search_time << endl;
    end_time = omp_get_wtime();
    tick = omp_get_wtick();
    printf("Время  %lf\n", end_time-start_time);
}

    int n;
    double ** A;
    double * B;
 
    int v;
    cout << "Прочитать из файла - 1, заполнить и записать в файл - 2): "
    cin << v;
    if (v == 1)
    {
        ifstream ifs("data.txt");
        ifs >> n;
        A=new double*[n];
        for(int i=0;i<n;i++)
        {
            A[i]=new double[n];
            for(int j=0;j<n;j++) ifs >> A[i][j];
        }
        B=new double[n];
        for(int i=0;i<n;i++) ifs >> B[i];
    }
    else if (v == 2)
    {
    ofstream ofs("data.txt")
    printf("Введите количество неизвестных\n");
    cin>>n;
    ofs << n << endl;
 
    A=new double*[n];
    for(int i=0;i<n;i++){
        A[i]=new double[n];
        for(int j=0;j<n;j++){
        A[i][j] = 0+rand()%1000;
        ofs << A[i][j] << endl;
        }
    }
    B=new double[n];
    for(int i=0;i<n;i++)
    {
        B[i] = 0+rand()%1000;
        ofs << B[i] << endl;
    }
 
    }