src/IntRank2/IntRank2A.cc
/* Given the three PAS Cartesian components of an irreducible rank 2 spatial
   tensor function AxAyAz will determine the isotropic, anisotropic, and 
   asymmetry terms.  These will be returned in a single coordinate.          */
Given the three PAS Cartesian components of an irreducible rank 2 spatial
tensor function AxAyAz will determine the isotropic, anisotropic, and
asymmetry terms.  These will be returned in a single coordinate.
Given the three PAS Cartesian components of an irreducible rank 2 spatial tensor function AxAyAz will determine the isotropic anisotropic and asymmetry terms . These will be returned in a single coordinate .
coord IntRank2A::AisoDelzEta(const coord& AxAyAz)	
  {
  double xx   = AxAyAz.x();			
  double yy   = AxAyAz.y();			
  double zz   = AxAyAz.z();			
  double AisoVal = (1.0/3.0)*(xx+yy+zz);	
  if(fabs(AisoVal) < 1.e-12) AisoVal = 0.0;	
  xx -= AisoVal;				
  yy -= AisoVal;				
  zz -= AisoVal;
  SortAxAyAz(xx, yy, zz);			
  double EtaVal = (xx-yy)/(zz);			
  if(fabs(EtaVal) < 1.e-7)      EtaVal = 0.0;	
  if(fabs(EtaVal - 1.) < 1.e-7) EtaVal = 1.0;	
  double PASDelzz = zz - 0.5*(xx+yy);		
  PASDelzz *= (2./3.);                          
  if(fabs(PASDelzz) < 1.e-12) PASDelzz = 0.0;   
  if(EtaVal == 1) PASDelzz = fabs(PASDelzz);    
  return coord(AisoVal,PASDelzz,EtaVal); 	
  }
coord IntRank2A :: AisoDelzEta ( const coord & AxAyAz ) { double xx = AxAyAz . x ( ) ; double yy = AxAyAz . y ( ) ; double zz = AxAyAz . z ( ) ; double AisoVal = ( 1.0 / 3.0 ) * ( xx + yy + zz ) ; if ( fabs ( AisoVal ) < 1.e-12 ) AisoVal = 0.0 ; xx -= AisoVal ; yy -= AisoVal ; zz -= AisoVal ; SortAxAyAz ( xx , yy , zz ) ; double EtaVal = ( xx - yy ) / ( zz ) ; if ( fabs ( EtaVal ) < 1.e-7 ) EtaVal = 0.0 ; if ( fabs ( EtaVal - 1. ) < 1.e-7 ) EtaVal = 1.0 ; double PASDelzz = zz - 0.5 * ( xx + yy ) ; PASDelzz *= ( 2. / 3. ) ; if ( fabs ( PASDelzz ) < 1.e-12 ) PASDelzz = 0.0 ; if ( EtaVal == 1 ) PASDelzz = fabs ( PASDelzz ) ; return coord ( AisoVal , PASDelzz , EtaVal ) ; }
Given the three PAS Cartesian components of an irreducible rank 2 spatial
tensor function AxAyAz will determine the isotropic, anisotropic, and
asymmetry terms.
Given the three PAS Cartesian components of an irreducible rank 2 spatial tensor function AxAyAz will determine the isotropic anisotropic and asymmetry terms .
// Static, not constant // Get Axx (PAS) // Get Ayy (PAS) // Get Azz (PAS) // Get isotropic component // In sure zero if below delcut // Remove isotropic component // before sorting elements // Keep |Azz| >= |Ayy| >= |Axx| // Calculate the asymmetry // -> Small values set zero // -> Set exactly 1 if close // Calculate delzz value //  THIS IS NOT GAMMA SCALED //  (Small values set zero) // Use + delzz if eta is 1 //   (where Ayy = - Azz)