DVector& operator*(const DVectorMatrix& mat, const DVector& vec) { DVector *vect; vect = new DVector(mat.vesize); for (vei_t j=0; j<mat.vesize; j++) vect->set_comp(vec * mat.vector[j], F64B_TYPE, j); return *vect; }