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

Timestamp:

May 23, 2005, 6:18:17 PM (20 years ago)

Author:

ansari

Message:

1/ Suppression du constructeur de copie avec specification d'organisation memoire (MemoryMapping) pour TMatrix<T> et TVector<T>
2/ Amelioration de l'impression des TArray/TMatrix avec la specification setw(largeur)
3/ Correction petit bug dans la lecture fichiers ASCII argument non transmis ds utilarr.cc

Reza 23 Mai 2005

File:

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

trunk/SophyaLib/TArray/tmatrix.cc

@@ -1 @@
-// $Id: tmatrix.cc,v 1.29 2004-09-10 09:55:07 cmv Exp $
+// $Id: tmatrix.cc,v 1.30 2005-05-23 16:18:17 ansari Exp $
 //                         C.Magneville          04/99
 #include "sopnamsp.h"
 #include "machdefs.h"
+#include <iostream>
+#include <iomanip>
 #include <stdio.h>
 #include <stdlib.h>
@@ -102 +104 @@
 
 //! Constructor of a matrix from a TArray \b a
-template <class T>
-TMatrix<T>::TMatrix(const TArray<T>& a)
-: TArray<T>(a)
-{
-  if (a.NbDimensions() > 2) 
-    throw SzMismatchError("TMatrix<T>::TMatrix(const TArray<T>& a) a.NbDimensions()>2 ");
-  if (a.NbDimensions() == 1) {
-    size_[1] = 1;
-    step_[1] = size_[0]*step_[0];
-    ndim_ = 2;
-  }
-  arrtype_ = 1;   // Type = Matrix
-  UpdateMemoryMapping(a, SameMemoryMapping);
-}
-
-//! Constructor of a matrix from a TArray \b a
 /*!
   \param a : TArray to be copied or shared
   \param share : if true, share data. If false copy data
-  \param mm : define the memory mapping type
- */
-template <class T>
-TMatrix<T>::TMatrix(const TArray<T>& a, bool share, short mm )
+ */
+template <class T>
+TMatrix<T>::TMatrix(const TArray<T>& a, bool share)
 : TArray<T>(a, share)
 {
@@ -135 +120 @@
   }
   arrtype_ = 1;   // Type = Matrix
-  UpdateMemoryMapping(a, mm);
+  UpdateMemoryMapping(a, SameMemoryMapping);
 }
 
@@ -264 +249 @@
 TMatrix<T> TMatrix<T>::SubMatrix(Range rline, Range rcol) const
 {
+  Range rx, ry;
   short mm = GetMemoryMapping();
-  Range rx, ry;
   if (mm == CMemoryMapping)  { rx = rcol;  ry = rline; }
   else { ry = rcol;  rx = rline; }
-  TMatrix sm(SubArray(rx, ry, Range(0), Range(0), Range(0)),true, mm);
+  TMatrix sm(SubArray(rx, ry, Range(0), Range(0), Range(0)),true);
   sm.UpdateMemoryMapping(mm);
   return(sm);
@@ -389 +374 @@
   if (maxprt < 0)  maxprt = max_nprt_;
   sa_size_t npr = 0;
+
+  // keep stream's io flags 
+  ios_base::fmtflags ioflg = os.flags();
+  os << right ;
+
   Show(os, si);
   if (ndim_ < 1)  return;
+  // Calcul de la largeur d'impression pour chaque element
+  int fprtw = os.precision()+7;
+  int prtw = 5;
+
+  if ( (typeid(T) == typeid( int_4 )) || (typeid(T) == typeid( uint_4 )) ) prtw = 8;
+  else if ( (typeid(T) == typeid( int_8 )) || (typeid(T) == typeid( uint_8 )) ) prtw = 11;
+  else if ( typeid(T) == typeid( r_4 ) ) prtw = fprtw;
+  else if ( typeid(T) == typeid( r_8 ) ) prtw = fprtw;
+  else if ( typeid(T) == typeid(complex<r_4>) ) prtw = fprtw;
+  else if ( typeid(T) == typeid(complex<r_8>) ) prtw = fprtw;
+
   sa_size_t kc,kr;  
   for(kr=0; kr<size_[marowi_]; kr++) {
@@ -396 +397 @@
     for(kc=0; kc<size_[macoli_]; kc++) {
       if(kc > 0) os << " ";  
-      os << (*this)(kr, kc);   npr++; 
+      os << setw(prtw) << (*this)(kr, kc);   npr++; 
       if (npr >= (sa_size_t) maxprt) {
         if (npr < totsize_)  os << "\n     .... " << endl; return;
@@ -404 +405 @@
   }
   os << endl;
+  os.flags(ioflg);  // reset stream io flags
 }