da_mesolve

Dimension adaptive algorithm for qutip.mesolve used to solve infinite-dimensional problems. Do NOT use this for finite-dimensional problems such as qubits.

How to use

The usage is similar to qutip.mesolve. However, the Hamiltonian H, the collapse operatos c_ops, and the expectation value operators e_ops must be passed in as functions of the Hilbert space dimension N.

---

---

CHANGELOG

[0.0.5]

• Fixed bug in code.da_mesolve where match_output_N does not work when e_ops_func is not specified.
• Fixed the logic for utils.highest_occupied_level.

[0.0.4]

• Added description to utils.highest_occupied_level.
• Removed N_init from core.da_mesolve. The algorithm starts from the input rho0.
• The first element in rho_t of core.da_mesolve is changed from rho0 to its adapted version. If the highest level of rho0 is occupied, an error is raised. If not, then a truncation is done according to N_keep.
• match_output_N of core.da_mesolve is now False by default.
• Added the args argument to core.da_mesolve.

[0.0.3]

• Changed **qutip_kwargs** to options to be compatible with qutip. args, meanwhile, is not yet compatible with this package.

[0.0.2]

• Moved match_N to be a user-level function adapt_rho in utils.py.
• Added error handling in utils.highest_occupied_level if the highest occupied level is the highest level in the density matrix.

[0.0.1]

• Added the algorithm.