00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012 #ifndef MTL_RANK_TWO_UPDATE_INCLUDE
00013 #define MTL_RANK_TWO_UPDATE_INCLUDE
00014
00015 #include <boost/numeric/mtl/operation/rank_one_update.hpp>
00016
00017 namespace mtl { namespace matrix {
00018
00020
00022 template <typename Matrix, typename VectorX, typename VectorY>
00023 inline void rank_two_update(Matrix& matrix, const VectorX& x, const VectorY& y)
00024 {
00025 rank_one_update(matrix, x, y);
00026 rank_one_update(matrix, y, x);
00027 }
00028
00029 }}
00030
00031 #endif // MTL_RANK_TWO_UPDATE_INCLUDE