Changeset 288 in Sophya for trunk/SophyaLib/NTools/matrix.cc

Timestamp:

May 3, 1999, 6:55:25 PM (26 years ago)

Author:

ansari

Message:

grosses modif avec refonte code dans tmatrix cmv 3/5/99

File:

• trunk/SophyaLib/NTools/matrix.cc (modified) (8 diffs)

trunk/SophyaLib/NTools/matrix.cc

@@ -1 @@
-// $Id: matrix.cc,v 1.2 1999-04-22 16:18:41 ansari Exp $
+// $Id: matrix.cc,v 1.3 1999-05-03 16:55:21 ansari Exp $
 
 #include "machdefs.h"
@@ -10 +10 @@
 #include "cvector.h"
 #include "generalfit.h"
+#include "tmatrix.h"
 
 // Sous Linux, LN_MAXDOUBLE LN_MINDOUBLE ne semblent pas etre definies
@@ -41 +42 @@
 //      Construit une matrice 1x1 (pour ppersist).
 //--
-: nr(1), nc(1), ndata(1), nalloc(1), data(new double[1])
+: nr(1), nc(1), ndata(1), nalloc(1), data(new double[1]), FromTmatrix(false)
 {
   // if (r<=0 || c<=0) THROW(rangeCheckErr);
@@ -53 +54 @@
 //      Construit une matrice de r lignes et c colonnes.
 //--
-: nr(r), nc(c), ndata(r*c), nalloc(r*c), data(new double[r*c])
+: nr(r), nc(c), ndata(r*c), nalloc(r*c), data(new double[r*c]), FromTmatrix(false)
 {
   // if (r<=0 || c<=0) THROW(rangeCheckErr);
@@ -66 +67 @@
 //      le tableau des valeurs : pas d'allocation.
 //--
-: nr(r), nc(c), ndata(r*c), nalloc(0), data(values)
+: nr(r), nc(c), ndata(r*c), nalloc(0), data(values), FromTmatrix(false)
 {
   END_CONSTRUCTOR
@@ -77 +78 @@
 //--
 : nr(a.nr), nc(a.nc), ndata(a.nr*a.nc), nalloc(a.nr*a.nc),
-  data(new double[a.nr*a.nc])
+  data(new double[a.nr*a.nc]), FromTmatrix(false)
 {
   memcpy(data, a.data, nalloc * sizeof(double));
@@ -83 +84 @@
 }
 
+//++
+Matrix::Matrix(const TMatrix<r_8>& a)
+//
+//      Constructeur par copie a partir d'une TMatrix<r_8>.
+//      Attention, les donnees sont partagees.
+//--
+: nr(a.NRows()), nc(a.NCols()), ndata(a.NRows()*a.NCols())
+, nalloc(a.NRows()*a.NCols())
+, data(const_cast<r_8 *>(a.Data()))
+, FromTmatrix(true)
+{
+  END_CONSTRUCTOR
+}
+
 
 Matrix::~Matrix()
 {
   DBASSERT(ndata == nr*nc);
-  if (nalloc) delete[] data;
+  if (nalloc && !FromTmatrix) delete[] data;
 }
 
@@ -102 +117 @@
 {
   DBASSERT(ndata == nr*nc);
-  for (int i=0; i<ndata; i++)
-    data[i] = 0;
+  for (int i=0; i<ndata; i++) data[i] = 0;
 }