Changeset 1081 in Sophya for trunk/SophyaLib/TArray/tarray.cc

Timestamp:

Jul 24, 2000, 2:51:27 PM (25 years ago)

Author:

ansari

Message:

Adaptation modifs MuTyV et services de copie entre tableaux de type different - Reza 24/7/2000

File:

• trunk/SophyaLib/TArray/tarray.cc (modified) (4 diffs)

trunk/SophyaLib/TArray/tarray.cc

@@ -131 +131 @@
 }
 
+//! Constructor with size and contents copied (after conversion) from a different type TArray
+template <class T>
+TArray<T>::TArray(const BaseArray& a)
+  : BaseArray() , mNDBlock() 
+{
+  string exmsg = "TArray<T>::TArray(const BaseArray&)";
+  if (!UpdateSizes(a, exmsg))  throw( ParmError(exmsg) );
+  mNDBlock.ReSize(totsize_);    
+  //  if (a.mInfo) mInfo = new DVList(*(a.mInfo));  - pb protected !
+  ConvertAndCopyElt(a);
+}
+
 //! Destructor 
 template <class T>
@@ -155 +167 @@
 }
 
+//! Set array elements equal to the \b a array elements, after conversion 
+template <class T>
+TArray<T>& TArray<T>::SetBA(const BaseArray& a)
+{
+  if (this == &a) return(*this);
+  if (a.NbDimensions() < 1) 
+    throw RangeCheckError("TArray<T>::SetBA(a ) - Array a not allocated ! ");
+  if (NbDimensions() < 1) {
+    string exmsg = "TArray<T>::SetBA(const BaseArray& a)";
+    if (!UpdateSizes(a, exmsg))  throw( ParmError(exmsg) );
+    mNDBlock.ReSize(totsize_);    
+  }
+  ConvertAndCopyElt(a);
+  return(*this);
+}
+
 //! Clone array \b a
 template <class T>
@@ -233 +261 @@
 }
 
-inline double _SqrtRz_(double x) { return sqrt(x); }  // Pb avec SGI-CC - $CHECK$ - Reza 04/2000
-
 //! Give value (in \b double) for element at position \b ip..
 template <class T>
-double TArray<T>::ValueAtPosition(uint_8 ip) const
+MuTyV & TArray<T>::ValueAtPosition(uint_8 ip) const
 {
 #ifdef SO_BOUNDCHECKING
   if (ip >= totsize_)  throw( ParmError("TArray<T>::ValueAtPosition(uint_8 ip) Out-of-bound Error") );
 #endif
-  return( (double)(*(mNDBlock.Begin()+Offset(ip))) );
-}
-
-//! Give value (in \b double) for element at position \b ip..
-/*!
-  For complex values, we return the module of the complex number
-*/
-double TArray< complex<r_4> >::ValueAtPosition(uint_8 ip) const
-{
-#ifdef SO_BOUNDCHECKING
-  if (ip >= totsize_)  throw( ParmError("TArray<T>::ValueAtPosition(uint_8 ip) Out-of-bound Error") );
-#endif
-  complex<r_4> c = *(mNDBlock.Begin()+Offset(ip));
-  double cr = (double)(c.real());
-  double ci = (double)(c.imag());
-  return( _SqrtRz_(cr*cr+ci*ci) );
-}
-
-//! Give value (in \b double) for element at position \b ip..
-/*!
-  For complex values, we return the module of the complex number
-*/
-double TArray< complex<r_8> >::ValueAtPosition(uint_8 ip) const
-{
-#ifdef SO_BOUNDCHECKING
-  if (ip >= totsize_)  throw( ParmError("TArray<T>::ValueAtPosition(uint_8 ip) Out-of-bound Error") );
-#endif
-  complex<r_8> c = *(mNDBlock.Begin()+Offset(ip));
-  double cr = (double)(c.real());
-  double ci = (double)(c.imag());
-  return( _SqrtRz_(cr*cr+ci*ci) );
+  my_mtv = *(mNDBlock.Begin()+Offset(ip));
+  return( my_mtv );
 }
 
@@ -754 +751 @@
       for(k=0, ka=0;  k<gpas;  k+=step, ka+=stepa)  pe[k] = pea[ka]; 
     }
+  }
+  return(*this);
+}
+
+//! Converts and Copy elements of \b a
+template <class T>
+TArray<T>& TArray<T>::ConvertAndCopyElt(const BaseArray& a)
+{
+  if (NbDimensions() < 1) 
+    throw RangeCheckError("TArray<T>::ConvertAndCopyElt(const TArray<T>& )  - Not Allocated Array ! ");
+  if (!CompareSizes(a)) 
+    throw(SzMismatchError("TArray<T>::ConvertAndCopyElt(const TArray<T>&) SizeMismatch")) ;
+
+  T * pe;
+  uint_8 j,k,ka;
+  uint_8 offa;
+  // Non regular data spacing ...
+  uint_4 ax = MaxSizeKA();
+  uint_8 step = Step(ax);
+  uint_8 stepa = a.Step(ax);
+  uint_8 gpas = Size(ax)*step;
+  uint_8 naxa = Size()/Size(ax);
+  for(j=0; j<naxa; j++)  {
+    pe = mNDBlock.Begin()+Offset(ax,j);
+    offa = a.Offset(ax,j);
+    for(k=0, ka=0;  k<gpas;  k+=step, ka+=stepa)  pe[k] = a.ValueAtPosition(offa+ka); 
   }
   return(*this);