Support multi-vector b

```julia
LLDL = lldl(A, memory = 10)
x = LLDL \ b
```
currently only vector b is supported. I am wondering if multi-vector (b as a matrix) could be supported (so that x will be 2d as well). A naive multi-threading seems to have too much overhead. 