zhsmatrix-calc.hpp File Reference

Go to the source code of this file.

Functions
_zhsmatrix	t (const zhsmatrix &mat)
_zhsmatrix	conj (const zhsmatrix &mat)
_zhsmatrix	conjt (const zhsmatrix &mat)
void	idamax (CPPL_INT &i, CPPL_INT &j, const zhsmatrix &mat)
comple	damax (const zhsmatrix &mat)

Function Documentation

_zhsmatrix t ( const zhsmatrix &  mat )

inline

return transposed zhsmatrix

Definition at line 3 of file zhsmatrix-calc.hpp.

References _(), conj(), and zhsmatrix::data.

{CPPL_VERBOSE_REPORT;
  zhsmatrix newmat(mat);
  
  const std::vector<zcomponent>::iterator newmat_data_end =newmat.data.end();
  for(std::vector<zcomponent>::iterator it=newmat.data.begin(); it!=newmat_data_end; it++){
    it->v =std::conj(it->v);
  }
  
  return _(newmat);
}

_zhsmatrix conj ( const zhsmatrix &  mat )

inline

return its conjugate matrix

Definition at line 21 of file zhsmatrix-calc.hpp.

References _(), and zhsmatrix::data.

Referenced by t().

{CPPL_VERBOSE_REPORT;
  zhsmatrix newmat(mat);
  
  const std::vector<zcomponent>::iterator newmat_data_end =newmat.data.end();
  for(std::vector<zcomponent>::iterator it=newmat.data.begin(); it!=newmat_data_end; it++){
    it->v =std::conj(it->v);
  }
  
  return _(newmat);
}

_zhsmatrix conjt ( const zhsmatrix &  mat )

inline

return its conjugate matrix

Definition at line 35 of file zhsmatrix-calc.hpp.

References _().

{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  WARNING_REPORT;
  std::cerr << "This function call has no effect since the matrix is Hermitian." << std::endl;
#endif//CPPL_DEBUG
  
  zhsmatrix newmat(mat);
  return _(newmat);
}

void idamax ( CPPL_INT &  i, CPPL_INT &  j, const zhsmatrix &  mat  )

inline

search the index of element having the largest absolute value in 0-based numbering system

Definition at line 53 of file zhsmatrix-calc.hpp.

References zhsmatrix::data.

{CPPL_VERBOSE_REPORT;
  std::vector<zcomponent>::const_iterator itx(mat.data.begin());
  double vmax =0.;
  
  const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    if( vmax < norm(it->v) ){
      vmax =norm(it->v);
      itx =it;
    }
  }
  
  i =itx->i;
  j =itx->j;
}

comple damax ( const zhsmatrix &  mat )

inline

return its largest absolute value

Definition at line 72 of file zhsmatrix-calc.hpp.

References zhsmatrix::data.

{CPPL_VERBOSE_REPORT;
  std::vector<zcomponent>::const_iterator itx(mat.data.begin());
  double vmax =0.;
  
  const std::vector<zcomponent>::const_iterator mat_data_end =mat.data.end();
  for(std::vector<zcomponent>::const_iterator it=mat.data.begin(); it!=mat_data_end; it++){
    if( vmax < norm(it->v) ){
      vmax =norm(it->v);
      itx =it;
    }
  }
  
  return itx->v;
}