vec.h

Go to the documentation of this file.

#ifndef VEC_H
#define VEC_H

#include <string>
#include <sstream>
using namespace std;

#include <typedefs.h>


template < typename A, typename T >
class vec
{
public:

  A ai;

  uintc N;

  vec(A ai_, uintc N_)
    : ai(ai_), N(N_) {}

  void operator = ( T const x );

  template< typename Z >
  void operator = ( Z const & zi );

  void operator = (vec<A,T> const & vb);

  template< typename Z >
  void operator *= ( Z const & zi );

  template< typename Z >
  void operator += ( Z const & zi );

  template< typename Z >
  void operator -= ( Z const & zi );

  void operator += ( T const x );

  void operator *= ( T const x );

  T const operator [] (uintc i) const
    { return ai[i]; }

  T & operator [] (uintc i)
    { return ai[i]; }

  
  void identity(uintc k, T const x=1);

  void difference();
 
  void sum();

  void binomial(uintc N);

  operator string () const;

};


//---------------------------------------------------------
//  Implementation

template < typename A, typename T >
vec<A,T>::operator string() const
{
  if (N==0)
    return "";

  string s;
  { stringstream ss;  ss << ai[0]; s+=ss.str(); }
  
  for (uint i=1; i<N; ++i)
  {
    s += " ";
    { stringstream ss;  ss << ai[i]; s+=ss.str(); }
  }
  
  return s;
}

template < typename A, typename T >
void vec<A,T>::binomial(uintc N)
{
  identity(0);
  for (uint i=0; i<N; ++i)
    sum();
}

template < typename A, typename T >
void vec<A,T>::identity(uintc k, T const x)
{
  for (uint i=0; i<N; ++i)
    ai[i] = 0;
  ai[k]=x;
}

template < typename A, typename T >
void vec<A,T>::difference()
{
  for (uint i=N-1; i>0; --i)
    ai[i] = ai[i] - ai[i-1];
}

template < typename A, typename T >
void vec<A,T>::sum()
{
  for (uint i=N-1; i>0; --i)
    ai[i] = ai[i] + ai[i-1];
}

template < typename A, typename T >
template< typename Z >
void vec<A,T>::operator = ( Z const & zi )
{
  for (uint i=0; i<N; ++i)
    ai[i] = zi[i];
}

template < typename A, typename T >
template< typename Z >
void vec<A,T>::operator *= ( Z const & zi )
{
  for (uint i=0; i<N; ++i)
    ai[i] *= zi[i];
}

template < typename A, typename T >
template< typename Z >
void vec<A,T>::operator += ( Z const & zi )
{
  for (uint i=0; i<N; ++i)
    ai[i] += zi[i];
}

template < typename A, typename T >
template< typename Z >
void vec<A,T>::operator -= ( Z const & zi )
{
  for (uint i=0; i<N; ++i)
    ai[i] -= zi[i];
}

template < typename A, typename T >
void vec<A,T>::operator = ( T x )
{
  for (uint i=0; i<N; ++i)
    ai[i] = x; 
}

template < typename A, typename T >
void vec<A,T>::operator *= ( T x )
{
  for (uint i=0; i<N; ++i)
    ai[i] *= x; 
}

template < typename A, typename T >
void vec<A,T>::operator += ( T x )
{
  for (uint i=0; i<N; ++i)
    ai[i] += x; 
}

template < typename A, typename T >
void vec<A,T>::operator = ( vec<A,T> const & bi )
{
  for (uint i=0; i<N; ++i)
    ai[i] = bi[i];
}


#endif