AnimaL

Tutorial

Documentation

---

Documentation

• Overview
• Modules
• Namespace List
• Class Hierarchy
• Alphabetical List
• Compound List
• File List
• Namespace Members
• Compound Members
• File Members
• Examples

---

vectorVec.h

Go to the documentation of this file.

#ifndef animal_vectorvec_hpp___________________________
#define animal_vectorvec_hpp___________________________

#include <animal/vector.h>

namespace animal {

  template<class Vec>
  struct VectorVec : public animal::vector<Vec> 
  {

VectorVec(){}

    VectorVec( std::size_t s )
      : animal::vector<Vec>(s)
    {}

    VectorVec( std::size_t s, Vec v )
      : animal::vector<Vec>(s,v)
    {}

  };

template <class V> inline
double v_dot( const VectorVec<V>& a, const VectorVec<V>& b )
{
  //cerr << "v_dot, a = " << animal::v_output(a) << endl;
  //cerr << "v_dot, b = " << animal::v_output(b) << endl;
  double    res = 0;
  typename animal::container_traits<VectorVec<V> const>::iterator ai,aend,bi;
  for( ai=animal::begin(a),aend=animal::end(a),bi=animal::begin(b); ai!=aend; ++ai, ++bi )
    {
      //res += animal::v_dot(*ai,*bi);
      res += animal::v_dot((*ai),(*bi));
      //cerr << "v_dot, ai = " << *ai <<", bi = "<< *bi << endl;
    }
  return res;
}

template <class V> inline
void v_eq ( VectorVec<V>& Y, const VectorVec<V>& X )
{
    assert( animal::size(Y) == animal::size(X) );
    typename animal::container_traits<VectorVec<V> >::iterator Yi = animal::begin(Y), Yend=animal::end(Y);
    typename animal::container_traits<VectorVec<V> const>::iterator Xi = animal::begin(X);
    for( ; Yi!=Yend; ++Yi, ++Xi )
        *(Yi) =  (*Xi);
}

template <class V> inline
void v_eq ( VectorVec<V>& Y, const VectorVec<V>& W, const double& a, const VectorVec<V>& X )
{
    assert( animal::size(Y) == animal::size(X) );
    assert( animal::size(Y) == animal::size(W) );
    typename animal::container_traits<VectorVec<V> >::iterator Yi = animal::begin(Y), Yend=animal::end(Y);
    typename animal::container_traits<VectorVec<V> const>::iterator Xi = animal::begin(X), Wi=animal::begin(W);
    for( ; Yi!=Yend; ++Yi, ++Xi, ++Wi )
        *(Yi) = (*Wi) + (*Xi)*a;
}

template <class V> inline
void v_peq ( VectorVec<V>& Y, double a, const VectorVec<V>& X )
{
    assert( animal::size(Y) == animal::size(X) );
    typename animal::container_traits<VectorVec<V> >::iterator Yi = animal::begin(Y), Yend=animal::end(Y);
    typename animal::container_traits<VectorVec<V> const>::iterator Xi = animal::begin(X);
    for( ; Yi!=Yend; ++Yi, ++Xi )
        *(Yi) += (*Xi)*a;
}

template <class V> inline
void v_meq ( VectorVec<V>& Y, const VectorVec<V>& X )
{
    assert( animal::size(Y) == animal::size(X) );
    typename animal::container_traits<VectorVec<V> >::iterator Yi = animal::begin(Y), Yend=animal::end(Y);
    typename animal::container_traits<VectorVec<V> const>::iterator Xi = animal::begin(X);
    for( ; Yi!=Yend; ++Yi, ++Xi )
        *(Yi) -= (*Xi);
}

template <class V> inline
void v_teq ( VectorVec<V>& v, const double& a )
{
    typename animal::container_traits<VectorVec<V> >::iterator Yi = animal::begin(v), Yend=animal::end(v);
    for( ; Yi!=Yend; ++Yi )
        *(Yi) *= a;
}

template <class V> inline
double v_norm( const VectorVec<V>& a ){ return sqrt(v_dot(a,a)); }


}//animal

#endif

---