Решение линейного уравнения методом Гаусса-Зейделя дает неверные результаты

Question

Использование R-пакета optR для решения линейного уравнения дает неверные результаты для метода гаусса-сейделя. Я пробовал решать с помощью разных методов, и это выглядит нормально. Например, решить или использовать optR, но с помощью метода gauss. Пример здесь. Или я делаю что-то неправильно?

library(optR) 
A<-matrix(c(5,-1,2, 3, 8,-2,1, 1,4), nrow=3,ncol=3, byrow = TRUE) 
b<-matrix(c(12,-25,6), nrow=3,ncol=1,byrow=TRUE) 

Z1<-optR(A, b,method="gaussseidel", iter=500, tol=1e-7) 
Z2 <-optR(A, b, method="gauss", iter=500, tol=1e-7) 
Z3 <- solve(A,b) 

print(Z1) 
print(Z2) 
print(Z3) 

A %*% Z1$beta - b 
A %*% Z2$beta - b 
A %*% Z3 - b 

Коэффициенты Z1 неверны:

> print(Z1) 
call: 
optR.default(x = A, y = b, method = "gaussseidel", iter = 500, 
    tol = 1e-07) 

Coefficients: 
      [,1] 
[1,] 4.183908 
[2,] -2.379310 
[3,] -1.781609 
> print(Z2) 
call: 
optR.default(x = A, y = b, method = "gauss", iter = 500, tol = 1e-07) 

Coefficients: 
    [,1] 
[1,] 1 
[2,] -3 
[3,] 2 
> print(Z3) 
    [,1] 
[1,] 1 
[2,] -3 
[3,] 2 
> 
> A %*% Z1$beta - b 
      [,1] 
[1,] 7.735632 
[2,] 22.080460 
[3,] -11.321839 
> A %*% Z2$beta - b 
       [,1] 
[1,] 0.000000e+00 
[2,] -3.552714e-15 
[3,] 0.000000e+00 
> A %*% Z3 - b 
       [,1] 
[1,] 0.000000e+00 
[2,] -3.552714e-15 
[3,] 0.000000e+00 

Спасибо заранее.

Answer 1

Так что это решение сейчас. Это была ошибка в пакете. Поскольку версия 1.2.5 optR решена и дает правильные результаты.