GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-08 17:59:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

KroneckerProduct: we have const_cast_derived so why not use it?

Browse Source
This commit is contained in:
Chen-Pang He 2012-10-14 01:38:38 +08:00
parent 20a984cd2e
commit f34db6578a
1 changed files with 4 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
unsupported/Eigen/src/KroneckerProduct/KroneckerTensorProduct.h
View File

@ -31,6 +31,8 @@ void kroneckerProduct_full(const Derived_A& A, const Derived_B& B, Derived_AB &
Ac = A.cols(), Ac = A.cols(),
Br = B.rows(), Br = B.rows(),
Bc = B.cols(); Bc = B.cols();
AB.resize(Ar*Br,Ac*Bc);
for (unsigned int i=0; i<Ar; ++i) for (unsigned int i=0; i<Ar; ++i)
for (unsigned int j=0; j<Ac; ++j) for (unsigned int j=0; j<Ac; ++j)
AB.block(i*Br,j*Bc,Br,Bc) = A(i,j)*B; AB.block(i*Br,j*Bc,Br,Bc) = A(i,j)*B;
@ -78,22 +80,6 @@ void kroneckerProduct_sparse(const Derived_A &A, const Derived_B &B, Derived_AB
} // end namespace internal } // end namespace internal
/*!
* Computes Kronecker tensor product of two dense matrices
*
* \param a Dense matrix a
* \param b Dense matrix b
* \param c Kronecker tensor product of a and b
*/
template<typename A,typename B,typename CScalar,int CRows,int CCols, int COptions, int CMaxRows, int CMaxCols>
void kroneckerProduct(const MatrixBase<A>& a, const MatrixBase<B>& b, Matrix<CScalar,CRows,CCols,COptions,CMaxRows,CMaxCols>& c)
{
c.resize(a.rows()*b.rows(),a.cols()*b.cols());
internal::kroneckerProduct_full(a.derived(), b.derived(), c);
}
/*! /*!
* Computes Kronecker tensor product of two dense matrices * Computes Kronecker tensor product of two dense matrices
* *
@ -107,10 +93,9 @@ void kroneckerProduct(const MatrixBase<A>& a, const MatrixBase<B>& b, Matrix<CSc
* \param c Kronecker tensor product of a and b * \param c Kronecker tensor product of a and b
*/ */
template<typename A,typename B,typename C> template<typename A,typename B,typename C>
void kroneckerProduct(const MatrixBase<A>& a, const MatrixBase<B>& b, MatrixBase<C> const & c_) void kroneckerProduct(const MatrixBase<A>& a, const MatrixBase<B>& b, const MatrixBase<C>& c)
{ {
MatrixBase<C>& c = const_cast<MatrixBase<C>& >(c_); internal::kroneckerProduct_full(a.derived(), b.derived(), c.const_cast_derived());
internal::kroneckerProduct_full(a.derived(), b.derived(), c.derived());
} }
/*! /*!