00001 #include <stdlib.h>
00002 #include "timecontr.h"
00003
00004 timecontr::timecontr (void)
00005 {
00006
00007 time=0.0;
00008
00009 start_time=0.0;
00010
00011 end_time=0.0;
00012
00013
00014 forwarddt=0.0;
00015
00016 bfdt=0.0;
00017
00018 bfdtu=0.0;
00019
00020 bsdt=0.0;
00021
00022 backwarddt=0.0;
00023
00024
00025 nit=0;
00026
00027 apit=0;
00028
00029 imptime=NULL;
00030
00031 it_step = -1.0;
00032
00033 iiit=0;
00034
00035
00036 tct=notct;
00037
00038 dtmin=0.0;
00039
00040 dtmax=0.0;
00041
00042
00043 dtminfun = NULL;
00044
00045
00046 dtmaxfun = NULL;
00047
00048 }
00049
00050 timecontr::~timecontr (void)
00051 {
00052 delete [] imptime;
00053 delete dtminfun;
00054 delete dtmaxfun;
00055 }
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065 void timecontr::read (XFILE *in)
00066 {
00067 long i;
00068
00069
00070
00071
00072
00073
00074 xfscanf (in,"%k%m","time_contr_type",&timecontrtype_kwdset,(int*)&tct);
00075
00076
00077 xfscanf (in,"%k%lf%k%lf","start_time",&start_time,"end_time",&end_time);
00078
00079
00080 xfscanf (in,"%k%ld","num_imp_times", &nit);
00081
00082
00083
00084 if (nit == -1){
00085 xfscanf (in,"%lf",&it_step);
00086 nit = long((end_time - start_time)/it_step);
00087 if (nit < 2)
00088 nit = 2;
00089
00090 imptime=new double [nit];
00091 imptime[0] = start_time + it_step;
00092
00093 for (i=1; i<nit-1; i++)
00094 imptime[i] = imptime[i-1] + it_step;
00095 imptime[nit-1] = end_time;
00096 }
00097 else{
00098 if (nit==0){
00099 imptime = new double [1];
00100 imptime[0]=end_time;
00101 }
00102 else{
00103 imptime=new double [nit];
00104 }
00105 for (i=0;i<nit;i++){
00106 xfscanf (in,"%lf",imptime+i);
00107 }
00108 }
00109
00110
00111 timefun.read (in);
00112
00113 switch (tct){
00114 case fixed:{
00115 dtmin = dtmax = initialtimeincr ();
00116 break;
00117 }
00118 case adaptive:{
00119 xfscanf (in,"%k%lf %k%lf","dtmin",&dtmin,"dtmax",&dtmax);
00120 break;
00121 }
00122 case adaptivemin:{
00123 xfscanf (in,"%k%lf","dtmin",&dtmin);
00124 dtmax = initialtimeincr ();
00125 bfdtu = initialtimeincr ();
00126 break;
00127 }
00128 case adaptivemax:{
00129 xfscanf (in,"%k%lf","dtmax",&dtmax);
00130 dtmin = initialtimeincr ();
00131 break;
00132 }
00133 case adaptive_minmax:{
00134 dtminfun = new gfunct();
00135 dtmaxfun = new gfunct();
00136 xfscanf(in, "%k", "dtmin_fun");
00137 dtminfun->read(in);
00138 xfscanf(in, "%k", "dtmax_fun");
00139 dtmaxfun->read(in);
00140 bfdtu = initialtimeincr ();
00141 break;
00142 }
00143 default:{
00144 print_err("unknown type of time controller is required",__FILE__,__LINE__,__func__);
00145 }
00146 }
00147
00148
00149
00150 time = start_time;
00151
00152 forwarddt = initialtimeincr ();
00153
00154 bfdt = forwarddt;
00155 }
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166 void timecontr::print (FILE *out)
00167 {
00168 long i;
00169
00170 fprintf (out,"\n%d\n", tct);
00171
00172 fprintf (out,"%e %e\n",start_time,end_time);
00173
00174
00175 if (it_step >= 0.0)
00176 fprintf (out,"-1 %e\n", it_step);
00177 else
00178 {
00179 fprintf (out,"%ld\n",nit);
00180
00181 for (i=0;i<nit;i++){
00182 fprintf (out," %e \n",imptime[i]);
00183 }
00184 }
00185
00186
00187 timefun.print (out);
00188
00189 switch (tct)
00190 {
00191 case fixed:
00192 break;
00193 case adaptive:
00194 fprintf(out,"%e %e\n", dtmin, dtmax);
00195 break;
00196 case adaptivemin:
00197 fprintf (out, "%e\n", dtmin);
00198 break;
00199 case adaptivemax:
00200 fprintf (out, "%e\n", dtmax);
00201 break;
00202 case adaptive_minmax:
00203 dtminfun->print(out);
00204 fprintf(out, "\n");
00205 dtmaxfun->print(out);
00206 fprintf(out, "\n");
00207 default:
00208 print_err("unknown type of time controller is required",__FILE__,__LINE__,__func__);
00209 }
00210 }
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221 double timecontr::newtime ()
00222 {
00223 double newtime;
00224
00225
00226 forwarddt=timefun.getval(time);
00227
00228
00229 newtime=time+forwarddt;
00230
00231
00232 iiit=0;
00233 if (apit<nit){
00234 if (imptime[apit]<=newtime){
00235 newtime=imptime[apit];
00236 iiit=1;
00237 apit++;
00238 }
00239
00240
00241 forwarddt=newtime-time;
00242
00243 time=newtime;
00244 }
00245 else{
00246
00247 time=newtime;
00248 }
00249
00250
00251 backwarddt=forwarddt;
00252 return time;
00253 }
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266 double timecontr::newtime (double &dt)
00267 {
00268 double newtime, forwarddtu;
00269
00270
00271
00272 forwarddtu=timefun.getval(time);
00273
00274 if (tct == fixed)
00275 {
00276
00277 if ((dt < forwarddtu) && (dt < bsdt))
00278 {
00279
00280
00281
00282 print_err("time increment %le required by solver is less\n"
00283 " than fixed time increment %le specified by user.\n"
00284 " Use adaptive time control instead of fixed one",
00285 __FILE__, __LINE__, __func__, dt, forwarddtu);
00286 abort();
00287 }
00288
00289 forwarddt = forwarddtu;
00290 }
00291 else
00292 {
00293
00294
00295 forwarddt=dt;
00296 if (iiit==1){
00297
00298
00299
00300
00301 if (dt >= bsdt)
00302 forwarddt=bfdt;
00303 }
00304
00305
00306 if (tct == adaptivemin)
00307 {
00308
00309 dtmax = forwarddtu;
00310
00311 if ((bfdtu < forwarddtu) && (dt >= bsdt))
00312 {
00313
00314
00315 forwarddt = forwarddtu;
00316 bfdtu = forwarddtu;
00317 }
00318
00319 }
00320
00321 if (tct == adaptivemax)
00322 dtmin = forwarddtu;
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332
00333
00334
00335
00336
00337
00338
00339
00340
00341
00342
00343 }
00344
00345 if (forwarddt>dtmax)
00346 forwarddt = dtmax;
00347
00348
00349 if (forwarddt<dtmin)
00350 {
00351
00352
00353
00354 print_err("time increment %le required by solver is less\n"
00355 " than time increment %le specified by user.\n"
00356 " Use adaptive time control instead of fixed one",
00357 __FILE__, __LINE__, __func__, dt, forwarddtu);
00358 abort();
00359 }
00360
00361
00362
00363 newtime=time+forwarddt;
00364
00365
00366 if (newtime > end_time)
00367 {
00368 newtime=end_time;
00369 forwarddt = newtime - time;
00370 }
00371
00372
00373
00374 bfdt=forwarddt;
00375
00376
00377 iiit=0;
00378 if (apit<nit){
00379 if (imptime[apit]<=newtime){
00380
00381 newtime=imptime[apit];
00382 iiit=1;
00383 apit++;
00384 }
00385
00386
00387 forwarddt=newtime-time;
00388
00389 time=newtime;
00390
00391 }
00392 else{
00393
00394 time=newtime;
00395 }
00396
00397
00398 backwarddt=forwarddt;
00399
00400
00401 dt=forwarddt;
00402
00403
00404 bsdt = dt;
00405
00406 return time;
00407 }
00408
00409
00410
00411
00412
00413
00414
00415
00416 void timecontr::oldtime ()
00417 {
00418 if (apit>0){
00419 if (imptime[apit-1]==time){
00420 apit--;
00421 }
00422 }
00423
00424
00425 time-=backwarddt;
00426
00427 if (imptime[apit]==time)
00428 iiit=1;
00429 else
00430 iiit=0;
00431
00432 }
00433
00434
00435
00436
00437
00438
00439
00440
00441 double timecontr::starttime ()
00442 {
00443 return start_time;
00444 }
00445
00446
00447
00448
00449
00450
00451
00452
00453 double timecontr::endtime ()
00454 {
00455 return end_time;
00456 }
00457
00458
00459
00460
00461
00462
00463
00464
00465 double timecontr::actualtime ()
00466 {
00467 return time;
00468 }
00469
00470
00471
00472
00473
00474
00475
00476
00477 double timecontr::actualforwtimeincr ()
00478 {
00479 return forwarddt;
00480 }
00481
00482
00483
00484
00485
00486
00487
00488
00489 double timecontr::actualbacktimeincr ()
00490 {
00491 return backwarddt;
00492 }
00493
00494
00495
00496
00497
00498
00499
00500
00501 double timecontr::initialtimeincr ()
00502 {
00503 double dt;
00504
00505 dt = timefun.getval(start_time);
00506
00507 if (nit==0){
00508 return dt;
00509 }
00510 else{
00511 if (start_time+dt<imptime[0]){
00512 return dt;
00513 }
00514 else{
00515 dt = imptime[0]-start_time;
00516 return dt;
00517 }
00518 }
00519 }
00520
00521
00522
00523
00524
00525
00526
00527
00528
00529
00530 void timecontr::take_values (timecontr &tc)
00531 {
00532 time = tc.time;
00533 backwarddt = tc.backwarddt;
00534 forwarddt = tc.forwarddt;
00535 start_time = tc.start_time;
00536 end_time = tc.end_time;
00537 }
00538
00539
00540
00541
00542
00543
00544
00545
00546 void timecontr::seconds_days ()
00547 {
00548 time = time/86400.0;
00549 backwarddt = backwarddt/86400.0;
00550 forwarddt = forwarddt/86400.0;
00551 start_time = start_time/86400.0;
00552 end_time = end_time/86400.0;
00553 }
00554
00555
00556
00557
00558
00559
00560
00561
00562
00563 long timecontr::isitimptime ()
00564 {
00565 return iiit;
00566 }
00567
00568
00569
00570
00571
00572
00573
00574
00575
00576
00577 void timecontr::initiate (timecontr &tc)
00578 {
00579 long i;
00580
00581 start_time = tc.start_time;
00582 end_time = tc.end_time;
00583
00584 timefun.initiate (tc.timefun);
00585
00586 nit = tc.nit;
00587 imptime = new double [nit];
00588 for (i=0;i<nit;i++){
00589 imptime[i]=tc.imptime[i];
00590 }
00591
00592 time = start_time;
00593
00594
00595 forwarddt = initialtimeincr ();
00596 }
00597
00598
00599
00600
00601
00602
00603
00604
00605
00606
00607
00608
00609
00610
00611
00612
00613
00614
00615
00616
00617
00618
00619
00620
00621
00622
00623
00624
00625
00626
00627
00628
00629
00630
00631
00632
00633
00634
00635
00636
00637
00638 void timecontr::be_copy_of (timecontr &tc)
00639 {
00640 time = tc.time;
00641 backwarddt = tc.backwarddt;
00642 forwarddt = tc.forwarddt;
00643 bfdt = tc.bfdt;
00644 apit = tc.apit;
00645 iiit = tc.iiit;
00646 }
00647
00648
00649
00650
00651
00652
00653
00654
00655
00656
00657
00658
00659
00660 void timecontr::save_txt (FILE *out, int prec)
00661 {
00662 fprintf(out,"%.*le\n",prec,time);
00663 fprintf(out,"%.*le\n",prec,backwarddt);
00664 fprintf(out,"%.*le\n",prec,forwarddt);
00665 fprintf(out,"%.*le\n",prec,bfdt);
00666 fprintf(out,"%ld\n",apit);
00667 fprintf(out,"%ld\n",iiit);
00668 }
00669
00670
00671
00672
00673
00674
00675
00676
00677
00678
00679
00680
00681 void timecontr::save_bin (FILE *out)
00682 {
00683 fwrite(&time, sizeof(time), 1, out);
00684 fwrite(&backwarddt, sizeof(backwarddt), 1, out);
00685 fwrite(&forwarddt, sizeof(forwarddt), 1, out);
00686 fwrite(&bfdt, sizeof(bfdt), 1, out);
00687 fwrite(&apit, sizeof(apit), 1, out);
00688 fwrite(&iiit, sizeof(iiit), 1, out);
00689 }
00690
00691
00692
00693
00694
00695
00696
00697
00698
00699
00700
00701
00702 void timecontr::restore_txt (FILE *in)
00703 {
00704 fscanf(in,"%le",&time);
00705 fscanf(in,"%le",&backwarddt);
00706 fscanf(in,"%le",&forwarddt);
00707 fscanf(in,"%le",&bfdt);
00708 fscanf(in,"%ld",&apit);
00709 fscanf(in,"%ld",&iiit);
00710 }
00711
00712
00713
00714
00715
00716
00717
00718
00719
00720
00721
00722
00723 void timecontr::restore_bin (FILE *in)
00724 {
00725 fread(&time, sizeof(time), 1, in);
00726 fread(&backwarddt, sizeof(backwarddt), 1, in);
00727 fread(&forwarddt, sizeof(forwarddt), 1, in);
00728 fread(&bfdt, sizeof(bfdt), 1, in);
00729 fread(&apit, sizeof(apit), 1, in);
00730 fread(&iiit, sizeof(iiit), 1, in);
00731 }
00732
00733
00734