Skip to content

Significantly speed up superoperators

Vicentini Filippo requested to merge pv/speedupmylindblad into master

Created by: PhilipVinc

Super-operators use the op.max_conn_size to estimate the maximum size of super operators.

Right now this size is over-estimated because the diagonal terms of all operators only contribute 1 connected element (since they don't flip any spin), not n_operators. I think originally this was not how it was implemented, but at some point I made it so that all diagonal contributions are stored in the same connected element but evidently forgot to update this max_conn_size accordingly.

(Tests already check that this value is correct, so if tests pass, it's good)

Merge request reports