Fixed bugs

Author: theogf

src/kernels/kernelproduct.jl

@@ -29,23 +29,26 @@ hadamard(x,y) = x.*y
 function kernelmatrix(
     κ::KernelProduct,
     X::AbstractMatrix;
-    obsdim::Int=defaultobs)
-    reduce(hadamard,kernelmatrix(κ.kernels[i],X,obsdim=obsdim) for i in 1:length(κ))
+    obsdim::Int=defaultobs,
+)
+    reduce(hadamard, kernelmatrix(κ.kernels[i], X, obsdim = obsdim) for i in 1:length(κ))
 end
 
 function kernelmatrix(
     κ::KernelProduct,
     X::AbstractMatrix,
     Y::AbstractMatrix;
-    obsdim::Int=defaultobs)
-    reduce(hadamard,_kernelmatrix(κ.kernels[i],X,Y,obsdim) for i in 1:length(κ))
+    obsdim::Int=defaultobs,
+)
+    reduce(hadamard, kernelmatrix(κ.kernels[i], X, Y, obsdim = obsdim) for i in 1:length(κ))
 end
 
 function kerneldiagmatrix(
     κ::KernelProduct,
     X::AbstractMatrix;
-    obsdim::Int=defaultobs) #TODO Add test
-    reduce(hadamard,kerneldiagmatrix(κ.kernels[i],X,obsdim=obsdim) for i in 1:length(κ))
+    obsdim::Int=defaultobs,
+) #TODO Add test
+    reduce(hadamard, kerneldiagmatrix(κ.kernels[i], X, obsdim = obsdim) for i in 1:length(κ))
 end
 
 function Base.show(io::IO, κ::KernelProduct)

src/kernels/kernelsum.jl

@@ -58,7 +58,7 @@ function kernelmatrix(
     Y::AbstractMatrix;
     obsdim::Int = defaultobs,
 )
-    sum(κ.weights[i] * _kernelmatrix(κ.kernels[i], X, Y, obsdim) for i in 1:length(κ))
+    sum(κ.weights[i] * kernelmatrix(κ.kernels[i], X, Y, obsdim = obsdim) for i in 1:length(κ))
 end
 
 function kerneldiagmatrix(

src/matrix/kernelmatrix.jl

@@ -24,15 +24,15 @@ function kernelmatrix!(
     K::AbstractMatrix,
     κ::Kernel,
     X::AbstractMatrix;
-    obsdim::Int = defaultobs
+    obsdim::Int = defaultobs,
 )
     return kernelmatrix!(K, κ, vec_of_vecs(X, obsdim = obsdim))
 end
 
 function kernelmatrix!(
     K::AbstractMatrix,
     κ::Kernel,
-    X::AbstractVector
+    X::AbstractVector,
     )
     if !check_dims(K, X, X)
         throw(DimensionMismatch("Dimensions of the target array K $(size(K)) are not consistent with X $(size(X))"))
@@ -60,7 +60,7 @@ function kernelmatrix!(
     κ::Kernel,
     X::AbstractMatrix,
     Y::AbstractMatrix;
-    obsdim::Int = defaultobs
+    obsdim::Int = defaultobs,
 )
     return kernelmatrix!(K, κ, vec_of_vecs(X, obsdim = obsdim), vec_of_vecs(Y, obsdim = obsdim))
 end
@@ -69,7 +69,7 @@ function kernelmatrix!(
     K::AbstractMatrix,
     κ::Kernel,
     X::AbstractVector,
-    Y::AbstractVector
+    Y::AbstractVector,
     )
     if !check_dims(K, X, Y)
         throw(DimensionMismatch("Dimensions of the target array K $(size(K)) are not consistent with X $(size(X)) and Y $(size(Y))"))
@@ -130,7 +130,7 @@ kerneldiagmatrix
 function kerneldiagmatrix(
     κ::Kernel,
     X::AbstractMatrix;
-    obsdim::Int = defaultobs
+    obsdim::Int = defaultobs,
     )
     return kerneldiagmatrix(κ, vec_of_vecs(X, obsdim = obsdim))
 end
@@ -146,7 +146,7 @@ function kerneldiagmatrix!(
     K::AbstractVector,
     κ::Kernel,
     X::AbstractMatrix;
-    obsdim::Int = defaultobs
+    obsdim::Int = defaultobs,
 )
     if length(K) != size(X,obsdim)
         throw(DimensionMismatch("Dimensions of the target array K $(size(K)) are not consistent with X $(size(X))"))
@@ -157,7 +157,7 @@ end
 function kerneldiagmatrix!(
     K::AbstractVector,
     κ::Kernel,
-    X::AbstractVector
+    X::AbstractVector,
 )
     if length(K) != length(X)
         throw(DimensionMismatch("Dimensions of the target array K $(size(K)) are not consistent with X $(length(X))"))

src/utils.jl

@@ -53,7 +53,7 @@ end
 
 Base.size(D::RowVecs) = (size(D.X, 1),)
 Base.getindex(D::RowVecs, i::Int) = view(D.X, i, :)
-Base.getindex(D::ColVecs, i::CartesianIndex{1}) = view(D.X, i, :)
+Base.getindex(D::RowVecs, i::CartesianIndex{1}) = view(D.X, i, :)
 Base.getindex(D::RowVecs, i) = RowVecs(view(D.X, i, :))
 
 # Take highest Float among possibilities