Fixes wrong dimensions in x = A * y.

Author: prj-

src/femlib/RNM.hpp

@@ -1019,9 +1019,9 @@ class KN :public KN_<R> { public:
 
 
    KN& operator =(const Mul_KNM_KN_<R> & u) 
-        { if(this->unset()) this->set(new R[u.b.N()],u.b.N());KN_<R>::operator=(u);return *this;}
+        { if(this->unset()) this->set(new R[u.A.N()],u.A.N());KN_<R>::operator=(u);return *this;}
     KN& operator =(const Mul_KNMh_KN_<R> & u)
-    { if(this->unset()) this->set(new R[u.b.N()],u.b.N());KN_<R>::operator=(u);return *this;}
+    { if(this->unset()) this->set(new R[u.A.M()],u.A.M());KN_<R>::operator=(u);return *this;}
 
 //   KN& operator =(const MatriceCreuseMulKN_<R> & Ax) 
 //       {if(this->unset()) set(new R[Ax.v.N()],Ax.v.N()); KN_<R>::operator=(Ax);return *this;}

src/femlib/RNM_op.hpp

@@ -94,7 +94,6 @@ template<class R>
 KN_<R>& KN_<R>::operator oper (const Mul_KNMh_KN_<R> & u)  {
     K_throwassert (SameShape(u.A.shapei) && !constant());
     R * l(v); KN_<const_R>  li(u.A(0,'.')); //  first line
-    std::cout << " Mul_KNMh_KN_ "<< endl;
     for (long i=0;i<n;i++,l += step,++li)
         *l oper (conj(li),u.b);
     return *this;}