CPPLapack Code

//=============================================================================
/*! clear vector */
inline void dcovector::clear()
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::clear()" << std::endl;
  std::cerr << "# [NOTE] dcovector::clear()" << " An array at " << Array << " is going to be cleared." << std::endl;
#endif//CPPL_VERBOSE

  L =0;
  delete [] Array;  Array =NULL;
}

//=============================================================================
/*! make vector into zero vector */
inline dcovector& dcovector::zero()
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::zero()" << std::endl;
#endif//CPPL_VERBOSE
  
  for(long i=0; i<L; i++){ Array[i] =0.0; }
  return *this;
}

//=============================================================================
/*! change sign(+/-) of the vector */
inline void dcovector::chsign()
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::chsign()" << std::endl;
#endif//CPPL_VERBOSE
  
  for(long i=0; i<L; i++){ Array[i] =-Array[i]; }
}

//=============================================================================
/*! make a deep copy of the dcovector */
inline void dcovector::copy(const dcovector& vec)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::copy(const dcovector&)" << std::endl;
  std::cerr << "# [NOTE] dcovector::copy(const dcovector&)" << " A dcovector at " << Array << " is going to be deleted. ";
#endif//CPPL_VERBOSE
  
  L =vec.L;
  CAP =vec.CAP;
  delete [] Array;
  Array =new double[CAP];
  dcopy_(vec.L, vec.Array, 1, Array, 1);
  
#ifdef  CPPL_VERBOSE
  std::cerr << "Then, a COPY of a dcovector has been cleated at " << Array << "." << std::endl;
#endif//CPPL_VERBOSE
}

//=============================================================================
/*! make a shallow copy of the vector\n
 This function is not desinged to be used in project codes. */
inline void dcovector::shallow_copy(const _dcovector& vec)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::shallow_copy(const _dcovector&)" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] dcovector::shallow_copy(const _dcovector&) " << "A dcovector at " << Array << " is going to be deleted " << "and point at " << vec.Array << " instead." << std::endl;
#endif//CPPL_VERBOSE
  
  L =vec.L;
  CAP =vec.CAP;
  delete [] Array;
  Array =vec.Array;
  
  vec.nullify();
}

//=============================================================================
/*! make an alias of the vector\n
 Be carefull to use this function not to cause double free. */
inline void dcovector::alias(const dcovector& vec)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::alias(const dcovector&)" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] dcovector::alias(const dcovector&) " << "A dcovector at " << Array << " is going to be deleted " << "and point at " << vec.Array << " instead." << std::endl;
#endif//CPPL_VERBOSE
  
  L =vec.L;
  CAP =vec.CAP;
  delete [] Array;
  Array =vec.Array;
}

//=============================================================================
/*! unalias the vector */
inline void dcovector::unalias()
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::unalias()" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_VERBOSE
  std::cerr << "# [NOTE] dcovector::unalias() " << "A dcovector at " << Array << " is going to be out of reference." << std::endl;
#endif//CPPL_VERBOSE
  
  L =0;
  CAP =0;
  Array =NULL;
}

//=============================================================================
/*! resize vector */
inline dcovector& dcovector::resize(const long& _l, const long margin)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::resize(const long&, const long&=0)" << std::endl;
#endif//CPPL_VERBOSE
  
#ifdef  CPPL_DEBUG
  if( _l<0 ){
    std::cerr << "[ERROR] dcovector::resize(const long&, const long&=0)" << std::endl
              << "Vector size must be positive integers." << std::endl
              << "Your input was (" << _l << ", " << margin << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  L =_l;
  CAP =L+margin;
  delete [] Array;
  Array =new double[CAP];
  
  return *this;
}

//=============================================================================
/*! stretch or shrink vector */
inline void dcovector::stretch(const long& dl)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dcovector::stretch(const long&)" << std::endl;
#endif//CPPL_VERBOSE
  
  if(dl==0){ return; }
  
#ifdef  CPPL_DEBUG
  if( l+dl<0 ){
    std::cerr << "[ERROR] dcovector::stretch(const long&)" << std::endl
              << "Vector size must be positive integers." << std::endl
              << "Your input was (" << dl << ")." << std::endl;
    exit(1);
  }
#endif//CPPL_DEBUG
  
  L +=dl;
  
  if(L>CAP){
    CAP += dl+CPPL_SECTOR_SIZE;
    double* newArray(new double[CAP]);
    dcopy_(L-dl, Array, 1, newArray, 1);
    delete [] Array;
    Array =newArray;
  }
}

//=============================================================================
/*! swap two vectors */
inline void swap(dcovector& u, dcovector& v)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] swap(dcovector&, dcovector&)" << std::endl;
#endif//CPPL_VERBOSE
  
  long u_cap(u.CAP), u_l(u.L);
  double* u_array(u.Array);
  u.L=v.L; u.CAP=v.CAP; u.Array=v.Array;
  v.L=u_l; v.CAP=u_cap; v.Array=u_array;
}

//=============================================================================
/*! convert user object to smart-temporary object */
inline _dcovector _(dcovector& vec)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] _(dcovector&)" << std::endl;
#endif//CPPL_VERBOSE

  _dcovector newvec;
  
  //////// shallow copy ////////
  newvec.L =vec.L;
  newvec.CAP =vec.CAP;
  newvec.Array =vec.Array;
  
  //////// nullify ////////
  vec.L =0;
  vec.CAP =0;
  vec.Array =NULL;
  
  return newvec;
}