solve_ols: Solve OLS
In shihhowie/hw2:

Description Usage Arguments Value Examples

Solve ordinary least square using iterative methods

1	solve_ols(A, b, method = "Jacobi", n_cores = 1, iter = 10)

`A`	design matrix
`b`	target vector
`method`	Jacobi(Jacobi), or Gauss-Seidal(GS)
`n_cores`	Number of parallel processes
`iter`	number of iteration for the iterative method

list(x, run_time)

A = diag(2, 10, 10)
b = 1:10
Jacobi_res = solve_ols(A,b,method="Jacobi",n_cores=2)
GS_res = solve_ols(A,b,method="GS")
print(Jacobi_res[1])
print(GS_res[1])