Is your feature request related to a problem? Please describe.
It would be great to allow for complex-valued states/operators/snapshot matrices.
Describe the solution you'd like
Allow integration of complex-valued ODEs + solve the least squares problem using complex-valued snapshot matrices. I believe this could be done easily by changing the dtype in certain definitions to allow complex128.