|
420 | 420 | @test compare(mul!, AT, rand(T, 2,2), rand(T, 2,1), f(rand(T, 2))) |
421 | 421 | end |
422 | 422 | end |
| 423 | + |
| 424 | + @testset "$T gemv zero-dim" for T in eltypes |
| 425 | + y, A, x = rand(T, 4), rand(T, 4, 0), rand(T, 0) |
| 426 | + |
| 427 | + @test compare(*, AT, A, x) |
| 428 | + @test compare(mul!, AT, y, A, x) |
| 429 | + |
| 430 | + y = rand(T, 4) |
| 431 | + @test compare(mul!, AT, y, A, x, Ref(T(4)), Ref(T(5))) |
| 432 | + end |
423 | 433 | end |
424 | 434 |
|
425 | 435 | @testsuite "linalg/mul!/matrix-matrix" (AT, eltypes)->begin |
|
434 | 444 | @test compare(mul!, AT, C, f(A), g(B), Ref(T(4)), Ref(T(5))) |
435 | 445 | @test typeof(AT(rand(T, 3, 3)) * AT(rand(T, 3, 3))) <: AbstractMatrix |
436 | 446 | end |
| 447 | + |
| 448 | + @testset "$T gemm zero-dim" for T in eltypes |
| 449 | + A, B, C = rand(T, 4, 0), rand(T, 0, 4), rand(T, 4, 4) |
| 450 | + |
| 451 | + @test compare(*, AT, A, B) |
| 452 | + @test compare(mul!, AT, C, A, B) |
| 453 | + |
| 454 | + C = rand(T, 4, 4) |
| 455 | + @test compare(mul!, AT, C, A, B, Ref(T(4)), Ref(T(5))) |
| 456 | + end |
437 | 457 | end |
438 | 458 |
|
439 | 459 | @testsuite "linalg/norm" (AT, eltypes)->begin |
|
0 commit comments