src/BWRRelax/relaxQCSA.cc
// ------------ CSA-Quadrupole, With An Applied RF-Field ----------------
CSA-Quadrupole, With An Applied RF-Field
CSA - Quadrupole With An Applied RF - Field
super_op RCQrf(const sys_dynamic& sys, gen_op& Heff, double Wrflab, int level)
   {
   int hs = sys.HS();				
   int ls = hs*hs;				
   Heff.set_EBR();				
   matrix mx(ls, ls, 0.0);			
   super_op LOp(mx, Heff.get_basis());
   double taus[5];				
   taust(taus, sys.taus());
   double chi = chit(sys.taus());		
   double *w;				
   w = new double[hs];
   if(abs(level) > 1)				
     Heff.eigvals(w);
   RCQrf(LOp,sys,Heff,w,Wrflab,taus,chi,level);	
   delete [] w;
   return LOp;
   }
super_op RCQrf ( const sys_dynamic & sys , gen_op & Heff , double Wrflab , int level ) { int hs = sys . HS ( ) ; int ls = hs * hs ; Heff . set_EBR ( ) ; matrix mx ( ls , ls , 0.0 ) ; super_op LOp ( mx , Heff . get_basis ( ) ) ; double taus [ 5 ] ; taust ( taus , sys . taus ( ) ) ; double chi = chit ( sys . taus ( ) ) ; double * w ; w = new double [ hs ] ; if ( abs ( level ) > 1 ) Heff . eigvals ( w ) ; RCQrf ( LOp , sys , Heff , w , Wrflab , taus , chi , level ) ; delete [ ] w ; return LOp ; }
CSA-Quadrupole, With An Applied RF-Field
CSA - Quadrupole With An Applied RF - Field
// Input		sys   : Dynamic spin system //			Heff  : Effective Hamiltonian (Hertz) //			Wrflab: RF-Field frequency, lab frame (Hertz) //			level : Relaxation treatment level // Output		LOp   : CSA-Quadrupole relaxation superoperator // Note			      :	Computed in the eigenbasis of Heff // Total system Hilbert space // Total system Liouville space // Insure Heff is in its eigenbasis (Hz) // Construct zero superoperator //		    Prepare the Spectral Density Components // Get 5 taus for spectral densities // Get the system chi value // Set up for system energy levels (LAB) // Needed for higher level computations //	         Get the CSA-Quadrupolar Relaxation Superoperator // Set LOp to CSA-Quad R superop