vector.hpp

Go to the documentation of this file.

#ifndef VECTOR_HPP
#define VECTOR_HPP
 
#include <iostream>
#include "debug.hpp"
 
using namespace std;
 
using Real = double;//float;
 
#ifdef DEBUG
struct Vector{
  size_t n;
  Real* data;
  Real& operator[](size_t i);
};
 
inline Real&
Vector::operator[](size_t i){
  assert(i<n);
  return data[i];
}
 
inline Vector
init_vector(size_t n){
  Vector v;
  v.n=n;
  v.data=new Real[n];
  return v;
}
 
inline void
delete_vector(Vector v){
  delete[] v.data;
}
 
inline bool
is_null(Vector v){
  return v.n==0;
}
 
static const Vector NullVector={0,nullptr};
 
#else
using Vector=Real*;
 
inline Vector
init_vector(size_t n){
  return new Real[n];
}
 
inline void
delete_vector(Vector v){
  delete[] v;
}
#endif
 
inline void
display(Vector v,size_t n){
  if(n==0){
    cout<<"[]"<<endl;
    return;
  }
  cout<<'['<<v[0];
  for(size_t i=1;i<n;++i){
    cout<<','<<v[i];
  }
  cout<<']'<<endl;
}
 
inline size_t
argmax(Vector v,size_t n){
  assert(n>0);
  size_t imax=0;
  Real vmax=v[0];
  for(size_t i=1;i<n;++i){
    if(v[i]>vmax){
      vmax=v[i];
      imax=i;
    }
  }
  return imax;
}
 
inline size_t
indice2(size_t i,size_t j,size_t nj){
  return i*nj+j;
}
 
inline size_t
indice3(size_t i,size_t j,size_t k,size_t nj,size_t nk){
  return (i*nj+j)*nk+k;
}
 
inline size_t
indice4(size_t i,size_t j,size_t k,size_t l,size_t nj,size_t nk,size_t nl){
  return ((i*nj+j)*nk+k)*nl+l;
}
 
#endif