Calculating atomic absorption spectra using steadystate

[Gal Winer]

Hi! 

I'm trying to use the steadystate function to obtain the imaginary part of the steady state density matrix for a range of parameters. 

My approach is as follows: 

1. A function, makeH(OmegaP,delta), returns the Hamiltonian for the range of parameters (e.g laser power and detuning) 

2. A function make calcAbsorption(H,Cops)  calls the steadystate function on passed H and collapse operator list, extracts the correct element from the density matrix and performs some algebraic operations on it. 

3. call these function for a grid of parameters. and use the output. 

I'm trying to speed this calculation up as much as possible as I will have to run it on large parameter grids. 

Is my approach reasonable? Do qutip object support symbolic expressions via Sympy? Any tricks I can use to speed things up? 

thanks for the help, 

Gal 

[Eric]

Hi

qutip does not support symbolic expression.

The best help I can give would be to test different methods.

Iterative method can take a "preconditioner". Maybe your Hamiltonian are sufficiently similar that you can reuse the preconditioner matrix.

http://qutip.org/docs/latest/apidoc/functions.html#module-qutip.steadystate

http://qutip.org/docs/latest/apidoc/functions.html#qutip.steadystate.build_preconditioner

Eric