PageRenderTime 61ms CodeModel.GetById 20ms RepoModel.GetById 0ms app.codeStats 1ms

/trunk/DOC/html/init__cond_8cc-source.html

https://bitbucket.org/mchandra/tarang-mpi-old
HTML | 1885 lines | 1884 code | 0 blank | 1 comment | 0 complexity | ef3171e1de4180d6bfb409ecb5e3b2da MD5 | raw file
Possible License(s): GPL-3.0

Large files files are truncated, but you can click here to view the full file

  1. <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
  2. <html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
  3. <title>TARANG-MPI: init_cond.cc Source File</title>
  4. <link href="doxygen.css" rel="stylesheet" type="text/css">
  5. <link href="tabs.css" rel="stylesheet" type="text/css">
  6. </head><body>
  7. <!-- Generated by Doxygen 1.5.6 -->
  8. <div class="navigation" id="top">
  9. <div class="tabs">
  10. <ul>
  11. <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
  12. <li><a href="annotated.html"><span>Data&nbsp;Structures</span></a></li>
  13. <li class="current"><a href="files.html"><span>Files</span></a></li>
  14. </ul>
  15. </div>
  16. <h1>init_cond.cc</h1><a href="init__cond_8cc.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
  17. <a name="l00002"></a>00002
  18. <a name="l00003"></a>00003 <span class="comment">// SPECTRAL Version 1.0</span>
  19. <a name="l00004"></a>00004 <span class="comment">// Date: 16 Sept. 2007</span>
  20. <a name="l00005"></a>00005 <span class="comment">// Author: M. K. Verma</span>
  21. <a name="l00006"></a>00006 <span class="comment">// Filename: init_cond.cc</span>
  22. <a name="l00007"></a>00007
  23. <a name="l00008"></a>00008
  24. <a name="l00009"></a>00009 <span class="comment">/*====================================================================================</span>
  25. <a name="l00010"></a>00010 <span class="comment"></span>
  26. <a name="l00011"></a>00011 <span class="comment"> Sets up initial conditions of the Vector and Scalar field</span>
  27. <a name="l00012"></a>00012 <span class="comment"></span>
  28. <a name="l00013"></a>00013 <span class="comment">======================================================================================*/</span>
  29. <a name="l00014"></a>00014
  30. <a name="l00015"></a>00015 <span class="preprocessor">#include "../IncFluid.h"</span>
  31. <a name="l00016"></a>00016
  32. <a name="l00017"></a>00017 <span class="keyword">extern</span> Uniform&lt;DP&gt; <a class="code" href="useful__fns_8cc.html#a8ff9eceb2e85b30b25e4cc3f18c34b3">SPECrand</a>;
  33. <a name="l00018"></a>00018
  34. <a name="l00019"></a><a class="code" href="init__cond_8cc.html#b63ee86ed8efc2de66921c1ae96c7696">00019</a> <span class="keyword">const</span> <span class="keywordtype">int</span> <a class="code" href="init__cond_8cc.html#b63ee86ed8efc2de66921c1ae96c7696">IC_MAX_NO_WAVENO</a> = 40;
  35. <a name="l00020"></a>00020
  36. <a name="l00021"></a>00021 <span class="keyword">extern</span> <span class="keywordtype">int</span> <a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a>; <span class="comment">// My process id</span>
  37. <a name="l00022"></a>00022 <span class="keyword">extern</span> <span class="keywordtype">int</span> <a class="code" href="basis__basicfn_8cc.html#31bfe068640b411bea3df0d69177ded5">numprocs</a>; <span class="comment">// No of processors</span>
  38. <a name="l00023"></a>00023 <span class="keyword">extern</span> <span class="keyword">const</span> <span class="keywordtype">int</span> <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>; <span class="comment">// Id of master proc</span>
  39. <a name="l00024"></a>00024 <span class="keyword">extern</span> ptrdiff_t <a class="code" href="basis__basicfn_8cc.html#5b280d17c0a165711d64dda77338d027">local_N1</a>, <a class="code" href="basis__basicfn_8cc.html#5c5df0dd795318288fed5e089341a3fc">local_N1_start</a>; <span class="comment">// N1 size and start of l1 in the currentproc</span>
  40. <a name="l00025"></a>00025 <span class="keyword">extern</span> Array&lt;complex&lt;double&gt;,3&gt; *<a class="code" href="init__cond_8cc.html#86d7b76b5f2b9e19340a08476fae4ee4">container</a>;
  41. <a name="l00026"></a>00026
  42. <a name="l00027"></a>00027
  43. <a name="l00028"></a>00028
  44. <a name="l00029"></a><a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">00029</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">IncFluid::Input_prefix</a>(ifstream&amp; field_in_file)
  45. <a name="l00030"></a>00030 {
  46. <a name="l00031"></a>00031 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>)
  47. <a name="l00032"></a>00032 {
  48. <a name="l00033"></a>00033 <span class="keywordflow">if</span> (! field_in_file.is_open())
  49. <a name="l00034"></a>00034 {
  50. <a name="l00035"></a>00035 cout &lt;&lt; <span class="stringliteral">"Unable to open field_in_file: Exiting Program: Exiting Program "</span> &lt;&lt; endl;
  51. <a name="l00036"></a>00036 exit(1);
  52. <a name="l00037"></a>00037 }
  53. <a name="l00038"></a>00038 <span class="keywordflow">else</span>
  54. <a name="l00039"></a>00039 cout &lt;&lt; <span class="stringliteral">"Reading field configurations from field_in_file "</span> &lt;&lt; endl ;
  55. <a name="l00040"></a>00040 }
  56. <a name="l00041"></a>00041 }
  57. <a name="l00042"></a>00042
  58. <a name="l00043"></a>00043
  59. <a name="l00044"></a>00044 <span class="comment">//</span>
  60. <a name="l00045"></a>00045 <span class="comment">//</span>
  61. <a name="l00046"></a><a class="code" href="classIncFluid.html#be3930f267a9faf356e1d737c66bf093">00046</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#be3930f267a9faf356e1d737c66bf093">IncFluid::Read_init_cond</a>(<span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  62. <a name="l00047"></a>00047 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  63. <a name="l00048"></a>00048 {
  64. <a name="l00049"></a>00049 <span class="keywordflow">switch</span> (field_input_proc) {
  65. <a name="l00050"></a>00050
  66. <a name="l00051"></a>00051 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  67. <a name="l00052"></a>00052
  68. <a name="l00053"></a>00053 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file with Nreduced D</span>
  69. <a name="l00054"></a>00054
  70. <a name="l00055"></a>00055 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, (Vy:3D),Theta</span>
  71. <a name="l00056"></a>00056
  72. <a name="l00057"></a>00057 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, vorticity, Theta</span>
  73. <a name="l00058"></a>00058
  74. <a name="l00059"></a>00059 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy spectrum</span>
  75. <a name="l00060"></a>00060
  76. <a name="l00061"></a>00061 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  77. <a name="l00062"></a>00062
  78. <a name="l00063"></a>00063 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(init_cond_para); <span class="keywordflow">break</span>;
  79. <a name="l00064"></a>00064
  80. <a name="l00065"></a>00065 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(init_cond_para); <span class="keywordflow">break</span>;
  81. <a name="l00066"></a>00066
  82. <a name="l00067"></a>00067 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  83. <a name="l00068"></a>00068
  84. <a name="l00069"></a>00069 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(N_in_reduced, no_Nreduced_procs, nos_input_field_mode); <span class="keywordflow">break</span>;
  85. <a name="l00070"></a>00070 <span class="comment">// read from field_in_file with Nreduced D: Not working very well for Nreduced.</span>
  86. <a name="l00071"></a>00071 }
  87. <a name="l00072"></a>00072 }
  88. <a name="l00073"></a>00073
  89. <a name="l00074"></a>00074 <span class="comment">//</span>
  90. <a name="l00075"></a>00075 <span class="comment">//</span>
  91. <a name="l00076"></a>00076
  92. <a name="l00077"></a><a class="code" href="classIncFluid.html#574d7bfbeff3e9057f07e4d48b795e56">00077</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#be3930f267a9faf356e1d737c66bf093">IncFluid::Read_init_cond</a>(<a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  93. <a name="l00078"></a>00078 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  94. <a name="l00079"></a>00079 {
  95. <a name="l00080"></a>00080 <span class="keywordflow">switch</span> (field_input_proc) {
  96. <a name="l00081"></a>00081
  97. <a name="l00082"></a>00082 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  98. <a name="l00083"></a>00083
  99. <a name="l00084"></a>00084 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(T, N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>;
  100. <a name="l00085"></a>00085 <span class="comment">// read from field_in_file with Nreduced D</span>
  101. <a name="l00086"></a>00086
  102. <a name="l00087"></a>00087 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, (Vy:3D),Theta</span>
  103. <a name="l00088"></a>00088
  104. <a name="l00089"></a>00089 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, vorticity, Theta</span>
  105. <a name="l00090"></a>00090
  106. <a name="l00091"></a>00091 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy spectrum</span>
  107. <a name="l00092"></a>00092
  108. <a name="l00093"></a>00093 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  109. <a name="l00094"></a>00094
  110. <a name="l00095"></a>00095 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize only V field</span>
  111. <a name="l00096"></a>00096
  112. <a name="l00097"></a>00097 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize only V field</span>
  113. <a name="l00098"></a>00098
  114. <a name="l00099"></a>00099 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  115. <a name="l00100"></a>00100
  116. <a name="l00101"></a>00101 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(T, N_in_reduced, no_Nreduced_procs, nos_input_field_mode); <span class="keywordflow">break</span>;
  117. <a name="l00102"></a>00102 <span class="comment">// read from field_in_file with Nreduced D</span>
  118. <a name="l00103"></a>00103 }
  119. <a name="l00104"></a>00104
  120. <a name="l00105"></a>00105 }
  121. <a name="l00106"></a>00106
  122. <a name="l00107"></a>00107 <span class="comment">//</span>
  123. <a name="l00108"></a>00108 <span class="comment">//</span>
  124. <a name="l00109"></a>00109
  125. <a name="l00110"></a><a class="code" href="classIncFluid.html#06bbc8a83da16ce339f8fb8e0e32dbf1">00110</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#be3930f267a9faf356e1d737c66bf093">IncFluid::Read_init_cond</a>(<a class="code" href="classIncVF.html">IncVF</a>&amp; W, <span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  126. <a name="l00111"></a>00111 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  127. <a name="l00112"></a>00112 {
  128. <a name="l00113"></a>00113 <span class="keywordflow">switch</span> (field_input_proc) {
  129. <a name="l00114"></a>00114
  130. <a name="l00115"></a>00115 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(W, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  131. <a name="l00116"></a>00116
  132. <a name="l00117"></a>00117 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(W, N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>;
  133. <a name="l00118"></a>00118 <span class="comment">// read from field_in_file with Nreduced D</span>
  134. <a name="l00119"></a>00119
  135. <a name="l00120"></a>00120 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(W); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, (Vy:3D),Theta</span>
  136. <a name="l00121"></a>00121
  137. <a name="l00122"></a>00122 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(W); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, vorticity, Theta</span>
  138. <a name="l00123"></a>00123
  139. <a name="l00124"></a>00124 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(W, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy spectrum</span>
  140. <a name="l00125"></a>00125
  141. <a name="l00126"></a>00126 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(W, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  142. <a name="l00127"></a>00127
  143. <a name="l00128"></a>00128 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(W, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V, W field</span>
  144. <a name="l00129"></a>00129
  145. <a name="l00130"></a>00130 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(W, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V,W field</span>
  146. <a name="l00131"></a>00131
  147. <a name="l00132"></a>00132 <span class="keywordflow">case</span> (11) : <a class="code" href="classIncFluid.html#6ad79b117329b4d33d6d0d2a48615b53">Init_cond_DYNAMO_SIX_MODE</a>(W, init_cond_para); <span class="keywordflow">break</span>;
  148. <a name="l00133"></a>00133
  149. <a name="l00134"></a>00134 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(W, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  150. <a name="l00135"></a>00135
  151. <a name="l00136"></a>00136 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(W, N_in_reduced, no_Nreduced_procs, nos_input_field_mode); <span class="keywordflow">break</span>;
  152. <a name="l00137"></a>00137 <span class="comment">// read from field_in_file with Nreduced D</span>
  153. <a name="l00138"></a>00138 }
  154. <a name="l00139"></a>00139
  155. <a name="l00140"></a>00140 }
  156. <a name="l00141"></a>00141
  157. <a name="l00142"></a>00142
  158. <a name="l00143"></a>00143 <span class="comment">//</span>
  159. <a name="l00144"></a>00144 <span class="comment">//</span>
  160. <a name="l00145"></a>00145
  161. <a name="l00146"></a><a class="code" href="classIncFluid.html#eebccc0335fd5ffab6aedfbb1dad5eaa">00146</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#be3930f267a9faf356e1d737c66bf093">IncFluid::Read_init_cond</a>(<a class="code" href="classIncVF.html">IncVF</a>&amp; W, <a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  162. <a name="l00147"></a>00147 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  163. <a name="l00148"></a>00148 {
  164. <a name="l00149"></a>00149 <span class="keywordflow">switch</span> (field_input_proc) {
  165. <a name="l00150"></a>00150
  166. <a name="l00151"></a>00151 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(W, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  167. <a name="l00152"></a>00152
  168. <a name="l00153"></a>00153 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(W, T, N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>;
  169. <a name="l00154"></a>00154 <span class="comment">// read from field_in_file with Nreduced D</span>
  170. <a name="l00155"></a>00155
  171. <a name="l00156"></a>00156 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(W, T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, (Vy:3D),Theta</span>
  172. <a name="l00157"></a>00157
  173. <a name="l00158"></a>00158 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(W, T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, vorticity, Theta</span>
  174. <a name="l00159"></a>00159
  175. <a name="l00160"></a>00160 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy spectrum</span>
  176. <a name="l00161"></a>00161
  177. <a name="l00162"></a>00162 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  178. <a name="l00163"></a>00163
  179. <a name="l00164"></a>00164 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V, W field</span>
  180. <a name="l00165"></a>00165
  181. <a name="l00166"></a>00166 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V, W field</span>
  182. <a name="l00167"></a>00167
  183. <a name="l00168"></a>00168 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(W, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  184. <a name="l00169"></a>00169
  185. <a name="l00170"></a>00170 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(W, T, N_in_reduced, no_Nreduced_procs, nos_input_field_mode); <span class="keywordflow">break</span>;
  186. <a name="l00171"></a>00171 <span class="comment">// read from field_in_file with Nreduced D</span>
  187. <a name="l00172"></a>00172 }
  188. <a name="l00173"></a>00173
  189. <a name="l00174"></a>00174 }
  190. <a name="l00175"></a>00175
  191. <a name="l00176"></a>00176
  192. <a name="l00177"></a>00177 <span class="comment">//</span>
  193. <a name="l00178"></a>00178 <span class="comment">//</span>
  194. <a name="l00179"></a>00179
  195. <a name="l00180"></a><a class="code" href="classIncFluid.html#cb5b5dec631611fcecc0d5cdddf5868d">00180</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#cb5b5dec631611fcecc0d5cdddf5868d">IncFluid::Read_init_cond_RB</a>(<span class="keywordtype">string</span> Pr_switch, <a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  196. <a name="l00181"></a>00181 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  197. <a name="l00182"></a>00182 {
  198. <a name="l00183"></a>00183 <span class="keywordflow">switch</span> (field_input_proc) {
  199. <a name="l00184"></a>00184
  200. <a name="l00185"></a>00185 <span class="keywordflow">case</span> (0) : <a class="code" href="classIncFluid.html#55d61b9682d43cb8ce9aafcbb8f9e9db">Init_cond_RB_Lorenz</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// first 3 entries are w101, th101, th200</span>
  201. <a name="l00186"></a>00186
  202. <a name="l00187"></a>00187 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(Pr_switch, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  203. <a name="l00188"></a>00188
  204. <a name="l00189"></a>00189 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(Pr_switch, T, N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>;
  205. <a name="l00190"></a>00190 <span class="comment">// read from field_in_file with Nreduced D</span>
  206. <a name="l00191"></a>00191
  207. <a name="l00192"></a>00192 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(Pr_switch, T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, (Vy:3D),Theta</span>
  208. <a name="l00193"></a>00193
  209. <a name="l00194"></a>00194 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(Pr_switch, T); <span class="keywordflow">break</span>; <span class="comment">// Modes - ki, Vx, vorticity, Theta</span>
  210. <a name="l00195"></a>00195
  211. <a name="l00196"></a>00196 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(Pr_switch, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy spectrum</span>
  212. <a name="l00197"></a>00197
  213. <a name="l00198"></a>00198 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(Pr_switch, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  214. <a name="l00199"></a>00199
  215. <a name="l00200"></a>00200 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(Pr_switch, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize only V field</span>
  216. <a name="l00201"></a>00201
  217. <a name="l00202"></a>00202 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(Pr_switch, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize only V field</span>
  218. <a name="l00203"></a>00203
  219. <a name="l00204"></a>00204 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(Pr_switch, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  220. <a name="l00205"></a>00205
  221. <a name="l00206"></a>00206 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(Pr_switch, T, N_in_reduced, no_Nreduced_procs, nos_input_field_mode); <span class="keywordflow">break</span>;
  222. <a name="l00207"></a>00207 <span class="comment">// read from field_in_file with Nreduced D</span>
  223. <a name="l00208"></a>00208 }
  224. <a name="l00209"></a>00209 }
  225. <a name="l00210"></a>00210
  226. <a name="l00211"></a><a class="code" href="classIncFluid.html#c94bd056e9c2b2c1288cdf9dc95ce479">00211</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#cb5b5dec631611fcecc0d5cdddf5868d">IncFluid::Read_init_cond_RB</a>(<span class="keywordtype">string</span> Pr_switch, <a class="code" href="classIncVF.html">IncVF</a>&amp; W, <a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">int</span> field_input_proc, Array&lt;DP,1&gt; init_cond_para,
  227. <a name="l00212"></a>00212 <span class="keywordtype">int</span> N_in_reduced[], <span class="keywordtype">int</span> no_Nreduced_procs, <span class="keywordtype">string</span> nos_input_field_mode)
  228. <a name="l00213"></a>00213 {
  229. <a name="l00214"></a>00214 <span class="keywordflow">switch</span> (field_input_proc) {
  230. <a name="l00215"></a>00215
  231. <a name="l00216"></a>00216 <span class="keywordflow">case</span> (0) : <a class="code" href="classIncFluid.html#55d61b9682d43cb8ce9aafcbb8f9e9db">Init_cond_RB_Lorenz</a>(T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// first 3 entries are w101, th101, th200 </span>
  232. <a name="l00217"></a>00217
  233. <a name="l00218"></a>00218 <span class="keywordflow">case</span> (1) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(Pr_switch, W, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  234. <a name="l00219"></a>00219
  235. <a name="l00220"></a>00220 <span class="keywordflow">case</span> (2) : <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">Init_cond</a>(Pr_switch, W, T, N_in_reduced, nos_input_field_mode); <span class="keywordflow">break</span>;
  236. <a name="l00221"></a>00221 <span class="comment">// read from field_in_file with Nreduced D</span>
  237. <a name="l00222"></a>00222
  238. <a name="l00223"></a>00223 <span class="keywordflow">case</span> (3) : <a class="code" href="classIncFluid.html#b81030f98a7ccac86253e6bdf2b2569d">Init_cond_modes_SIMPLE</a>(Pr_switch, W, T); <span class="keywordflow">break</span>; <span class="comment">// Triad - ki, Vx, (Vy:3D),Theta</span>
  239. <a name="l00224"></a>00224
  240. <a name="l00225"></a>00225 <span class="keywordflow">case</span> (4) : <a class="code" href="classIncFluid.html#212c6db89a6a788b7022e3911c2d5ffb">Init_cond_modes_VORTICITY</a>(Pr_switch, W, T); <span class="keywordflow">break</span>; <span class="comment">// Triad - ki, Vx, vorticity, Theta</span>
  241. <a name="l00226"></a>00226
  242. <a name="l00227"></a>00227 <span class="keywordflow">case</span> (5) : <a class="code" href="classIncFluid.html#f051aceda92b40d5acb1e643da67b3da">Init_cond_energy_spectrum</a>(Pr_switch, W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given spectrum</span>
  243. <a name="l00228"></a>00228
  244. <a name="l00229"></a>00229 <span class="keywordflow">case</span> (6) : <a class="code" href="classIncFluid.html#902e6e061b97dff5aaa085be10004d65">Init_cond_energy_helicity_spectrum</a>(Pr_switch, W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// given energy and hel spectrum</span>
  245. <a name="l00230"></a>00230
  246. <a name="l00231"></a>00231 <span class="keywordflow">case</span> (7) : <a class="code" href="classIncFluid.html#e5dd5720bc72f0ea698764622b55af3e">Init_cond_Taylor_Green</a>(Pr_switch, W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V, W field</span>
  247. <a name="l00232"></a>00232
  248. <a name="l00233"></a>00233 <span class="keywordflow">case</span> (8) : <a class="code" href="classIncFluid.html#3d791dbda9966f3268e01785b046ab65">Init_cond_ABC</a>(Pr_switch, W, T, init_cond_para); <span class="keywordflow">break</span>; <span class="comment">// initialize V, W field</span>
  249. <a name="l00234"></a>00234
  250. <a name="l00235"></a>00235 <span class="keywordflow">case</span> (101) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(Pr_switch, W, T, nos_input_field_mode); <span class="keywordflow">break</span>; <span class="comment">// read from field_in_file</span>
  251. <a name="l00236"></a>00236
  252. <a name="l00237"></a>00237 <span class="keywordflow">case</span> (102) : <a class="code" href="classIncFluid.html#2bdeccfa4f5b792b6a6caac6bbc6a4ed">Init_cond_split_arrays</a>(Pr_switch, W, T, N_in_reduced, no_Nreduced_procs, nos_input_field_mode);
  253. <a name="l00238"></a>00238 <span class="keywordflow">break</span>;
  254. <a name="l00239"></a>00239 <span class="comment">// read from field_in_file with Nreduced D</span>
  255. <a name="l00240"></a>00240 }
  256. <a name="l00241"></a>00241
  257. <a name="l00242"></a>00242 }
  258. <a name="l00243"></a>00243
  259. <a name="l00244"></a>00244
  260. <a name="l00245"></a>00245 <span class="comment">/*====================================================================================</span>
  261. <a name="l00246"></a>00246 <span class="comment"></span>
  262. <a name="l00247"></a>00247 <span class="comment"> Input from a file: field_in_file </span>
  263. <a name="l00248"></a>00248 <span class="comment"> Arrays are split </span>
  264. <a name="l00249"></a>00249 <span class="comment"></span>
  265. <a name="l00250"></a>00250 <span class="comment">======================================================================================*/</span>
  266. <a name="l00251"></a>00251
  267. <a name="l00252"></a><a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">00252</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<span class="keywordtype">string</span> nos_input_field_mode)
  268. <a name="l00253"></a>00253 {
  269. <a name="l00254"></a>00254 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  270. <a name="l00255"></a>00255
  271. <a name="l00256"></a>00256 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>); <span class="comment">// *VF_temp is termporary array useful for input</span>
  272. <a name="l00257"></a>00257
  273. <a name="l00258"></a>00258 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>)
  274. <a name="l00259"></a>00259 cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  275. <a name="l00260"></a>00260 }
  276. <a name="l00261"></a>00261
  277. <a name="l00262"></a>00262
  278. <a name="l00263"></a><a class="code" href="classIncFluid.html#e4819a6a9bce7830c119ba326d4262b5">00263</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">string</span> nos_input_field_mode)
  279. <a name="l00264"></a>00264 {
  280. <a name="l00265"></a>00265 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  281. <a name="l00266"></a>00266
  282. <a name="l00267"></a>00267 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  283. <a name="l00268"></a>00268 T.<a class="code" href="classCSF.html#d24b8acaffa60136ec0e3d4376cec1c8">CS_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  284. <a name="l00269"></a>00269
  285. <a name="l00270"></a>00270 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>)
  286. <a name="l00271"></a>00271 cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  287. <a name="l00272"></a>00272 }
  288. <a name="l00273"></a>00273
  289. <a name="l00274"></a>00274 <span class="comment">//</span>
  290. <a name="l00275"></a>00275 <span class="comment">//</span>
  291. <a name="l00276"></a>00276
  292. <a name="l00277"></a><a class="code" href="classIncFluid.html#dd6e7aada9433ac059c165516f755b29">00277</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<a class="code" href="classIncVF.html">IncVF</a>&amp; W, <span class="keywordtype">string</span> nos_input_field_mode)
  293. <a name="l00278"></a>00278 {
  294. <a name="l00279"></a>00279 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  295. <a name="l00280"></a>00280
  296. <a name="l00281"></a>00281 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  297. <a name="l00282"></a>00282 W.<a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  298. <a name="l00283"></a>00283
  299. <a name="l00284"></a>00284 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>)
  300. <a name="l00285"></a>00285 cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  301. <a name="l00286"></a>00286 }
  302. <a name="l00287"></a>00287
  303. <a name="l00288"></a>00288 <span class="comment">//</span>
  304. <a name="l00289"></a>00289
  305. <a name="l00290"></a><a class="code" href="classIncFluid.html#89e2c395efe4c200aea501495d7a7159">00290</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<a class="code" href="classIncVF.html">IncVF</a>&amp; W, <a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">string</span> nos_input_field_mode)
  306. <a name="l00291"></a>00291 {
  307. <a name="l00292"></a>00292 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  308. <a name="l00293"></a>00293
  309. <a name="l00294"></a>00294 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  310. <a name="l00295"></a>00295 W.<a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  311. <a name="l00296"></a>00296 T.<a class="code" href="classCSF.html#d24b8acaffa60136ec0e3d4376cec1c8">CS_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  312. <a name="l00297"></a>00297
  313. <a name="l00298"></a>00298 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>)
  314. <a name="l00299"></a>00299 cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  315. <a name="l00300"></a>00300 }
  316. <a name="l00301"></a>00301 <span class="comment">/*==================================================================================== </span>
  317. <a name="l00302"></a>00302 <span class="comment"> Input from a file: field_in_file(Nreduced[]) </span>
  318. <a name="l00303"></a>00303 <span class="comment">======================================================================================*/</span>
  319. <a name="l00304"></a>00304
  320. <a name="l00305"></a>00305 <span class="comment">// Fluid</span>
  321. <a name="l00306"></a><a class="code" href="classIncFluid.html#ab174c94b72bec24036c357b749eb2e4">00306</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<span class="keywordtype">int</span> Nreduced[], <span class="keywordtype">string</span> nos_input_field_mode)
  322. <a name="l00307"></a>00307 {
  323. <a name="l00308"></a>00308 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  324. <a name="l00309"></a>00309
  325. <a name="l00310"></a>00310 (*V1) = 0.0; (*V2) = 0.0; (*V3) = 0.0;
  326. <a name="l00311"></a>00311
  327. <a name="l00312"></a>00312 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, Nreduced, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  328. <a name="l00313"></a>00313
  329. <a name="l00314"></a>00314 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>) cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  330. <a name="l00315"></a>00315 }
  331. <a name="l00316"></a>00316
  332. <a name="l00317"></a>00317
  333. <a name="l00318"></a>00318 <span class="comment">// Passive scalar + RB convection</span>
  334. <a name="l00319"></a>00319
  335. <a name="l00320"></a><a class="code" href="classIncFluid.html#00c57aede998adea7dd676f67b4d150c">00320</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<a class="code" href="classIncSF.html">IncSF</a>&amp; T, <span class="keywordtype">int</span> Nreduced[], <span class="keywordtype">string</span> nos_input_field_mode)
  336. <a name="l00321"></a>00321 {
  337. <a name="l00322"></a>00322 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  338. <a name="l00323"></a>00323
  339. <a name="l00324"></a>00324
  340. <a name="l00325"></a>00325 (*V1) = 0.0; (*V2) = 0.0; (*V3) = 0.0;
  341. <a name="l00326"></a>00326 (*T.<a class="code" href="classCSF.html#3b75b0fe8368d57759f1a67a050e6367">F</a>) = 0.0;
  342. <a name="l00327"></a>00327
  343. <a name="l00328"></a>00328 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, Nreduced, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  344. <a name="l00329"></a>00329 T.<a class="code" href="classCSF.html#d24b8acaffa60136ec0e3d4376cec1c8">CS_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, Nreduced, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  345. <a name="l00330"></a>00330
  346. <a name="l00331"></a>00331 <span class="keywordflow">if</span> (<a class="code" href="basis__basicfn_8cc.html#f9650afe4b09491cf2583300c86c81fb">my_id</a> == <a class="code" href="basis__basicfn_8cc.html#effbd5f4f8929fc0d3c5fc83e500abbb">master_id</a>) cout &lt;&lt; <span class="stringliteral">"Reading of field configurations ended successfully"</span> &lt;&lt; endl;
  347. <a name="l00332"></a>00332 }
  348. <a name="l00333"></a>00333
  349. <a name="l00334"></a>00334 <span class="comment">// MHD</span>
  350. <a name="l00335"></a>00335
  351. <a name="l00336"></a><a class="code" href="classIncFluid.html#0a7fbe275a5dd7d0209ca5e153bb5450">00336</a> <span class="keywordtype">void</span> <a class="code" href="classIncFluid.html#2cc82cf45ab178363a3883dd13c77f40">IncFluid::Init_cond</a>(<a class="code" href="classIncVF.html">IncVF</a>&amp; W, <span class="keywordtype">int</span> Nreduced[], <span class="keywordtype">string</span> nos_input_field_mode)
  352. <a name="l00337"></a>00337 {
  353. <a name="l00338"></a>00338 <a class="code" href="classIncFluid.html#1145d2815506aae853374229e13eddda">Input_prefix</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>);
  354. <a name="l00339"></a>00339
  355. <a name="l00340"></a>00340 (*V1) = 0.0; (*V2) = 0.0; (*V3) = 0.0;
  356. <a name="l00341"></a>00341 (*W.<a class="code" href="classCVF.html#77e2105a2d7d2a35cc65987f28951841">V1</a>) = 0.0; (*W.<a class="code" href="classCVF.html#baeccade011a2b615777dc89f2710929">V2</a>) = 0.0; (*W.<a class="code" href="classCVF.html#9f69a238438f20339b5c7e05faa3d1f6">V3</a>) = 0.0;
  357. <a name="l00342"></a>00342
  358. <a name="l00343"></a>00343 <a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="code" href="classIncFluid.html#6a2ec36a5aedecf16533d7926ddf855b">field_in_file</a>, Nreduced, *<a class="code" href="classIncVF.html#4c9070634e826d4ac9aba1dd9ed3444e">VF_temp</a>);
  359. <a name="l00344"></a>00344 W.<a class="code" href="classCVF.html#77508ab009c7f35dc391693901caa50f">CV_input</a>(<a class="c…

Large files files are truncated, but you can click here to view the full file