You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
23 lines
677 B
23 lines
677 B
2 years ago
|
#include <Eigen/Dense>
|
||
|
|
#include <iostream>
|
||
|
|
|
||
|
|
using namespace Eigen;
|
||
|
|
|
||
|
|
template <typename Derived1, typename Derived2>
|
||
|
|
void copyUpperTriangularPart(MatrixBase<Derived1>& dst, const MatrixBase<Derived2>& src)
|
||
|
|
{
|
||
|
|
/* Note the 'template' keywords in the following line! */
|
||
|
|
dst.template triangularView<Upper>() = src.template triangularView<Upper>();
|
||
|
|
}
|
||
|
|
|
||
|
|
int main()
|
||
|
|
{
|
||
|
|
MatrixXi m1 = MatrixXi::Ones(5,5);
|
||
|
|
MatrixXi m2 = MatrixXi::Random(4,4);
|
||
|
|
std::cout << "m2 before copy:" << std::endl;
|
||
|
|
std::cout << m2 << std::endl << std::endl;
|
||
|
|
copyUpperTriangularPart(m2, m1.topLeftCorner(4,4));
|
||
|
|
std::cout << "m2 after copy:" << std::endl;
|
||
|
|
std::cout << m2 << std::endl << std::endl;
|
||
|
|
}
|