vector by transposed vector multiplication

include/Geometry/matrix.h

@@ -161,11 +161,15 @@ class Matrix
 
 // Vector / Matrix multiplication
 template <unsigned int M, unsigned int N, typename T>
-Vector<M,T> operator*(const Vector<N,T>& v, const Matrix<M,N,T>& m);
+Vector<M,T> operator*(const Vector<N,T>& v, const Matrix<M,N,T>& m) ;
 
 // Scalar / Matrix multiplication
 template <unsigned int M, unsigned int N, typename T>
-Matrix<M,N,T> operator*(T s, const Matrix<M,N,T>& m);
+Matrix<M,N,T> operator*(T s, const Matrix<M,N,T>& m) ;
+
+// Vector / Transposed vector multiplication
+template <unsigned int M, unsigned int N, typename T>
+Matrix<M,N,T> transposed_vectors_mult(const Vector<M,T>& v1, const Vector<N,T>& v2) ;
 
 
 /**********************************************/

include/Geometry/matrix.hpp

@@ -405,9 +405,11 @@ bool Matrix<M,N,T>::operator==(const Matrix<M,N,T>& m) const {
 template <unsigned int M, unsigned int N, typename T>
 std::ostream& operator<<(std::ostream& out, const Matrix<M,N,T>& m)
 {
-	for(unsigned int i = 0; i < M; ++i)
+	for(unsigned int i = 0; i < M; ++i) {
 		for(unsigned int j = 0; j < N; ++j)
 			out << m(i,j) << " " ;
+		std::cout << std::endl ;
+	}
 	return out ;
 }
 
@@ -446,6 +448,15 @@ Matrix<M,N,T> operator*(T s, const Matrix<M,N,T>& m)
 	return res ;
 }
 
+template <unsigned int M, unsigned int N, typename T>
+Matrix<M,N,T> transposed_vectors_mult(const Vector<M,T>& v1, const Vector<N,T>& v2) {
+	Matrix<M,N,T> res ;
+	for(unsigned int i = 0; i < M; ++i)
+		for(unsigned int j = 0; j < N; ++j)
+			res(i,j) = v1[i] * v2[j] ;
+	return res ;
+}
+
 } // namespace Geom
 
 } // namespace CGoGN