src/LSLib/LSprop.cc
/* These propagators may or may not have equivalent Hilbert space counterparts.
   Density operator evolution occurs according either to

                            -iLt                        
                           e     |sigma>  ==> G |sigma> 

   if the Liouvillian acts directly on the density operator or as


   if the inputLouvillian acts on the difference density operator.
   Note that in these functions the Liouvillian should have the bulk of its
   frequency components as imaginary terms and the bulk of its relaxation and
   exchange components are real terms, e.g. L = -i[H, ] + R + X.  See the
   documentation (theory section) for specifics as to how the superoparator
   propagator is formed when the system evolves to a steady state.

                                Leff	: Active Liouvillian (1/sec)
	  			sigmass	: Steady state density operator
                                tevol	: Evolution time (seconds)
           Output               LSprop	: Propagator for input Liouvillian  */
These propagators may or may not have equivalent Hilbert space counterparts.
Density operator evolution occurs according either to

iLt
e     |sigma>  ==> G |sigma>

if the Liouvillian acts directly on the density operator or as

if the inputLouvillian acts on the difference density operator.
Note that in these functions the Liouvillian should have the bulk of its
frequency components as imaginary terms and the bulk of its relaxation and
exchange components are real terms, e.g.

Leff    : Active Liouvillian (1/sec)
sigmass : Steady state density operator
tevol   : Evolution time (seconds)
Output               LSprop  : Propagator for input Liouvillian
These propagators may or may not have equivalent Hilbert space counterparts . Density operator evolution occurs according either to iLt e |sigma > == > G |sigma > if the Liouvillian acts directly on the density operator or as if the inputLouvillian acts on the difference density operator . Note that in these functions the Liouvillian should have the bulk of its frequency components as imaginary terms and the bulk of its relaxation and exchange components are real terms e . g . Leff : Active Liouvillian ( 1 / sec ) sigmass : Steady state density operator tevol : Evolution time ( seconds ) Output LSprop : Propagator for input Liouvillian
LSprop::LSprop(const super_op& Leff, double tevol)
  {
  GOp =	exp(Leff, -tevol);		
  Gt  = tevol;				
  }
LSprop :: LSprop ( const super_op & Leff , double tevol ) { GOp = exp ( Leff , - tevol ) ; Gt = tevol ; }
These propagators may or may not have equivalent Hilbert space counterparts.
These propagators may or may not have equivalent Hilbert space counterparts .
// Exponential Liouvillian // Set our evolution time