|
996 | 996 | @test A * D ≈ mul!(S, A, D) ≈ M * D |
997 | 997 | @test D * A ≈ mul!(S, D, A) ≈ D * M |
998 | 998 | @test mul!(copy(S), D, A, 2, 2) ≈ D * M * 2 + S * 2 |
| 999 | + @test mul!(copy(S), D, A, 0, 2) ≈ D * M * 0 + S * 2 |
999 | 1000 | @test mul!(copy(S), A, D, 2, 2) ≈ M * D * 2 + S * 2 |
| 1001 | + @test mul!(copy(S), A, D, 0, 2) ≈ M * D * 0 + S * 2 |
1000 | 1002 |
|
1001 | 1003 | A2 = Bidiagonal(dv, zero(ev), uplo) |
1002 | 1004 | M2 = Array(A2) |
@@ -1074,10 +1076,12 @@ end |
1074 | 1076 | @test B * v ≈ M * v |
1075 | 1077 | @test mul!(similar(v), B, v) ≈ M * v |
1076 | 1078 | @test mul!(ones(size(v)), B, v, 2, 3) ≈ M * v * 2 .+ 3 |
| 1079 | + @test mul!(ones(size(v)), B, v, 0, 3) ≈ M * v * 0 .+ 3 |
1077 | 1080 |
|
1078 | 1081 | @test B * B ≈ M * M |
1079 | 1082 | @test mul!(similar(B, size(B)), B, B) ≈ M * M |
1080 | 1083 | @test mul!(ones(size(B)), B, B, 2, 4) ≈ M * M * 2 .+ 4 |
| 1084 | + @test mul!(ones(size(B)), B, B, 0, 4) ≈ M * M * 0 .+ 4 |
1081 | 1085 |
|
1082 | 1086 | for m in 0:6 |
1083 | 1087 | AL = rand(m,n) |
|
0 commit comments