Fix segfault
Add direct solvers test; failures for vector ranges and slices to investigate