src/Pulses/PulComposite.cc
// ____________________________________________________________________________
// vii             CLASS COMPOSITE PULSE STEADY-STATE HANDLING
// ____________________________________________________________________________
vii             CLASS COMPOSITE PULSE STEADY-STATE HANDLING
vii CLASS COMPOSITE PULSE STEADY - STATE HANDLING
void PulComposite::SetSSs()
  {
  if(!WFsteps) return;			
  if(!CheckH(10,15)) CPulfatality(16);	
  if(!R.size()) return;			
  if(!SigSSs)                           
    SigSSs = new densop[WFsteps];	
  super_op Leff;
  for(int i=0; i<WFsteps; i++)
    {
    if(i == Hindex[i])			
      {					
      Leff = complexi*Lsteps[i] + R;	
      SigSSs[i]=densop(Leff,R,sigmaeq);	
      }
    }
  }
void PulComposite :: SetSSs ( ) { if ( ! WFsteps ) return ; if ( ! CheckH ( 10 , 15 ) ) CPulfatality ( 16 ) ; if ( ! R . size ( ) ) return ; if ( ! SigSSs ) SigSSs = new densop [ WFsteps ] ; super_op Leff ; for ( int i = 0 ; i < WFsteps ; i ++ ) { if ( i == Hindex [ i ] ) { Leff = complexi * Lsteps [ i ] + R ; SigSSs [ i ] = densop ( Leff , R , sigmaeq ) ; } } }
____________________________________________________________________________
vii             CLASS COMPOSITE PULSE STEADY-STATE HANDLING
____________________________________________________________________________
____________________________________________________________________________ vii CLASS COMPOSITE PULSE STEADY - STATE HANDLING ____________________________________________________________________________
// Input                CPul	: A composite pulse (this) // Output               void	: Steady-state operators for each //				  step of the composite pulse are //			   	  calculated // Exit if NO steps // Insure Hamiltonians exist // No can do if no relaxation set // If propagators don't yet exist // allocate some space for them // If unique Hamiltonian  // then calculate unique ss // Full Liouvillian, ith step // Steady state dens.op. ith step