solve_ols: Solves Ax=b using Gauss-Seidel or Jacobi iteration.

In DavidJamesKent/stsci6520_hw2: STSCI6520 HW2

Description

Solves Ax=b using Gauss-Seidel or Jacobi iteration.

Usage

solve_ols(A, b, x0 = NULL, niter = 100, method = "Gauss-Seidel",
  ncores = 1, verbose = T)

Arguments

A	n x n matrix.
b	Length n vector.
x0	(Optional) length n initial value for iterations. All zeroes if not provided.
method	(Optional) method, one of "Gauss-Seidel" or "Jacobi".
ncores	(Optional) number of cores to use. Ignored if method is "Gauss-Seidel".
verbose	(Optional). If TRUE, print information about progress to the screen.
iter	(Optional) number of iterations to perform.

Value

Value of x after iter iterations of the method.

DavidJamesKent/stsci6520_hw2 documentation built on May 26, 2019, 12:30 a.m.