diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-06-24 15:13:00 +0000 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-06-24 15:13:00 +0000 |
commit | 3b94436d2fc7d8f6716b056ed544acecf256c4fe (patch) | |
tree | 0736b130a5ea62a21be73a58a6eba1f90dd025b7 /doc/snippets/MatrixBase_isUnitary.cpp | |
parent | c9560df4a0c274eb5011f0596682a3cf3274363e (diff) |
* vectorize dot product, copying code from sum.
* make the conj functor vectorizable: it is just identity in real case,
and complex doesn't use the vectorized path anyway.
* fix bug in Block: a 3x1 block in a 4x4 matrix (all fixed-size)
should not be vectorizable, since in fixed-size we are assuming
the size to be a multiple of packet size. (Or would you prefer
Vector3d to be flagged "packetaccess" even though no packet access
is possible on vectors of that type?)
* rename:
isOrtho for vectors ---> isOrthogonal
isOrtho for matrices ---> isUnitary
* add normalize()
* reimplement normalized with quotient1 functor
Diffstat (limited to 'doc/snippets/MatrixBase_isUnitary.cpp')
-rw-r--r-- | doc/snippets/MatrixBase_isUnitary.cpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/doc/snippets/MatrixBase_isUnitary.cpp b/doc/snippets/MatrixBase_isUnitary.cpp new file mode 100644 index 000000000..d29aa5f71 --- /dev/null +++ b/doc/snippets/MatrixBase_isUnitary.cpp @@ -0,0 +1,5 @@ +Matrix3d m = Matrix3d::identity(); +m(0,2) = 1e-4; +cout << "Here's the matrix m:" << endl << m << endl; +cout << "m.isUnitary() returns: " << m.isUnitary() << endl; +cout << "m.isUnitary(1e-3) returns: " << m.isUnitary(1e-3) << endl; |