Hbar in qutip

[Real Solace]

Hello, I think there might have been a few issues on this in the past however in this case I would like to know whether it is possible to use the hbar constant in its full form and not equal to 1?

I am seeing the dynamics which I want to see when using qutip but the values of my position and momentum expectation are up to ten orders of magnitude higher than they are supposed to be, and if i used hbar = 6.7e-16 the dynamics dissapear.

I have tried qutip's autotidyup = False, but if theres any other way to handle the way qutip calculates very small values I would be grateful/

Kind Regards

Sahra

[Andrew M. C. Dawes]

The fundamental issue is representing small numbers in a digital format (float, double etc). These limits barely cover values of hbar, let alone anything smaller, so it’s not really optional to do anything other than re-scale so hbar=1. That said, you should also be able to rescale your expected results in order to compare to the numerical result. If you post more details here, we may be able to help more specifically.

Andy

--
You received this message because you are subscribed to the Google Groups "QuTiP: Quantum Toolbox in Python" group.
To unsubscribe from this group and stop receiving emails from it, send an email to qutip+un...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/qutip/24cd4e1c-17a4-4e9f-b322-29c759fe564fn%40googlegroups.com.

[Real Solace]

Thanks, I've re-scaled and they look as expected so far. I will post more details once I decide exactly what I want.