Implement MPSMatrixSolve

[vboussange]

Hey there,

using Metal
A = Metal.randn(100,100)
b = Metal.randn(100,1)
A\b

throws

ERROR: MethodError: no method matching getrs!(::Char, ::MtlMatrix{Float32, Private}, ::MtlVector{UInt32, Private}, ::Matrix{Float32})

Any idea on what is happening?

[tgymnich]

Solving of linear equations is not yet implemented in Metal.jl (see #145). However you can already use LU factorisation in the meantime.

using LinearAlgebra
A = Metal.randn(100,100)
lua = lu(A)
# ...

[cortner]

How would you use the lu factorisation if back-substitution is not implemented? I would also be extremely grateful if the full Metal linear algebra suite (small as it is) could be interfaced.