src/Level1/SpinT.cc
// Output		SOp   : Spin Operator for SphT Component T
//								  0m 
Output         SOp   : Spin Operator for SphT Component T
0m
Output SOp : Spin Operator for SphT Component T 0m
spin_op T20SS(const spin_sys& sys, int spin, const coord& vect, int m, int rev)
  {
  spin_op SOp;
  complex vp,vm;
  double vz;
  vp = complex(vect.x(),vect.y());	
  vm = complex(vect.x(),-vect.y());	
  vz = vect.z();			
  switch (m)
    {
    case 0:				
      SOp = -1.0/sqrt(3.0)		
	     *(Iz(sys,spin)*vect.z() 	
              + 0.5*(Im(sys,spin)*vm
               + Ip(sys,spin)*vp));
      break;
    default:
      spin_T_error(3);
      std::cout << " 0," << m << "\n";
      break;
    }
  return SOp;				
  rev = 0;				
  }
spin_op T20SS ( const spin_sys & sys , int spin , const coord & vect , int m , int rev ) { spin_op SOp ; complex vp , vm ; double vz ; vp = complex ( vect . x ( ) , vect . y ( ) ) ; vm = complex ( vect . x ( ) , - vect . y ( ) ) ; vz = vect . z ( ) ; switch ( m ) { case 0 : SOp = -1.0 / sqrt ( 3.0 ) * ( Iz ( sys , spin ) * vect . z ( ) + 0.5 * ( Im ( sys , spin ) * vm + Ip ( sys , spin ) * vp ) ) ; break ; default : spin_T_error ( 3 ) ; std :: cout << " " << m << " \n " ; break ; } return SOp ; rev = 0 ; }
Output         SOp   : Spin Operator for SphT Component T
0m
Output SOp : Spin Operator for SphT Component T 0m
// u+ = ux + iuy // u- = ux - iuy // uz //  (2)            [1]       1 [          ] // T   (i) = - sqrt|-| I u + - |I u + I u | //  00             [3]  z z  2 [ - -   + +] // Return SOp // Compiler likes this to be used