std_vec_ops.h

Go to the documentation of this file.

#ifndef _STD_VEC_OPS_H_
#define _STD_VEC_OPS_H_

#include <math.h>
#include <vector>
#include <basic.h>

namespace StdVecOps {

  inline
  long sum( const std::vector<int>& a_iv  )
  {
    long retval(0);
    for (int i=0; i<a_iv.size(); i++) retval += a_iv[i];
    return retval;
  }

  inline
  long product( const std::vector<int>& a_iv  )
  {
    long retval(1);
    for (int i=0; i<a_iv.size(); i++) retval *= a_iv[i];
    return retval;
  }

  inline
  void copyFrom(  std::vector<int>& a_iv,       
                  const int* const  a_iv_carr,  
                  int               a_n         
               )
  {
    a_iv.resize(a_n);
    for (int i=0; i<a_n; i++) a_iv[i] = a_iv_carr[i];
    return;
  }

  inline
  void add( std::vector<int>& a_iv, 
            const int         a_a    )
  {
    for (int i=0; i<a_iv.size(); i++) a_iv[i] += a_a;
    return;
  }

  inline
  void copyFrom(  std::vector<double>& a_iv,       
                  const double* const  a_iv_carr,  
                  int                  a_n         
               )
  {
    a_iv.resize(a_n);
    for (int i=0; i<a_n; i++) a_iv[i] = a_iv_carr[i];
    return;
  }

  inline
  void createNormalVector(std::vector<double>&  a_normal_vec, 
                          const int*            a_vec,        
                          const int             a_size        )
  {
    a_normal_vec = std::vector<double>(a_size, 0.0);
    double magn = 0.0;
    for (int i=0; i<a_size; i++) {
      magn += (double)(a_vec[i]*a_vec[i]);
    }
    magn = sqrt(magn);
    if (magn > _MACHINE_ZERO_) {
      for (int i=0; i<a_size; i++) {
        a_normal_vec[i] = a_vec[i] / magn;
      }
    }
    return;
  }

  inline
  void createNormalVector(std::vector<double>&    a_normal_vec, 
                          const std::vector<int>& a_vec         )
  {
    createNormalVector( a_normal_vec,
                        a_vec.data(),
                        a_vec.size() );
    return;
  }
                          
  inline
  double compute2Norm(  const std::vector<double>& a_a, 
                        const std::vector<double>& a_b   )
  {
    double retval = 0.0;

    for (int i = 0; i < std::min(a_a.size(),a_b.size()); i++) {
      retval += ( (a_a[i]-a_b[i]) * (a_a[i]-a_b[i]) );
    }
    retval = sqrt(retval);

    return retval;
  }
                          
}

#endif