kadir-gunel / xles.jl Goto Github PK

function mapOrthogonal(X::T, Y::T; λ::Float32=Float32(1)) where {T}
    F = CUDA.CUBLAS.svd(X * Y')
    W = permutedims(F.U * F.Vt * cuinv((X * X') + λ .* CuMatrix{Float32}(I, 300, 300)))
    return W, F.S
end

The above code snippet is kind of misleading. Thinking in theory lead me to think that the matrix multiplication inside the decomposition operation X * Y' is wrong. Because what we do is that we rotate the source space towards the target space where as this operation leads to opposite direction!

We should calculate the rotation operation based on X and not Y!

To do so :

1. Just transpose the operation : X * Y' to Y * X'.
2. Probably, we need to shift the Vt and U operations.
3. Debug for verification.