// array.c

#ifndef array_C
#define array_C

#include "trfutil/TRFMath.h"

#ifndef DEFECT_CMATH_NOT_STD
using std::abs;
#endif

//**********************************************************************
// Member functions.
//**********************************************************************

// Return the smallest element.
template <class T>
T array<T>::min( ) const {
  T t = data[0];
  for ( int i=1; i<size; ++i )
    if ( data[i] < t ) t = data[i];
  return t;
}

//**********************************************************************

// Return the largest element.
template <class T>
T array<T>::max( ) const {
  T t = data[0];
  for ( int i=1; i<size; ++i )
    if ( data[i] > t ) t = data[i];
  return t;
}

//**********************************************************************

// Return the absolute smallest element.
template <class T>
T array<T>::amin( ) const {
  T t = abs(data[0]);
  for ( int i=1; i<size; ++i )
    if ( abs(data[i]) < t ) t = abs(data[i]);
  return t;
}

//**********************************************************************

// Return the absolute largest element.
// Useful for checking "equality".
template <class T>
T array<T>::amax( ) const {
  T t = abs(data[0]);
  for ( int i=1; i<size; ++i )
    if ( abs(data[i]) > t ) t = abs(data[i]);
  return t;
}

//**********************************************************************

// difference

template <class T>
array<T> array<T>::operator-( ) {
  array<T> arr(size);
  for ( int i=0; i<arr.size; i++ ) arr.data[i] = -data[i];
  return arr;
}

//**********************************************************************
// Friends.
//**********************************************************************

// output stream

template <class T>
std::ostream& operator<<(std::ostream& stream, const array<T>& arr)
{
  arr.ostr(stream);
  return stream;
}

//**********************************************************************

// equality

template <class T>
bool operator==(const array<T>& arr1, const array<T>& arr2) {
  if ( arr1.size != arr2.size ) return false;
  for ( int i=0; i<arr1.size; i++ ) 
    if ( arr1(i) != arr2(i) ) return false;
  return true;
}

//**********************************************************************

// inequality

template <class T>
bool operator!=(const array<T>& arr1, const array<T>& arr2) {
  return ! (arr1 == arr2);
}

//**********************************************************************

// equality with tolerance

template <class T>
bool is_equal(const array<T>& arr1, const array<T>& arr2) {
  if ( arr1.size != arr2.size ) return false;
  for ( int i=0; i<arr1.size; i++ ) 
    if ( ! is_equal( arr1(i), arr2(i) ) ) return false;
  return true;
}

//**********************************************************************

#endif