How to make division for SparceMatrices?

In my below MWE, I am trying to use \ for two Sparce matrices, but I am having an error. Any way to solve it?

using SparseArrays;
using LinearAlgebra;
using .Threads
using BenchmarkTools
using StaticArrays
A = [1 3 3; 6 9 7;-3 3 4]; As = sparse([1 3 3; 6 9 7;-3 3 4]); b = [1,2,3]; bs = sparse([1,2,3]);
E = As \ bs;

ERROR: LoadError: MethodError: no method matching ldiv!(::SuiteSparse.UMFPACK.UmfpackLU{Float64, Int64}, ::SparseVector{Float64, Int64})
Closest candidates are:
  ldiv!(::LU{T, Tridiagonal{T, V}}, ::AbstractVecOrMat{T} where T) where {T, V} at C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\LinearAlgebra\src\lu.jl:583
  ldiv!(::UpperHessenberg, ::AbstractVecOrMat{T} where T; shift) at C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\LinearAlgebra\src\hessenberg.jl:127    
  ldiv!(::LDLt{var"#s832", var"#s831"} where {var"#s832", var"#s831"<:SymTridiagonal}, ::AbstractVecOrMat{T} where T) at C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\LinearAlgebra\src\ldlt.jl:169

Only the matrix should be sparse, not the vector.

Is there any reason for that?

Is there any trick to change the vector b to matrix of single dimension?