src/Matrix/h_matrix.cc
// sosiz - the above routine is a kludge, but it worked OK. A permutation 
//         on hmx will produce hmx but the code to do so is more difficult
//          within h_matrix.  So I just spawned it to normal matrix then
//          cast it back to hermitian.  This uses extra arrays ......
//         I've since fixed the routines to just stay Hermitian and the
//         above can be deleted when I'm convinced they work fine.
the above routine is a kludge, but it worked OK. A permutation
on hmx will produce hmx but the code to do so is more difficult
within h_matrix.  So I just spawned it to normal matrix then
cast it back to hermitian.  This uses extra arrays
I've since fixed the routines to just stay Hermitian and the
above can be deleted when I'm convinced they work fine.
the above routine is a kludge but it worked OK . A permutation on hmx will produce hmx but the code to do so is more difficult within h_matrix . So I just spawned it to normal matrix then cast it back to hermitian . This uses extra arrays I ' ve since fixed the routines to just stay Hermitian and the above can be deleted when I ' m convinced they work fine .
double h_matrix::maxRe() const
  { 
  double maxval=-HUGE_VAL;
  for(int i=0; i<size; i++) { maxval = gmax(Re(data[i]), maxval); }
  return (size)?maxval:0;  
  }
double h_matrix :: maxRe ( ) const { double maxval = - HUGE_VAL ; for ( int i = 0 ; i < size ; i ++ ) { maxval = gmax ( Re ( data [ i ] ) , maxval ) ; } return ( size ) ? maxval : 0 ; }
sosiz - the above routine is a kludge, but it worked OK. A permutation
on hmx will produce hmx but the code to do so is more difficult
within h_matrix.
sosiz - the above routine is a kludge but it worked OK . A permutation on hmx will produce hmx but the code to do so is more difficult within h_matrix .