src/BWRRelax/relaxQuad.cc
// --------------- Quad-Quad, With An Applied RF-Field ------------------
Quad-Quad, With An Applied RF-Field
Quad - Quad With An Applied RF - Field
void RQQrf(super_op& LOp, const sys_dynamic& sys, gen_op& Heff, double*w,
                 double Wrflab, double* taus, double chi, int type, int level)
   {
   matrix Xis = xiQ(sys);			
   int ns = sys.spins();			
   spin_T *T;					
   T = new spin_T[ns];				
   int i;
   for(i=0; i<ns; i++)				
     if(Re(Xis.get(i,i)))
       T[i] = T_Q(sys,i);
   space_T *A; 					
   A = new space_T[ns];				
   for(i=0; i<ns; i++)				
     if(Re(Xis.get(i,i)))
       A[i] = sys.TQ(i);
   Rrfij(LOp, sys, Heff, w, Wrflab, Xis, Xis, A, A,	
                         T, T, taus, chi, type, level);
   return;
   }
void RQQrf ( super_op & LOp , const sys_dynamic & sys , gen_op & Heff , double * w , double Wrflab , double * taus , double chi , int type , int level ) { matrix Xis = xiQ ( sys ) ; int ns = sys . spins ( ) ; spin_T * T ; T = new spin_T [ ns ] ; int i ; for ( i = 0 ; i < ns ; i ++ ) if ( Re ( Xis . get ( i , i ) ) ) T [ i ] = T_Q ( sys , i ) ; space_T * A ; A = new space_T [ ns ] ; for ( i = 0 ; i < ns ; i ++ ) if ( Re ( Xis . get ( i , i ) ) ) A [ i ] = sys . TQ ( i ) ; Rrfij ( LOp , sys , Heff , w , Wrflab , Xis , Xis , A , A , T , T , taus , chi , type , level ) ; return ; }
Quad-Quad, With An Applied RF-Field
Quad - Quad With An Applied RF - Field
// Input                LOp   : Superoperator // 			sys   : Dynamic spin system //			Heff  : Effective Hamiltonian //			Wrflab: Field freqency in the lab frame (Hz) //			w     : Transition frequencies //			taus  : Sytem 5 effective correlation times //			chi   : Spin system chi value // 			type  : Relaxation type to compute //				   0 = Quad-Quad AC & CC //				   + = Quad-Quad AC //				   - = Quad-Quad CC //			level : Relaxation treatment level // Output		LOp   : Quad relaxation superoperator // Note			      :	Computed in the eigenbasis of Ho //		    Prepare the Interaction Constants // Quadrupolar interaction constants //			Prepare the Spin Tensors // Total number of spins // Spin tensors for each spin // Compiler didn't like spin_T T[ns] // Set them all to Quad spin tensors //			Prepare the Spatial Tensors // Spatial tensors for each spin // Compiler didn't like space_T A[ns] // Set them all to Quad space tensors //		    Perform Computation Over All Spins // All Quad parameters