Changeset 813 in Sophya for trunk/SophyaLib/TArray/tmatrix.h

Timestamp:

Apr 5, 2000, 5:44:19 PM (25 years ago)

Author:

ansari

Message:

Correction bug/amelioarions TArray,TMatrix,TVector - Reza 5/4/2000

File:

• trunk/SophyaLib/TArray/tmatrix.h (modified) (4 diffs)

trunk/SophyaLib/TArray/tmatrix.h

@@ -34 +34 @@
   void Realloc(uint_4 r,uint_4 c, short mm=SameMemoryMapping, bool force=false);
 
-  // Sub-matrix extraction $CHECK$ Reza 03/2000  Doit-on declarer cette methode const ?
-  TMatrix<T> operator () (Range rline, Range rcol) const ;
+  // Sub-matrix extraction $CHECK$ Reza 03/2000  Doit-on declarer ces methode const ?
+  TMatrix<T> SubMatrix(Range rline, Range rcol) const ;
+  inline TMatrix<T> operator () (Range rline, Range rcol) const 
+                    { return SubMatrix(rline, rcol); }
+  // Lignes et colonnes de la matrice
+  inline TMatrix<T> Row(uint_4 ir) const
+                    { return SubMatrix(Range(ir,ir), Range(0,NCols()-1)); }
+  inline TMatrix<T> Column(uint_4 ic) const
+                    { return SubMatrix(Range(0,NRows()-1), Range(ic,ic)); }
 
   // Inline element acces methods 
@@ -53 +60 @@
   inline  TMatrix<T>& operator = (IdentityMatrix imx) { return SetIdentity(imx); }
 
+  inline  TMatrix<T>&  operator = (Sequence seq)    { SetSeq(seq); return(*this); }
+
   // Operations diverses  avec une constante
-  inline  TMatrix<T>&  operator = (T x)             { Set(x); return(*this); }
+  inline  TMatrix<T>&  operator = (T x)             { SetT(x); return(*this); }
   inline  TMatrix<T>&  operator += (T x)            { Add(x); return(*this); }
   inline  TMatrix<T>&  operator -= (T x)            { Sub(x); return(*this); }
@@ -61 +70 @@
 
   //  operations avec matrices 
-  inline  TMatrix<T>&  operator += (const TMatrix<T>& a)  { return AddElt(a); }
-  inline  TMatrix<T>&  operator -= (const TMatrix<T>& a)  { return SubElt(a); }
+  inline  TMatrix<T>&  operator += (const TMatrix<T>& a)  { AddElt(a); return(*this); }
+  inline  TMatrix<T>&  operator -= (const TMatrix<T>& a)  { SubElt(a); return(*this); }
+  // Produit matriciel Multiply : C = (*this)*B
+  TMatrix<T>  Multiply(const TMatrix<T>& b, short mm=SameMemoryMapping) const;
+  inline  TMatrix<T>&  operator *= (const TMatrix<T>& b)
+          { this->Set(Multiply(b));  return(*this); }
 
-  // Input-Output
+  // I/O print, ...
+  virtual string InfoString() const;
   virtual void  Print(ostream& os, int_4 maxprt=-1, bool si=false) const ;
-
-  // Produit matriciel 
-  TMatrix<T>  Multiply(const TMatrix<T>& b, short mm=SameMemoryMapping) const;
-  //  inline TMatrix<T>& operator *= (const TMatrix<T>& b) 
 
 protected:
@@ -98 +108 @@
 
 
+// Surcharge d'operateurs C = A (+,-) B
+// $CHECK$ Reza 3/4/2000 Pas necessaire  de redefinir les operateurs
+// Defini au niveau de TArray<T> - Pour ameliorer l'efficacite
+// Doit-on le faire aussi pour les constantes ? - Fin de $CHECK$ Reza 3/4/2000
+
+template <class T>
+inline TMatrix<T> operator + (const TMatrix<T>& a,const TMatrix<T>& b)
+    {TMatrix<T> result(a); result.SetTemp(true); result.AddElt(b); return result;}
+
+template <class T>
+inline TMatrix<T> operator - (const TMatrix<T>& a,const TMatrix<T>& b)
+    {TMatrix<T> result(a); result.SetTemp(true); result.SubElt(b); return result;}
+
 // Surcharge d'operateurs C = A * B
 
 template <class T> inline TMatrix<T> operator * (const TMatrix<T>& a, const TMatrix<T>& b)
-{ TMatrix<T> result(a); result.SetTemp(true); result.Multiply(b);  return result;}
+{ TMatrix<T> result(a); result.SetTemp(true); return(result.Multiply(b)); }
 
 typedef TMatrix<r_8> Matrix;