Why the current interface in MatCreateVecs requires that a matrix has been already either allocated or setup ? Shouldn’t be enough to check if the layout has been already setup? Stefano