File tree Expand file tree Collapse file tree 1 file changed +8
-4
lines changed Expand file tree Collapse file tree 1 file changed +8
-4
lines changed Original file line number Diff line number Diff line change @@ -104,19 +104,23 @@ A = poisson(1000)
104
104
A = float .(A)
105
105
ml = ruge_stuben (A)
106
106
x = solve (ml, A * ones (1000 ))
107
- @test sum (abs2, x - ones (1000 )) < 1e-10
107
+ @test sum (abs2, x - ones (1000 )) < 1e-8
108
108
109
109
A = load (" randlap.jld" )[" G" ]
110
- ml = ruge_stuben (A)
110
+ smoother = GaussSeidel (ForwardSweep ())
111
+ ml = ruge_stuben (A, presmoother = smoother,
112
+ postsmoother = smoother)
111
113
x = solve (ml, A * ones (100 ))
112
- @test sum (abs2, x - zeros (100 )) < 1e-10
114
+ @test sum (abs2, x - zeros (100 )) < 1e-8
113
115
114
116
end
115
117
116
118
@testset " Preconditioning" begin
117
119
A = load (" thing.jld" )[" G" ]
118
120
n = size (A, 1 )
119
- ml = ruge_stuben (A)
121
+ smoother = GaussSeidel (ForwardSweep ())
122
+ ml = ruge_stuben (A, presmoother = smoother,
123
+ postsmoother = smoother)
120
124
p = aspreconditioner (ml)
121
125
b = zeros (n)
122
126
b[1 ] = 1
You can’t perform that action at this time.
0 commit comments