dgematrix-misc.hpp

Go to the documentation of this file.

//=============================================================================
/*! clear all the matrix data and set the sizes 0 */
inline void dgematrix::clear()
{
#ifdef  CPPSL_DEBUG
  std::cerr << "# [NOTE] dgematrix::clear() "
            << "An array at " << array << " is going to be cleared."
            << std::endl;
#endif//CPPSL_DEBUG
  
  M =N =0;
  Ml =Nl =0;
  for(int i=0; i<9; i++){ Desc[i]=0; }
  delete [] Array;
  Array =new double[0];
}

//=============================================================================
/*! change the matrix into a zero matrix */
inline void dgematrix::zero()
{
  for(long i=0; i<Ml*Nl; i++){ Array[i] =0.0; }
}

//=============================================================================
/*! change the matrix into an identity matrix */
inline void dgematrix::identity()
{
#ifdef  CPPSL_DEBUG
  if(M!=N){
    std::cerr << "[ERROR] dgematrix::identity()" << std::endl
              << "Only square matrix can be a identity matrix." << std::endl
              << "The matrix size was " << M << "x" << N << "." << std::endl;
    exit(1);
  }
#endif//CPPSL_DEBUG
  
  for(long i=0; i<Ml*Nl; i++){ Array[i] =0.0; }
  for(long i=0; i<M; i++){ operator()(i,i) =1.0; }
}

//=============================================================================
/*! change sign(+/-) of the matrix */
inline void dgematrix::chsign()
{
  for(long i=0; i<Ml*Nl; i++){ Array[i] =-Array[i]; }
}

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

//=============================================================================
/*! make a shallow copy of the matrix\n
  This function is not designed to be used in project codes. */
inline void dgematrix::shallow_copy(const _dgematrix& mat)
{
#ifdef  CPPSL_DEBUG
  std::cerr << "# [NOTE] dgematrix:shallow_copy(const _dgematrix&) "
            << "A dgematrix at " << Array << " is going to be deleted, "
            << "and point to " << mat.Array << " instead." << std::endl;
#endif//CPPSL_DEBUG

  delete [] Array;
  M =mat.M; N =mat.N;
  Ml =mat.Ml; Nl =mat.Nl;
  for(int i=0; i<9; i++){ Desc[i] =mat.Desc[i]; }
  Array =mat.Array;
}

//=============================================================================
/*! resize the matrix */
inline void dgematrix::resize(const long& _m, const long& _n)
{
#ifdef  CPPSL_DEBUG
  if( _m<0 || _n<0 ){
    std::cerr << "[ERROR] dgematrix::resize(const long&, const long&)"
              << std::endl
              << "Matrix sizes must be positive integers." << std::endl
              << "Your input was (" << _m << "," << _n << ")." << std::endl;
    exit(1);
  }
#endif//CPPSL_DEBUG
  
  delete [] Array;
  M =_m; N =_n;
  Ml =numroc_( M, mb, myrow, 0, nprow );
  Nl =numroc_( N, nb, mycol, 0, npcol );
  desc[0]=1; desc[1]=icontxt; desc[2]=M; desc[3]=N; desc[4]=mb; desc[5]=nb;
  desc[6]=0; desc[7]=0; desc[8]=max(1,Ml);
  Array =new double[Ml*Nl];
}

//=============================================================================
/*! swap two matrices */
inline void swap(dgematrix& A, dgematrix& B)
{
  long A_M(A.M), A_N(A.N), A_Ml(A.Ml), A_Nl(A.Nl), A_Desc[9];
  double* A_Array(A.Array);
  A.M=B.M; A.N=B.N; A.Ml=B.Ml; A.Nl=B.Nl;  A.Array=B.Array;
  B.M=A_M; B.N=A_N; B.Ml=A_Ml; B.Nl=A_Nl;  B.Array=A_Array;
  for(int i=0; i<9; i++){
    A_Desc[i] =A.Desc[i];
    A.Desc[i] =B.Desc[i];
    B.Desc[i] =A_Desc[i];
  }
}

---