include("getDivGrad.jl") include("ssor.jl") A = getDivGrad(32,32,32) rhs = randn(size(A,1)) tol = 1e-9 x1,flag1,relres1,iter1 = ssor(A,rhs,tol=tol,maxIter=1000,out=2)
