#include "main.h" namespace Eigen {template <typename Lhs,typename Rhs>const Product<Lhs,Rhs>const Lhs& lhs, const Rhs& rhs)return Product<Lhs,Rhs>(lhs,rhs);template <typename Lhs,typename Rhs>const Product<Lhs,Rhs,LazyProduct>const Lhs& lhs, const Rhs& rhs)return Product<Lhs,Rhs,LazyProduct>(lhs,rhs);template <typename DstXprType, typename SrcXprType>const EigenBase<DstXprType> &dst, const SrcXprType &src)typename DstXprType::Scalar,typename SrcXprType::Scalar>());return dst.const_cast_derived();template <typename DstXprType, template <typename > class StorageBase, typename SrcXprType>const DstXprType& copy_using_evaluator(const NoAlias<DstXprType, StorageBase>& dst, const SrcXprType &src)typename DstXprType::Scalar,typename return dst.expression();template <typename DstXprType, typename SrcXprType>const PlainObjectBase<DstXprType> &dst, const SrcXprType &src)#ifdef EIGEN_NO_AUTOMATIC_RESIZING0 || (IsVectorAtCompileTime ? (dst.size() == src.size())"Size mismatch. Automatic resizing is disabled because EIGEN_NO_AUTOMATIC_RESIZING is defined" );#else #endif typename DstXprType::Scalar,typename SrcXprType::Scalar>());return dst.const_cast_derived();template <typename DstXprType, typename SrcXprType>void add_assign_using_evaluator(const DstXprType& dst, const SrcXprType& src)typedef typename DstXprType::Scalar Scalar;const_cast <DstXprType&>(dst), src.derived(), internal::add_assign_op<Scalar,typename SrcXprType::Scalar>());template <typename DstXprType, typename SrcXprType>void subtract_assign_using_evaluator(const DstXprType& dst, const SrcXprType& src)typedef typename DstXprType::Scalar Scalar;const_cast <DstXprType&>(dst), src.derived(), internal::sub_assign_op<Scalar,typename SrcXprType::Scalar>());template <typename DstXprType, typename SrcXprType>void multiply_assign_using_evaluator(const DstXprType& dst, const SrcXprType& src)typedef typename DstXprType::Scalar Scalar;alar,typename SrcXprType::Scalar>());template <typename DstXprType, typename SrcXprType>void divide_assign_using_evaluator(const DstXprType& dst, const SrcXprType& src)typedef typename DstXprType::Scalar Scalar;alar,typename SrcXprType::Scalar>());template <typename DstXprType, typename SrcXprType>void swap_using_evaluator(const DstXprType& dst, const SrcXprType& src)typedef typename DstXprType::Scalar Scalar;_assign_op<Scalar>());namespace internal {template <typename Dst, template <typename > class StorageBase, typename Src, typename Func>void call_assignment(const NoAlias<Dst,StorageBase>& dst, const Src& src, const Func& func)template <typename Dst, template <typename > class StorageBase, typename Src, typename Func>void call_restricted_packet_assignment(const NoAlias<Dst,StorageBase>& dst, const Src& src, const Func& func)template <typename XprType> long get_cost(const XprType& ) { return Eigen::internal::evaluator<XprType>::CoeffReadCost; }using namespace std;#define VERIFY_IS_APPROX_EVALUATOR(DEST,EXPR) VERIFY_IS_APPROX(copy_using_evaluator(DEST,(EXPR)), (EXPR).eval());#define VERIFY_IS_APPROX_EVALUATOR2(DEST,EXPR,REF) VERIFY_IS_APPROX(copy_using_evaluator(DEST,(EXPR)), (REF).eval());// Testing Matrix evaluator and Transpose const Vector2d v_const(v);// Testing Transpose // Transpose as rvalue // Transpose as lvalue // Testing Array evaluator 2 ,3 );3 ,2 );1 ,2 ,3 , 4 ,5 ,6 ;const ArrayXXf a_const(a);// Testing CwiseNullaryOp evaluator 1 ).all() && (w.array() <= 1 ).all()); // not easy to test ... 3 ));// mix CwiseNullaryOp and transpose // test product expressions int s = internal::random<int >(1 ,100 );// check that prod works even with aliasing present // check compound assignment of products // test product with all possible sizes int s = internal::random<int >(1 ,100 );float , 1 , 1 > m11, res11; m11.setRandom(1 ,1 );float , 1 , 4 > m14, res14; m14.setRandom(1 ,4 );float , 1 ,Dynamic> m1X, res1X; m1X.setRandom(1 ,s);float , 4 , 1 > m41, res41; m41.setRandom(4 ,1 );float , 4 , 4 > m44, res44; m44.setRandom(4 ,4 );float , 4 ,Dynamic> m4X, res4X; m4X.setRandom(4 ,s);float ,Dynamic, 1 > mX1, resX1; mX1.setRandom(s,1 );float ,Dynamic, 4 > mX4, resX4; mX4.setRandom(s,4 );float ,Dynamic,Dynamic> mXX, resXX; mXX.setRandom(s,s);2 ,3 );3 ,2 );1 ,2 ,3 , 4 ,5 ,6 ;const ArrayXXf a_const(a);// this does not work because Random is eval-before-nested: // copy_using_evaluator(w, Vector2d::Random().transpose()); // test CwiseUnaryOp 3 * v);3 * v).transpose());3 ).transpose());2 * a_const + 3 ).transpose());// test CwiseBinaryOp ctor2d::Constant(3 )));// dynamic matrices and arrays 6 ,6 ), mat2(6 ,6 );6 ,6 ));6 ,6 ), arr2(6 ,6 );6 ,6 , 3 .0 ));// test automatic resizing 3 ,3 );9 ,9 );// test direct traversal // matrix, nullary // TODO: find a way to test direct traversal with array // transpose 2 * Matrix3f::Identity()); // unary // binary 0 ,0 ,2 ,2 ), Matrix3f::Identity().block(1 ,1 ,2 ,2 )); // block // test linear traversal // matrix, nullary // array // transpose 2 * Matrix3f::Zero()); // unary // binary // test inner vectorization // matrix // array // transpose // TODO: find out why Matrix4f::Zero() does not allow inner vectorization 2 * m4src); // unary // binary // test linear vectorization 6 ,6 ), mXsrc = MatrixXf::Random(6 ,6 );6 ,6 ), aXsrc = ArrayXXf::Random(6 ,6 );// matrix // array // transpose 6 ,6 )); // nullary 2 * mXsrc); // unary // binary // test blocks and slice vectorization 4 ,4 >(1 ,0 )));10 , 10 , 3 .0 ).block(2 , 3 , 6 , 6 ));1 , 1 , 2 , 3 ), m3.bottomRows(2 ));1 , 1 , 2 , 3 ) = m3.bottomRows(2 );20 ,20 );20 ,20 );9 ,12 );4 , 4 , 9 , 12 ), mXsrc);4 , 4 , 9 , 12 ) = mXsrc;// test Map const float raw[3 ] = {1 ,2 ,3 };float buffer[3 ] = {0 ,0 ,0 };const Vector3f>(raw));const Array3f>(raw));2 *v3;0 ] == 2 );1 ] == 4 );2 ] == 6 );// test CwiseUnaryView 6 ,6 ), matXcd_ref(6 ,6 );// test Select 0 ).select(aXsrc, -aXsrc));// test Replicate 6 , 6 );6 );6 , 6 );12 , 12 );2 ,2 ));2 ,2 >()));// test partial reductions 6 );// test MatrixWrapper and ArrayWrapper 6 ,6 );6 ,6 );2 ).matrix());2 );6 ,6 );6 ,6 ).array());// test Reverse // test Diagonal 5 );1 ));1 >());1 ), vec1);1 ) = vec1;1 >(), mat1.diagonal(1 ));1 >() = mat2.diagonal(1 );// test swapping 6 , 6 );2 * mat1 + MatrixXd::Identity(6 , 6 );0 , 0 , 3 , 3 ), mat2.block(3 , 3 , 3 , 3 ));0 , 0 , 3 , 3 ).swap(mat2ref.block(3 , 3 , 3 , 3 ));2 ), mat2.col(3 ).transpose());2 ).swap(mat2.col(3 ).transpose());// test compound assignment const Matrix4d mat_const = Matrix4d::Random(); 1 ), 2 *mat.row(2 ));1 ) -= 2 *mat_ref.row(2 );const ArrayXXf arr_const = ArrayXXf::Random(5 ,3 ); 5 , 3 , 0 .5 );1 ), arr.row(2 ) + 1 );1 ) /= (arr_ref.row(2 ) + 1 );// test triangular shapes 6 ,6 ), B(6 ,6 ), C(6 ,6 ), D(6 ,6 );pper>()));iew<UnitLower>()));iew<UnitUpper>()));"copy_using_evaluator(B.triangularView<Upper>(), A)" );"copy_using_evaluator(B.triangularView<Lower>(), A.triangularView<Lower>())" ););"copy_using_evaluator(B.triangularView<Lower>(), A.triangularView<Lower>().transpose())" );"swap_using_evaluator(B.triangularView<Upper>(), A.triangularView<Upper>())" );larView<Upper>()*A));jointView<Upper>()*A));// test diagonal shapes 6 );6 ,6 ), B(6 ,6 );*A));()));// test CoeffReadCost 1 );3 );2 *a+b), 4 );1 );15 );1 );15 );1 );15 );// regression test for PR 544 and bug 1622 (introduced in #71609c4) // test restricted_packet_assignment with an unaligned destination const size_t M = 2 ;const size_t K = 2 ;const size_t N = 5 ;float *destMem = new float [(M*N) + 1 ];float *dest = (internal::UIntPtr(destMem)%EIGEN_MAX_ALIGN_BYTES) == 0 ? destMem+1 : destMem;const Matrix<float , Dynamic, Dynamic, RowMajor> a = Matrix<float , Dynamic, Dynamic, RowMajor>::Random(M, K);const Matrix<float , Dynamic, Dynamic, RowMajor> b = Matrix<float , Dynamic, Dynamic, RowMajor>::Random(K, N);float , Dynamic, Dynamic, RowMajor> > z(dest, M, N);;float , Dynamic, Dynamic, RowMajor>, Matrix<float , Dynamic, Dynamic, RowMajor>, LazyProduct> tmp(a,b);assign_op<float , float >());delete [] destMem;
Messung V0.5 in Prozent C=94 H=92 G=92
¤ Dauer der Verarbeitung: 0.10 Sekunden
(vorverarbeitet am 2026-06-06)
¤ *© Formatika GbR, Deutschland
