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

Timestamp:

Jul 29, 2004, 2:31:16 PM (21 years ago)

Author:

ansari

Message:

1/ Remplacement des methodes Add/Sub/Mul/DivElt(a) par

Add/Sub/Mul/DivElt(TArray a, TArray res)

2/ Operateurs += -= A+B A-B TArray et TMatrix/TVecteur modifies en consequence
3/ Ajout methode TArray::ScalarProduct()
4/ Methode TArray::SetT renomme en SetCst() SetT garde en alias
5/ Ajout parametre bool fzero (mise a zero) ajoute ds constructeur et

ReSize() de TMatrix et TVecteur.

Reza 29/07/2004

File:

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

trunk/SophyaLib/TArray/tmatrix.h

@@ -15 +15 @@
   // Creation / destruction 
   TMatrix();
-  TMatrix(sa_size_t r,sa_size_t c, short mm=BaseArray::AutoMemoryMapping);
+  TMatrix(sa_size_t r,sa_size_t c, short mm=BaseArray::AutoMemoryMapping, bool fzero=true);
   TMatrix(const TMatrix<T>& a);
   TMatrix(const TMatrix<T>& a, bool share);
@@ -48 +48 @@
   inline sa_size_t NCol()  const {return size_[macoli_]; } // back-compat Peida
 
-  void ReSize(sa_size_t r,sa_size_t c, short mm=BaseArray::SameMemoryMapping);
+  void ReSize(sa_size_t r,sa_size_t c, short mm=BaseArray::SameMemoryMapping, bool fzero=true);
   //! a synonym (alias) for method ReSize(sa_size_t, sa_size_t, short)
-  inline void SetSize(sa_size_t r,sa_size_t c, short mm=BaseArray::SameMemoryMapping)
-                 { ReSize(r, c, mm); }
+  inline void SetSize(sa_size_t r,sa_size_t c, short mm=BaseArray::SameMemoryMapping, bool fzero=true)
+                 { ReSize(r, c, mm, fzero); }
   // Reallocation de place
   void Realloc(sa_size_t r,sa_size_t c, short mm=BaseArray::SameMemoryMapping, bool force=false);
@@ -103 +103 @@
   //  operations avec matrices 
   //! += : add a matrix
-  inline  TMatrix<T>&  operator += (const TMatrix<T>& a)  { AddElt(a); return(*this); }
+  inline  TMatrix<T>&  operator += (const TMatrix<T>& a)  { AddElt(a,*this); return(*this); }
   //! -= : substract a matrix
-  inline  TMatrix<T>&  operator -= (const TMatrix<T>& a)  { SubElt(a); return(*this); }
+  inline  TMatrix<T>&  operator -= (const TMatrix<T>& a)  { SubElt(a,*this); return(*this); }
+
   TMatrix<T>  Multiply(const TMatrix<T>& b, short mm=BaseArray::SameMemoryMapping) const;
-  //! *= : matrix product : C = (*this)*B
-  inline  TMatrix<T>&  operator *= (const TMatrix<T>& b)
-          { this->Set(Multiply(b));  return(*this); }
+  //A supprimer ? Reza Juillet 2004 ! *= : matrix product : C = (*this)*B
+  //  inline  TMatrix<T>&  operator *= (const TMatrix<T>& b)
+  //       { this->Set(Multiply(b));  return(*this); }
 
   // I/O print, ...
@@ -201 +202 @@
   the original matrix elements.  */
 template <class T> inline TMatrix<T> operator - (const TMatrix<T>& a)
-    {TMatrix<T> result; result.CloneOrShare(a); result.SetTemp(true); 
-    result.NegateElt(); return result;}
+    {TMatrix<T> result; result.SetTemp(true); 
+    a.NegateElt(result); return result;}
 
 
@@ -215 +216 @@
 inline TMatrix<T> operator + (const TMatrix<T>& a,const TMatrix<T>& b)
     {TMatrix<T> result; result.SetTemp(true); 
-    if (b.IsTemp())  { result.Share(b); result.AddElt(a); }
-    else { result.CloneOrShare(a); result.AddElt(b); }
-    return result; }
+    a.AddElt(b, result);     return result; }
 
 
@@ -225 +224 @@
 inline TMatrix<T> operator - (const TMatrix<T>& a,const TMatrix<T>& b)
     {TMatrix<T> result; result.SetTemp(true); 
-    if (b.IsTemp())  { result.Share(b); result.SubElt(a, true); }
-    else { result.CloneOrShare(a); result.SubElt(b); }
-    return result; }
+    a.SubElt(b, result);     return result; }
+    
 
 // Surcharge d'operateurs C = A * B