Chapter 2

f_02_16.jl

BICGSTAB scheme (right preconditioned)

Other versions: j_02_16.jl, nm216.jl


In [1]:
using NMfE

In [2]:
A=Float64[16 4 8;4 5 -4;8 -4 22];

In [3]:
b=Float64[4, 2, 5];

In [4]:
c = A\b


Out[4]:
3-element Array{Float64,1}:
 -0.25
  1.0 
  0.5 

In [5]:
round(A * c, 14) == b


Out[5]:
true