need to improve the matrix packing
strides for packing, no mallocs in solver loop