Changeset 4035 in Sophya for trunk/SophyaLib/TArray/tmatrix.cc

Timestamp:

Nov 14, 2011, 5:28:25 PM (14 years ago)

Author:

ansari

Message:

1/ modif mineure ds TArray<T>::::ReadASCII() au print level global de BaseArray
2/ Correction bug gestion memoire au niveau des constructeurs de copie TArray/TMatrix/TVector
avec un BaseArray en argument. Ajout argument optionnel bool pack a ces constructeurs
3/ On autorise desormais la creation des objets TArray/TMatrix/TVector par constructeur de
copie sur des objets non alloues

Reza, 14/11/2011

File:

• trunk/SophyaLib/TArray/tmatrix.cc (modified) (3 diffs)

trunk/SophyaLib/TArray/tmatrix.cc

@@ -1 @@
-// $Id: tmatrix.cc,v 1.40 2010-08-05 12:16:19 ansari Exp $
+// $Id: tmatrix.cc,v 1.41 2011-11-14 16:28:25 ansari Exp $
 //                         C.Magneville          04/99
 #include "sopnamsp.h"
@@ -127 +127 @@
 : TArray<T>(a, share)
 {
+  arrtype_ = 1;   // Type = Matrix
+  if (a.NbDimensions() == 0)  {   // Reza-Nov 2011: we allow copy contrsuctor on non allocated arrays  
+    UpdateMemoryMapping(a, SameMemoryMapping);
+    return;
+  }
   if (a.NbDimensions() > 2) 
     throw SzMismatchError("TMatrix<T>::TMatrix(const TArray<T>& a, ...) a.NbDimensions()>2");
@@ -134 +139 @@
     ndim_ = 2;
   }
+  UpdateMemoryMapping(a, SameMemoryMapping);
+}
+
+//! Constructor of a matrix from a TArray \b a , with different data type
+/*!
+  Matrix size and memory layout are copied from the array \b a, or a packed matrix is created if \b pack==true.
+  \param a : original array, to copy sizes and data from
+  \param pack : if \b true , create a packed matrix, else same memory layout as \b a. 
+*/
+template <class T>
+TMatrix<T>::TMatrix(const BaseArray& a, bool pack)
+  : TArray<T>()
+{
+  // On ne peut pas passer par TArray<T>(const BaseArray&), car il faut initialiser  arrtype_ d'abord !
   arrtype_ = 1;   // Type = Matrix
+  if (a.NbDimensions() == 0)  {   // Reza-Nov 2011: we allow copy contrsuctor on non allocated arrays  
+    UpdateMemoryMapping(a, SameMemoryMapping);
+    return;
+  }
+  if (a.NbDimensions() > 2) 
+    throw SzMismatchError("TMatrix<T>::TMatrix((const BaseArray& a, bool) a.NbDimensions()>2");
+  string exmsg = "TMatrix<T>::TMatrix(const BaseArray&,bool)";
+  TArray<T>::ReSize(a,pack,false);
   UpdateMemoryMapping(a, SameMemoryMapping);
-}
-
-//! Constructor of a matrix from a TArray \b a , with a different data type
-template <class T>
-TMatrix<T>::TMatrix(const BaseArray& a)
-: TArray<T>()
-{
-  arrtype_ = 1;   // Type = Matrix
-  UpdateMemoryMapping(a, SameMemoryMapping);
-  SetBA(a);
+  TArray<T>::ConvertAndCopyElt(a);
 }