At this point it is more flexible using a second template parameter. We can choose float* or DenseVector or ... The only thing we have to implement for other vector types is the triangular solver things.