dcovector-misc.hpp

Go to the documentation of this file.

//=============================================================================
/*! clear vector */
inline void dcovector::clear()
{
#ifdef  CPPSL_DEBUG
  std::cerr << "# [NOTE] dcovector::clear() "
            << "An array at " << Array << " is going to be cleared." << std::endl;
#endif//CPPSL_DEBUG
  
  L =0;
  Ll =0;
  Desc[0]=1; Desc[1]=icontxt; Desc[2]=0; Desc[3]=0; Desc[4]=mb; Desc[5]=nb;
  Desc[6]=0; Desc[7]=0; Desc[8]=max(1,Ll);
  delete [] Array;
  Array =new double[0];
}

//=============================================================================
/*! make vector into zero vector */
inline void dcovector::zero()
{
  for(long i=0; i<Ll; i++){ Array[i] =0.0; }
}

//=============================================================================
/*! change sign(+/-) of the vector */
inline void dcovector::chsign()
{
  for(long i=0; i<Ll; i++){ Array[i] =-Array[i]; }
}

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

//=============================================================================
/*! make a shallow copy of the vector */
inline void dcovector::shallow_copy(const _dcovector& vec)
{
#ifdef  CPPSL_DEBUG
  std::cerr << "# [NOTE] dcovector::shallow_copy(const dcovector&) "
            << "A dcovector at " << Array << " is going to be deleted "
            << "and point to " << vec.Array << " instead." << std::endl;
#endif//CPPSL_DEBUG
  
  delete [] Array;
  L =vec.L;
  Ll =vec.Ll;
  for(int i=0; i<9; i++){ Desc[i] =vec.Desc[i]; }
  Array =vec.Array;
}

//=============================================================================
/*! resize vector */
inline void dcovector::resize(const long& _l)
{
#ifdef  CPPSL_DEBUG
  if( _l<0 ){
    std::cerr << "[ERROR] dcovector::resize(const long&)" << std::endl
         << "Vector size must be positive integers." << std::endl
         << "Your input was (" << _l << ")." << std::endl;
    exit(1);
  }
#endif//CPPSL_DEBUG
  
  delete [] Array;
  L =_l;
  long ml(numroc_(L, mb, myrow, 0, nprow));
  if(mycol==0){ Ll=ml; } else{ Ll=0; }
  Desc[0]=1; Desc[1]=icontxt; Desc[2]=L; Desc[3]=1; Desc[4]=mb; Desc[5]=nb;
  Desc[6]=0; Desc[7]=0; Desc[8]=max(1,Ll);
  Array =new double[Ll];
}

//=============================================================================
/*! swap two vectors */
inline void swap(dcovector& u, dcovector& v)
{
  long u_L(u.L), u_Ll(u.Ll), u_Desc[9];
  double* u_Array(u.Array);
  u.L=v.L; u.Ll=v.Ll;  u.Array=v.Array;
  v.L=u_L; v.Ll=u_Ll;  v.Array=u_Array;
  for(int i=0; i<9; i++){
    u_Desc[i] =u.Desc[i];
    u.Desc[i] =v.Desc[i];
    v.Desc[i] =u_Desc[i];
  }
}

---