205
205
x = rand (n,n)
206
206
y = rand (n,n)
207
207
208
- Pl, Pr = LinearSolve. scaling_preconditioner (s)
208
+ Pl, Pr = LinearSolve. scaling_preconditioner (1 / s)
209
209
210
210
mul! (y, Pl, x); @test y ≈ s * x
211
211
mul! (y, Pr, x); @test y ≈ s \ x
@@ -219,16 +219,19 @@ end
219
219
220
220
@testset " vector scaling_preconditioner" begin
221
221
s = rand (n)
222
- Pl, Pr = LinearSolve. scaling_preconditioner (s)
222
+ Pl, Pr = LinearSolve. scaling_preconditioner (1 ./ s)
223
223
224
- mul! (y, Pl, x); @test y ≈ s * x
225
- mul! (y, Pr, x); @test y ≈ s \ x
224
+ x = rand (n,n)
225
+ y = rand (n,n)
226
226
227
- y . = x; ldiv! ( Pl, x); @test x ≈ s \ y
228
- y . = x; ldiv! ( Pr, x); @test x ≈ s * y
227
+ mul! (y, Pl, x); @test y ≈ s .* x
228
+ mul! (y, Pr, x); @test y ≈ s .\ x
229
229
230
- ldiv! (y, Pl, x); @test y ≈ s \ x
231
- ldiv! (y, Pr, x); @test y ≈ s * x
230
+ y .= x; ldiv! (Pl, x); @test x ≈ s .\ y
231
+ y .= x; ldiv! (Pr, x); @test x ≈ s .* y
232
+
233
+ ldiv! (y, Pl, x); @test y ≈ s .\ x
234
+ ldiv! (y, Pr, x); @test y ≈ s .* x
232
235
end
233
236
234
237
@testset " ComposePreconditioenr" begin
0 commit comments