CPPLapack Code

//=============================================================================
/*!  */
template<long n>
inline const dsymatrix_small<n>& operator+(const dsymatrix_small<n>& A)
{
  return A;
}

//=============================================================================
/*!  */
template<long n>
inline dsymatrix_small<n> operator-(const dsymatrix_small<n>& A)
{
  dsymatrix_small<n> X;
  for(long i=0; i<n; i++){
    for(long j=0; j<=i; j++){
      X(i,j)=-A(i,j);
    }
  }
  return X;
}