Changeset 2575 in Sophya for trunk/SophyaLib/TArray/tvector.cc

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/tvector.cc (modified) (5 diffs)

trunk/SophyaLib/TArray/tvector.cc

@@ -1 @@
-// $Id: tvector.cc,v 1.15 2002-11-15 09:35:44 ansari Exp $
+// $Id: tvector.cc,v 1.16 2004-07-29 12:31:16 ansari Exp $
 //                         C.Magneville          04/99
 #include "machdefs.h"
@@ -52 +52 @@
   \param lcv : line or column vector ?
   \param mm : memory mapping type
+  \param fzero : if \b true , set vector elements to zero 
   \sa SelectVectorType
  */
 template <class T>
-TVector<T>::TVector(sa_size_t n, short lcv, short mm)
+TVector<T>::TVector(sa_size_t n, short lcv, short mm, bool fzero)
 // Constructeur 
-  : TMatrix<T>(1,1,mm)
+  : TMatrix<T>(1,1,mm,false)
 {
   arrtype_ = 2;   // Type = Vector
   lcv = SelectVectorType(lcv);
-  ReSize(n,lcv);
+  ReSize(n,lcv,fzero);
 }
 
@@ -145 +146 @@
  */
 template <class T>
-void TVector<T>::ReSize(sa_size_t n, short lcv)
+void TVector<T>::ReSize(sa_size_t n, short lcv, bool fzero)
 {
   if( n == 0 ) 
@@ -154 +155 @@
   if (lcv == ColumnVector) { r = n;  c = 1; }
   else { c = n; r = 1; }
-  TMatrix<T>::ReSize(r,c);
+  TMatrix<T>::ReSize(r,c,BaseArray::SameMemoryMapping,fzero);
   veceli_ = (lcv ==  ColumnVector ) ?  marowi_ : macoli_;
 }
@@ -191 +192 @@
   TVector sv( mtx(rr, cr) , true, GetVectorType(), GetMemoryMapping() );
   return(sv);  
-}
-
-//! Return the norm \f$ \sum{V(i)^2} \f$
-template <class T>
-T TVector<T>::Norm2() const
-{
-  T ret = 0;
-  for(sa_size_t k=0; k<Size(); k++)   ret += (*this)(k)*(*this)(k);
-  return ret;
 }