Changeset 2575 in Sophya for trunk/SophyaLib/TArray/tarray.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/tarray.h (modified) (5 diffs)

trunk/SophyaLib/TArray/tarray.h

@@ -125 +125 @@
   //! Fill TArray with Sequence \b seq
   inline  TArray<T>&  operator = (Sequence const & seq)    { return SetSeq(seq); }
+
 // A = x (tous les elements a x)
-  virtual TArray<T>&  SetT(T x);
+  virtual TArray<T>&  SetCst(T x);
+  //! Fill an array with a constant value \b x ( alias for \b SetCst() method )
+  inline  TArray<T>&  SetT(T x) { return SetCst(x); }
   //! Fill TArray with all elements equal to \b x
   inline  TArray<T>&  operator = (T x)             { return SetT(x); }
@@ -138 +141 @@
 
 // A += -= *= /= x (ajoute, soustrait, ... x a tous les elements)
+  // Methodes Add/Sub/Mul/Div() sont la pour compatibilite avant V=2 (1.818)
+  // Faut-il les garder ? Reza, Juillet 2004
   inline  TArray<T>&  Add(T x)                     { return AddCst(x, *this); }
   inline  TArray<T>&  Sub(T x, bool fginv=false)   { return SubCst(x, *this, fginv); }
@@ -158 +163 @@
 
 // A += -=  (ajoute, soustrait element par element les deux tableaux )
-  virtual TArray<T>&  AddElt(const TArray<T>& a);
-  //! Operator TArray += TArray
-  inline  TArray<T>&  operator += (const TArray<T>& a)  { return AddElt(a); }
-  virtual TArray<T>&  SubElt(const TArray<T>& a, bool fginv=false);
-  //! Operator TArray -= TArray
-  inline  TArray<T>&  operator -= (const TArray<T>& a)  { return SubElt(a); }
+  virtual TArray<T>&  AddElt(const TArray<T>& a, TArray<T>& res) const ;
+  virtual TArray<T>&  SubElt(const TArray<T>& a, TArray<T>& res, bool fginv=false) const ;
 // Multiplication, division element par element les deux tableaux 
-  virtual TArray<T>&  MulElt(const TArray<T>& a);
-  virtual TArray<T>&  DivElt(const TArray<T>& a, bool fginv=false, bool divzero=false);
+  virtual TArray<T>&  MulElt(const TArray<T>& a, TArray<T>& res) const ;
+  virtual TArray<T>&  DivElt(const TArray<T>& a, TArray<T>& res, bool fginv=false, bool divzero=false) const ;
+
+  //! Operator TArray += TArray (element by element addition in place)
+  inline  TArray<T>&  operator += (const TArray<T>& a)  { return AddElt(a, *this); }
+  //! Operator TArray -= TArray (element by element subtraction in place)
+  inline  TArray<T>&  operator -= (const TArray<T>& a)  { return SubElt(a, *this); }
+
 // Recopie des valeurs, element par element
   virtual TArray<T>&  CopyElt(const TArray<T>& a);
 // Recopie des valeurs avec conversion prealable, element par element
   virtual TArray<T>&  ConvertAndCopyElt(const BaseArray& a);
+
+// Calcul du produit scalaire ( Somme_i (*this)(i)*a(i) )
+  virtual T ScalarProduct(const TArray<T>& a) const ;
+// Norme(^2) 
+  //! Returns the squarred of the array norm, defined as  Sum_k (*this)(k)*(*this)(k)
+  inline T Norm2() const { return ScalarProduct(*this); }
 
 // Somme et produit des elements
@@ -283 +296 @@
 inline TArray<T> operator + (const TArray<T>& a,const TArray<T>& b)
     { TArray<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; }
 
 /*! \ingroup TArray \fn operator-(const TArray<T>&,const TArray<T>&)
@@ -292 +303 @@
 inline TArray<T> operator - (const TArray<T>& a,const TArray<T>& b)
     { TArray<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; }