src/IntRank2/CartMx2A.cc
/* Given a possible row in a 3x3 Jacobian matrix, <i|a|0>, <i|a|1>, & <i|a|2>,
   this function does a very simple check to insure the row will not produce
   a singlular Jacobian matrix.  It does so by checking that 1.) row i is not
   zero and 2.) row j where j<i is not a constant multiple of row i. The
   latter check is done for all rows j where j<i.                            */
Given a possible row in a 3x3 Jacobian matrix, , , & ,
this function does a very simple check to insure the row will not produce
a singlular Jacobian matrix.  It does so by checking that 1.) row i is not
zero and 2.) row j where j<i is not a constant multiple of row i. The
latter check is done for all rows j where j<i.
Given a possible row in a 3x3 Jacobian matrix & this function does a very simple check to insure the row will not produce a singlular Jacobian matrix . It does so by checking that 1 . ) row i is not zero and 2 . ) row j where j<i is not a constant multiple of row i . The latter check is done for all rows j where j<i .
bool CartMx2A::CheckSing(double ai0, double ai1, double ai2, int i)
  {
  if(ai0==0 && J.getRe(i,0) != 0) return false;
  if(ai0!=0 && J.getRe(i,0) == 0) return false;
  double sf = ai0/J.getRe(i,0);
  if(fabs(ai1/J.getRe(i,1) - sf) > 1.e-10) return false;
  if(fabs(ai2/J.getRe(i,2) - sf) > 1.e-10) return false;
  return true;
  }
bool CartMx2A :: CheckSing ( double ai0 , double ai1 , double ai2 , int i ) { if ( ai0 == 0 && J . getRe ( i , 0 ) != 0 ) return false ; if ( ai0 != 0 && J . getRe ( i , 0 ) == 0 ) return false ; double sf = ai0 / J . getRe ( i , 0 ) ; if ( fabs ( ai1 / J . getRe ( i , 1 ) - sf ) > 1.e-10 ) return false ; if ( fabs ( ai2 / J . getRe ( i , 2 ) - sf ) > 1.e-10 ) return false ; return true ; }
Given a possible row in a 3x3 Jacobian matrix, <i|a|0>, <i|a|1>, & <i|a|2>,
this function does a very simple check to insure the row will not produce
a singlular Jacobian matrix.
Given a possible row in a 3x3 Jacobian matrix <i|a|0 > <i|a|1 > & <i|a|2 > this function does a very simple check to insure the row will not produce a singlular Jacobian matrix .