Changeset 4054 in Sophya for trunk/Eval/JET/smtx.cc

Timestamp:

Mar 17, 2012, 11:13:26 PM (14 years ago)

Author:

ansari

Message:

ajout programme telacc.cc pour test rapidite acces aux elements, Reza 17/03/2012

File:

• trunk/Eval/JET/smtx.cc (modified) (7 diffs)

trunk/Eval/JET/smtx.cc

@@ -17 +17 @@
 template <class T>
 SimpleMatrix<T>::SimpleMatrix()
-  : MtxSize(), JETExpression<T, JETExpMtx<T> >(), dbv_(&data_, this) 
+  : MtxSize(), JETExpression<T, JETExpMtx<T> >(), dbv_(&data_, this),
+    data_p_(NULL)
 {
   this->Expression().JETExpMtxSet(&data_, this);
@@ -24 +25 @@
 template <class T>
 SimpleMatrix<T>::SimpleMatrix(sa_size_t nr, sa_size_t nc)
-  : MtxSize(nr, nc), JETExpression<T, JETExpMtx<T> >(), data_(nr*nc), dbv_(&data_, this)
+  : MtxSize(nr, nc), JETExpression<T, JETExpMtx<T> >(), data_(nr*nc), dbv_(&data_, this),
+    data_p_(data_.Begin())
 {
   this->Expression().JETExpMtxSet(&data_, this);
@@ -31 +33 @@
 template <class T>
 SimpleMatrix<T>::SimpleMatrix(const SimpleMatrix & m)
-  : MtxSize(m), JETExpression<T, JETExpMtx<T> >(), data_(m.data_), dbv_(&data_, this)
+  : MtxSize(m), JETExpression<T, JETExpMtx<T> >(), data_(m.data_), dbv_(&data_, this),
+    data_p_(data_.Begin())    
 {
   this->Expression().JETExpMtxSet(&data_, this);
@@ -38 +41 @@
 template <class T>
 SimpleMatrix<T>::SimpleMatrix(const SimpleMatrix & m, bool share)
-  : MtxSize(m), JETExpression<T, JETExpMtx<T> >(dbv_), data_(m.data_, share), dbv_(&data_, this)
+  : MtxSize(m), JETExpression<T, JETExpMtx<T> >(dbv_), data_(m.data_, share), dbv_(&data_, this),
+    data_p_(data_.Begin())
 {
 }
@@ -55 +59 @@
   MtxSize::ReSize(a);
   data_.Realloc(nrow_*ncol_);
+  data_p_=data_.Begin();
   return;
 }
@@ -70 +75 @@
     throw(SzMismatchError("SimpleMatrix::Set() Size(a) != Size(b)"));
   data_ = b.data_;
+  data_p_=data_.Begin();  
   return(*this);
 }
@@ -88 +94 @@
     throw RangeCheckError("SimpleMatrix<T>::Init(T c. T s) Size=0");
   for(sa_size_t k=0; k<nrow_*ncol_; k++) {
-    data_(k) = c; c += step;
+    data_p_[k] = c; c += step;
   }
   return(*this);