Changeset 1314 in Sophya for trunk/SophyaLib/TArray/basarr.cc

Timestamp:

Nov 9, 2000, 6:45:46 PM (25 years ago)

Author:

ansari

Message:

Ajout methode IndexAtPosition ds BaseArray - Reza 9/11/2000

File:

• trunk/SophyaLib/TArray/basarr.cc (modified) (2 diffs)

trunk/SophyaLib/TArray/basarr.cc

@@ -417 +417 @@
   //#ifdef SO_BOUNDCHECKING
   if (rest != 0) 
-    throw PError("BaseArray::Offset(sa_size_t ip) GUG !!! rest != 0");
+    throw PError("BaseArray::Offset(sa_size_t ip) BUG !!! rest != 0");
   //#endif
 //   if (rest != 0) cerr << " BUG ---- BaseArray::Offset( " << ip << " )" << rest << endl;
@@ -426 +426 @@
   for(k=0; k<ndim_; k++)  off += idx[k]*step_[k];
   return (off);
+}
+//! return index of element \b ip, along the five array axes 
+void BaseArray::IndexAtPosition(sa_size_t ip, sa_size_t & ix, sa_size_t & iy, 
+                                sa_size_t & iz, sa_size_t & it, sa_size_t & iu) const
+{
+  ix = iy = iz = it = iu = 0;
+  if ( (ndim_ < 1) || (ip == 0) )  return;
+  if (ip >= totsize_)
+    throw RangeCheckError("BaseArray::IndexAtPosition(...) Out of range index ip");
+  sa_size_t idx[BASEARRAY_MAXNDIMS];
+  int_4 k;
+  sa_size_t rest = ip;
+  for(k=0; k<ndim_; k++) {
+    idx[k] = rest%size_[k];   rest /= size_[k];
+    if (rest == 0)  break;
+  }
+  if (rest != 0) 
+    throw PError("BaseArray::IndexAtPosition(...) BUG !!! rest != 0");
+  ix = idx[0];
+  iy = idx[1];
+  iz = idx[2];
+  it = idx[3];
+  iu = idx[4];  
+  return;
 }