|
| 1 | +### Needs to be kept out for julia 1.0 |
| 2 | + |
| 3 | +struct baseSE <: KernelFunctions.BaseKernel end |
| 4 | +(k::baseSE)(x, y) = exp(-evaluate(SqEuclidean(), x, y)) |
| 5 | + |
1 | 6 | @testset "kernelmatrix" begin
|
2 | 7 |
|
3 | 8 | rng = MersenneTwister(123456)
|
|
17 | 22 | Kdiag = [zeros(dims[1]),zeros(dims[2])]
|
18 | 23 | s = rand(rng)
|
19 | 24 | k = SqExponentialKernel()
|
20 |
| - struct baseSE <: KernelFunctions.BaseKernel end |
21 |
| - (k::baseSE)(x, y) = exp(-evaluate(SqEuclidean(), x, y)) |
22 | 25 | newk = baseSE()
|
23 | 26 | kt = transform(SqExponentialKernel(),s)
|
24 | 27 |
|
|
29 | 32 | @test kernelmatrix!(K[2], k, vA) ≈ kernelmatrix(k, A) atol = 1e-5
|
30 | 33 | @test kernelmatrix!(K[2], k, vA, vB) ≈ kernelmatrix(k, A, B) atol = 1e-5
|
31 | 34 | for obsdim in [1,2]
|
32 |
| - @show obsdim |
33 | 35 | @test kernelmatrix!(K[obsdim], k, A, B, obsdim = obsdim) == kernelmatrix(k, A, B, obsdim = obsdim)
|
34 | 36 | @test kernelmatrix!(K[obsdim], k, A, obsdim = obsdim) == kernelmatrix(k, A, obsdim = obsdim)
|
35 | 37 | @test kerneldiagmatrix!(Kdiag[obsdim], k, A, obsdim = obsdim) == kerneldiagmatrix(k, A, obsdim = obsdim)
|
|
54 | 56 | @test k(A,obsdim=obsdim) == kernelmatrix(k,A,obsdim=obsdim)
|
55 | 57 | # @test KernelFunctions._kernel(k,1.0,2.0) == KernelFunctions._kernel(k,[1.0],[2.0])
|
56 | 58 | @test_throws DimensionMismatch kernelmatrix(k,A,C,obsdim=obsdim)
|
57 |
| - @test kernelmatrix!(K[obsdim], newk, A, B, obsdim = obsdim) ≈ kernelmatrix(k, A, B, obsdim = obsdim) |
58 |
| - @test kernelmatrix!(K[obsdim], newk, A, obsdim = obsdim) ≈ kernelmatrix(k, A, obsdim = obsdim) |
59 |
| - @test kerneldiagmatrix!(Kdiag[obsdim], newk, A, obsdim = obsdim) ≈ kerneldiagmatrix(k, A, obsdim = obsdim) |
| 59 | + @test kernelmatrix(newk, A, B, obsdim = obsdim) ≈ kernelmatrix(k, A, B, obsdim = obsdim) |
| 60 | + @test kernelmatrix(newk, A, obsdim = obsdim) ≈ kernelmatrix(k, A, obsdim = obsdim) |
| 61 | + @test kerneldiagmatrix(newk, A, obsdim = obsdim) ≈ kerneldiagmatrix(k, A, obsdim = obsdim) |
60 | 62 | end
|
61 | 63 | end
|
62 | 64 | @testset "Transformed Kernel Matrix Operations" begin
|
|
0 commit comments