Update open systems simulator to use ndarray::linalg::kron.

[cgranade]

Please describe what you would like the maintenance work to be.
The open systems simulator uses its own implementation of the flattened tensor product (the Kronecker product), but ndarray has added a new ndarray::linalg::kron function as of 0.15.4. Using that function in favor of our own implementation decreases the API surface that needs to be tested and improves performance.

Additional context
See rust-ndarray/ndarray#1105 for further discussion of the ndarray::linalg::kron feature.