 |
tode.cpp.html - numeric - C++ library with numerical algorithms |
 |
git clone git://src.adamsgaard.dk/numeric |
 |
Log |
|
Files |
|
Refs |
|
LICENSE |
|
--- |
|
tode.cpp.html (17180B) |
|
--- |
|
1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/… |
|
2 <html> |
|
3 <head> |
|
4 <meta http-equiv="content-type" content="text/html; charset=UTF-8"> |
|
5 <title>~/code/numeric/exam/ode.cpp.html</title> |
|
6 <meta name="Generator" content="Vim/7.4"> |
|
7 <meta name="plugin-version" content="vim7.4_v1"> |
|
8 <meta name="syntax" content="cpp"> |
|
9 <meta name="settings" content="number_lines,use_css,pre_wrap,no_foldcolu… |
|
10 <meta name="colorscheme" content="desert"> |
|
11 <style type="text/css"> |
|
12 <!-- |
|
13 pre { white-space: pre-wrap; font-family: monospace; color: #ffffff; bac… |
|
14 body { font-family: monospace; color: #ffffff; background-color: #000000… |
|
15 * { font-size: 1em; } |
|
16 .Type { color: #008000; } |
|
17 .Statement { color: #804000; } |
|
18 .LineNr { color: #804000; } |
|
19 .PreProc { color: #c000c0; } |
|
20 .Constant { color: #af5f00; } |
|
21 .Special { color: #c000c0; } |
|
22 .Comment { color: #008080; } |
|
23 --> |
|
24 </style> |
|
25 |
|
26 <script type='text/javascript'> |
|
27 <!-- |
|
28 |
|
29 /* function to open any folds containing a jumped-to line before jumping… |
|
30 function JumpToLine() |
|
31 { |
|
32 var lineNum; |
|
33 lineNum = window.location.hash; |
|
34 lineNum = lineNum.substr(1); /* strip off '#' */ |
|
35 |
|
36 if (lineNum.indexOf('L') == -1) { |
|
37 lineNum = 'L'+lineNum; |
|
38 } |
|
39 lineElem = document.getElementById(lineNum); |
|
40 /* Always jump to new location even if the line was hidden inside a fo… |
|
41 * we corrected the raw number to a line ID. |
|
42 */ |
|
43 if (lineElem) { |
|
44 lineElem.scrollIntoView(true); |
|
45 } |
|
46 return true; |
|
47 } |
|
48 if ('onhashchange' in window) { |
|
49 window.onhashchange = JumpToLine; |
|
50 } |
|
51 |
|
52 --> |
|
53 </script> |
|
54 </head> |
|
55 <body onload='JumpToLine();'> |
|
56 <pre id='vimCodeElement'> |
|
57 <span id="L1" class="LineNr"> 1 </span><span class="PreProc">#include <… |
|
58 <span id="L2" class="LineNr"> 2 </span><span class="PreProc">#include <… |
|
59 <span id="L3" class="LineNr"> 3 </span><span class="PreProc">#include <… |
|
60 <span id="L4" class="LineNr"> 4 </span><span class="PreProc">#include <… |
|
61 <span id="L5" class="LineNr"> 5 </span><span class="PreProc">#include <… |
|
62 <span id="L6" class="LineNr"> 6 </span><span class="PreProc">#include <… |
|
63 <span id="L7" class="LineNr"> 7 </span><span class="PreProc">#include <… |
|
64 <span id="L8" class="LineNr"> 8 </span> |
|
65 <span id="L9" class="LineNr"> 9 </span><span class="Comment">// Constru… |
|
66 <span id="L10" class="LineNr"> 10 </span>ODE::ODE(std::vector<std::<s… |
|
67 <span id="L11" class="LineNr"> 11 </span> (*f_in)(<span c… |
|
68 <span id="L12" class="LineNr"> 12 </span> <span c… |
|
69 <span id="L13" class="LineNr"> 13 </span> <span class="Type">con… |
|
70 <span id="L14" class="LineNr"> 14 </span> <span class="Type">con… |
|
71 <span id="L15" class="LineNr"> 15 </span> <span class="Type">con… |
|
72 <span id="L16" class="LineNr"> 16 </span> <span class="Type">con… |
|
73 <span id="L17" class="LineNr"> 17 </span> <span class="Type">con… |
|
74 <span id="L18" class="LineNr"> 18 </span> <span class="Type">con… |
|
75 <span id="L19" class="LineNr"> 19 </span> <span class="Type">con… |
|
76 <span id="L20" class="LineNr"> 20 </span> <span class="Type">con… |
|
77 <span id="L21" class="LineNr"> 21 </span> <span class="Type">con… |
|
78 <span id="L22" class="LineNr"> 22 </span> : f(f_in), |
|
79 <span id="L23" class="LineNr"> 23 </span> a(a_in), b(b_in), |
|
80 <span id="L24" class="LineNr"> 24 </span> h((b_in-a_in)*h_start), |
|
81 <span id="L25" class="LineNr"> 25 </span> n_max(max_steps), |
|
82 <span id="L26" class="LineNr"> 26 </span> delta(delta_in), epsilon(ep… |
|
83 <span id="L27" class="LineNr"> 27 </span> power(power_in), safety(saf… |
|
84 <span id="L28" class="LineNr"> 28 </span>{ |
|
85 <span id="L29" class="LineNr"> 29 </span> x_list.push_back(a); |
|
86 <span id="L30" class="LineNr"> 30 </span> y_list.push_back(y_start); |
|
87 <span id="L31" class="LineNr"> 31 </span> |
|
88 <span id="L32" class="LineNr"> 32 </span> <span class="Comment">// Perf… |
|
89 <span id="L33" class="LineNr"> 33 </span> rkdriver(); |
|
90 <span id="L34" class="LineNr"> 34 </span>} |
|
91 <span id="L35" class="LineNr"> 35 </span> |
|
92 <span id="L36" class="LineNr"> 36 </span> |
|
93 <span id="L37" class="LineNr"> 37 </span><span class="Comment">// Estima… |
|
94 <span id="L38" class="LineNr"> 38 </span>Floattype ODE::tau(<span class=… |
|
95 <span id="L39" class="LineNr"> 39 </span> <span class=… |
|
96 <span id="L40" class="LineNr"> 40 </span>{ |
|
97 <span id="L41" class="LineNr"> 41 </span> <span class="Statement">retur… |
|
98 <span id="L42" class="LineNr"> 42 </span>} |
|
99 <span id="L43" class="LineNr"> 43 </span> |
|
100 <span id="L44" class="LineNr"> 44 </span><span class="Comment">// Runge-… |
|
101 <span id="L45" class="LineNr"> 45 </span><span class="Type">void</span> … |
|
102 <span id="L46" class="LineNr"> 46 </span> <span class=… |
|
103 <span id="L47" class="LineNr"> 47 </span> std::vector&… |
|
104 <span id="L48" class="LineNr"> 48 </span> std::vector&… |
|
105 <span id="L49" class="LineNr"> 49 </span>{ |
|
106 <span id="L50" class="LineNr"> 50 </span> <span class="Comment">// Deno… |
|
107 <span id="L51" class="LineNr"> 51 </span> <span class="Type">const</spa… |
|
108 <span id="L52" class="LineNr"> 52 </span> |
|
109 <span id="L53" class="LineNr"> 53 </span> <span class="Comment">// Eval… |
|
110 <span id="L54" class="LineNr"> 54 </span> (<span class="Type">void</spa… |
|
111 <span id="L55" class="LineNr"> 55 </span> <span class="Type">const</spa… |
|
112 <span id="L56" class="LineNr"> 56 </span> <span class="Type">const</spa… |
|
113 <span id="L57" class="LineNr"> 57 </span> |
|
114 <span id="L58" class="LineNr"> 58 </span> <span class="Comment">// Writ… |
|
115 <span id="L59" class="LineNr"> 59 </span> y1 = y0 + k12*h; |
|
116 <span id="L60" class="LineNr"> 60 </span> dy = (k0 - k12) * h/den2; |
|
117 <span id="L61" class="LineNr"> 61 </span>} |
|
118 <span id="L62" class="LineNr"> 62 </span> |
|
119 <span id="L63" class="LineNr"> 63 </span> |
|
120 <span id="L64" class="LineNr"> 64 </span><span class="Comment">// ODE dr… |
|
121 <span id="L65" class="LineNr"> 65 </span><span class="Type">void</span> … |
|
122 <span id="L66" class="LineNr"> 66 </span>{ |
|
123 <span id="L67" class="LineNr"> 67 </span> std::vector<std::<span cla… |
|
124 <span id="L68" class="LineNr"> 68 </span> std::vector<std::<span cla… |
|
125 <span id="L69" class="LineNr"> 69 </span> |
|
126 <span id="L70" class="LineNr"> 70 </span> std::<span class="Type">compl… |
|
127 <span id="L71" class="LineNr"> 71 </span> Floattype err, tol; |
|
128 <span id="L72" class="LineNr"> 72 </span> std::vector<std::<span cla… |
|
129 <span id="L73" class="LineNr"> 73 </span> |
|
130 <span id="L74" class="LineNr"> 74 </span> <span class="Statement">while… |
|
131 <span id="L75" class="LineNr"> 75 </span> || x_list.back().imag() &… |
|
132 <span id="L76" class="LineNr"> 76 </span> { |
|
133 <span id="L77" class="LineNr"> 77 </span> <span class="Comment">// Ge… |
|
134 <span id="L78" class="LineNr"> 78 </span> x = x_list.back(); |
|
135 <span id="L79" class="LineNr"> 79 </span> y = y_list.back(); |
|
136 <span id="L80" class="LineNr"> 80 </span> |
|
137 <span id="L81" class="LineNr"> 81 </span> <span class="Comment">// Ma… |
|
138 <span id="L82" class="LineNr"> 82 </span> <span class="Statement">if<… |
|
139 <span id="L83" class="LineNr"> 83 </span> || (x + h).imag() > b.i… |
|
140 <span id="L84" class="LineNr"> 84 </span> h = b - x; |
|
141 <span id="L85" class="LineNr"> 85 </span> |
|
142 <span id="L86" class="LineNr"> 86 </span> <span class="Comment">// Ru… |
|
143 <span id="L87" class="LineNr"> 87 </span> rkstep12(x, y, y1, dy); |
|
144 <span id="L88" class="LineNr"> 88 </span> |
|
145 <span id="L89" class="LineNr"> 89 </span> <span class="Comment">// De… |
|
146 <span id="L90" class="LineNr"> 90 </span> err = cnorm(dy); <span clas… |
|
147 <span id="L91" class="LineNr"> 91 </span> tol = tau(y, h); <span clas… |
|
148 <span id="L92" class="LineNr"> 92 </span> <span class="Statement">if<… |
|
149 <span id="L93" class="LineNr"> 93 </span> x_list.push_back(x+h); |
|
150 <span id="L94" class="LineNr"> 94 </span> y_list.push_back(y1); |
|
151 <span id="L95" class="LineNr"> 95 </span> } |
|
152 <span id="L96" class="LineNr"> 96 </span> |
|
153 <span id="L97" class="LineNr"> 97 </span> <span class="Comment">// Ch… |
|
154 <span id="L98" class="LineNr"> 98 </span> <span class="Statement">if<… |
|
155 <span id="L99" class="LineNr"> 99 </span> std::cerr << <span … |
|
156 <span id="L100" class="LineNr">100 </span> << <span… |
|
157 <span id="L101" class="LineNr">101 </span> << <span… |
|
158 <span id="L102" class="LineNr">102 </span> << <span… |
|
159 <span id="L103" class="LineNr">103 </span> << <span… |
|
160 <span id="L104" class="LineNr">104 </span> <span class="Statement">… |
|
161 <span id="L105" class="LineNr">105 </span> } |
|
162 <span id="L106" class="LineNr">106 </span> |
|
163 <span id="L107" class="LineNr">107 </span> <span class="Comment">// D… |
|
164 <span id="L108" class="LineNr">108 </span> std::<span class="Type">co… |
|
165 <span id="L109" class="LineNr">109 </span> <span class="Statement">if… |
|
166 <span id="L110" class="LineNr">110 </span> h = h*pow(tol/err, power… |
|
167 <span id="L111" class="LineNr">111 </span> <span class="Statement">el… |
|
168 <span id="L112" class="LineNr">112 </span> h = multiplicator*h; |
|
169 <span id="L113" class="LineNr">113 </span> } |
|
170 <span id="L114" class="LineNr">114 </span>} |
|
171 <span id="L115" class="LineNr">115 </span> |
|
172 <span id="L116" class="LineNr">116 </span> |
|
173 <span id="L117" class="LineNr">117 </span><span class="Comment">// Retur… |
|
174 <span id="L118" class="LineNr">118 </span>Inttype ODE::steps() |
|
175 <span id="L119" class="LineNr">119 </span>{ |
|
176 <span id="L120" class="LineNr">120 </span> <span class="Statement">retu… |
|
177 <span id="L121" class="LineNr">121 </span>} |
|
178 <span id="L122" class="LineNr">122 </span> |
|
179 <span id="L123" class="LineNr">123 </span><span class="Type">void</span>… |
|
180 <span id="L124" class="LineNr">124 </span>{ |
|
181 <span id="L125" class="LineNr">125 </span> <span class="Statement">for<… |
|
182 <span id="L126" class="LineNr">126 </span> std::cout << x_list[… |
|
183 <span id="L127" class="LineNr">127 </span> <span class="Statement">fo… |
|
184 <span id="L128" class="LineNr">128 </span> std::cout << y_lis… |
|
185 <span id="L129" class="LineNr">129 </span> std::cout << <span c… |
|
186 <span id="L130" class="LineNr">130 </span> } |
|
187 <span id="L131" class="LineNr">131 </span>} |
|
188 <span id="L132" class="LineNr">132 </span> |
|
189 <span id="L133" class="LineNr">133 </span><span class="Comment">// Write… |
|
190 <span id="L134" class="LineNr">134 </span><span class="Type">void</span>… |
|
191 <span id="L135" class="LineNr">135 </span>{ |
|
192 <span id="L136" class="LineNr">136 </span> std::ofstream outstream; |
|
193 <span id="L137" class="LineNr">137 </span> |
|
194 <span id="L138" class="LineNr">138 </span> <span class="Comment">// Ope… |
|
195 <span id="L139" class="LineNr">139 </span> outstream.open(filename); |
|
196 <span id="L140" class="LineNr">140 </span> <span class="Statement">if</… |
|
197 <span id="L141" class="LineNr">141 </span> std::cerr << <span c… |
|
198 <span id="L142" class="LineNr">142 </span> << filenam… |
|
199 <span id="L143" class="LineNr">143 </span> <span class="Statement">re… |
|
200 <span id="L144" class="LineNr">144 </span> } |
|
201 <span id="L145" class="LineNr">145 </span> |
|
202 <span id="L146" class="LineNr">146 </span> <span class="Comment">// Wri… |
|
203 <span id="L147" class="LineNr">147 </span> <span class="Statement">for<… |
|
204 <span id="L148" class="LineNr">148 </span> outstream << x_list[… |
|
205 <span id="L149" class="LineNr">149 </span> outstream << x_list[… |
|
206 <span id="L150" class="LineNr">150 </span> <span class="Statement">fo… |
|
207 <span id="L151" class="LineNr">151 </span> outstream << y_lis… |
|
208 <span id="L152" class="LineNr">152 </span> outstream << y_lis… |
|
209 <span id="L153" class="LineNr">153 </span> } |
|
210 <span id="L154" class="LineNr">154 </span> outstream << <span c… |
|
211 <span id="L155" class="LineNr">155 </span> } |
|
212 <span id="L156" class="LineNr">156 </span> |
|
213 <span id="L157" class="LineNr">157 </span> <span class="Comment">// Clo… |
|
214 <span id="L158" class="LineNr">158 </span> outstream.close(); |
|
215 <span id="L159" class="LineNr">159 </span> |
|
216 <span id="L160" class="LineNr">160 </span> <span class="Statement">if</… |
|
217 <span id="L161" class="LineNr">161 </span> std::cout << <span c… |
|
218 <span id="L162" class="LineNr">162 </span>} |
|
219 <span id="L163" class="LineNr">163 </span> |
|
220 </pre> |
|
221 </body> |
|
222 </html> |
|
223 <!-- vim: set foldmethod=manual : --> |
