gaussSeidel.h

Go to the documentation of this file.

 
#pragma once
#include "../../Matrix.h"
#include "../../matrix_utils.h"
#include "LU.h"
#include "backwardSub.h"
#include "forwardSub.h"
#include <vector>
 
Matrix<double> gaussSeidel(const Matrix<double>& A, const Matrix<double>& b) {
  auto LR = LU(A);
 
  auto bCopy = b;
  for(size_t i = 0; i < b.rows(); i++) { bCopy.SetRow(i, b(LR.second[i])); }
 
  auto c = forwardSub(LR.first, bCopy);
  return backwardSub(LR.first, c);
}