| 1 | #include "datatable.h" | 
|---|
| 2 | #include "sopnamsp.h" | 
|---|
| 3 | #include "strutil.h" | 
|---|
| 4 | #include "pexceptions.h" | 
|---|
| 5 | #include "fiosegdb.h" | 
|---|
| 6 |  | 
|---|
| 7 |  | 
|---|
| 8 | /*! | 
|---|
| 9 | \class SOPHYA::DataTable | 
|---|
| 10 | \ingroup HiStats | 
|---|
| 11 | This class can be used to organize data in table (row-column) form. | 
|---|
| 12 | Each column holds homogeneous data (same data type), while different | 
|---|
| 13 | columns can be used for different data types | 
|---|
| 14 | (integer, float, string ...). | 
|---|
| 15 | The whole data set is kept in memory. | 
|---|
| 16 | \sa SOPHYA::MuTyV | 
|---|
| 17 | \sa SOPHYA::DataTableRow | 
|---|
| 18 | \sa SOPHYA::BaseDataTable | 
|---|
| 19 | \sa SOPHYA::SegDataBlock | 
|---|
| 20 |  | 
|---|
| 21 | \code | 
|---|
| 22 | #include "datatable.h" | 
|---|
| 23 | // ... | 
|---|
| 24 | DataTable dt(64); | 
|---|
| 25 | dt.AddFloatColumn("X0_f"); | 
|---|
| 26 | dt.AddFloatColumn("X1_f"); | 
|---|
| 27 | dt.AddDoubleColumn("X0X0pX1X1_d"); | 
|---|
| 28 | double x[5]; | 
|---|
| 29 | for(int i=0; i<63; i++) { | 
|---|
| 30 | x[0] = (i%9)-4.;  x[1] = (i/9)-3.;  x[2] = x[0]*x[0]+x[1]*x[1]; | 
|---|
| 31 | dt.AddRow(x); | 
|---|
| 32 | } | 
|---|
| 33 | // Printing table info | 
|---|
| 34 | cout << dt ; | 
|---|
| 35 | // Saving object into a PPF file | 
|---|
| 36 | POutPersist po("dtable.ppf"); | 
|---|
| 37 | po << dt ; | 
|---|
| 38 | \endcode | 
|---|
| 39 |  | 
|---|
| 40 | DataTableRow objects can also be used to fill the table, or access the data, as | 
|---|
| 41 | shown in the example below: | 
|---|
| 42 | \code | 
|---|
| 43 | // We use the same table as in the example above: | 
|---|
| 44 | DataTableRow row = dt.EmptyRow(); | 
|---|
| 45 | row[0] = 83.3; row[1] = 20.;   row[2] = 83.3*83.3+20.*20.; | 
|---|
| 46 | dt.AddRow(row); | 
|---|
| 47 | row["X0_f"] = 7.5; row["X1_f"] = 6.;  row["X0X0pX1X1_d"] = 7.5*7.5+6.*6.; | 
|---|
| 48 | dt.AddRow(row); | 
|---|
| 49 | \endcode | 
|---|
| 50 | */ | 
|---|
| 51 |  | 
|---|
| 52 |  | 
|---|
| 53 | //! Default constructor - optional specification of block (or segment) size | 
|---|
| 54 | DataTable::DataTable(sa_size_t segsz) | 
|---|
| 55 | : BaseDataTable(segsz) | 
|---|
| 56 | { | 
|---|
| 57 | } | 
|---|
| 58 |  | 
|---|
| 59 | //! copy constructor | 
|---|
| 60 | /*! | 
|---|
| 61 | The copy constructur shares the data if \b share=true. | 
|---|
| 62 | Otherwise, the Clone() method is called to make a complete copy. | 
|---|
| 63 | */ | 
|---|
| 64 | DataTable::DataTable(DataTable const & a, bool share) | 
|---|
| 65 | : BaseDataTable(a.SegmentSize()) | 
|---|
| 66 | { | 
|---|
| 67 | if (share) Share(a); | 
|---|
| 68 | else Clone(a); | 
|---|
| 69 | mNEnt = a.mNEnt; | 
|---|
| 70 | mNSeg = a.mNSeg; | 
|---|
| 71 | if (a.mInfo)  mInfo = new DVList(*(a.mInfo)); | 
|---|
| 72 | } | 
|---|
| 73 | //! Copy the table structure from \b a and shares the data (columns content) | 
|---|
| 74 | void DataTable::Share(DataTable const & a) | 
|---|
| 75 | { | 
|---|
| 76 | // On copie la structure de table | 
|---|
| 77 | CopyStructure(a); | 
|---|
| 78 | // Et on partage les donnees des colonnes | 
|---|
| 79 | for (size_t kk=0; kk<mNames.size(); kk++) { | 
|---|
| 80 | sa_size_t sk = mNames[kk].ser; | 
|---|
| 81 | sa_size_t ska = a.mNames[kk].ser; | 
|---|
| 82 | switch (mNames[kk].type) { | 
|---|
| 83 | case IntegerField : | 
|---|
| 84 | mICols[sk].Share(a.mICols[ska]); | 
|---|
| 85 | break; | 
|---|
| 86 | case LongField : | 
|---|
| 87 | mLCols[sk].Share(a.mLCols[ska]); | 
|---|
| 88 | break; | 
|---|
| 89 | case FloatField : | 
|---|
| 90 | mFCols[sk].Share(a.mFCols[ska]); | 
|---|
| 91 | break; | 
|---|
| 92 | case DoubleField : | 
|---|
| 93 | case DateTimeField : | 
|---|
| 94 | mDCols[sk].Share(a.mDCols[ska]); | 
|---|
| 95 | break; | 
|---|
| 96 | case ComplexField : | 
|---|
| 97 | mYCols[sk].Share(a.mYCols[ska]); | 
|---|
| 98 | break; | 
|---|
| 99 | case DoubleComplexField : | 
|---|
| 100 | mZCols[sk].Share(a.mZCols[ska]); | 
|---|
| 101 | break; | 
|---|
| 102 | case StringField : | 
|---|
| 103 | mSCols[sk].Share(a.mSCols[ska]); | 
|---|
| 104 | break; | 
|---|
| 105 | default: | 
|---|
| 106 | throw ForbiddenError("DataTable::Share() : unknown column type "); | 
|---|
| 107 | break; | 
|---|
| 108 | } | 
|---|
| 109 | } | 
|---|
| 110 | } | 
|---|
| 111 |  | 
|---|
| 112 | //! Copy the table structure from \b a and duplicate (copy) the data (columns content) | 
|---|
| 113 | void DataTable::Clone(DataTable const & a) | 
|---|
| 114 | { | 
|---|
| 115 | // On copie la structure de table | 
|---|
| 116 | CopyStructure(a); | 
|---|
| 117 | // Et on partage les donnees des colonnes | 
|---|
| 118 | for (size_t kk=0; kk<mNames.size(); kk++) { | 
|---|
| 119 | sa_size_t sk = mNames[kk].ser; | 
|---|
| 120 | sa_size_t ska = a.mNames[kk].ser; | 
|---|
| 121 | switch (mNames[kk].type) { | 
|---|
| 122 | case IntegerField : | 
|---|
| 123 | mICols[sk].Clone(a.mICols[ska], true); | 
|---|
| 124 | break; | 
|---|
| 125 | case LongField : | 
|---|
| 126 | mLCols[sk].Clone(a.mLCols[ska], true); | 
|---|
| 127 | break; | 
|---|
| 128 | case FloatField : | 
|---|
| 129 | mFCols[sk].Clone(a.mFCols[ska], true); | 
|---|
| 130 | break; | 
|---|
| 131 | case DoubleField : | 
|---|
| 132 | case DateTimeField : | 
|---|
| 133 | mDCols[sk].Clone(a.mDCols[ska], true); | 
|---|
| 134 | break; | 
|---|
| 135 | case ComplexField : | 
|---|
| 136 | mYCols[sk].Clone(a.mYCols[ska], true); | 
|---|
| 137 | break; | 
|---|
| 138 | case DoubleComplexField : | 
|---|
| 139 | mZCols[sk].Clone(a.mZCols[ska], true); | 
|---|
| 140 | break; | 
|---|
| 141 | case StringField : | 
|---|
| 142 | mSCols[sk].Clone(a.mSCols[ska], true); | 
|---|
| 143 | break; | 
|---|
| 144 | default: | 
|---|
| 145 | throw ForbiddenError("DataTable::Clone() : unknown column type "); | 
|---|
| 146 | break; | 
|---|
| 147 | } | 
|---|
| 148 | } | 
|---|
| 149 | } | 
|---|
| 150 | //! Reset (/clear) the table content and structure | 
|---|
| 151 | void DataTable::Clear() | 
|---|
| 152 | { | 
|---|
| 153 | if ( (NVar() == 0) && (NEntry() == 0)) return; | 
|---|
| 154 | mNEnt = 0; | 
|---|
| 155 | mNSeg = 0; | 
|---|
| 156 | if (mVarD) delete[] mVarD; | 
|---|
| 157 | mVarD = NULL; | 
|---|
| 158 | if (mVarMTV) delete[] mVarMTV; | 
|---|
| 159 | mVarMTV = NULL; | 
|---|
| 160 | mNames.clear(); | 
|---|
| 161 | if (mInfo) delete mInfo; | 
|---|
| 162 | mInfo = NULL; | 
|---|
| 163 | mMin.clear(); | 
|---|
| 164 | mMax.clear(); | 
|---|
| 165 | mMinMaxNEnt.clear(); | 
|---|
| 166 | mIColsP.clear(); | 
|---|
| 167 | mLColsP.clear(); | 
|---|
| 168 | mFColsP.clear(); | 
|---|
| 169 | mDColsP.clear(); | 
|---|
| 170 | mYColsP.clear(); | 
|---|
| 171 | mZColsP.clear(); | 
|---|
| 172 | mSColsP.clear(); | 
|---|
| 173 |  | 
|---|
| 174 | mIColIdx.clear(); | 
|---|
| 175 | mLColIdx.clear(); | 
|---|
| 176 | mFColIdx.clear(); | 
|---|
| 177 | mDColIdx.clear(); | 
|---|
| 178 | mYColIdx.clear(); | 
|---|
| 179 | mZColIdx.clear(); | 
|---|
| 180 | mSColIdx.clear(); | 
|---|
| 181 |  | 
|---|
| 182 | mICols.clear(); | 
|---|
| 183 | mLCols.clear(); | 
|---|
| 184 | mFCols.clear(); | 
|---|
| 185 | mDCols.clear(); | 
|---|
| 186 | mYCols.clear(); | 
|---|
| 187 | mZCols.clear(); | 
|---|
| 188 | mSCols.clear(); | 
|---|
| 189 | } | 
|---|
| 190 |  | 
|---|
| 191 |  | 
|---|
| 192 | /*! | 
|---|
| 193 | Implements the action defined in the BaseDataTable interface. | 
|---|
| 194 | In the current implementation, throws an exception (ParmError) | 
|---|
| 195 | if the table contains some data already. | 
|---|
| 196 | */ | 
|---|
| 197 | sa_size_t DataTable::AddColumn(FieldType ft, string const & cnom) | 
|---|
| 198 | { | 
|---|
| 199 | if (NEntry() > 0) | 
|---|
| 200 | throw ParmError("DataTable::AddColumn() Table contains already data "); | 
|---|
| 201 | CheckColName(cnom); | 
|---|
| 202 | sa_size_t ser; | 
|---|
| 203 | sa_size_t idx = NVar(); | 
|---|
| 204 | switch (ft) { | 
|---|
| 205 | case IntegerField : | 
|---|
| 206 | ser = mICols.size(); | 
|---|
| 207 | mICols.push_back(SegDataBlock<int_4>(mSegSz)); | 
|---|
| 208 | mIColIdx.push_back(idx); | 
|---|
| 209 | mIColsP.push_back(NULL); | 
|---|
| 210 | for(sa_size_t kk=0; kk<mICols.size(); kk++) | 
|---|
| 211 | mIColsP[kk] = &(mICols[kk]); | 
|---|
| 212 | break; | 
|---|
| 213 | case LongField : | 
|---|
| 214 | ser = mLCols.size(); | 
|---|
| 215 | mLCols.push_back(SegDataBlock<int_8>(mSegSz)); | 
|---|
| 216 | mLColIdx.push_back(idx); | 
|---|
| 217 | mLColsP.push_back(NULL); | 
|---|
| 218 | for(sa_size_t kk=0; kk<mLCols.size(); kk++) | 
|---|
| 219 | mLColsP[kk] = &(mLCols[kk]); | 
|---|
| 220 | break; | 
|---|
| 221 | case FloatField : | 
|---|
| 222 | ser = mFCols.size(); | 
|---|
| 223 | mFCols.push_back(SegDataBlock<r_4>(mSegSz)); | 
|---|
| 224 | mFColIdx.push_back(idx); | 
|---|
| 225 | mFColsP.push_back(NULL); | 
|---|
| 226 | for(sa_size_t kk=0; kk<mFCols.size(); kk++) | 
|---|
| 227 | mFColsP[kk] = &(mFCols[kk]); | 
|---|
| 228 | break; | 
|---|
| 229 | case DoubleField : | 
|---|
| 230 | case DateTimeField : | 
|---|
| 231 | ser = mDCols.size(); | 
|---|
| 232 | mDCols.push_back(SegDataBlock<r_8>(mSegSz)); | 
|---|
| 233 | mDColIdx.push_back(idx); | 
|---|
| 234 | mDColsP.push_back(NULL); | 
|---|
| 235 | for(sa_size_t kk=0; kk<mDCols.size(); kk++) | 
|---|
| 236 | mDColsP[kk] = &(mDCols[kk]); | 
|---|
| 237 | break; | 
|---|
| 238 | case ComplexField : | 
|---|
| 239 | ser = mYCols.size(); | 
|---|
| 240 | mYCols.push_back(SegDataBlock< complex<r_4> >(mSegSz)); | 
|---|
| 241 | mYColIdx.push_back(idx); | 
|---|
| 242 | mYColsP.push_back(NULL); | 
|---|
| 243 | for(sa_size_t kk=0; kk<mYCols.size(); kk++) | 
|---|
| 244 | mYColsP[kk] = &(mYCols[kk]); | 
|---|
| 245 | break; | 
|---|
| 246 | case DoubleComplexField : | 
|---|
| 247 | ser = mZCols.size(); | 
|---|
| 248 | mZCols.push_back(SegDataBlock< complex<r_8> >(mSegSz)); | 
|---|
| 249 | mZColIdx.push_back(idx); | 
|---|
| 250 | mZColsP.push_back(NULL); | 
|---|
| 251 | for(sa_size_t kk=0; kk<mZCols.size(); kk++) | 
|---|
| 252 | mZColsP[kk] = &(mZCols[kk]); | 
|---|
| 253 | break; | 
|---|
| 254 | case StringField : | 
|---|
| 255 | ser = mSCols.size(); | 
|---|
| 256 | mSCols.push_back(SegDataBlock<string>(mSegSz)); | 
|---|
| 257 | mSColIdx.push_back(idx); | 
|---|
| 258 | mSColsP.push_back(NULL); | 
|---|
| 259 | for(sa_size_t kk=0; kk<mSCols.size(); kk++) | 
|---|
| 260 | mSColsP[kk] = &(mSCols[kk]); | 
|---|
| 261 | break; | 
|---|
| 262 | default: | 
|---|
| 263 | throw ParmError("DataTable::AddColumn() unknown field type "); | 
|---|
| 264 | break; | 
|---|
| 265 | } | 
|---|
| 266 | colst col; | 
|---|
| 267 | col.nom = cnom; | 
|---|
| 268 | col.type = ft; | 
|---|
| 269 | col.ser = ser; | 
|---|
| 270 | mNames.push_back(col); | 
|---|
| 271 |  | 
|---|
| 272 | // Tableaux de calcul min-max | 
|---|
| 273 | mMin.push_back(0.); | 
|---|
| 274 | mMax.push_back(0.); | 
|---|
| 275 | mMinMaxNEnt.push_back(0); | 
|---|
| 276 |  | 
|---|
| 277 | return NVar(); | 
|---|
| 278 | } | 
|---|
| 279 |  | 
|---|