LifeV
fsi/examples/challenge_VPH/ud_functions.cpp
Go to the documentation of this file.
1 /* -*- Mode : c++; c-tab-always-indent: t; indent-tabs-mode: nil; -*-
2 
3  <short description here>
4 
5  Gilles Fourestey gilles.fourestey@epfl.ch
6 
7 */
8 /** \file ud_functions.cpp
9 */
10 
11 #include "ud_functions.hpp"
12 //#include "flowConditions.hpp"
13 #include "lifev/core/array/VectorEpetra.hpp"
14 
15 //#define ANEURISM100170
16 #define AORTA
17 namespace LifeV
18 {
19 
20 Real abdominalAorta (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
21 {
22  return 1e2 * (std::cos (6.283 * t) - 1) / 2;
23 }
24 
25 Real aortaFluxIn (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
26 {
27  if (t <= 0.00 + 0.01)
28  {
29  return 0.0000e+00;
30  }
31  if (t <= 0.01 + 0.01)
32  {
33  return 0.0000e+00;
34  }
35  if (t <= 0.02 + 0.01)
36  {
37  return 0.0000e+00;
38  }
39  if (t <= 0.03 + 0.01)
40  {
41  return 0.0000e+00;
42  }
43  if (t <= 0.04 + 0.01)
44  {
45  return -9.1759e-06 * 1e6;
46  }
47  if (t <= 0.05 + 0.01)
48  {
49  return -3.0930e-05 * 1e6;
50  }
51  if (t <= 0.06 + 0.01)
52  {
53  return -6.2639e-05 * 1e6;
54  }
55  if (t <= 0.07 + 0.01)
56  {
57  return -1.0212e-04 * 1e6;
58  }
59  if (t <= 0.08 + 0.01)
60  {
61  return -1.4760e-04 * 1e6;
62  }
63  if (t <= 0.09 + 0.01)
64  {
65  return -1.9726e-04 * 1e6;
66  }
67  if (t <= 0.10 + 0.01)
68  {
69  return -2.4980e-04 * 1e6;
70  }
71  if (t <= 0.11 + 0.01)
72  {
73  return -2.9526e-04 * 1e6;
74  }
75  if (t <= 0.12 + 0.01)
76  {
77  return -3.2956e-04 * 1e6;
78  }
79  if (t <= 0.13 + 0.01)
80  {
81  return -3.5469e-04 * 1e6;
82  }
83  if (t <= 0.14 + 0.01)
84  {
85  return -3.7250e-04 * 1e6;
86  }
87  if (t <= 0.15 + 0.01)
88  {
89  return -3.8429e-04 * 1e6;
90  }
91  if (t <= 0.16 + 0.01)
92  {
93  return -3.9123e-04 * 1e6;
94  }
95  if (t <= 0.17 + 0.01)
96  {
97  return -3.9431e-04 * 1e6;
98  }
99  if (t <= 0.18 + 0.01)
100  {
101  return -3.9349e-04 * 1e6;
102  }
103  if (t <= 0.19 + 0.01)
104  {
105  return -3.8858e-04 * 1e6;
106  }
107  if (t <= 0.20 + 0.01)
108  {
109  return -3.7985e-04 * 1e6;
110  }
111  if (t <= 0.21 + 0.01)
112  {
113  return -3.6756e-04 * 1e6;
114  }
115  if (t <= 0.22 + 0.01)
116  {
117  return -3.5207e-04 * 1e6;
118  }
119  if (t <= 0.23 + 0.01)
120  {
121  return -3.3408e-04 * 1e6;
122  }
123  if (t <= 0.24 + 0.01)
124  {
125  return -3.1402e-04 * 1e6;
126  }
127  if (t <= 0.25 + 0.01)
128  {
129  return -2.9288e-04 * 1e6;
130  }
131  if (t <= 0.26 + 0.01)
132  {
133  return -2.7154e-04 * 1e6;
134  }
135  if (t <= 0.27 + 0.01)
136  {
137  return -2.5054e-04 * 1e6;
138  }
139  if (t <= 0.28 + 0.01)
140  {
141  return -2.2979e-04 * 1e6;
142  }
143  if (t <= 0.29 + 0.01)
144  {
145  return -2.0904e-04 * 1e6;
146  }
147  if (t <= 0.30 + 0.01)
148  {
149  return -1.8880e-04 * 1e6;
150  }
151  if (t <= 0.31 + 0.01)
152  {
153  return -1.6899e-04 * 1e6;
154  }
155  if (t <= 0.32 + 0.01)
156  {
157  return -1.4864e-04 * 1e6;
158  }
159  if (t <= 0.33 + 0.01)
160  {
161  return -1.2730e-04 * 1e6;
162  }
163  if (t <= 0.34 + 0.01)
164  {
165  return -1.0400e-04 * 1e6;
166  }
167  if (t <= 0.35 + 0.01)
168  {
169  return -7.9755e-05 * 1e6;
170  }
171  if (t <= 0.36 + 0.01)
172  {
173  return -5.8719e-05 * 1e6;
174  }
175  if (t <= 0.37 + 0.01)
176  {
177  return -4.0345e-05 * 1e6;
178  }
179  if (t <= 0.38 + 0.01)
180  {
181  return -2.4596e-05 * 1e6;
182  }
183  if (t <= 0.39 + 0.01)
184  {
185  return -1.2259e-05 * 1e6;
186  }
187  if (t <= 0.40 + 0.01)
188  {
189  return -3.8110e-06 * 1e6;
190  }
191  if (t <= 0.41 + 0.01)
192  {
193  return 0.0000e+00;
194  }
195  if (t <= 0.42 + 0.01)
196  {
197  return 0.0000e+00;
198  }
199  if (t <= 0.43 + 0.01)
200  {
201  return 0.0000e+00;
202  }
203  if (t <= 0.44 + 0.01)
204  {
205  return 0.0000e+00;
206  }
207  if (t <= 0.45 + 0.01)
208  {
209  return 0.0000e+00;
210  }
211  if (t <= 0.46 + 0.01)
212  {
213  return 0.0000e+00;
214  }
215  if (t <= 0.47 + 0.01)
216  {
217  return 0.0000e+00;
218  }
219  if (t <= 0.48 + 0.01)
220  {
221  return 0.0000e+00;
222  }
223  if (t <= 0.49 + 0.01)
224  {
225  return 0.0000e+00;
226  }
227  if (t <= 0.50 + 0.01)
228  {
229  return 0.0000e+00;
230  }
231  if (t <= 0.51 + 0.01)
232  {
233  return 0.0000e+00;
234  }
235  if (t <= 0.52 + 0.01)
236  {
237  return 0.0000e+00;
238  }
239  if (t <= 0.53 + 0.01)
240  {
241  return 0.0000e+00;
242  }
243  if (t <= 0.54 + 0.01)
244  {
245  return 0.0000e+00;
246  }
247  if (t <= 0.55 + 0.01)
248  {
249  return 0.0000e+00;
250  }
251  if (t <= 0.56 + 0.01)
252  {
253  return 0.0000e+00;
254  }
255  if (t <= 0.57 + 0.01)
256  {
257  return 0.0000e+00;
258  }
259  if (t <= 0.58 + 0.01)
260  {
261  return 0.0000e+00;
262  }
263  if (t <= 0.59 + 0.01)
264  {
265  return 0.0000e+00;
266  }
267  if (t <= 0.60 + 0.01)
268  {
269  return 0.0000e+00;
270  }
271  if (t <= 0.61 + 0.01)
272  {
273  return 0.0000e+00;
274  }
275  if (t <= 0.62 + 0.01)
276  {
277  return 0.0000e+00;
278  }
279  if (t <= 0.63 + 0.01)
280  {
281  return 0.0000e+00;
282  }
283  if (t <= 0.64 + 0.01)
284  {
285  return 0.0000e+00;
286  }
287  if (t <= 0.65 + 0.01)
288  {
289  return 0.0000e+00;
290  }
291  if (t <= 0.66 + 0.01)
292  {
293  return 0.0000e+00;
294  }
295  if (t <= 0.67 + 0.01)
296  {
297  return 0.0000e+00;
298  }
299  if (t <= 0.68 + 0.01)
300  {
301  return 0.0000e+00;
302  }
303  if (t <= 0.69 + 0.01)
304  {
305  return 0.0000e+00;
306  }
307  if (t <= 0.70 + 0.01)
308  {
309  return 0.0000e+00;
310  }
311  if (t <= 0.71 + 0.01)
312  {
313  return 0.0000e+00;
314  }
315  if (t <= 0.72 + 0.01)
316  {
317  return 0.0000e+00;
318  }
319  if (t <= 0.73 + 0.01)
320  {
321  return 0.0000e+00;
322  }
323  if (t <= 0.74 + 0.01)
324  {
325  return 0.0000e+00;
326  }
327  if (t <= 0.75 + 0.01)
328  {
329  return 0.0000e+00;
330  }
331  if (t <= 0.76 + 0.01)
332  {
333  return 0.0000e+00;
334  }
335  if (t <= 0.77 + 0.01)
336  {
337  return 0.0000e+00;
338  }
339  if (t <= 0.78 + 0.01)
340  {
341  return 0.0000e+00;
342  }
343  if (t <= 0.79 + 0.01)
344  {
345  return 0.0000e+00;
346  }
347 }
348 
349 Real aortaPhisPress (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
350 {
351  /*switch(i) {
352  case 1:
353  return 0.0;
354  break;*/
355  // case 2:
356  if (t <= 0.00)
357  {
358  return -1.e1 * 11017;
359  }
360  if (t <= 0.01)
361  {
362  return -1.e1 * 10954;
363  }
364  if (t <= 0.02)
365  {
366  return -1.e1 * 10893;
367  }
368  if (t <= 0.03)
369  {
370  return -1.e1 * 10832;
371  }
372  if (t <= 0.04)
373  {
374  return -1.e1 * 10771;
375  }
376  if (t <= 0.05)
377  {
378  return -1.e1 * 10712;
379  }
380  if (t <= 0.06)
381  {
382  return -1.e1 * 10653;
383  }
384  if (t <= 0.07)
385  {
386  return -1.e1 * 11113;
387  }
388  if (t <= 0.08)
389  {
390  return -1.e1 * 11544;
391  }
392  if (t <= 0.09)
393  {
394  return -1.e1 * 11869;
395  }
396  if (t <= 0.10)
397  {
398  return -1.e1 * 12146;
399  }
400  if (t <= 0.11)
401  {
402  return -1.e1 * 12394;
403  }
404  if (t <= 0.12)
405  {
406  return -1.e1 * 12635;
407  }
408  if (t <= 0.13)
409  {
410  return -1.e1 * 12889;
411  }
412  if (t <= 0.14)
413  {
414  return -1.e1 * 13151;
415  }
416  if (t <= 0.15)
417  {
418  return -1.e1 * 13398;
419  }
420  if (t <= 0.16)
421  {
422  return -1.e1 * 13620;
423  }
424  if (t <= 0.17)
425  {
426  return -1.e1 * 13833;
427  }
428  if (t <= 0.18)
429  {
430  return -1.e1 * 14035;
431  }
432  if (t <= 0.19)
433  {
434  return -1.e1 * 14229;
435  }
436  if (t <= 0.20)
437  {
438  return -1.e1 * 14436;
439  }
440  if (t <= 0.21)
441  {
442  return -1.e1 * 14613;
443  }
444  if (t <= 0.22)
445  {
446  return -1.e1 * 14753;
447  }
448  if (t <= 0.23)
449  {
450  return -1.e1 * 14878;
451  }
452  if (t <= 0.24)
453  {
454  return -1.e1 * 14974;
455  }
456  if (t <= 0.25)
457  {
458  return -1.e1 * 15032;
459  }
460  if (t <= 0.26)
461  {
462  return -1.e1 * 15047;
463  }
464  if (t <= 0.27)
465  {
466  return -1.e1 * 15025;
467  }
468  if (t <= 0.28)
469  {
470  return -1.e1 * 14975;
471  }
472  if (t <= 0.29)
473  {
474  return -1.e1 * 14899;
475  }
476  if (t <= 0.30)
477  {
478  return -1.e1 * 14822;
479  }
480  if (t <= 0.31)
481  {
482  return -1.e1 * 14721;
483  }
484  if (t <= 0.32)
485  {
486  return -1.e1 * 14594;
487  }
488  if (t <= 0.33)
489  {
490  return -1.e1 * 14496;
491  }
492  if (t <= 0.34)
493  {
494  return -1.e1 * 14375;
495  }
496  if (t <= 0.35)
497  {
498  return -1.e1 * 14198;
499  }
500  if (t <= 0.36)
501  {
502  return -1.e1 * 13990;
503  }
504  if (t <= 0.37)
505  {
506  return -1.e1 * 13726;
507  }
508  if (t <= 0.38)
509  {
510  return -1.e1 * 13397;
511  }
512  if (t <= 0.39)
513  {
514  return -1.e1 * 13167;
515  }
516  if (t <= 0.40)
517  {
518  return -1.e1 * 13132;
519  }
520  if (t <= 0.41)
521  {
522  return -1.e1 * 13315;
523  }
524  if (t <= 0.42)
525  {
526  return -1.e1 * 13271;
527  }
528  if (t <= 0.43)
529  {
530  return -1.e1 * 13157;
531  }
532  if (t <= 0.44)
533  {
534  return -1.e1 * 13028;
535  }
536  if (t <= 0.45)
537  {
538  return -1.e1 * 12975;
539  }
540  if (t <= 0.46)
541  {
542  return -1.e1 * 12933;
543  }
544  if (t <= 0.47)
545  {
546  return -1.e1 * 12891;
547  }
548  if (t <= 0.48)
549  {
550  return -1.e1 * 12836;
551  }
552  if (t <= 0.49)
553  {
554  return -1.e1 * 12768;
555  }
556  if (t <= 0.50)
557  {
558  return -1.e1 * 12700;
559  }
560  if (t <= 0.51)
561  {
562  return -1.e1 * 12641;
563  }
564  if (t <= 0.52)
565  {
566  return -1.e1 * 12592;
567  }
568  if (t <= 0.53)
569  {
570  return -1.e1 * 12548;
571  }
572  if (t <= 0.54)
573  {
574  return -1.e1 * 12504;
575  }
576  if (t <= 0.55)
577  {
578  return -1.e1 * 12456;
579  }
580  if (t <= 0.56)
581  {
582  return -1.e1 * 12405;
583  }
584  if (t <= 0.57)
585  {
586  return -1.e1 * 12353;
587  }
588  if (t <= 0.58)
589  {
590  return -1.e1 * 12300;
591  }
592  if (t <= 0.59)
593  {
594  return -1.e1 * 12244;
595  }
596  if (t <= 0.60)
597  {
598  return -1.e1 * 12184;
599  }
600  if (t <= 0.61)
601  {
602  return -1.e1 * 12122;
603  }
604  if (t <= 0.62)
605  {
606  return -1.e1 * 12058;
607  }
608  if (t <= 0.63)
609  {
610  return -1.e1 * 11995;
611  }
612  if (t <= 0.64)
613  {
614  return -1.e1 * 11933;
615  }
616  if (t <= 0.65)
617  {
618  return -1.e1 * 11871;
619  }
620  if (t <= 0.66)
621  {
622  return -1.e1 * 11810;
623  }
624  if (t <= 0.67)
625  {
626  return -1.e1 * 11747;
627  }
628  if (t <= 0.68)
629  {
630  return -1.e1 * 11684;
631  }
632  if (t <= 0.69)
633  {
634  return -1.e1 * 11620;
635  }
636  if (t <= 0.70)
637  {
638  return -1.e1 * 11556;
639  }
640  if (t <= 0.71)
641  {
642  return -1.e1 * 11492;
643  }
644  if (t <= 0.72)
645  {
646  return -1.e1 * 11428;
647  }
648  if (t <= 0.73)
649  {
650  return -1.e1 * 11365;
651  }
652  if (t <= 0.74)
653  {
654  return -1.e1 * 11302;
655  }
656  if (t <= 0.75)
657  {
658  return -1.e1 * 11240;
659  }
660  if (t <= 0.76)
661  {
662  return -1.e1 * 11179;
663  }
664  if (t <= 0.77)
665  {
666  return -1.e1 * 11120;
667  }
668  if (t <= 0.78)
669  {
670  return -1.e1 * 11062;
671  }
672  if (t <= 0.79)
673  {
674  return -1.e1 * 11006;
675  }
676  // break;
677  /* case 3:
678  return 0.0;
679  break;}
680  return 0.;*/
681 }
682 
683 
684 Real aortaFlux3 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
685 {
686  if (t <= 0.00)
687  {
688  return 3.350E-05 * 1e6;
689  }
690  if (t <= 0.01)
691  {
692  return 3.373E-05 * 1e6;
693  }
694  if (t <= 0.02)
695  {
696  return 3.402E-05 * 1e6;
697  }
698  if (t <= 0.03)
699  {
700  return 3.434E-05 * 1e6;
701  }
702  if (t <= 0.04)
703  {
704  return 3.466E-05 * 1e6;
705  }
706  if (t <= 0.05)
707  {
708  return 3.495E-05 * 1e6;
709  }
710  if (t <= 0.06)
711  {
712  return 3.519E-05 * 1e6;
713  }
714  if (t <= 0.07)
715  {
716  return 3.539E-05 * 1e6;
717  }
718  if (t <= 0.08)
719  {
720  return 3.564E-05 * 1e6;
721  }
722  if (t <= 0.09)
723  {
724  return 3.617E-05 * 1e6;
725  }
726  if (t <= 0.10)
727  {
728  return 3.773E-05 * 1e6;
729  }
730  if (t <= 0.11)
731  {
732  return 4.176E-05 * 1e6;
733  }
734  if (t <= 0.12)
735  {
736  return 5.037E-05 * 1e6;
737  }
738  if (t <= 0.13)
739  {
740  return 6.546E-05 * 1e6;
741  }
742  if (t <= 0.14)
743  {
744  return 8.701E-05 * 1e6;
745  }
746  if (t <= 0.15)
747  {
748  return 1.117E-04 * 1e6;
749  }
750  if (t <= 0.16)
751  {
752  return 1.345E-04 * 1e6;
753  }
754  if (t <= 0.17)
755  {
756  return 1.519E-04 * 1e6;
757  }
758  if (t <= 0.18)
759  {
760  return 1.642E-04 * 1e6;
761  }
762  if (t <= 0.19)
763  {
764  return 1.737E-04 * 1e6;
765  }
766  if (t <= 0.20)
767  {
768  return 1.821E-04 * 1e6;
769  }
770  if (t <= 0.21)
771  {
772  return 1.897E-04 * 1e6;
773  }
774  if (t <= 0.22)
775  {
776  return 1.958E-04 * 1e6;
777  }
778  if (t <= 0.23)
779  {
780  return 1.999E-04 * 1e6;
781  }
782  if (t <= 0.24)
783  {
784  return 2.019E-04 * 1e6;
785  }
786  if (t <= 0.25)
787  {
788  return 2.020E-04 * 1e6;
789  }
790  if (t <= 0.26)
791  {
792  return 2.004E-04 * 1e6;
793  }
794  if (t <= 0.27)
795  {
796  return 1.972E-04 * 1e6;
797  }
798  if (t <= 0.28)
799  {
800  return 1.926E-04 * 1e6;
801  }
802  if (t <= 0.29)
803  {
804  return 1.868E-04 * 1e6;
805  }
806  if (t <= 0.30)
807  {
808  return 1.798E-04 * 1e6;
809  }
810  if (t <= 0.31)
811  {
812  return 1.719E-04 * 1e6;
813  }
814  if (t <= 0.32)
815  {
816  return 1.632E-04 * 1e6;
817  }
818  if (t <= 0.33)
819  {
820  return 1.540E-04 * 1e6;
821  }
822  if (t <= 0.34)
823  {
824  return 1.446E-04 * 1e6;
825  }
826  if (t <= 0.35)
827  {
828  return 1.350E-04 * 1e6;
829  }
830  if (t <= 0.36)
831  {
832  return 1.254E-04 * 1e6;
833  }
834  if (t <= 0.37)
835  {
836  return 1.158E-04 * 1e6;
837  }
838  if (t <= 0.38)
839  {
840  return 1.062E-04 * 1e6;
841  }
842  if (t <= 0.39)
843  {
844  return 9.651E-05 * 1e6;
845  }
846  if (t <= 0.40)
847  {
848  return 8.634E-05 * 1e6;
849  }
850  if (t <= 0.41)
851  {
852  return 7.558E-05 * 1e6;
853  }
854  if (t <= 0.42)
855  {
856  return 6.447E-05 * 1e6;
857  }
858  if (t <= 0.43)
859  {
860  return 5.382E-05 * 1e6;
861  }
862  if (t <= 0.44)
863  {
864  return 4.484E-05 * 1e6;
865  }
866  if (t <= 0.45)
867  {
868  return 3.865E-05 * 1e6;
869  }
870  if (t <= 0.46)
871  {
872  return 3.556E-05 * 1e6;
873  }
874  if (t <= 0.47)
875  {
876  return 3.473E-05 * 1e6;
877  }
878  if (t <= 0.48)
879  {
880  return 3.457E-05 * 1e6;
881  }
882  if (t <= 0.49)
883  {
884  return 3.373E-05 * 1e6;
885  }
886  if (t <= 0.50)
887  {
888  return 3.191E-05 * 1e6;
889  }
890  if (t <= 0.51)
891  {
892  return 2.975E-05 * 1e6;
893  }
894  if (t <= 0.52)
895  {
896  return 2.809E-05 * 1e6;
897  }
898  if (t <= 0.53)
899  {
900  return 2.730E-05 * 1e6;
901  }
902  if (t <= 0.54)
903  {
904  return 2.718E-05 * 1e6;
905  }
906  if (t <= 0.55)
907  {
908  return 2.732E-05 * 1e6;
909  }
910  if (t <= 0.56)
911  {
912  return 2.744E-05 * 1e6;
913  }
914  if (t <= 0.57)
915  {
916  return 2.753E-05 * 1e6;
917  }
918  if (t <= 0.58)
919  {
920  return 2.772E-05 * 1e6;
921  }
922  if (t <= 0.59)
923  {
924  return 2.811E-05 * 1e6;
925  }
926  if (t <= 0.60)
927  {
928  return 2.866E-05 * 1e6;
929  }
930  if (t <= 0.61)
931  {
932  return 2.929E-05 * 1e6;
933  }
934  if (t <= 0.62)
935  {
936  return 2.990E-05 * 1e6;
937  }
938  if (t <= 0.63)
939  {
940  return 3.044E-05 * 1e6;
941  }
942  if (t <= 0.64)
943  {
944  return 3.091E-05 * 1e6;
945  }
946  if (t <= 0.65)
947  {
948  return 3.132E-05 * 1e6;
949  }
950  if (t <= 0.66)
951  {
952  return 3.168E-05 * 1e6;
953  }
954  if (t <= 0.67)
955  {
956  return 3.199E-05 * 1e6;
957  }
958  if (t <= 0.68)
959  {
960  return 3.224E-05 * 1e6;
961  }
962  if (t <= 0.69)
963  {
964  return 3.244E-05 * 1e6;
965  }
966  if (t <= 0.70)
967  {
968  return 3.259E-05 * 1e6;
969  }
970  if (t <= 0.71)
971  {
972  return 3.270E-05 * 1e6;
973  }
974  if (t <= 0.72)
975  {
976  return 3.277E-05 * 1e6;
977  }
978  if (t <= 0.73)
979  {
980  return 3.281E-05 * 1e6;
981  }
982  if (t <= 0.74)
983  {
984  return 3.282E-05 * 1e6;
985  }
986  if (t <= 0.75)
987  {
988  return 3.283E-05 * 1e6;
989  }
990  if (t <= 0.76)
991  {
992  return 3.283E-05 * 1e6;
993  }
994  if (t <= 0.77)
995  {
996  return 3.286E-05 * 1e6;
997  }
998  if (t <= 0.78)
999  {
1000  return 3.291E-05 * 1e6;
1001  }
1002  if (t <= 0.79)
1003  {
1004  return 3.300E-05 * 1e6;
1005  }
1006 }//thoracic aorta,
1007 
1008 Real aortaFlux5 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
1009 {
1010  if (t <= 0.00)
1011  {
1012  return 3.033E-06 * 1e6;
1013  }
1014  if (t <= 0.01)
1015  {
1016  return 3.041E-06 * 1e6;
1017  }
1018  if (t <= 0.02)
1019  {
1020  return 3.062E-06 * 1e6;
1021  }
1022  if (t <= 0.03)
1023  {
1024  return 3.094E-06 * 1e6;
1025  }
1026  if (t <= 0.04)
1027  {
1028  return 3.127E-06 * 1e6;
1029  }
1030  if (t <= 0.05)
1031  {
1032  return 3.150E-06 * 1e6;
1033  }
1034  if (t <= 0.06)
1035  {
1036  return 3.152E-06 * 1e6;
1037  }
1038  if (t <= 0.07)
1039  {
1040  return 3.141E-06 * 1e6;
1041  }
1042  if (t <= 0.08)
1043  {
1044  return 3.196E-06 * 1e6;
1045  }
1046  if (t <= 0.09)
1047  {
1048  return 3.574E-06 * 1e6;
1049  }
1050  if (t <= 0.10)
1051  {
1052  return 4.778E-06 * 1e6;
1053  }
1054  if (t <= 0.11)
1055  {
1056  return 7.387E-06 * 1e6;
1057  }
1058  if (t <= 0.12)
1059  {
1060  return 1.150E-05 * 1e6;
1061  }
1062  if (t <= 0.13)
1063  {
1064  return 1.609E-05 * 1e6;
1065  }
1066  if (t <= 0.14)
1067  {
1068  return 1.933E-05 * 1e6;
1069  }
1070  if (t <= 0.15)
1071  {
1072  return 2.007E-05 * 1e6;
1073  }
1074  if (t <= 0.16)
1075  {
1076  return 1.885E-05 * 1e6;
1077  }
1078  if (t <= 0.17)
1079  {
1080  return 1.706E-05 * 1e6;
1081  }
1082  if (t <= 0.18)
1083  {
1084  return 1.569E-05 * 1e6;
1085  }
1086  if (t <= 0.19)
1087  {
1088  return 1.481E-05 * 1e6;
1089  }
1090  if (t <= 0.20)
1091  {
1092  return 1.401E-05 * 1e6;
1093  }
1094  if (t <= 0.21)
1095  {
1096  return 1.294E-05 * 1e6;
1097  }
1098  if (t <= 0.22)
1099  {
1100  return 1.160E-05 * 1e6;
1101  }
1102  if (t <= 0.23)
1103  {
1104  return 1.018E-05 * 1e6;
1105  }
1106  if (t <= 0.24)
1107  {
1108  return 8.832E-06 * 1e6;
1109  }
1110  if (t <= 0.25)
1111  {
1112  return 7.609E-06 * 1e6;
1113  }
1114  if (t <= 0.26)
1115  {
1116  return 6.578E-06 * 1e6;
1117  }
1118  if (t <= 0.27)
1119  {
1120  return 5.843E-06 * 1e6;
1121  }
1122  if (t <= 0.28)
1123  {
1124  return 5.472E-06 * 1e6;
1125  }
1126  if (t <= 0.29)
1127  {
1128  return 5.412E-06 * 1e6;
1129  }
1130  if (t <= 0.30)
1131  {
1132  return 5.491E-06 * 1e6;
1133  }
1134  if (t <= 0.31)
1135  {
1136  return 5.527E-06 * 1e6;
1137  }
1138  if (t <= 0.32)
1139  {
1140  return 5.420E-06 * 1e6;
1141  }
1142  if (t <= 0.33)
1143  {
1144  return 5.169E-06 * 1e6;
1145  }
1146  if (t <= 0.34)
1147  {
1148  return 4.829E-06 * 1e6;
1149  }
1150  if (t <= 0.35)
1151  {
1152  return 4.465E-06 * 1e6;
1153  }
1154  if (t <= 0.36)
1155  {
1156  return 4.111E-06 * 1e6;
1157  }
1158  if (t <= 0.37)
1159  {
1160  return 3.750E-06 * 1e6;
1161  }
1162  if (t <= 0.38)
1163  {
1164  return 3.304E-06 * 1e6;
1165  }
1166  if (t <= 0.39)
1167  {
1168  return 2.668E-06 * 1e6;
1169  }
1170  if (t <= 0.40)
1171  {
1172  return 1.800E-06 * 1e6;
1173  }
1174  if (t <= 0.41)
1175  {
1176  return 8.269E-07 * 1e6;
1177  }
1178  if (t <= 0.42)
1179  {
1180  return 9.760E-08 * 1e6;
1181  }
1182  if (t <= 0.43)
1183  {
1184  return 7.311E-08 * 1e6;
1185  }
1186  if (t <= 0.44)
1187  {
1188  return 1.041E-06 * 1e6;
1189  }
1190  if (t <= 0.45)
1191  {
1192  return 2.783E-06 * 1e6;
1193  }
1194  if (t <= 0.46)
1195  {
1196  return 4.537E-06 * 1e6;
1197  }
1198  if (t <= 0.47)
1199  {
1200  return 5.488E-06 * 1e6;
1201  }
1202  if (t <= 0.48)
1203  {
1204  return 5.431E-06 * 1e6;
1205  }
1206  if (t <= 0.49)
1207  {
1208  return 4.863E-06 * 1e6;
1209  }
1210  if (t <= 0.50)
1211  {
1212  return 4.452E-06 * 1e6;
1213  }
1214  if (t <= 0.51)
1215  {
1216  return 4.499E-06 * 1e6;
1217  }
1218  if (t <= 0.52)
1219  {
1220  return 4.824E-06 * 1e6;
1221  }
1222  if (t <= 0.53)
1223  {
1224  return 5.059E-06 * 1e6;
1225  }
1226  if (t <= 0.54)
1227  {
1228  return 4.989E-06 * 1e6;
1229  }
1230  if (t <= 0.55)
1231  {
1232  return 4.671E-06 * 1e6;
1233  }
1234  if (t <= 0.56)
1235  {
1236  return 4.292E-06 * 1e6;
1237  }
1238  if (t <= 0.57)
1239  {
1240  return 3.981E-06 * 1e6;
1241  }
1242  if (t <= 0.58)
1243  {
1244  return 3.749E-06 * 1e6;
1245  }
1246  if (t <= 0.59)
1247  {
1248  return 3.553E-06 * 1e6;
1249  }
1250  if (t <= 0.60)
1251  {
1252  return 3.377E-06 * 1e6;
1253  }
1254  if (t <= 0.61)
1255  {
1256  return 3.255E-06 * 1e6;
1257  }
1258  if (t <= 0.62)
1259  {
1260  return 3.224E-06 * 1e6;
1261  }
1262  if (t <= 0.63)
1263  {
1264  return 3.281E-06 * 1e6;
1265  }
1266  if (t <= 0.64)
1267  {
1268  return 3.377E-06 * 1e6;
1269  }
1270  if (t <= 0.65)
1271  {
1272  return 3.452E-06 * 1e6;
1273  }
1274  if (t <= 0.66)
1275  {
1276  return 3.472E-06 * 1e6;
1277  }
1278  if (t <= 0.67)
1279  {
1280  return 3.441E-06 * 1e6;
1281  }
1282  if (t <= 0.68)
1283  {
1284  return 3.389E-06 * 1e6;
1285  }
1286  if (t <= 0.69)
1287  {
1288  return 3.343E-06 * 1e6;
1289  }
1290  if (t <= 0.70)
1291  {
1292  return 3.312E-06 * 1e6;
1293  }
1294  if (t <= 0.71)
1295  {
1296  return 3.289E-06 * 1e6;
1297  }
1298  if (t <= 0.72)
1299  {
1300  return 3.262E-06 * 1e6;
1301  }
1302  if (t <= 0.73)
1303  {
1304  return 3.223E-06 * 1e6;
1305  }
1306  if (t <= 0.74)
1307  {
1308  return 3.177E-06 * 1e6;
1309  }
1310  if (t <= 0.75)
1311  {
1312  return 3.132E-06 * 1e6;
1313  }
1314  if (t <= 0.76)
1315  {
1316  return 3.094E-06 * 1e6;
1317  }
1318  if (t <= 0.77)
1319  {
1320  return 3.065E-06 * 1e6;
1321  }
1322  if (t <= 0.78)
1323  {
1324  return 3.040E-06 * 1e6;
1325  }
1326  if (t <= 0.79)
1327  {
1328  return 3.016E-06 * 1e6;
1329  }
1330 }//first branchstd::placeholders::_1,
1331 
1332 Real aortaFlux6 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
1333 {
1334  if (t <= 0.00)
1335  {
1336  return 7.817E-07 * 1e6;
1337  }
1338  if (t <= 0.01)
1339  {
1340  return 7.879E-07 * 1e6;
1341  }
1342  if (t <= 0.02)
1343  {
1344  return 7.977E-07 * 1e6;
1345  }
1346  if (t <= 0.03)
1347  {
1348  return 8.077E-07 * 1e6;
1349  }
1350  if (t <= 0.04)
1351  {
1352  return 8.144E-07 * 1e6;
1353  }
1354  if (t <= 0.05)
1355  {
1356  return 8.164E-07 * 1e6;
1357  }
1358  if (t <= 0.06)
1359  {
1360  return 8.136E-07 * 1e6;
1361  }
1362  if (t <= 0.07)
1363  {
1364  return 8.129E-07 * 1e6;
1365  }
1366  if (t <= 0.08)
1367  {
1368  return 8.462E-07 * 1e6;
1369  }
1370  if (t <= 0.09)
1371  {
1372  return 9.858E-07 * 1e6;
1373  }
1374  if (t <= 0.10)
1375  {
1376  return 1.317E-06 * 1e6;
1377  }
1378  if (t <= 0.11)
1379  {
1380  return 1.854E-06 * 1e6;
1381  }
1382  if (t <= 0.12)
1383  {
1384  return 2.476E-06 * 1e6;
1385  }
1386  if (t <= 0.13)
1387  {
1388  return 3.005E-06 * 1e6;
1389  }
1390  if (t <= 0.14)
1391  {
1392  return 3.366E-06 * 1e6;
1393  }
1394  if (t <= 0.15)
1395  {
1396  return 3.626E-06 * 1e6;
1397  }
1398  if (t <= 0.16)
1399  {
1400  return 3.860E-06 * 1e6;
1401  }
1402  if (t <= 0.17)
1403  {
1404  return 4.044E-06 * 1e6;
1405  }
1406  if (t <= 0.18)
1407  {
1408  return 4.065E-06 * 1e6;
1409  }
1410  if (t <= 0.19)
1411  {
1412  return 3.824E-06 * 1e6;
1413  }
1414  if (t <= 0.20)
1415  {
1416  return 3.320E-06 * 1e6;
1417  }
1418  if (t <= 0.21)
1419  {
1420  return 2.659E-06 * 1e6;
1421  }
1422  if (t <= 0.22)
1423  {
1424  return 2.006E-06 * 1e6;
1425  }
1426  if (t <= 0.23)
1427  {
1428  return 1.503E-06 * 1e6;
1429  }
1430  if (t <= 0.24)
1431  {
1432  return 1.215E-06 * 1e6;
1433  }
1434  if (t <= 0.25)
1435  {
1436  return 1.117E-06 * 1e6;
1437  }
1438  if (t <= 0.26)
1439  {
1440  return 1.136E-06 * 1e6;
1441  }
1442  if (t <= 0.27)
1443  {
1444  return 1.190E-06 * 1e6;
1445  }
1446  if (t <= 0.28)
1447  {
1448  return 1.222E-06 * 1e6;
1449  }
1450  if (t <= 0.29)
1451  {
1452  return 1.219E-06 * 1e6;
1453  }
1454  if (t <= 0.30)
1455  {
1456  return 1.197E-06 * 1e6;
1457  }
1458  if (t <= 0.31)
1459  {
1460  return 1.179E-06 * 1e6;
1461  }
1462  if (t <= 0.32)
1463  {
1464  return 1.175E-06 * 1e6;
1465  }
1466  if (t <= 0.33)
1467  {
1468  return 1.176E-06 * 1e6;
1469  }
1470  if (t <= 0.34)
1471  {
1472  return 1.164E-06 * 1e6;
1473  }
1474  if (t <= 0.35)
1475  {
1476  return 1.129E-06 * 1e6;
1477  }
1478  if (t <= 0.36)
1479  {
1480  return 1.063E-06 * 1e6;
1481  }
1482  if (t <= 0.37)
1483  {
1484  return 9.647E-07 * 1e6;
1485  }
1486  if (t <= 0.38)
1487  {
1488  return 8.310E-07 * 1e6;
1489  }
1490  if (t <= 0.39)
1491  {
1492  return 6.635E-07 * 1e6;
1493  }
1494  if (t <= 0.40)
1495  {
1496  return 4.774E-07 * 1e6;
1497  }
1498  if (t <= 0.41)
1499  {
1500  return 3.116E-07 * 1e6;
1501  }
1502  if (t <= 0.42)
1503  {
1504  return 2.251E-07 * 1e6;
1505  }
1506  if (t <= 0.43)
1507  {
1508  return 2.627E-07 * 1e6;
1509  }
1510  if (t <= 0.44)
1511  {
1512  return 4.099E-07 * 1e6;
1513  }
1514  if (t <= 0.45)
1515  {
1516  return 5.913E-07 * 1e6;
1517  }
1518  if (t <= 0.46)
1519  {
1520  return 7.359E-07 * 1e6;
1521  }
1522  if (t <= 0.47)
1523  {
1524  return 8.403E-07 * 1e6;
1525  }
1526  if (t <= 0.48)
1527  {
1528  return 9.515E-07 * 1e6;
1529  }
1530  if (t <= 0.49)
1531  {
1532  return 1.097E-06 * 1e6;
1533  }
1534  if (t <= 0.50)
1535  {
1536  return 1.245E-06 * 1e6;
1537  }
1538  if (t <= 0.51)
1539  {
1540  return 1.331E-06 * 1e6;
1541  }
1542  if (t <= 0.52)
1543  {
1544  return 1.314E-06 * 1e6;
1545  }
1546  if (t <= 0.53)
1547  {
1548  return 1.204E-06 * 1e6;
1549  }
1550  if (t <= 0.54)
1551  {
1552  return 1.048E-06 * 1e6;
1553  }
1554  if (t <= 0.55)
1555  {
1556  return 9.004E-07 * 1e6;
1557  }
1558  if (t <= 0.56)
1559  {
1560  return 7.937E-07 * 1e6;
1561  }
1562  if (t <= 0.57)
1563  {
1564  return 7.358E-07 * 1e6;
1565  }
1566  if (t <= 0.58)
1567  {
1568  return 7.163E-07 * 1e6;
1569  }
1570  if (t <= 0.59)
1571  {
1572  return 7.186E-07 * 1e6;
1573  }
1574  if (t <= 0.60)
1575  {
1576  return 7.280E-07 * 1e6;
1577  }
1578  if (t <= 0.61)
1579  {
1580  return 7.359E-07 * 1e6;
1581  }
1582  if (t <= 0.62)
1583  {
1584  return 7.400E-07 * 1e6;
1585  }
1586  if (t <= 0.63)
1587  {
1588  return 7.423E-07 * 1e6;
1589  }
1590  if (t <= 0.64)
1591  {
1592  return 7.465E-07 * 1e6;
1593  }
1594  if (t <= 0.65)
1595  {
1596  return 7.556E-07 * 1e6;
1597  }
1598  if (t <= 0.66)
1599  {
1600  return 7.700E-07 * 1e6;
1601  }
1602  if (t <= 0.67)
1603  {
1604  return 7.871E-07 * 1e6;
1605  }
1606  if (t <= 0.68)
1607  {
1608  return 8.028E-07 * 1e6;
1609  }
1610  if (t <= 0.69)
1611  {
1612  return 8.129E-07 * 1e6;
1613  }
1614  if (t <= 0.70)
1615  {
1616  return 8.151E-07 * 1e6;
1617  }
1618  if (t <= 0.71)
1619  {
1620  return 8.101E-07 * 1e6;
1621  }
1622  if (t <= 0.72)
1623  {
1624  return 8.006E-07 * 1e6;
1625  }
1626  if (t <= 0.73)
1627  {
1628  return 7.905E-07 * 1e6;
1629  }
1630  if (t <= 0.74)
1631  {
1632  return 7.830E-07 * 1e6;
1633  }
1634  if (t <= 0.75)
1635  {
1636  return 7.792E-07 * 1e6;
1637  }
1638  if (t <= 0.76)
1639  {
1640  return 7.782E-07 * 1e6;
1641  }
1642  if (t <= 0.77)
1643  {
1644  return 7.782E-07 * 1e6;
1645  }
1646  if (t <= 0.78)
1647  {
1648  return 7.775E-07 * 1e6;
1649  }
1650  if (t <= 0.79)
1651  {
1652  return 7.754E-07 * 1e6;
1653  }
1654 }//branch 1_2 smallest
1655 
1656 Real aortaFlux7 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
1657 {
1658  if (t <= 0.00)
1659  {
1660  return 1.930E-06 * 1e6;
1661  }
1662  if (t <= 0.01)
1663  {
1664  return 1.710E-06 * 1e6;
1665  }
1666  if (t <= 0.02)
1667  {
1668  return 1.495E-06 * 1e6;
1669  }
1670  if (t <= 0.03)
1671  {
1672  return 1.289E-06 * 1e6;
1673  }
1674  if (t <= 0.04)
1675  {
1676  return 1.096E-06 * 1e6;
1677  }
1678  if (t <= 0.05)
1679  {
1680  return 9.184E-07 * 1e6;
1681  }
1682  if (t <= 0.06)
1683  {
1684  return 7.568E-07 * 1e6;
1685  }
1686  if (t <= 0.07)
1687  {
1688  return 6.240E-07 * 1e6;
1689  }
1690  if (t <= 0.08)
1691  {
1692  return 5.854E-07 * 1e6;
1693  }
1694  if (t <= 0.09)
1695  {
1696  return 8.113E-07 * 1e6;
1697  }
1698  if (t <= 0.10)
1699  {
1700  return 1.559E-06 * 1e6;
1701  }
1702  if (t <= 0.11)
1703  {
1704  return 3.032E-06 * 1e6;
1705  }
1706  if (t <= 0.12)
1707  {
1708  return 5.221E-06 * 1e6;
1709  }
1710  if (t <= 0.13)
1711  {
1712  return 7.903E-06 * 1e6;
1713  }
1714  if (t <= 0.14)
1715  {
1716  return 1.081E-05 * 1e6;
1717  }
1718  if (t <= 0.15)
1719  {
1720  return 1.376E-05 * 1e6;
1721  }
1722  if (t <= 0.16)
1723  {
1724  return 1.661E-05 * 1e6;
1725  }
1726  if (t <= 0.17)
1727  {
1728  return 1.919E-05 * 1e6;
1729  }
1730  if (t <= 0.18)
1731  {
1732  return 2.134E-05 * 1e6;
1733  }
1734  if (t <= 0.19)
1735  {
1736  return 2.292E-05 * 1e6;
1737  }
1738  if (t <= 0.20)
1739  {
1740  return 2.389E-05 * 1e6;
1741  }
1742  if (t <= 0.21)
1743  {
1744  return 2.433E-05 * 1e6;
1745  }
1746  if (t <= 0.22)
1747  {
1748  return 2.435E-05 * 1e6;
1749  }
1750  if (t <= 0.23)
1751  {
1752  return 2.406E-05 * 1e6;
1753  }
1754  if (t <= 0.24)
1755  {
1756  return 2.353E-05 * 1e6;
1757  }
1758  if (t <= 0.25)
1759  {
1760  return 2.280E-05 * 1e6;
1761  }
1762  if (t <= 0.26)
1763  {
1764  return 2.188E-05 * 1e6;
1765  }
1766  if (t <= 0.27)
1767  {
1768  return 2.078E-05 * 1e6;
1769  }
1770  if (t <= 0.28)
1771  {
1772  return 1.947E-05 * 1e6;
1773  }
1774  if (t <= 0.29)
1775  {
1776  return 1.797E-05 * 1e6;
1777  }
1778  if (t <= 0.30)
1779  {
1780  return 1.628E-05 * 1e6;
1781  }
1782  if (t <= 0.31)
1783  {
1784  return 1.445E-05 * 1e6;
1785  }
1786  if (t <= 0.32)
1787  {
1788  return 1.254E-05 * 1e6;
1789  }
1790  if (t <= 0.33)
1791  {
1792  return 1.060E-05 * 1e6;
1793  }
1794  if (t <= 0.34)
1795  {
1796  return 8.684E-06 * 1e6;
1797  }
1798  if (t <= 0.35)
1799  {
1800  return 6.838E-06 * 1e6;
1801  }
1802  if (t <= 0.36)
1803  {
1804  return 5.084E-06 * 1e6;
1805  }
1806  if (t <= 0.37)
1807  {
1808  return 3.412E-06 * 1e6;
1809  }
1810  if (t <= 0.38)
1811  {
1812  return 1.784E-06 * 1e6;
1813  }
1814  if (t <= 0.39)
1815  {
1816  return 1.534E-07 * 1e6;
1817  }
1818  if (t <= 0.40)
1819  {
1820  return -1.494E-06 * 1e6;
1821  }
1822  if (t <= 0.41)
1823  {
1824  return -3.093E-06 * 1e6;
1825  }
1826  if (t <= 0.42)
1827  {
1828  return -4.495E-06 * 1e6;
1829  }
1830  if (t <= 0.43)
1831  {
1832  return -5.521E-06 * 1e6;
1833  }
1834  if (t <= 0.44)
1835  {
1836  return -6.086E-06 * 1e6;
1837  }
1838  if (t <= 0.45)
1839  {
1840  return -6.252E-06 * 1e6;
1841  }
1842  if (t <= 0.46)
1843  {
1844  return -6.160E-06 * 1e6;
1845  }
1846  if (t <= 0.47)
1847  {
1848  return -5.908E-06 * 1e6;
1849  }
1850  if (t <= 0.48)
1851  {
1852  return -5.508E-06 * 1e6;
1853  }
1854  if (t <= 0.49)
1855  {
1856  return -4.944E-06 * 1e6;
1857  }
1858  if (t <= 0.50)
1859  {
1860  return -4.234E-06 * 1e6;
1861  }
1862  if (t <= 0.51)
1863  {
1864  return -3.441E-06 * 1e6;
1865  }
1866  if (t <= 0.52)
1867  {
1868  return -2.639E-06 * 1e6;
1869  }
1870  if (t <= 0.53)
1871  {
1872  return -1.873E-06 * 1e6;
1873  }
1874  if (t <= 0.54)
1875  {
1876  return -1.160E-06 * 1e6;
1877  }
1878  if (t <= 0.55)
1879  {
1880  return -5.019E-07 * 1e6;
1881  }
1882  if (t <= 0.56)
1883  {
1884  return 1.024E-07 * 1e6;
1885  }
1886  if (t <= 0.57)
1887  {
1888  return 6.543E-07 * 1e6;
1889  }
1890  if (t <= 0.58)
1891  {
1892  return 1.156E-06 * 1e6;
1893  }
1894  if (t <= 0.59)
1895  {
1896  return 1.614E-06 * 1e6;
1897  }
1898  if (t <= 0.60)
1899  {
1900  return 2.034E-06 * 1e6;
1901  }
1902  if (t <= 0.61)
1903  {
1904  return 2.419E-06 * 1e6;
1905  }
1906  if (t <= 0.62)
1907  {
1908  return 2.765E-06 * 1e6;
1909  }
1910  if (t <= 0.63)
1911  {
1912  return 3.065E-06 * 1e6;
1913  }
1914  if (t <= 0.64)
1915  {
1916  return 3.307E-06 * 1e6;
1917  }
1918  if (t <= 0.65)
1919  {
1920  return 3.486E-06 * 1e6;
1921  }
1922  if (t <= 0.66)
1923  {
1924  return 3.600E-06 * 1e6;
1925  }
1926  if (t <= 0.67)
1927  {
1928  return 3.654E-06 * 1e6;
1929  }
1930  if (t <= 0.68)
1931  {
1932  return 3.657E-06 * 1e6;
1933  }
1934  if (t <= 0.69)
1935  {
1936  return 3.621E-06 * 1e6;
1937  }
1938  if (t <= 0.70)
1939  {
1940  return 3.554E-06 * 1e6;
1941  }
1942  if (t <= 0.71)
1943  {
1944  return 3.465E-06 * 1e6;
1945  }
1946  if (t <= 0.72)
1947  {
1948  return 3.357E-06 * 1e6;
1949  }
1950  if (t <= 0.73)
1951  {
1952  return 3.233E-06 * 1e6;
1953  }
1954  if (t <= 0.74)
1955  {
1956  return 3.094E-06 * 1e6;
1957  }
1958  if (t <= 0.75)
1959  {
1960  return 2.941E-06 * 1e6;
1961  }
1962  if (t <= 0.76)
1963  {
1964  return 2.774E-06 * 1e6;
1965  }
1966  if (t <= 0.77)
1967  {
1968  return 2.591E-06 * 1e6;
1969  }
1970  if (t <= 0.78)
1971  {
1972  return 2.395E-06 * 1e6;
1973  }
1974  if (t <= 0.79)
1975  {
1976  return 2.185E-06 * 1e6;
1977  }
1978 }//R. Brachia, branch 1_3
1979 
1980 Real aortaFlux8 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
1981 {
1982  if (t <= 0.00)
1983  {
1984  return 2.445E-06 * 1e6;
1985  }
1986  if (t <= 0.01)
1987  {
1988  return 2.470E-06 * 1e6;
1989  }
1990  if (t <= 0.02)
1991  {
1992  return 2.519E-06 * 1e6;
1993  }
1994  if (t <= 0.03)
1995  {
1996  return 2.578E-06 * 1e6;
1997  }
1998  if (t <= 0.04)
1999  {
2000  return 2.631E-06 * 1e6;
2001  }
2002  if (t <= 0.05)
2003  {
2004  return 2.665E-06 * 1e6;
2005  }
2006  if (t <= 0.06)
2007  {
2008  return 2.674E-06 * 1e6;
2009  }
2010  if (t <= 0.07)
2011  {
2012  return 2.674E-06 * 1e6;
2013  }
2014  if (t <= 0.08)
2015  {
2016  return 2.777E-06 * 1e6;
2017  }
2018  if (t <= 0.09)
2019  {
2020  return 3.260E-06 * 1e6;
2021  }
2022  if (t <= 0.10)
2023  {
2024  return 4.539E-06 * 1e6;
2025  }
2026  if (t <= 0.11)
2027  {
2028  return 6.956E-06 * 1e6;
2029  }
2030  if (t <= 0.12)
2031  {
2032  return 1.045E-05 * 1e6;
2033  }
2034  if (t <= 0.13)
2035  {
2036  return 1.438E-05 * 1e6;
2037  }
2038  if (t <= 0.14)
2039  {
2040  return 1.769E-05 * 1e6;
2041  }
2042  if (t <= 0.15)
2043  {
2044  return 1.958E-05 * 1e6;
2045  }
2046  if (t <= 0.16)
2047  {
2048  return 1.991E-05 * 1e6;
2049  }
2050  if (t <= 0.17)
2051  {
2052  return 1.910E-05 * 1e6;
2053  }
2054  if (t <= 0.18)
2055  {
2056  return 1.775E-05 * 1e6;
2057  }
2058  if (t <= 0.19)
2059  {
2060  return 1.628E-05 * 1e6;
2061  }
2062  if (t <= 0.20)
2063  {
2064  return 1.479E-05 * 1e6;
2065  }
2066  if (t <= 0.21)
2067  {
2068  return 1.326E-05 * 1e6;
2069  }
2070  if (t <= 0.22)
2071  {
2072  return 1.165E-05 * 1e6;
2073  }
2074  if (t <= 0.23)
2075  {
2076  return 1.001E-05 * 1e6;
2077  }
2078  if (t <= 0.24)
2079  {
2080  return 8.448E-06 * 1e6;
2081  }
2082  if (t <= 0.25)
2083  {
2084  return 7.065E-06 * 1e6;
2085  }
2086  if (t <= 0.26)
2087  {
2088  return 5.930E-06 * 1e6;
2089  }
2090  if (t <= 0.27)
2091  {
2092  return 5.089E-06 * 1e6;
2093  }
2094  if (t <= 0.28)
2095  {
2096  return 4.566E-06 * 1e6;
2097  }
2098  if (t <= 0.29)
2099  {
2100  return 4.350E-06 * 1e6;
2101  }
2102  if (t <= 0.30)
2103  {
2104  return 4.382E-06 * 1e6;
2105  }
2106  if (t <= 0.31)
2107  {
2108  return 4.552E-06 * 1e6;
2109  }
2110  if (t <= 0.32)
2111  {
2112  return 4.729E-06 * 1e6;
2113  }
2114  if (t <= 0.33)
2115  {
2116  return 4.798E-06 * 1e6;
2117  }
2118  if (t <= 0.34)
2119  {
2120  return 4.698E-06 * 1e6;
2121  }
2122  if (t <= 0.35)
2123  {
2124  return 4.428E-06 * 1e6;
2125  }
2126  if (t <= 0.36)
2127  {
2128  return 4.026E-06 * 1e6;
2129  }
2130  if (t <= 0.37)
2131  {
2132  return 3.513E-06 * 1e6;
2133  }
2134  if (t <= 0.38)
2135  {
2136  return 2.875E-06 * 1e6;
2137  }
2138  if (t <= 0.39)
2139  {
2140  return 2.076E-06 * 1e6;
2141  }
2142  if (t <= 0.40)
2143  {
2144  return 1.131E-06 * 1e6;
2145  }
2146  if (t <= 0.41)
2147  {
2148  return 1.861E-07 * 1e6;
2149  }
2150  if (t <= 0.42)
2151  {
2152  return -4.577E-0 * 1e6;
2153  }
2154  if (t <= 0.43)
2155  {
2156  return -4.657E-0 * 1e6;
2157  }
2158  if (t <= 0.44)
2159  {
2160  return 3.124E-07 * 1e6;
2161  }
2162  if (t <= 0.45)
2163  {
2164  return 1.684E-06 * 1e6;
2165  }
2166  if (t <= 0.46)
2167  {
2168  return 3.174E-06 * 1e6;
2169  }
2170  if (t <= 0.47)
2171  {
2172  return 4.306E-06 * 1e6;
2173  }
2174  if (t <= 0.48)
2175  {
2176  return 4.873E-06 * 1e6;
2177  }
2178  if (t <= 0.49)
2179  {
2180  return 4.980E-06 * 1e6;
2181  }
2182  if (t <= 0.50)
2183  {
2184  return 4.876E-06 * 1e6;
2185  }
2186  if (t <= 0.51)
2187  {
2188  return 4.757E-06 * 1e6;
2189  }
2190  if (t <= 0.52)
2191  {
2192  return 4.681E-06 * 1e6;
2193  }
2194  if (t <= 0.53)
2195  {
2196  return 4.603E-06 * 1e6;
2197  }
2198  if (t <= 0.54)
2199  {
2200  return 4.458E-06 * 1e6;
2201  }
2202  if (t <= 0.55)
2203  {
2204  return 4.224E-06 * 1e6;
2205  }
2206  if (t <= 0.56)
2207  {
2208  return 3.928E-06 * 1e6;
2209  }
2210  if (t <= 0.57)
2211  {
2212  return 3.618E-06 * 1e6;
2213  }
2214  if (t <= 0.58)
2215  {
2216  return 3.335E-06 * 1e6;
2217  }
2218  if (t <= 0.59)
2219  {
2220  return 3.098E-06 * 1e6;
2221  }
2222  if (t <= 0.60)
2223  {
2224  return 2.913E-06 * 1e6;
2225  }
2226  if (t <= 0.61)
2227  {
2228  return 2.778E-06 * 1e6;
2229  }
2230  if (t <= 0.62)
2231  {
2232  return 2.697E-06 * 1e6;
2233  }
2234  if (t <= 0.63)
2235  {
2236  return 2.670E-06 * 1e6;
2237  }
2238  if (t <= 0.64)
2239  {
2240  return 2.692E-06 * 1e6;
2241  }
2242  if (t <= 0.65)
2243  {
2244  return 2.746E-06 * 1e6;
2245  }
2246  if (t <= 0.66)
2247  {
2248  return 2.811E-06 * 1e6;
2249  }
2250  if (t <= 0.67)
2251  {
2252  return 2.868E-06 * 1e6;
2253  }
2254  if (t <= 0.68)
2255  {
2256  return 2.904E-06 * 1e6;
2257  }
2258  if (t <= 0.69)
2259  {
2260  return 2.915E-06 * 1e6;
2261  }
2262  if (t <= 0.70)
2263  {
2264  return 2.903E-06 * 1e6;
2265  }
2266  if (t <= 0.71)
2267  {
2268  return 2.872E-06 * 1e6;
2269  }
2270  if (t <= 0.72)
2271  {
2272  return 2.826E-06 * 1e6;
2273  }
2274  if (t <= 0.73)
2275  {
2276  return 2.769E-06 * 1e6;
2277  }
2278  if (t <= 0.74)
2279  {
2280  return 2.704E-06 * 1e6;
2281  }
2282  if (t <= 0.75)
2283  {
2284  return 2.636E-06 * 1e6;
2285  }
2286  if (t <= 0.76)
2287  {
2288  return 2.569E-06 * 1e6;
2289  }
2290  if (t <= 0.77)
2291  {
2292  return 2.511E-06 * 1e6;
2293  }
2294  if (t <= 0.78)
2295  {
2296  return 2.465E-06 * 1e6;
2297  }
2298  if (t <= 0.79)
2299  {
2300  return 2.433E-06 * 1e6;
2301  }
2302 }// 15, LCCA, branch 2
2303 
2304 Real aortaFlux9 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2305 {
2306  if (t <= 0.00)
2307  {
2308  return 7.854E-07 * 1e6;
2309  }
2310  if (t <= 0.01)
2311  {
2312  return 7.900E-07 * 1e6;
2313  }
2314  if (t <= 0.02)
2315  {
2316  return 7.988E-07 * 1e6;
2317  }
2318  if (t <= 0.03)
2319  {
2320  return 8.077E-07 * 1e6;
2321  }
2322  if (t <= 0.04)
2323  {
2324  return 8.131E-07 * 1e6;
2325  }
2326  if (t <= 0.05)
2327  {
2328  return 8.141E-07 * 1e6;
2329  }
2330  if (t <= 0.06)
2331  {
2332  return 8.120E-07 * 1e6;
2333  }
2334  if (t <= 0.07)
2335  {
2336  return 8.085E-07 * 1e6;
2337  }
2338  if (t <= 0.08)
2339  {
2340  return 8.082E-07 * 1e6;
2341  }
2342  if (t <= 0.09)
2343  {
2344  return 8.255E-07 * 1e6;
2345  }
2346  if (t <= 0.10)
2347  {
2348  return 8.927E-07 * 1e6;
2349  }
2350  if (t <= 0.11)
2351  {
2352  return 1.055E-06 * 1e6;
2353  }
2354  if (t <= 0.12)
2355  {
2356  return 1.343E-06 * 1e6;
2357  }
2358  if (t <= 0.13)
2359  {
2360  return 1.728E-06 * 1e6;
2361  }
2362  if (t <= 0.14)
2363  {
2364  return 2.118E-06 * 1e6;
2365  }
2366  if (t <= 0.15)
2367  {
2368  return 2.419E-06 * 1e6;
2369  }
2370  if (t <= 0.16)
2371  {
2372  return 2.598E-06 * 1e6;
2373  }
2374  if (t <= 0.17)
2375  {
2376  return 2.673E-06 * 1e6;
2377  }
2378  if (t <= 0.18)
2379  {
2380  return 2.642E-06 * 1e6;
2381  }
2382  if (t <= 0.19)
2383  {
2384  return 2.465E-06 * 1e6;
2385  }
2386  if (t <= 0.20)
2387  {
2388  return 2.118E-06 * 1e6;
2389  }
2390  if (t <= 0.21)
2391  {
2392  return 1.660E-06 * 1e6;
2393  }
2394  if (t <= 0.22)
2395  {
2396  return 1.210E-06 * 1e6;
2397  }
2398  if (t <= 0.23)
2399  {
2400  return 8.723E-07 * 1e6;
2401  }
2402  if (t <= 0.24)
2403  {
2404  return 6.812E-07 * 1e6;
2405  }
2406  if (t <= 0.25)
2407  {
2408  return 6.155E-07 * 1e6;
2409  }
2410  if (t <= 0.26)
2411  {
2412  return 6.323E-07 * 1e6;
2413  }
2414  if (t <= 0.27)
2415  {
2416  return 6.927E-07 * 1e6;
2417  }
2418  if (t <= 0.28)
2419  {
2420  return 7.708E-07 * 1e6;
2421  }
2422  if (t <= 0.29)
2423  {
2424  return 8.538E-07 * 1e6;
2425  }
2426  if (t <= 0.30)
2427  {
2428  return 9.397E-07 * 1e6;
2429  }
2430  if (t <= 0.31)
2431  {
2432  return 1.030E-06 * 1e6;
2433  }
2434  if (t <= 0.32)
2435  {
2436  return 1.120E-06 * 1e6;
2437  }
2438  if (t <= 0.33)
2439  {
2440  return 1.203E-06 * 1e6;
2441  }
2442  if (t <= 0.34)
2443  {
2444  return 1.263E-06 * 1e6;
2445  }
2446  if (t <= 0.35)
2447  {
2448  return 1.292E-06 * 1e6;
2449  }
2450  if (t <= 0.36)
2451  {
2452  return 1.288E-06 * 1e6;
2453  }
2454  if (t <= 0.37)
2455  {
2456  return 1.258E-06 * 1e6;
2457  }
2458  if (t <= 0.38)
2459  {
2460  return 1.207E-06 * 1e6;
2461  }
2462  if (t <= 0.39)
2463  {
2464  return 1.142E-06 * 1e6;
2465  }
2466  if (t <= 0.40)
2467  {
2468  return 1.064E-06 * 1e6;
2469  }
2470  if (t <= 0.41)
2471  {
2472  return 9.816E-07 * 1e6;
2473  }
2474  if (t <= 0.42)
2475  {
2476  return 9.133E-07 * 1e6;
2477  }
2478  if (t <= 0.43)
2479  {
2480  return 8.870E-07 * 1e6;
2481  }
2482  if (t <= 0.44)
2483  {
2484  return 9.268E-07 * 1e6;
2485  }
2486  if (t <= 0.45)
2487  {
2488  return 1.035E-06 * 1e6;
2489  }
2490  if (t <= 0.46)
2491  {
2492  return 1.183E-06 * 1e6;
2493  }
2494  if (t <= 0.47)
2495  {
2496  return 1.329E-06 * 1e6;
2497  }
2498  if (t <= 0.48)
2499  {
2500  return 1.445E-06 * 1e6;
2501  }
2502  if (t <= 0.49)
2503  {
2504  return 1.529E-06 * 1e6;
2505  }
2506  if (t <= 0.50)
2507  {
2508  return 1.587E-06 * 1e6;
2509  }
2510  if (t <= 0.51)
2511  {
2512  return 1.612E-06 * 1e6;
2513  }
2514  if (t <= 0.52)
2515  {
2516  return 1.587E-06 * 1e6;
2517  }
2518  if (t <= 0.53)
2519  {
2520  return 1.501E-06 * 1e6;
2521  }
2522  if (t <= 0.54)
2523  {
2524  return 1.370E-06 * 1e6;
2525  }
2526  if (t <= 0.55)
2527  {
2528  return 1.230E-06 * 1e6;
2529  }
2530  if (t <= 0.56)
2531  {
2532  return 1.112E-06 * 1e6;
2533  }
2534  if (t <= 0.57)
2535  {
2536  return 1.029E-06 * 1e6;
2537  }
2538  if (t <= 0.58)
2539  {
2540  return 9.776E-07 * 1e6;
2541  }
2542  if (t <= 0.59)
2543  {
2544  return 9.460E-07 * 1e6;
2545  }
2546  if (t <= 0.60)
2547  {
2548  return 9.239E-07 * 1e6;
2549  }
2550  if (t <= 0.61)
2551  {
2552  return 9.061E-07 * 1e6;
2553  }
2554  if (t <= 0.62)
2555  {
2556  return 8.917E-07 * 1e6;
2557  }
2558  if (t <= 0.63)
2559  {
2560  return 8.821E-07 * 1e6;
2561  }
2562  if (t <= 0.64)
2563  {
2564  return 8.786E-07 * 1e6;
2565  }
2566  if (t <= 0.65)
2567  {
2568  return 8.809E-07 * 1e6;
2569  }
2570  if (t <= 0.66)
2571  {
2572  return 8.871E-07 * 1e6;
2573  }
2574  if (t <= 0.67)
2575  {
2576  return 8.940E-07 * 1e6;
2577  }
2578  if (t <= 0.68)
2579  {
2580  return 8.978E-07 * 1e6;
2581  }
2582  if (t <= 0.69)
2583  {
2584  return 8.959E-07 * 1e6;
2585  }
2586  if (t <= 0.70)
2587  {
2588  return 8.874E-07 * 1e6;
2589  }
2590  if (t <= 0.71)
2591  {
2592  return 8.734E-07 * 1e6;
2593  }
2594  if (t <= 0.72)
2595  {
2596  return 8.561E-07 * 1e6;
2597  }
2598  if (t <= 0.73)
2599  {
2600  return 8.383E-07 * 1e6;
2601  }
2602  if (t <= 0.74)
2603  {
2604  return 8.221E-07 * 1e6;
2605  }
2606  if (t <= 0.75)
2607  {
2608  return 8.088E-07 * 1e6;
2609  }
2610  if (t <= 0.76)
2611  {
2612  return 7.986E-07 * 1e6;
2613  }
2614  if (t <= 0.77)
2615  {
2616  return 7.909E-07 * 1e6;
2617  }
2618  if (t <= 0.78)
2619  {
2620  return 7.852E-07 * 1e6;
2621  }
2622  if (t <= 0.79)
2623  {
2624  return 7.807E-07 * 1e6;
2625  }
2626 }// 20 LVA branch 3_1
2627 
2628 Real aortaFlux4 (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2629 {
2630  if (t <= 0.00)
2631  {
2632  return 2.581E-06 * 1e6;
2633  }
2634  if (t <= 0.01)
2635  {
2636  return 2.499E-06 * 1e6;
2637  }
2638  if (t <= 0.02)
2639  {
2640  return 2.399E-06 * 1e6;
2641  }
2642  if (t <= 0.03)
2643  {
2644  return 2.281E-06 * 1e6;
2645  }
2646  if (t <= 0.04)
2647  {
2648  return 2.147E-06 * 1e6;
2649  }
2650  if (t <= 0.05)
2651  {
2652  return 2.002E-06 * 1e6;
2653  }
2654  if (t <= 0.06)
2655  {
2656  return 1.849E-06 * 1e6;
2657  }
2658  if (t <= 0.07)
2659  {
2660  return 1.695E-06 * 1e6;
2661  }
2662  if (t <= 0.08)
2663  {
2664  return 1.565E-06 * 1e6;
2665  }
2666  if (t <= 0.09)
2667  {
2668  return 1.537E-06 * 1e6;
2669  }
2670  if (t <= 0.10)
2671  {
2672  return 1.766E-06 * 1e6;
2673  }
2674  if (t <= 0.11)
2675  {
2676  return 2.460E-06 * 1e6;
2677  }
2678  if (t <= 0.12)
2679  {
2680  return 3.773E-06 * 1e6;
2681  }
2682  if (t <= 0.13)
2683  {
2684  return 5.709E-06 * 1e6;
2685  }
2686  if (t <= 0.14)
2687  {
2688  return 8.131E-06 * 1e6;
2689  }
2690  if (t <= 0.15)
2691  {
2692  return 1.087E-05 * 1e6;
2693  }
2694  if (t <= 0.16)
2695  {
2696  return 1.379E-05 * 1e6;
2697  }
2698  if (t <= 0.17)
2699  {
2700  return 1.675E-05 * 1e6;
2701  }
2702  if (t <= 0.18)
2703  {
2704  return 1.957E-05 * 1e6;
2705  }
2706  if (t <= 0.19)
2707  {
2708  return 2.205E-05 * 1e6;
2709  }
2710  if (t <= 0.20)
2711  {
2712  return 2.403E-05 * 1e6;
2713  }
2714  if (t <= 0.21)
2715  {
2716  return 2.538E-05 * 1e6;
2717  }
2718  if (t <= 0.22)
2719  {
2720  return 2.606E-05 * 1e6;
2721  }
2722  if (t <= 0.23)
2723  {
2724  return 2.611E-05 * 1e6;
2725  }
2726  if (t <= 0.24)
2727  {
2728  return 2.562E-05 * 1e6;
2729  }
2730  if (t <= 0.25)
2731  {
2732  return 2.474E-05 * 1e6;
2733  }
2734  if (t <= 0.26)
2735  {
2736  return 2.362E-05 * 1e6;
2737  }
2738  if (t <= 0.27)
2739  {
2740  return 2.238E-05 * 1e6;
2741  }
2742  if (t <= 0.28)
2743  {
2744  return 2.111E-05 * 1e6;
2745  }
2746  if (t <= 0.29)
2747  {
2748  return 1.981E-05 * 1e6;
2749  }
2750  if (t <= 0.30)
2751  {
2752  return 1.850E-05 * 1e6;
2753  }
2754  if (t <= 0.31)
2755  {
2756  return 1.715E-05 * 1e6;
2757  }
2758  if (t <= 0.32)
2759  {
2760  return 1.576E-05 * 1e6;
2761  }
2762  if (t <= 0.33)
2763  {
2764  return 1.432E-05 * 1e6;
2765  }
2766  if (t <= 0.34)
2767  {
2768  return 1.284E-05 * 1e6;
2769  }
2770  if (t <= 0.35)
2771  {
2772  return 1.132E-05 * 1e6;
2773  }
2774  if (t <= 0.36)
2775  {
2776  return 9.768E-06 * 1e6;
2777  }
2778  if (t <= 0.37)
2779  {
2780  return 8.180E-06 * 1e6;
2781  }
2782  if (t <= 0.38)
2783  {
2784  return 6.543E-06 * 1e6;
2785  }
2786  if (t <= 0.39)
2787  {
2788  return 4.831E-06 * 1e6;
2789  }
2790  if (t <= 0.40)
2791  {
2792  return 3.030E-06 * 1e6;
2793  }
2794  if (t <= 0.41)
2795  {
2796  return 1.163E-06 * 1e6;
2797  }
2798  if (t <= 0.42)
2799  {
2800  return -6.817E-07 * 1e6;
2801  }
2802  if (t <= 0.43)
2803  {
2804  return -2.362E-06 * 1e6;
2805  }
2806  if (t <= 0.44)
2807  {
2808  return -3.738E-06 * 1e6;
2809  }
2810  if (t <= 0.45)
2811  {
2812  return -4.742E-06 * 1e6;
2813  }
2814  if (t <= 0.46)
2815  {
2816  return -5.400E-06 * 1e6;
2817  }
2818  if (t <= 0.47)
2819  {
2820  return -5.785E-06 * 1e6;
2821  }
2822  if (t <= 0.48)
2823  {
2824  return -5.956E-06 * 1e6;
2825  }
2826  if (t <= 0.49)
2827  {
2828  return -5.932E-06 * 1e6;
2829  }
2830  if (t <= 0.50)
2831  {
2832  return -5.723E-06 * 1e6;
2833  }
2834  if (t <= 0.51)
2835  {
2836  return -5.358E-06 * 1e6;
2837  }
2838  if (t <= 0.52)
2839  {
2840  return -4.889E-06 * 1e6;
2841  }
2842  if (t <= 0.53)
2843  {
2844  return -4.370E-06 * 1e6;
2845  }
2846  if (t <= 0.54)
2847  {
2848  return -3.846E-06 * 1e6;
2849  }
2850  if (t <= 0.55)
2851  {
2852  return -3.341E-06 * 1e6;
2853  }
2854  if (t <= 0.56)
2855  {
2856  return -2.866E-06 * 1e6;
2857  }
2858  if (t <= 0.57)
2859  {
2860  return -2.422E-06 * 1e6;
2861  }
2862  if (t <= 0.58)
2863  {
2864  return -2.004E-06 * 1e6;
2865  }
2866  if (t <= 0.59)
2867  {
2868  return -1.601E-06 * 1e6;
2869  }
2870  if (t <= 0.60)
2871  {
2872  return -1.206E-06 * 1e6;
2873  }
2874  if (t <= 0.61)
2875  {
2876  return -8.123E-07 * 1e6;
2877  }
2878  if (t <= 0.62)
2879  {
2880  return -4.195E-07 * 1e6;
2881  }
2882  if (t <= 0.63)
2883  {
2884  return -3.198E-08 * 1e6;
2885  }
2886  if (t <= 0.64)
2887  {
2888  return 3.429E-07 * 1e6;
2889  }
2890  if (t <= 0.65)
2891  {
2892  return 6.970E-07 * 1e6;
2893  }
2894  if (t <= 0.66)
2895  {
2896  return 1.024E-06 * 1e6;
2897  }
2898  if (t <= 0.67)
2899  {
2900  return 1.320E-06 * 1e6;
2901  }
2902  if (t <= 0.68)
2903  {
2904  return 1.583E-06 * 1e6;
2905  }
2906  if (t <= 0.69)
2907  {
2908  return 1.813E-06 * 1e6;
2909  }
2910  if (t <= 0.70)
2911  {
2912  return 2.011E-06 * 1e6;
2913  }
2914  if (t <= 0.71)
2915  {
2916  return 2.180E-06 * 1e6;
2917  }
2918  if (t <= 0.72)
2919  {
2920  return 2.323E-06 * 1e6;
2921  }
2922  if (t <= 0.73)
2923  {
2924  return 2.442E-06 * 1e6;
2925  }
2926  if (t <= 0.74)
2927  {
2928  return 2.539E-06 * 1e6;
2929  }
2930  if (t <= 0.75)
2931  {
2932  return 2.613E-06 * 1e6;
2933  }
2934  if (t <= 0.76)
2935  {
2936  return 2.665E-06 * 1e6;
2937  }
2938  if (t <= 0.77)
2939  {
2940  return 2.691E-06 * 1e6;
2941  }
2942  if (t <= 0.78)
2943  {
2944  return 2.690E-06 * 1e6;
2945  }
2946  if (t <= 0.79)
2947  {
2948  return 2.661E-06 * 1e6;
2949  }
2950 }//21, L. Brachia, bhanch 3_2
2951 
2952 
2953 
2954 
2955 
2956 
2957 
2958 Real f (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2959 {
2960  return 0.;
2961 }
2962 
2963 Real u1 (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2964 {
2965  return 0.0;
2966 }
2967 
2968 Real fZero (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2969 {
2970  return 0.0;
2971 }
2972 
2973 // Initial velocity
2974 Real u0 (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2975 {
2976  return 0.0;
2977 }
2978 
2979 Real p0 (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2980 {
2981  return 0.0;
2982 }
2983 
2984 
2985 Real E (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
2986 {
2987  return 1000;
2988  // if(z<5)
2989  // return /*-*/10*1e3;//-29;//*e-1*5; // about [(110-60)*(133.332*10)]/[10*(2.08-1.8)] ( 10 because of mm, 133... because of mmHg)
2990  // // (see paper by Liu, Dang, etc). in their plot x and y are inverted.
2991  // if(z>5 && z<6)
2992  // return 10+20*(z-5)*1e3;
2993  // if(z>=6)
2994  // return 30*1e3;
2995  // if ( z<0 && z>=-3 )
2996  // return 10-3*(z)*1e3;
2997  // if(z<-3)
2998  // return 19*1e3;
2999 }
3000 
3001 
3002 // Initial displacement and velocity
3003 Real d0 (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& i)
3004 {
3005  switch (i)
3006  {
3007  case 1:
3008  return 0.;
3009  break;
3010  case 2:
3011  return 0.;
3012  break;
3013  case 3:
3014  return 0.;
3015  break;
3016  default:
3017  ERROR_MSG ("This entrie is not allowed: ud_functions.hpp");
3018  break;
3019  }
3020 }
3021 
3022 Real w0 (const Real& /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& i)
3023 {
3024 
3025  switch (i)
3026  {
3027  case 1:
3028  return 0.0;
3029  break;
3030  case 2:
3031  return 0.0;
3032  break;
3033  case 3:
3034  return 0.0;
3035  break;
3036  default:
3037  ERROR_MSG ("This entrie is not allowed: ud_functions.hpp");
3038  break;
3039  }
3040 }
3041 
3042 
3043 
3044 Real linearPress2 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3045 {
3046  if (t > 0.8)
3047  {
3048  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3049  }
3050 
3051  Real ti = floor (t * 1000) / 1000;
3052  Real tii = ti + 0.001;
3053  return (aortaPhisPress (tii) - aortaPhisPress (ti) ) / (0.001) * (t - (ti) ) + aortaPhisPress (ti) + 115000;
3054 }
3055 
3056 
3057 Real linearFlux3 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3058 {
3059  if (t > 0.8)
3060  {
3061  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3062  }
3063 
3064  Real ti = floor (t * 1000) / 1000;
3065  Real tii = ti + 0.001;
3066  return (aortaFlux3 (tii) - aortaFlux3 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux3 (ti);
3067 }
3068 
3069 
3070 Real linearFluxIn (Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3071 {
3072  if (t > 0.8)
3073  {
3074  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3075  }
3076 
3077  Real ti = floor (t * 1000) / 1000;
3078  Real tii = ti + 0.001;
3079 
3080  return (aortaFluxIn (tii) - aortaFluxIn (ti) ) / (0.001) * (t - (ti) ) + aortaFluxIn (ti);
3081 }
3082 
3083 
3084 Real linearFlux4 (Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3085 {
3086  if (t > 0.8)
3087  {
3088  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3089  }
3090 
3091  Real ti = floor (t * 1000) / 1000;
3092  Real tii = ti + 0.001;
3093  return (aortaFlux4 (tii) - aortaFlux4 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux4 (ti);
3094 }
3095 
3096 
3097 
3098 Real linearFlux5 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3099 {
3100  if (t > 0.8)
3101  {
3102  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3103  }
3104 
3105  Real ti = floor (t * 1000) / 1000;
3106  Real tii = ti + 0.001;
3107  return (aortaFlux5 (tii) - aortaFlux5 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux5 (ti);
3108 }
3109 
3110 
3111 
3112 Real linearFlux6 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3113 {
3114  if (t > 0.8)
3115  {
3116  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3117  }
3118 
3119  Real ti = floor (t * 1000) / 1000;
3120  Real tii = ti + 0.001;
3121  return (aortaFlux6 (tii) - aortaFlux6 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux6 (ti);
3122 }
3123 
3124 
3125 Real linearFlux7 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3126 {
3127  if (t > 0.8)
3128  {
3129  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3130  }
3131 
3132  Real ti = floor (t * 1000) / 1000;
3133  Real tii = ti + 0.001;
3134  return (aortaFlux7 (tii) - aortaFlux7 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux7 (ti);
3135 }
3136 
3137 
3138 Real linearFlux8 (Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3139 {
3140  if (t > 0.8)
3141  {
3142  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3143  }
3144 
3145  Real ti = floor (t * 1000) / 1000;
3146  Real tii = ti + 0.001;
3147  return (aortaFlux8 (tii) - aortaFlux8 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux8 (ti);
3148 }
3149 
3150 Real linearFlux9 ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3151 {
3152  if (t > 0.8)
3153  {
3154  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3155  }
3156 
3157  Real ti = floor (t * 1000) / 1000;
3158  Real tii = ti + 0.001;
3159  return (aortaFlux9 (tii) - aortaFlux9 (ti) ) / (0.001) * (t - (ti) ) + aortaFlux9 (ti);
3160 }
3161 
3162 
3163 Real linearFlux3_ (Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3164 {
3165  if (t > 0.8)
3166  {
3167  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3168  }
3169 
3170  Real ti = floor (t * 1000) / 1000;
3171  Real tii = ti + 0.001;
3172  return (aortaFlux3_ (tii) - aortaFlux3_ (ti) ) / (0.001) * (t - (ti) ) + aortaFlux3_ (ti);
3173 }
3174 
3175 
3176 Real linearFlux6_ ( Real t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3177 {
3178  if (t > 0.8)
3179  {
3180  t = ( ( (int) floor (t * 1000) ) % 800) / 1000;
3181  }
3182 
3183  Real ti = floor (t * 1000) / 1000;
3184  Real tii = ti + 0.001;
3185  return (aortaFlux6_ (tii) - aortaFlux6_ (ti) ) / (0.001) * (t - (ti) ) + aortaFlux6_ (ti);
3186 }
3187 
3188 Real aortaFlux3_ (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3189 {
3190  if (t <= 0.00 + 0.01)
3191  {
3192  return 33.6 ;
3193  }
3194  if (t <= 0.01 + 0.01)
3195  {
3196  return 33.83 ;
3197  }
3198  if (t <= 0.02 + 0.01)
3199  {
3200  return 34.12 ;
3201  }
3202  if (t <= 0.03 + 0.01)
3203  {
3204  return 34.44 ;
3205  }
3206  if (t <= 0.04 + 0.01)
3207  {
3208  return 34.76 ;
3209  }
3210  if (t <= 0.05 + 0.01)
3211  {
3212  return 35.05 ;
3213  }
3214  if (t <= 0.06 + 0.01)
3215  {
3216  return 35.29 ;
3217  }
3218  if (t <= 0.07 + 0.01)
3219  {
3220  return 35.49 ;
3221  }
3222  if (t <= 0.08 + 0.01)
3223  {
3224  return 35.74 ;
3225  }
3226  if (t <= 0.09 + 0.01)
3227  {
3228  return 36.27 ;
3229  }
3230  if (t <= 0.10 + 0.01)
3231  {
3232  return 37.83 ;
3233  }
3234  if (t <= 0.11 + 0.01)
3235  {
3236  return 41.86 ;
3237  }
3238  if (t <= 0.12 + 0.01)
3239  {
3240  return 50.47 ;
3241  }
3242  if (t <= 0.13 + 0.01)
3243  {
3244  return 65.55999999999999;
3245  }
3246  if (t <= 0.14 + 0.01)
3247  {
3248  return 87.10999999999999;
3249  }
3250  if (t <= 0.15 + 0.01)
3251  {
3252  return 111.8 ;
3253  }
3254  if (t <= 0.16 + 0.01)
3255  {
3256  return 134.6 ;
3257  }
3258  if (t <= 0.17 + 0.01)
3259  {
3260  return 152 ;
3261  }
3262  if (t <= 0.18 + 0.01)
3263  {
3264  return 164.3 ;
3265  }
3266  if (t <= 0.19 + 0.01)
3267  {
3268  return 173.8 ;
3269  }
3270  if (t <= 0.20 + 0.01)
3271  {
3272  return 182.2 ;
3273  }
3274  if (t <= 0.21 + 0.01)
3275  {
3276  return 189.8 ;
3277  }
3278  if (t <= 0.22 + 0.01)
3279  {
3280  return 195.9 ;
3281  }
3282  if (t <= 0.23 + 0.01)
3283  {
3284  return 200 ;
3285  }
3286  if (t <= 0.24 + 0.01)
3287  {
3288  return 202 ;
3289  }
3290  if (t <= 0.25 + 0.01)
3291  {
3292  return 202.1 ;
3293  }
3294  if (t <= 0.26 + 0.01)
3295  {
3296  return 200.5 ;
3297  }
3298  if (t <= 0.27 + 0.01)
3299  {
3300  return 197.3 ;
3301  }
3302  if (t <= 0.28 + 0.01)
3303  {
3304  return 192.7 ;
3305  }
3306  if (t <= 0.29 + 0.01)
3307  {
3308  return 186.9 ;
3309  }
3310  if (t <= 0.30 + 0.01)
3311  {
3312  return 179.9 ;
3313  }
3314  if (t <= 0.31 + 0.01)
3315  {
3316  return 172 ;
3317  }
3318  if (t <= 0.32 + 0.01)
3319  {
3320  return 163.3 ;
3321  }
3322  if (t <= 0.33 + 0.01)
3323  {
3324  return 154.1 ;
3325  }
3326  if (t <= 0.34 + 0.01)
3327  {
3328  return 144.7 ;
3329  }
3330  if (t <= 0.35 + 0.01)
3331  {
3332  return 135.1 ;
3333  }
3334  if (t <= 0.36 + 0.01)
3335  {
3336  return 125.5 ;
3337  }
3338  if (t <= 0.37 + 0.01)
3339  {
3340  return 115.9 ;
3341  }
3342  if (t <= 0.38 + 0.01)
3343  {
3344  return 106.3 ;
3345  }
3346  if (t <= 0.39 + 0.01)
3347  {
3348  return 96.60999999999999;
3349  }
3350  if (t <= 0.40 + 0.01)
3351  {
3352  return 86.44 ;
3353  }
3354  if (t <= 0.41 + 0.01)
3355  {
3356  return 75.67999999999999;
3357  }
3358  if (t <= 0.42 + 0.01)
3359  {
3360  return 64.56999999999999;
3361  }
3362  if (t <= 0.43 + 0.01)
3363  {
3364  return 53.92 ;
3365  }
3366  if (t <= 0.44 + 0.01)
3367  {
3368  return 44.94 ;
3369  }
3370  if (t <= 0.45 + 0.01)
3371  {
3372  return 38.75 ;
3373  }
3374  if (t <= 0.46 + 0.01)
3375  {
3376  return 35.66 ;
3377  }
3378  if (t <= 0.47 + 0.01)
3379  {
3380  return 34.83000000000001;
3381  }
3382  if (t <= 0.48 + 0.01)
3383  {
3384  return 34.67 ;
3385  }
3386  if (t <= 0.49 + 0.01)
3387  {
3388  return 33.83 ;
3389  }
3390  if (t <= 0.50 + 0.01)
3391  {
3392  return 32.01 ;
3393  }
3394  if (t <= 0.51 + 0.01)
3395  {
3396  return 29.85 ;
3397  }
3398  if (t <= 0.52 + 0.01)
3399  {
3400  return 28.19 ;
3401  }
3402  if (t <= 0.53 + 0.01)
3403  {
3404  return 27.4 ;
3405  }
3406  if (t <= 0.54 + 0.01)
3407  {
3408  return 27.28 ;
3409  }
3410  if (t <= 0.55 + 0.01)
3411  {
3412  return 27.42 ;
3413  }
3414  if (t <= 0.56 + 0.01)
3415  {
3416  return 27.54 ;
3417  }
3418  if (t <= 0.57 + 0.01)
3419  {
3420  return 27.63 ;
3421  }
3422  if (t <= 0.58 + 0.01)
3423  {
3424  return 27.82 ;
3425  }
3426  if (t <= 0.59 + 0.01)
3427  {
3428  return 28.21 ;
3429  }
3430  if (t <= 0.60 + 0.01)
3431  {
3432  return 28.76 ;
3433  }
3434  if (t <= 0.61 + 0.01)
3435  {
3436  return 29.39 ;
3437  }
3438  if (t <= 0.62 + 0.01)
3439  {
3440  return 30 ;
3441  }
3442  if (t <= 0.63 + 0.01)
3443  {
3444  return 30.54 ;
3445  }
3446  if (t <= 0.64 + 0.01)
3447  {
3448  return 31.01 ;
3449  }
3450  if (t <= 0.65 + 0.01)
3451  {
3452  return 31.42 ;
3453  }
3454  if (t <= 0.66 + 0.01)
3455  {
3456  return 31.78 ;
3457  }
3458  if (t <= 0.67 + 0.01)
3459  {
3460  return 32.09 ;
3461  }
3462  if (t <= 0.68 + 0.01)
3463  {
3464  return 32.34 ;
3465  }
3466  if (t <= 0.69 + 0.01)
3467  {
3468  return 32.54 ;
3469  }
3470  if (t <= 0.70 + 0.01)
3471  {
3472  return 32.69 ;
3473  }
3474  if (t <= 0.71 + 0.01)
3475  {
3476  return 32.8 ;
3477  }
3478  if (t <= 0.72 + 0.01)
3479  {
3480  return 32.87 ;
3481  }
3482  if (t <= 0.73 + 0.01)
3483  {
3484  return 32.91 ;
3485  }
3486  if (t <= 0.74 + 0.01)
3487  {
3488  return 32.92 ;
3489  }
3490  if (t <= 0.75 + 0.01)
3491  {
3492  return 32.93000000000001;
3493  }
3494  if (t <= 0.76 + 0.01)
3495  {
3496  return 32.93000000000001;
3497  }
3498  if (t <= 0.77 + 0.01)
3499  {
3500  return 32.96 ;
3501  }
3502  if (t <= 0.78 + 0.01)
3503  {
3504  return 33.01000000000001;
3505  }
3506  if (t <= 0.79 + 0.01)
3507  {
3508  return 33.1 ;
3509  }
3510  RETURN_UNDEFINED;
3511 }
3512 
3513 Real aortaFlux6_ (const Real& t, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3514 {
3515  if (t <= 0.00 + 0.01)
3516  {
3517  return 0.6737071250000001;
3518  }
3519  if (t <= 0.01 + 0.01)
3520  {
3521  return 0.6799071250000001;
3522  }
3523  if (t <= 0.02 + 0.01)
3524  {
3525  return 0.6897071250000001;
3526  }
3527  if (t <= 0.03 + 0.01)
3528  {
3529  return 0.6997071250000001;
3530  }
3531  if (t <= 0.04 + 0.01)
3532  {
3533  return 0.706407125 ;
3534  }
3535  if (t <= 0.05 + 0.01)
3536  {
3537  return 0.7084071250000001;
3538  }
3539  if (t <= 0.06 + 0.01)
3540  {
3541  return 0.7056071250000001;
3542  }
3543  if (t <= 0.07 + 0.01)
3544  {
3545  return 0.7049071250000001;
3546  }
3547  if (t <= 0.08 + 0.01)
3548  {
3549  return 0.7382071250000001;
3550  }
3551  if (t <= 0.09 + 0.01)
3552  {
3553  return 0.8778071249999999;
3554  }
3555  if (t <= 0.10 + 0.01)
3556  {
3557  return 1.209007125 ;
3558  }
3559  if (t <= 0.11 + 0.01)
3560  {
3561  return 1.746007125 ;
3562  }
3563  if (t <= 0.12 + 0.01)
3564  {
3565  return 2.368007125 ;
3566  }
3567  if (t <= 0.13 + 0.01)
3568  {
3569  return 2.897007125 ;
3570  }
3571  if (t <= 0.14 + 0.01)
3572  {
3573  return 3.258007125 ;
3574  }
3575  if (t <= 0.15 + 0.01)
3576  {
3577  return 3.518007125 ;
3578  }
3579  if (t <= 0.16 + 0.01)
3580  {
3581  return 3.752007125 ;
3582  }
3583  if (t <= 0.17 + 0.01)
3584  {
3585  return 3.936007124999999 ;
3586  }
3587  if (t <= 0.18 + 0.01)
3588  {
3589  return 3.957007125 ;
3590  }
3591  if (t <= 0.19 + 0.01)
3592  {
3593  return 3.716007125 ;
3594  }
3595  if (t <= 0.20 + 0.01)
3596  {
3597  return 3.212007125 ;
3598  }
3599  if (t <= 0.21 + 0.01)
3600  {
3601  return 2.551007125 ;
3602  }
3603  if (t <= 0.22 + 0.01)
3604  {
3605  return 1.898007125 ;
3606  }
3607  if (t <= 0.23 + 0.01)
3608  {
3609  return 1.395007125 ;
3610  }
3611  if (t <= 0.24 + 0.01)
3612  {
3613  return 1.107007125 ;
3614  }
3615  if (t <= 0.25 + 0.01)
3616  {
3617  return 1.009007125 ;
3618  }
3619  if (t <= 0.26 + 0.01)
3620  {
3621  return 1.028007125 ;
3622  }
3623  if (t <= 0.27 + 0.01)
3624  {
3625  return 1.082007125 ;
3626  }
3627  if (t <= 0.28 + 0.01)
3628  {
3629  return 1.114007125 ;
3630  }
3631  if (t <= 0.29 + 0.01)
3632  {
3633  return 1.111007125 ;
3634  }
3635  if (t <= 0.30 + 0.01)
3636  {
3637  return 1.089007125 ;
3638  }
3639  if (t <= 0.31 + 0.01)
3640  {
3641  return 1.071007125 ;
3642  }
3643  if (t <= 0.32 + 0.01)
3644  {
3645  return 1.067007125 ;
3646  }
3647  if (t <= 0.33 + 0.01)
3648  {
3649  return 1.068007125 ;
3650  }
3651  if (t <= 0.34 + 0.01)
3652  {
3653  return 1.056007125 ;
3654  }
3655  if (t <= 0.35 + 0.01)
3656  {
3657  return 1.021007125 ;
3658  }
3659  if (t <= 0.36 + 0.01)
3660  {
3661  return 0.955007125 ;
3662  }
3663  if (t <= 0.37 + 0.01)
3664  {
3665  return 0.856707125 ;
3666  }
3667  if (t <= 0.38 + 0.01)
3668  {
3669  return 0.7230071250000001;
3670  }
3671  if (t <= 0.39 + 0.01)
3672  {
3673  return 0.5555071250000001;
3674  }
3675  if (t <= 0.40 + 0.01)
3676  {
3677  return 0.3694071249999999;
3678  }
3679  if (t <= 0.41 + 0.01)
3680  {
3681  return 0.203607125 ;
3682  }
3683  if (t <= 0.42 + 0.01)
3684  {
3685  return 0.117107125 ;
3686  }
3687  if (t <= 0.43 + 0.01)
3688  {
3689  return 0.154707125 ;
3690  }
3691  if (t <= 0.44 + 0.01)
3692  {
3693  return 0.3019071249999999;
3694  }
3695  if (t <= 0.45 + 0.01)
3696  {
3697  return 0.4833071250000001;
3698  }
3699  if (t <= 0.46 + 0.01)
3700  {
3701  return 0.6279071250000001;
3702  }
3703  if (t <= 0.47 + 0.01)
3704  {
3705  return 0.7323071250000001;
3706  }
3707  if (t <= 0.48 + 0.01)
3708  {
3709  return 0.8435071250000001;
3710  }
3711  if (t <= 0.49 + 0.01)
3712  {
3713  return 0.989007125 ;
3714  }
3715  if (t <= 0.50 + 0.01)
3716  {
3717  return 1.137007125 ;
3718  }
3719  if (t <= 0.51 + 0.01)
3720  {
3721  return 1.223007125 ;
3722  }
3723  if (t <= 0.52 + 0.01)
3724  {
3725  return 1.206007125 ;
3726  }
3727  if (t <= 0.53 + 0.01)
3728  {
3729  return 1.096007125 ;
3730  }
3731  if (t <= 0.54 + 0.01)
3732  {
3733  return 0.9400071249999998;
3734  }
3735  if (t <= 0.55 + 0.01)
3736  {
3737  return 0.7924071250000001;
3738  }
3739  if (t <= 0.56 + 0.01)
3740  {
3741  return 0.6857071250000001;
3742  }
3743  if (t <= 0.57 + 0.01)
3744  {
3745  return 0.6278071250000001;
3746  }
3747  if (t <= 0.58 + 0.01)
3748  {
3749  return 0.6083071250000001;
3750  }
3751  if (t <= 0.59 + 0.01)
3752  {
3753  return 0.610607125 ;
3754  }
3755  if (t <= 0.60 + 0.01)
3756  {
3757  return 0.620007125 ;
3758  }
3759  if (t <= 0.61 + 0.01)
3760  {
3761  return 0.6279071250000001;
3762  }
3763  if (t <= 0.62 + 0.01)
3764  {
3765  return 0.6320071250000001;
3766  }
3767  if (t <= 0.63 + 0.01)
3768  {
3769  return 0.6343071250000001;
3770  }
3771  if (t <= 0.64 + 0.01)
3772  {
3773  return 0.6385071250000001;
3774  }
3775  if (t <= 0.65 + 0.01)
3776  {
3777  return 0.6476071250000001;
3778  }
3779  if (t <= 0.66 + 0.01)
3780  {
3781  return 0.6620071250000001;
3782  }
3783  if (t <= 0.67 + 0.01)
3784  {
3785  return 0.679107125 ;
3786  }
3787  if (t <= 0.68 + 0.01)
3788  {
3789  return 0.6948071250000001;
3790  }
3791  if (t <= 0.69 + 0.01)
3792  {
3793  return 0.7049071250000001;
3794  }
3795  if (t <= 0.70 + 0.01)
3796  {
3797  return 0.7071071250000001;
3798  }
3799  if (t <= 0.71 + 0.01)
3800  {
3801  return 0.7021071250000001;
3802  }
3803  if (t <= 0.72 + 0.01)
3804  {
3805  return 0.6926071250000001;
3806  }
3807  if (t <= 0.73 + 0.01)
3808  {
3809  return 0.6825071250000001;
3810  }
3811  if (t <= 0.74 + 0.01)
3812  {
3813  return 0.675007125 ;
3814  }
3815  if (t <= 0.75 + 0.01)
3816  {
3817  return 0.671207125 ;
3818  }
3819  if (t <= 0.76 + 0.01)
3820  {
3821  return 0.6702071250000001;
3822  }
3823  if (t <= 0.77 + 0.01)
3824  {
3825  return 0.6702071250000001;
3826  }
3827  if (t <= 0.78 + 0.01)
3828  {
3829  return 0.669507125 ;
3830  }
3831  if (t <= 0.79 + 0.01)
3832  {
3833  return 0.6674071250000001;
3834  }
3835  RETURN_UNDEFINED;
3836 }
3837 
3838 Real u2 (Real /*t*/, const Real& /*x*/, const Real& /*y*/, const Real& /*z*/, const ID& /*i*/)
3839 {
3840  return -1e4;
3841 }
3842 
3843 }
fZero
Real fZero(const Real &, const Real &, const Real &, const Real &, const ID &)
Definition: cavity_stokes.cpp:87
LifeV::d0
Real d0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3017
LifeV::aortaFluxIn
Real aortaFluxIn(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:20
LifeV::aortaFlux6_
Real aortaFlux6_(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3529
LifeV::E
Real E(const Real &, const Real &, const Real &, const Real &, const ID &)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2989
RETURN_UNDEFINED
#define RETURN_UNDEFINED
Definition: LifeAssert.hpp:70
LifeV::aortaFlux6
Real aortaFlux6(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:1331
LifeV::linearFlux8
Real linearFlux8(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3154
LifeV::p0
Real p0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2983
LifeV::linearFlux6_
Real linearFlux6_(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3192
LifeV::linearPress2
Real linearPress2(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3060
ETCurrentFE::updateInverseJacobian
void updateInverseJacobian(const UInt &iQuadPt)
Definition: ETCurrentFE.cpp:405
LifeV::aortaFlux3_
Real aortaFlux3_(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3204
LifeV::w0
Real w0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3037
ERROR_MSG
#define ERROR_MSG(A)
Definition: LifeAssert.hpp:69
LifeV::ID
uint32_type ID
IDs.
Definition: LifeV.hpp:194
LifeV::linearFlux3_
Real linearFlux3_(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3179
LifeV::linearFlux7
Real linearFlux7(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3141
LifeV::f
Real f(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2962
LifeV::aortaFlux3
Real aortaFlux3(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:681
LifeV::linearFlux5
Real linearFlux5(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3114
LifeV::aortaPhisPress
Real aortaPhisPress(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:345
LifeV::Real
double Real
Generic real data.
Definition: LifeV.hpp:175
LifeV::abdominalAorta
Real abdominalAorta(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/challenge_VPH/ud_functions.cpp:20
LifeV::u1
Real u1(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2967
LifeV::linearFlux9
Real linearFlux9(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3166
LifeV::aortaFlux8
Real aortaFlux8(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:1981
LifeV::aortaFlux4
Real aortaFlux4(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2631
LifeV::u0
Real u0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2978
LifeV::linearFlux3
Real linearFlux3(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3073
LifeV::aortaFlux9
Real aortaFlux9(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:2306
LifeV::aortaFlux7
Real aortaFlux7(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:1656
LifeV::u2
Real u2(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3854
LifeV::linearFlux6
Real linearFlux6(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3128
LifeV::aortaFlux5
Real aortaFlux5(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:1006
LifeV::linearFlux4
Real linearFlux4(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3100
LifeV::linearFluxIn
Real linearFluxIn(Real t, const Real &, const Real &, const Real &, const ID &i)
Definition: fsi/examples/application_aortaFSI/ud_functions.cpp:3086