mjp98

Call `domain(sp)` rather than `domain(cfs)` in PeriodicLine hilbert

1

This addresses https://github.com/JuliaApproximation/SingularIntegralEquations.jl/issues/152 by correcting calls of `domain(f::AbstractVector)` to `domain(S::Space)` in `hilbert(S::Space{1/(1+z^2),PeriodicLine()) hilbertF = z-> -z/(1+z^2) z = pi @test hilbert(F,z) ≈ hilbertF(z) ```

`hilbert(S::Space{<:PeriodicLine},f,z)` calls domain`(f::AbstractVector)`

1

On line 14 of [periodicline.jl](https://github.com/JuliaApproximation/SingularIntegralEquations.jl/blob/master/src/periodicline.jl) there appears to be a typo, leading to a call to domain(f::AbstractVector) which is not defined. ``` function hilbert(S::Space{

It looks like there is no rule for `LinRange` ``` julia> using ChainRules julia> isnothing(ChainRules.rrule(LinRange, 1, 2, 10)) true ```

Currently (Distances v0.10.7, Julia 1.8.5) testing equality between metrics parametrised by vectors returns `false` when comparing two instances, e.g. ``` julia> m = PeriodicEuclidean(ones(1)) PeriodicEuclidean{Vector{Float64}}([1.0]) julia> m == deepcopy(m) false...

Independent mokernel failing numerical equality test randomly

2

Testing locally on 1.8.5, I find that https://github.com/JuliaGaussianProcesses/KernelFunctions.jl/blob/master/test/mokernels/independent.jl#L17 occasionally fails, ```Julia independent: Test Failed at .../KernelFunctions.jl/test/mokernels/independent.jl:17 Expression: kernelmatrix(k, x, y) == kernelmatrix(k, collect(x), collect(y)) Evaluated: [0.6308251776147069 0.7710805544766418 … 0.0 0.0;...

Define `isequal` and `==` more widely between `Kernel`s

16

It would be nice to define `isequal` and/or `==` as widely as possible between `Kernel`s and `Transform`s. At the moment some basic cases seem to be missing, such as ```Julia...

Error in `LinearAlgebra.copy_oftype` on addition of symmetric `KeyedArray` and `UniformScaling`

2

In Julia 1.8.2 ```Julia julia> using AxisKeys, LinearAlgebra julia> A = Symmetric(rand(2,2)) 2×2 Symmetric{Float64, Matrix{Float64}}: 0.358751 0.612696 0.612696 0.236462 julia> K = KeyedArray(A, ([:a,:b], [:c, :d])) 2-dimensional KeyedArray(...) with keys:...

Closes #133

Currently `tr` is computed block by block, which works when all blocks are square. For consistency with `LinearAlgebra` this should probably work as long as the total matrix is square....