source: Sophya/trunk/SophyaProg/Tests/tnt.cc@ 3399

Last change on this file since 3399 was 2852, checked in by ansari, 20 years ago

Amelioration prog test DataTable - Reza 21/11/2005

File size: 14.6 KB
Line 
1#include <stdio.h>
2#include <stdlib.h>
3#include <math.h>
4#include <iostream>
5#include <fstream>
6
7#include "sopnamsp.h"
8#include "histinit.h"
9#include "dvlist.h"
10#include "ntuple.h"
11#include "xntuple.h"
12#include "datatable.h"
13#include "swppfdtable.h"
14
15/* Programme test des classes DVList, NTuple,, XNTuple */
16/* DataTable et les handlers PPF */
17/* SOPHYA - R. Ansari (LAL) - 2000 - 2005 */
18
19void test_dvl();
20void test_ntuple();
21void test_xntuple() ;
22void test_Xntuple() ;
23void test_DataTable() ;
24void test_SwPPFDataTable() ;
25
26int main(int narg, char *arg[])
27{
28 SophyaInit();
29 if (narg < 2) {
30 cout << " tnt/Erreur arg - Usage: tnt d/n/x/X/DT/SWDT \n"
31 << " d:DVList n:NTuple x/X:XNTuple DT,SWDT: SwPPF/DataTable test \n" << endl;
32 exit(0);
33 }
34 int rc = 0;
35 try {
36 if (*arg[1] == 'd') test_dvl();
37 else if (*arg[1] == 'n') test_ntuple();
38 else if (*arg[1] == 'x') test_xntuple();
39 else if (*arg[1] == 'X') test_Xntuple();
40 else if (strcmp(arg[1],"DT") == 0) test_DataTable();
41 else if (strcmp(arg[1],"SWDT") == 0) test_SwPPFDataTable();
42 }
43 catch(PThrowable exc ) {
44 cerr << "tnt-main() , Catched exception: \n" << exc.Msg() << endl;
45 rc = 97;
46 }
47 catch(std::exception ex) {
48 cerr << "tnt-main() , Catched exception ! " << (string)(ex.what()) << endl;
49 rc = 98;
50 }
51 catch(...) {
52 cerr << "tnt-main() , Catched ... ! " << endl;
53 rc = 99;
54 }
55 return rc;
56}
57
58/* ***** Test de NTuple simple ***** */
59void test_ntuple()
60{
61 char * names[3] = {"XPos", "YPos", "Val"};
62 int i,j, k;
63 float xnt[3];
64
65 cout << "======= test_ntuple: Testing NTuple ======= " << endl;
66
67 cout << "Creation Ntuple avec X,Y,Val ... " << endl;
68 NTuple nt1(3, names, 20);
69
70 k = 0;
71 for(j=0; j<8; j++)
72 for(i=0; i<12; i++)
73 { xnt[0] = i+0.5; xnt[1] = j+0.5; xnt[2] = k;
74 nt1.Fill(xnt); k++; }
75
76 nt1.Info().Comment() = "NTuple de Test - Cree par tnt.cc";
77 nt1.Info()["Version"] = SophyaVersion();
78 nt1.Show(cout);
79 nt1.Print(0, 5);
80 nt1.Print(18, 5);
81 nt1.Print(94, 5);
82
83 string fn = "nt.ppf";
84 {
85 cout << "Ecriture NTuple ds nt.ppf ... \n" << endl;
86 ObjFileIO<NTuple> fio(&nt1);
87 fio.Write(fn);
88 }
89
90 {
91 cout << "Lecture NTuple (nt2) ds nt.ppf ... \n" << endl;
92 PInPersist fsi(fn);
93 NTuple nt2;
94 fsi >> nt2;
95 nt2.Show(cout);
96 nt2.Print(0, 5);
97 nt2.Print(18, 5);
98 nt2.Print(94, 5);
99 }
100
101
102}
103
104/* ***** Test de dvlist ***** */
105void test_dvl()
106{
107DVList dvl;
108
109cout << "\n ======= test_dvl: Testing MuTyV ======= " << endl;
110MuTyV zvs = " ( 1.41 -2.71) ";
111MuTyV dvs = "434.898";
112MuTyV fvu = 314.1596;
113MuTyV ivu = 7654321;
114cout << " float->string: fvu= " << fvu << " (string)fvu=" << (string)fvu << endl;
115cout << " int->string: ivu= " << ivu << " (string)ivu=" << (string)ivu << endl;
116 complex<double> zzd = zvs;
117cout << "String->complex<double>: zvs = " << zvs
118 << " (complex<double>)zvs= " << zzd << endl;
119cout << "String->double: dvs = " << dvs
120 << " (double)dvs= " << (double)dvs << endl;
121TimeStamp ts;
122MuTyV tv = ts;
123cout << "MuTyV=TimeStamp->String: tv = " << tv
124 << " ->double= " << (double)tv << " ->TimeStamp:"
125 << TimeStamp((double)tv).ToString() << endl;
126tv = "1995-09-23T15:20:32";
127string tts = (string)tv;
128cout << "MuTyV=TimeStamp tv=" << tv << " ->TimeStamp(tv)="
129 << TimeStamp(tts) << endl;
130
131cout << "\n\n ======= test_dvl: Testing DVList ======= " << endl;
132dvl.SetI("Var XXX", 12345);
133dvl.SetI("IV1-80", 80);
134dvl.SetI("IV2-330", 330);
135
136dvl.SetD("DV1-0.2", 0.2);
137dvl.SetD("DV2-4.5", 4.5);
138dvl.SetZ("ZV", complex<r_8>(2.0,-1.44));
139
140dvl.SetI("IVV3-783", 783);
141dvl("IVV3-783-O") = 7783;
142// dvl["Avec[]"] = "operateur [] !";
143dvl.SetD("DVV3", 3.141592652141592652);
144dvl.SetComment("DVV3", "Comment for DVV3, r_8 type variable");
145dvl("DVV3avec()") = 44.555e-8;
146
147dvl.SetS("Blanc_White", "123.456Ma premiere chaine");
148dvl.SetI("IntegerValue", 55777);
149dvl.SetComment("IntegerValue", "This variable has a comment");
150dvl.Comment() = "This is a test DVL produced by the program tdvl.cc \n Using SOPHYA , Feb 2000";
151
152dvl["Sinf"] = "inf 0985";
153dvl["ToDay"] = ts;
154dvl.Print();
155
156
157double d = dvl("DV2-4.5");
158float f = dvl("DVV3");
159int i = dvl("IVV3-783");
160complex<r_8> z = dvl("ZV");
161
162printf("\n \n Essai1 (IVV3 DVV3 DV2= ) %d %.20g %g \n", i, f, (float)d);
163printf("\n \n Essai ZV= (%.2g %g I) \n", z.real(), z.imag());
164cout << "Test Comment/IntegerValue: " << dvl.GetComment("IntegerValue") << endl;
165cout << "Test Comment/DVV3: " << dvl.GetComment("DVV3") << endl;
166
167cout << "Test string recup " << (string)(dvl["Blanc_White"]) << endl;
168cout << "Test string recup(int=80) " << (string)(dvl["IV1-80"]) << endl;
169dvl("DVV3") = (double)3.141592652141592652;
170cout << "Test string recup(double=Pi..i) " << (string)(dvl["DVV3"]) << endl;
171
172 {
173 cout << " Writing DVList in file PPF dvl.ppf " << endl;
174 POutPersist os("dvl.ppf");
175 os << dvl ;
176 }
177
178 cout << "-------------------------- \n\n"
179 << " reading DVList from file dvl.ppf ... \n" << endl;
180
181 DVList dvlr("dvl.ppf");
182 double df1 = (double)( dvlr["DVV3"] );
183 double df2 = (double)3.141592652141592652 - df1;
184 cout << " Test Precision : Pi-Pi= " << df2 << "DVV3= " << df1 << endl;
185
186 cout << dvlr;
187
188}
189void test_DataTable()
190{
191 cout << "======= test_DataTable: simple DataTable test ======= " << endl;
192 DataTable dt(64);
193 dt.AddIntegerColumn("line");
194 dt.AddStringColumn("sline");
195 dt.AddDoubleColumn("x");
196 dt.AddFloatColumn("f_sin");
197 dt.AddFloatColumn("f_cos");
198 dt.AddDoubleColumn("f_sinxcos");
199 dt.AddDoubleColumn("f_x2");
200 dt.AddComplexColumn("cmplx_cos_sin");
201 dt.AddDoubleComplexColumn("dcmplx_cos_sin");
202 dt.AddDateTimeColumn("datime");
203 MuTyV rec[10];
204 TimeStamp ts("01/3/2005","0:0:0");
205 TimeStamp ts2(0.);
206
207 cout << " 1/ First 1000 lines .... " << endl;
208 for(int k = 0; k<1000; k++) {
209 rec[0] = k;
210 string sline = "L-";
211 sline += (string)rec[0];
212 rec[1] = sline;
213 double x = M_PI*k/100.;
214 double fx = sin(x)*cos(x);
215 double sx = sin(x);
216 double cx = cos(x);
217 rec[2] = x;
218 rec[3] = sx;
219 rec[4] = cx;
220 rec[5] = fx;
221 rec[6] = x*x;
222 rec[7] = complex<r_4>(cx, sx);
223 rec[8] = complex<r_8>(cx, sx);
224 ts2.Set(ts.ToDays()+(double)k/24.);
225 rec[9] = ts2;
226 dt.AddLine(rec);
227 }
228 cout << "1.b/ Use of = operator to make copy of dt " << endl;
229 DataTable dtc1;
230 dtc1 = dt;
231
232 dt.Show();
233 cout << " 2/ Lines 1000-2000 .... " << endl;
234 for(int k = 1000; k<2000; k++) {
235 rec[0] = k;
236 string sline = "L-";
237 sline += (string)rec[0];
238 rec[1] = sline;
239 double x = M_PI*k/100.;
240 double fx = sin(x)*cos(x);
241 double sx = sin(x);
242 double cx = cos(x);
243 rec[2] = x;
244 rec[3] = sx;
245 rec[4] = cx;
246 rec[5] = fx;
247 rec[6] = x*x;
248 rec[7] = complex<r_4>(cx, sx);
249 rec[8] = complex<r_8>(cx, sx);
250 ts2.Set(ts.ToDays()+(double)k/24.);
251 rec[9] = ts2;
252 dt.AddLine(rec);
253 }
254 cout << "2.b/ dt.Show(); : " << endl;
255 dt.Show();
256 cout << "2.c/ dtc1.Show(); : " << endl;
257 dtc1.Show();
258 cout << "2.d/ dt.LineHeaderToString() dt.LineToString(k) : " << endl;
259 cout << dt.LineHeaderToString() << endl;
260 for(int k = 0; k<1500; k+=75)
261 cout << "Line[" << k << "] " << dt.LineToString(k) << endl ;
262 {
263 cout << "3/ Writing DataTable dt to PPF stream dtable.ppf " << endl;
264 POutPersist po("dtable.ppf");
265 po << dt;
266 }
267 {
268 cout << "4/ Reading DataTable dtr from PPF stream dtable.ppf " << endl;
269 PInPersist pi("dtable.ppf");
270 DataTable dtr;
271 pi >> dtr;
272 cout << "4.b/ cout << dtr; " << endl;
273 cout << dtr;
274 cout << "4.c/ dtr.LineHeaderToString() dtr.LineToString(k) : " << endl;
275 cout << dtr.LineHeaderToString() << endl;
276 for(int k = 0; k<1500; k+=75)
277 cout << "Line[" << k << "] " << dtr.LineToString(k) << endl ;
278
279 }
280 {
281 cout << "5/ Reading DataTable dtr from text file dtable.txt " << endl;
282 {
283 cout << "5.a Writing lines to file dtable.txt " << endl;
284 ofstream os("dtable.txt");
285 dt.Print(os,0,450);
286 }
287 cout << "5.b CopyStructure + FillFromASCIIFile from dtable.txt " << endl;
288 ifstream is("dtable.txt");
289 DataTable dtr;
290 dtr.CopyStructure(dt);
291 dtr.FillFromASCIIFile(is);
292 cout << dtr;
293 cout << "5.c/ dtr.LineHeaderToString() dtr.LineToString(k) : " << endl;
294 cout << dtr.LineHeaderToString() << endl;
295 for(int k = 0; k<10; k++)
296 cout << "Line[" << k << "] " << dtr.LineToString(k) << endl ;
297 for(int k = 50; k<300; k+=25)
298 cout << "Line[" << k << "] " << dtr.LineToString(k) << endl ;
299
300 }
301 cout << "6./Test with DataTableRow object " << endl;
302 DataTable dtrow(64);
303 dtrow.AddDoubleColumn("x");
304 dtrow.AddFloatColumn("f_cos");
305 dtrow.AddStringColumn("sline");
306 dtrow.AddIntegerColumn("line");
307 dtrow.AddDateTimeColumn("datime");
308
309 DataTableRow row = dtrow.EmptyRow();
310 for(int k = 0; k<25; k++) {
311 MuTyV mtk = k;
312 string sline = "L-";
313 sline += (string)mtk;
314 row["sline"] = sline;
315 row["line"] = k;
316 double x = M_PI*k/25.;
317 double cx = cos(x);
318 row[0] = x;
319 row[1] = cx;
320 ts2.Set(ts.ToDays()+(double)k);
321 row["datime"] = ts2;
322 dtrow.AddRow(row);
323 }
324 cout << dtrow << endl;
325 cout << dtrow.LineHeaderToString() << endl;
326 for(int k = 0; k<25; k+=3) {
327 dtrow.GetRow(k, row);
328 cout << "k=" << k << " GetRow(): " << row << endl;
329 cout << " ...LineToString(): " << dtrow.LineToString(k) << endl;
330 }
331 cout << endl;
332 cout << "============ FIN test_DataTable ============== " << endl;
333}
334
335void test_SwPPFDataTable()
336{
337 cout << "======= test_SwPPFDataTable(): simple SwPPFDataTable test ======= " << endl;
338 {
339 POutPersist po("swdtable.ppf");
340 SwPPFDataTable dt(po, 64);
341 dt.AddIntegerColumn("line");
342 dt.AddStringColumn("sline");
343 dt.AddDoubleColumn("x");
344 dt.AddFloatColumn("f_sin");
345 dt.AddFloatColumn("f_cos");
346 dt.AddDoubleColumn("f_sinxcos");
347 dt.AddDoubleColumn("f_x2");
348 dt.AddComplexColumn("cmplx_cos_sin");
349 dt.AddDoubleComplexColumn("dcmplx_cos_sin");
350 dt.AddDateTimeColumn("datime");
351 MuTyV rec[10];
352 TimeStamp ts("01/3/2005","0:0:0");
353 TimeStamp ts2(0.);
354 cout << " 1/ First 1000 lines .... " << endl;
355 for(int k = 0; k<1000; k++) {
356 rec[0] = k;
357 string sline = "L-";
358 sline += (string)rec[0];
359 rec[1] = sline;
360 double x = M_PI*k/100.;
361 double fx = sin(x)*cos(x);
362 double sx = sin(x);
363 double cx = cos(x);
364 rec[2] = x;
365 rec[3] = sx;
366 rec[4] = cx;
367 rec[5] = fx;
368 rec[6] = x*x;
369 rec[7] = complex<r_4>(cx, sx);
370 rec[8] = complex<r_8>(cx, sx);
371 ts2.Set(ts.ToDays()+(double)k/24.);
372 rec[9] = ts2;
373 dt.AddLine(rec);
374 }
375 cout << "1.b/ Use of = operator to make copy of dt " << endl;
376 SwPPFDataTable dtc1;
377 dtc1 = dt;
378 dt.Show();
379 cout << " 2/ Lines 1000-2000 .... " << endl;
380 for(int k = 1000; k<2000; k++) {
381 rec[0] = k;
382 string sline = "L-";
383 sline += (string)rec[0];
384 rec[1] = sline;
385 double x = M_PI*k/100.;
386 double fx = sin(x)*cos(x);
387 double sx = sin(x);
388 double cx = cos(x);
389 rec[2] = x;
390 rec[3] = sx;
391 rec[4] = cx;
392 rec[5] = fx;
393 rec[6] = x*x;
394 rec[7] = complex<r_4>(cx, sx);
395 rec[8] = complex<r_8>(cx, sx);
396 ts2.Set(ts.ToDays()+(double)k/24.);
397 rec[9] = ts2;
398 dt.AddLine(rec);
399 }
400 cout << "2.b/ dt.Show(); : " << endl;
401 dt.Show();
402 cout << "2.c/ dtc1.Show(); : " << endl;
403 dtc1.Show();
404 cout << "3/ Writing SwPPFDataTable dt to PPF stream swdtable.ppf " << endl;
405 po << dt;
406 }
407 {
408 cout << "4/ Reading SwPPFDataTable dtr from PPF stream swdtable.ppf " << endl;
409 PInPersist pi("swdtable.ppf");
410 SwPPFDataTable dtr;
411 pi >> dtr;
412 cout << "4.b/ cout << dtr; " << endl;
413 cout << dtr;
414 cout << "4.c/ dtr.LineHeaderToString() dtr.LineToString(k) : " << endl;
415 cout << dtr.LineHeaderToString() << endl;
416 for(int k = 0; k<1500; k+=75)
417 cout << "Line[" << k << "] " << dtr.LineToString(k) << endl ;
418 }
419
420 cout << "============ FIN test_SwPPFDataTable() ======== ======= " << endl;
421}
422
423
424
425
426void test_Xntuple()
427{
428 char* names[] = {"str1", "str2", "str3", "str4", "str5"} ;
429 XNTuple nt(0, 0, 0, 5, names) ;
430 char** ce = new char*[5] ;
431 int i;
432 for(i = 0 ; i < 5 ; i++)
433 ce[i] = new char[20] ;
434
435 cout << "======= test_Xtuple: simple XNTuple test ======= " << endl;
436
437 strncpy(ce[1], "toto a une auto", 20) ;
438 strncpy(ce[2], "titi a une iti", 20) ;
439 strncpy(ce[3], "tutu a une utu", 20) ;
440 strncpy(ce[4], "tata a une ata", 20) ;
441 for(i = 0 ; i < 100000 ; i++) {
442 sprintf(ce[0], "%d", i) ;
443 nt.Fill(NULL, NULL, NULL, ce) ;
444 }
445
446 nt.Show() ;
447 cout << nt.LineHeaderToString() ;
448 cout << nt.LineToString(5027) << endl ;
449
450 char* names2[] = {"d0", "d1", "f0", "f1", "f2", "i0", "str0", "str1"} ;
451 XNTuple nt2(2, 3, 1, 2, names2) ;
452 double de[2] ; float fe[3] ; int ie ;
453 char** ce2 = new char*[2] ;
454 for(i = 0 ; i < 2 ; i++) ce2[i] = new char[20] ;
455 strncpy(ce2[1], "glop glop", 20) ;
456
457 for(i = 0 ; i < 100000 ; i++) {
458 de[0] = i ;
459 de[1] = sin((double)i) ;
460 fe[0] = i ;
461 fe[1] = i * cos((double)i) ;
462 fe[2] = 2*i ;
463 ie = -i;
464 sprintf(ce[0], "%d", i) ;
465 nt2.Fill(de, fe, &ie, ce) ;
466 }
467 nt2.Show() ;
468 nt2.LineHeaderToString() ;
469 // nt2.LineToString(20) ;
470}
471
472
473void test_xntuple()
474{
475 char* names[] = {"dblval", "floval", "intval", "strval"} ;
476 XNTuple nt(1, 1, 1, 1, names) ;
477 double de ;
478 float fe ;
479 int ie ;
480 char* ce = new char[22] ;
481
482
483 cout << "======= test_Xtuple: XNTuple test ======= " << endl;
484 int i;
485 for(i = 0 ; i < nt.NVar() ; i++)
486 printf(" +++ %s <--> %d \n",
487 nt.NomIndex(i).c_str(), nt.IndexNom(nt.NomIndex(i).c_str())) ;
488
489 for(i = 0 ; i < 100000 ; i++) {
490 de = fe = ie = i ;
491 sprintf(ce, "%d", i) ;
492 nt.Fill(&de, &fe, &ie, &ce) ;
493 }
494 nt.Show() ;
495 cout << nt.VarList_C("toto") ;
496 cout << nt.LineHeaderToString() ;
497 cout << nt.LineToString(20) << endl << endl ;
498
499
500 XNTuple nt2 ;
501 nt2.SwapPath() = "/tmp/sop/" ;
502 nt2.Copy(nt) ;
503 nt2.Show() ;
504 for(i = 0 ; i < 100000 ; i++) {
505 de = fe = ie = i ;
506 sprintf(ce, "%d", i) ;
507 nt2.Fill(&de, &fe, &ie, &ce) ;
508 }
509 nt2.Show() ;
510
511 {
512 POutPersist os("xnt.ppf");
513 os << nt2 ;
514 }
515
516 XNTuple::SetSwapPath("/tmp/sop/") ;
517 XNTuple nt3("xnt.ppf") ;
518 nt3.Show() ;
519
520 for(i = 0 ; i < nt3.NEntry() ; i+= 1000)
521 printf("%f %f %d %s\n",
522 nt3.GetDVal(i,0),
523 nt3.GetFVal(i,1),
524 nt3.GetIVal(i,2),
525 nt3.GetSVal(i,3).c_str()) ;
526
527 double min, max ;
528 for(i = 0 ; i < nt3.NVar() ; i++) {
529 nt3.GetMinMax(i, min, max) ;
530 printf("GetMinMax(%s) : %f/%f\n",
531 nt3.NomIndex(i).c_str(), min, max) ;
532 }
533 // nt3.Show() ;
534}
535
Note: See TracBrowser for help on using the repository browser.