CPPLapack Code

//=============================================================================
template<long n>
inline void idamax(long& I, long& J, const dsymatrix_small<n>& A)
{
  double max(-1.);
  for(int i=0; i<n; i++){
    for(int j=0; j<=i; j++){
      if( max<fabs(A(i,j)) ){
        I=i;
        J=j;
        max =fabs(A(i,j));
      }
    }
  }
  return;  
}

//=============================================================================
template<long n>
inline double damax(const dsymatrix_small<n>& A)
{
  long i(0),j(0);
  idamax(i,j,A);
  return A(i,j);
}