You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
When setting thermalmode(n,mode), the Hilbert space of the mode is usually truncated when it no longer meaningfully contributes.
However, since the function is defined to have tr(thermalmode(mode,\bar{n})) = 1 for any Hilbert space dimension, rounding error that runs against machine precision begins to occur for large N Hilbert space. There should be a built-in mechanism that catches this and removes rounding error and sets higher mode number to zero.
The following plot represents the RMSE of the solution to timeevolution.schrodinger_dynamic between iterations for a density matrix with thermal mode n = 0.5 and increasing Hilbert space dimension.
The text was updated successfully, but these errors were encountered:
I should note that convergence for say n = 5 under the same setup (only a different thermal mode) had a clear convergence pattern. It is possible that the code isn't properly handling specifically non-integer values.
When setting
thermalmode(n,mode)
, the Hilbert space of the mode is usually truncated when it no longer meaningfully contributes.However, since the function is defined to have
tr(thermalmode(mode,\bar{n})) = 1
for any Hilbert space dimension, rounding error that runs against machine precision begins to occur for large N Hilbert space. There should be a built-in mechanism that catches this and removes rounding error and sets higher mode number to zero.The following plot represents the RMSE of the solution to
timeevolution.schrodinger_dynamic
between iterations for a density matrix with thermal mode n = 0.5 and increasing Hilbert space dimension.The text was updated successfully, but these errors were encountered: