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