Метод Гаусса-Зейделя и методы связанных градиентов для матрицы - PullRequest
0 голосов
/ 28 декабря 2018

У меня есть матрица

, и я должен написать код, используя методы Гаусса-Зейделя и связанный градиент, взяв структуру матрицы.

Ax = e, где A - матрица, а e - вектор со значениями 1

. Я не знаю, как писать код, используя методы связанного градиента и мой алгоритм Гаусса-Зейделя.нет основной части, где я добавляю это все думает

//part of gauss-seidel method    

const int N = 128;             //size of array
const int no_of_iter = 128;    //iterations

int main() {
  double result[N];            //array for result
  double result_pom[N];        //temporary result array
  double sum = 0.0;

  int x, y;                            
  for (int i = 0; i < no_of_iter; i++) {
    for (y = 0; y < N; y++) {
       result_pom[y] = result[y]; //set values of result to result_pom
    }

    for (x = 0; x < N; x++) {
      sum = 0.0;
      //for functions where I add (x,y) el

      result[x] = 0.25 * (1 - sum); //because 4 is dominant el of matrix and
                                    //1 is value of vector e
    }
  }
}
...