jiri/LA_library
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

smat - implemented rank2 update

Browse Source
This commit is contained in:
jiri
2026-10-23 16:35:22 +02:00
parent 338381e160
commit c4055403f6
2 changed files with 40 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files
smat.cc
+39
View File
@@ -677,6 +677,45 @@ const double NRSMat< std::complex<double> >::norm(const std::complex<double> sca
}
/***************************************************************************//**
* symmetric (hermitian) rank-2 update of this matrix stored in packed form
* \f[S \leftarrow \alpha x y^H + \alpha^* y x^H + S\f]
* generic version, BLAS-based specializations for double and complex<double> follow
******************************************************************************/
template <typename T>
void NRSMat<T>::rank2update(const T alpha, const NRVec<T> &x, const NRVec<T> &y) {
#ifdef DEBUG
if(nn != x.size() || nn != y.size()) laerror("incompatible dimensions in NRSMat<T>::rank2update");
#endif
NOT_GPU(*this); NOT_GPU(x); NOT_GPU(y);
copyonwrite();
const T alphac = LA_traits<T>::conjugate(alpha);
T *p = v;
for(int i=0; i<nn; ++i) for(int j=0; j<=i; ++j)
*p++ += alpha*x[i]*LA_traits<T>::conjugate(y[j]) + alphac*y[i]*LA_traits<T>::conjugate(x[j]);
}
template<>
void NRSMat<double>::rank2update(const double alpha, const NRVec<double> &x, const NRVec<double> &y) {
#ifdef DEBUG
if(nn != x.size() || nn != y.size()) laerror("incompatible dimensions in NRSMat<double>::rank2update");
#endif
NOT_GPU(*this); NOT_GPU(x); NOT_GPU(y);
copyonwrite();
cblas_dspr2(CblasRowMajor, CblasLower, nn, alpha, x, 1, y, 1, v);
}
template<>
void NRSMat<std::complex<double> >::rank2update(const std::complex<double> alpha, const NRVec<std::complex<double> > &x, const NRVec<std::complex<double> > &y) {
#ifdef DEBUG
if(nn != x.size() || nn != y.size()) laerror("incompatible dimensions in NRSMat<std::complex<double> >::rank2update");
#endif
NOT_GPU(*this); NOT_GPU(x); NOT_GPU(y);
copyonwrite();
cblas_zhpr2(CblasRowMajor, CblasLower, nn, &alpha, x, 1, y, 1, v);
}
/***************************************************************************//**
* for this real double-precision symmetric matrix \f$S\f$ stored in packed form,
* real scalar value \f$\alpha\f$ and real double-precision symmetric matrix \f$T\f$, compute