Changeset 2667 in Sophya for trunk/SophyaLib/HiStats/ntuple.cc

Timestamp:

Apr 15, 2005, 1:51:38 PM (20 years ago)

Author:

ansari

Message:

Correction bug ds NTuple::GetVecD() (else manquant) suite a passage a NTuple double/ou/float - Reza 15/4/2005

File:

• trunk/SophyaLib/HiStats/ntuple.cc (modified) (3 diffs)

trunk/SophyaLib/HiStats/ntuple.cc

@@ -15 +15 @@
    \class SOPHYA::NTuple
    \ingroup HiStats
-   Simple NTuple class (a Table or 2-D data set) with floating value
-   columns. 
+   Simple NTuple class (a Table or 2-D data set) with double or 
+   single precision floating point value columns. 
    \sa SOPHYA::ObjFileIO<NTuple>
 
@@ -23 +23 @@
    // ...
    char * names[3] = {"XPos", "YPos", "Val"};
-   // NTuple (Table) creation with 3 columns
+   // NTuple (Table) creation with 3 columns (double precision)
    NTuple  nt(3, names);
    // Filling the NTuple
-   r_4 x[3];
+   r_8 x[3];
    for(int i=0; i<63; i++) {
      x[0] = (i%9)-4.;  x[1] = (i/9)-3.;  x[2] = x[0]*x[0]+x[1]*x[1];
@@ -385 +385 @@
 int numb = n/mBlk;
 int offb = n-numb*mBlk;
-if (mFgDouble) 
-  memcpy(ret, (mPtrD[numb]+offb*mNVar), mNVar*sizeof(r_8));
-for(i=0; i<mNVar; i++) ret[i] = (mPtr[numb]+offb*mNVar)[i];
+if (mFgDouble) memcpy(ret, (mPtrD[numb]+offb*mNVar), mNVar*sizeof(r_8));
+else for(i=0; i<mNVar; i++) ret[i] = (mPtr[numb]+offb*mNVar)[i];
 return(ret);
 }