1
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
2
<html xmlns="http://www.w3.org/1999/xhtml">
4
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
5
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
6
<title>tesseract: /usr/local/google/home/jbreiden/tesseract-ocr-read-only/ccstruct/quspline.cpp Source File</title>
8
<link href="tabs.css" rel="stylesheet" type="text/css"/>
9
<link href="doxygen.css" rel="stylesheet" type="text/css" />
10
<link href="navtree.css" rel="stylesheet" type="text/css"/>
11
<script type="text/javascript" src="jquery.js"></script>
12
<script type="text/javascript" src="resize.js"></script>
13
<script type="text/javascript" src="navtree.js"></script>
14
<script type="text/javascript">
15
$(document).ready(initResizable);
17
<link href="search/search.css" rel="stylesheet" type="text/css"/>
18
<script type="text/javascript" src="search/search.js"></script>
19
<script type="text/javascript">
20
$(document).ready(function() { searchBox.OnSelectItem(0); });
25
<div id="top"><!-- do not remove this div! -->
29
<table cellspacing="0" cellpadding="0">
31
<tr style="height: 56px;">
34
<td style="padding-left: 0.5em;">
35
<div id="projectname">tesseract
36
 <span id="projectnumber">3.03</span>
48
<!-- Generated by Doxygen 1.7.6.1 -->
49
<script type="text/javascript">
50
var searchBox = new SearchBox("searchBox", "search",false,'Search');
52
<div id="navrow1" class="tabs">
54
<li><a href="index.html"><span>Main Page</span></a></li>
55
<li><a href="pages.html"><span>Related Pages</span></a></li>
56
<li><a href="modules.html"><span>Modules</span></a></li>
57
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
58
<li><a href="annotated.html"><span>Classes</span></a></li>
59
<li class="current"><a href="files.html"><span>Files</span></a></li>
61
<div id="MSearchBox" class="MSearchBoxInactive">
63
<img id="MSearchSelect" src="search/mag_sel.png"
64
onmouseover="return searchBox.OnSearchSelectShow()"
65
onmouseout="return searchBox.OnSearchSelectHide()"
67
<input type="text" id="MSearchField" value="Search" accesskey="S"
68
onfocus="searchBox.OnSearchFieldFocus(true)"
69
onblur="searchBox.OnSearchFieldFocus(false)"
70
onkeyup="searchBox.OnSearchFieldChange(event)"/>
71
</span><span class="right">
72
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
78
<div id="navrow2" class="tabs2">
80
<li><a href="files.html"><span>File List</span></a></li>
81
<li><a href="globals.html"><span>File Members</span></a></li>
85
<div id="side-nav" class="ui-resizable side-nav-resizable">
87
<div id="nav-tree-contents">
90
<div id="splitbar" style="-moz-user-select:none;"
91
class="ui-resizable-handle">
94
<script type="text/javascript">
95
initNavTree('a00785.html','');
97
<div id="doc-content">
99
<div class="headertitle">
100
<div class="title">/usr/local/google/home/jbreiden/tesseract-ocr-read-only/ccstruct/quspline.cpp</div> </div>
102
<div class="contents">
103
<a href="a00785.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/**********************************************************************</span>
104
<a name="l00002"></a>00002 <span class="comment"> * File: quspline.cpp (Formerly qspline.c)</span>
105
<a name="l00003"></a>00003 <span class="comment"> * Description: Code for the QSPLINE class.</span>
106
<a name="l00004"></a>00004 <span class="comment"> * Author: Ray Smith</span>
107
<a name="l00005"></a>00005 <span class="comment"> * Created: Tue Oct 08 17:16:12 BST 1991</span>
108
<a name="l00006"></a>00006 <span class="comment"> *</span>
109
<a name="l00007"></a>00007 <span class="comment"> * (C) Copyright 1991, Hewlett-Packard Ltd.</span>
110
<a name="l00008"></a>00008 <span class="comment"> ** Licensed under the Apache License, Version 2.0 (the "License");</span>
111
<a name="l00009"></a>00009 <span class="comment"> ** you may not use this file except in compliance with the License.</span>
112
<a name="l00010"></a>00010 <span class="comment"> ** You may obtain a copy of the License at</span>
113
<a name="l00011"></a>00011 <span class="comment"> ** http://www.apache.org/licenses/LICENSE-2.0</span>
114
<a name="l00012"></a>00012 <span class="comment"> ** Unless required by applicable law or agreed to in writing, software</span>
115
<a name="l00013"></a>00013 <span class="comment"> ** distributed under the License is distributed on an "AS IS" BASIS,</span>
116
<a name="l00014"></a>00014 <span class="comment"> ** WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</span>
117
<a name="l00015"></a>00015 <span class="comment"> ** See the License for the specific language governing permissions and</span>
118
<a name="l00016"></a>00016 <span class="comment"> ** limitations under the License.</span>
119
<a name="l00017"></a>00017 <span class="comment"> *</span>
120
<a name="l00018"></a>00018 <span class="comment"> **********************************************************************/</span>
121
<a name="l00019"></a>00019
122
<a name="l00020"></a>00020 <span class="preprocessor">#include "allheaders.h"</span>
123
<a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="a00838.html">memry.h</a>"</span>
124
<a name="l00022"></a>00022 <span class="preprocessor">#include "<a class="code" href="a00783.html">quadlsq.h</a>"</span>
125
<a name="l00023"></a>00023 <span class="preprocessor">#include "<a class="code" href="a00786.html">quspline.h</a>"</span>
126
<a name="l00024"></a>00024
127
<a name="l00025"></a>00025 <span class="comment">// Include automatically generated configuration file if running autoconf.</span>
128
<a name="l00026"></a>00026 <span class="preprocessor">#ifdef HAVE_CONFIG_H</span>
129
<a name="l00027"></a>00027 <span class="preprocessor"></span><span class="preprocessor">#include "<a class="code" href="a00950.html">config_auto.h</a>"</span>
130
<a name="l00028"></a>00028 <span class="preprocessor">#endif</span>
131
<a name="l00029"></a>00029 <span class="preprocessor"></span>
132
<a name="l00030"></a><a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">00030</a> <span class="preprocessor">#define QSPLINE_PRECISION 16 //no of steps to draw</span>
133
<a name="l00031"></a>00031 <span class="preprocessor"></span>
134
<a name="l00032"></a>00032 <span class="comment">/**********************************************************************</span>
135
<a name="l00033"></a>00033 <span class="comment"> * QSPLINE::QSPLINE</span>
136
<a name="l00034"></a>00034 <span class="comment"> *</span>
137
<a name="l00035"></a>00035 <span class="comment"> * Constructor to build a QSPLINE given the components used in the old code.</span>
138
<a name="l00036"></a>00036 <span class="comment"> **********************************************************************/</span>
139
<a name="l00037"></a>00037
140
<a name="l00038"></a><a class="code" href="a00513.html#ac9ee7256ceae55a1fdbba01fb6103f54">00038</a> <a class="code" href="a00513.html#a550dbc1b50898c27be5237d72812d1c7">QSPLINE::QSPLINE</a>( <span class="comment">//constructor</span>
141
<a name="l00039"></a>00039 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> <a class="code" href="a01042.html#acd4ea858d7133e40518e0832ff9d94b2">count</a>, <span class="comment">//no of segments</span>
142
<a name="l00040"></a>00040 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *xstarts, <span class="comment">//start coords</span>
143
<a name="l00041"></a>00041 <span class="keywordtype">double</span> *coeffs <span class="comment">//coefficients</span>
144
<a name="l00042"></a>00042 ) {
145
<a name="l00043"></a>00043 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> index; <span class="comment">//segment index</span>
146
<a name="l00044"></a>00044
147
<a name="l00045"></a>00045 <span class="comment">//get memory</span>
148
<a name="l00046"></a>00046 xcoords = (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((count + 1) * <span class="keyword">sizeof</span> (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
149
<a name="l00047"></a>00047 quadratics = (<a class="code" href="a00514.html">QUAD_COEFFS</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> (count * <span class="keyword">sizeof</span> (<a class="code" href="a00514.html">QUAD_COEFFS</a>));
150
<a name="l00048"></a>00048 segments = <a class="code" href="a01042.html#acd4ea858d7133e40518e0832ff9d94b2">count</a>;
151
<a name="l00049"></a>00049 <span class="keywordflow">for</span> (index = 0; index < segments; index++) {
152
<a name="l00050"></a>00050 <span class="comment">//copy them</span>
153
<a name="l00051"></a>00051 xcoords[index] = xstarts[index];
154
<a name="l00052"></a>00052 quadratics[index] = <a class="code" href="a00514.html">QUAD_COEFFS</a> (coeffs[index * 3],
155
<a name="l00053"></a>00053 coeffs[index * 3 + 1],
156
<a name="l00054"></a>00054 coeffs[index * 3 + 2]);
157
<a name="l00055"></a>00055 }
158
<a name="l00056"></a>00056 <span class="comment">//right edge</span>
159
<a name="l00057"></a>00057 xcoords[index] = xstarts[index];
160
<a name="l00058"></a>00058 }
161
<a name="l00059"></a>00059
162
<a name="l00060"></a>00060
163
<a name="l00061"></a>00061 <span class="comment">/**********************************************************************</span>
164
<a name="l00062"></a>00062 <span class="comment"> * QSPLINE::QSPLINE</span>
165
<a name="l00063"></a>00063 <span class="comment"> *</span>
166
<a name="l00064"></a>00064 <span class="comment"> * Constructor to build a QSPLINE by appproximation of points.</span>
167
<a name="l00065"></a>00065 <span class="comment"> **********************************************************************/</span>
168
<a name="l00066"></a>00066
169
<a name="l00067"></a><a class="code" href="a00513.html#af79e3ad5c225e40126555bd88566d2bc">00067</a> <a class="code" href="a00513.html#a550dbc1b50898c27be5237d72812d1c7">QSPLINE::QSPLINE</a> ( <span class="comment">//constructor</span>
170
<a name="l00068"></a>00068 <span class="keywordtype">int</span> xstarts[], <span class="comment">//spline boundaries</span>
171
<a name="l00069"></a>00069 <span class="keywordtype">int</span> segcount, <span class="comment">//no of segments</span>
172
<a name="l00070"></a>00070 <span class="keywordtype">int</span> xpts[], <span class="comment">//points to fit</span>
173
<a name="l00071"></a>00071 <span class="keywordtype">int</span> ypts[], <span class="keywordtype">int</span> pointcount, <span class="comment">//no of pts</span>
174
<a name="l00072"></a>00072 <span class="keywordtype">int</span> degree <span class="comment">//fit required</span>
175
<a name="l00073"></a>00073 ) {
176
<a name="l00074"></a>00074 <span class="keyword">register</span> <span class="keywordtype">int</span> pointindex; <span class="comment">/*no along text line */</span>
177
<a name="l00075"></a>00075 <span class="keyword">register</span> <span class="keywordtype">int</span> segment; <span class="comment">/*segment no */</span>
178
<a name="l00076"></a>00076 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *ptcounts; <span class="comment">//no in each segment</span>
179
<a name="l00077"></a>00077 <a class="code" href="a00511.html">QLSQ</a> qlsq; <span class="comment">/*accumulator */</span>
180
<a name="l00078"></a>00078
181
<a name="l00079"></a>00079 segments = segcount;
182
<a name="l00080"></a>00080 xcoords = (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((segcount + 1) * <span class="keyword">sizeof</span> (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
183
<a name="l00081"></a>00081 ptcounts = (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((segcount + 1) * <span class="keyword">sizeof</span> (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
184
<a name="l00082"></a>00082 quadratics = (<a class="code" href="a00514.html">QUAD_COEFFS</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> (segcount * <span class="keyword">sizeof</span> (<a class="code" href="a00514.html">QUAD_COEFFS</a>));
185
<a name="l00083"></a>00083 memmove (xcoords, xstarts, (segcount + 1) * <span class="keyword">sizeof</span> (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
186
<a name="l00084"></a>00084 ptcounts[0] = 0; <span class="comment">/*none in any yet */</span>
187
<a name="l00085"></a>00085 <span class="keywordflow">for</span> (segment = 0, pointindex = 0; pointindex < pointcount; pointindex++) {
188
<a name="l00086"></a>00086 <span class="keywordflow">while</span> (segment < segcount && xpts[pointindex] >= xstarts[segment]) {
189
<a name="l00087"></a>00087 segment++; <span class="comment">/*try next segment */</span>
190
<a name="l00088"></a>00088 <span class="comment">/*cumulative counts */</span>
191
<a name="l00089"></a>00089 ptcounts[segment] = ptcounts[segment - 1];
192
<a name="l00090"></a>00090 }
193
<a name="l00091"></a>00091 ptcounts[segment]++; <span class="comment">/*no in previous partition */</span>
194
<a name="l00092"></a>00092 }
195
<a name="l00093"></a>00093 <span class="keywordflow">while</span> (segment < segcount) {
196
<a name="l00094"></a>00094 segment++;
197
<a name="l00095"></a>00095 <span class="comment">/*zero the rest */</span>
198
<a name="l00096"></a>00096 ptcounts[segment] = ptcounts[segment - 1];
199
<a name="l00097"></a>00097 }
200
<a name="l00098"></a>00098
201
<a name="l00099"></a>00099 <span class="keywordflow">for</span> (segment = 0; segment < segcount; segment++) {
202
<a name="l00100"></a>00100 qlsq.<a class="code" href="a00511.html#a9f6ee88370df57093e248a2918ab20b2">clear</a> ();
203
<a name="l00101"></a>00101 <span class="comment">/*first blob */</span>
204
<a name="l00102"></a>00102 pointindex = ptcounts[segment];
205
<a name="l00103"></a>00103 <span class="keywordflow">if</span> (pointindex > 0
206
<a name="l00104"></a>00104 && xpts[pointindex] != xpts[pointindex - 1]
207
<a name="l00105"></a>00105 && xpts[pointindex] != xstarts[segment])
208
<a name="l00106"></a>00106 qlsq.<a class="code" href="a00511.html#adbbaa8f4d49392cc26b46fabc476cc88">add</a> (xstarts[segment],
209
<a name="l00107"></a>00107 ypts[pointindex - 1]
210
<a name="l00108"></a>00108 + (ypts[pointindex] - ypts[pointindex - 1])
211
<a name="l00109"></a>00109 * (xstarts[segment] - xpts[pointindex - 1])
212
<a name="l00110"></a>00110 / (xpts[pointindex] - xpts[pointindex - 1]));
213
<a name="l00111"></a>00111 <span class="keywordflow">for</span> (; pointindex < ptcounts[segment + 1]; pointindex++) {
214
<a name="l00112"></a>00112 qlsq.<a class="code" href="a00511.html#adbbaa8f4d49392cc26b46fabc476cc88">add</a> (xpts[pointindex], ypts[pointindex]);
215
<a name="l00113"></a>00113 }
216
<a name="l00114"></a>00114 <span class="keywordflow">if</span> (pointindex > 0 && pointindex < pointcount
217
<a name="l00115"></a>00115 && xpts[pointindex] != xstarts[segment + 1])
218
<a name="l00116"></a>00116 qlsq.<a class="code" href="a00511.html#adbbaa8f4d49392cc26b46fabc476cc88">add</a> (xstarts[segment + 1],
219
<a name="l00117"></a>00117 ypts[pointindex - 1]
220
<a name="l00118"></a>00118 + (ypts[pointindex] - ypts[pointindex - 1])
221
<a name="l00119"></a>00119 * (xstarts[segment + 1] - xpts[pointindex - 1])
222
<a name="l00120"></a>00120 / (xpts[pointindex] - xpts[pointindex - 1]));
223
<a name="l00121"></a>00121 qlsq.<a class="code" href="a00511.html#a847b2fa0f918e971a56e8db0ce8ed18f">fit</a> (degree);
224
<a name="l00122"></a>00122 quadratics[segment].<a class="code" href="a00514.html#aa9ab5e8947242fa1de50d9863165e032">a</a> = qlsq.<a class="code" href="a00511.html#a748cebc68a35dcb80970e2cfc1cace2d">get_a</a> ();
225
<a name="l00123"></a>00123 quadratics[segment].<a class="code" href="a00514.html#a58eae9f9c9aa7b738c5493b8da560747">b</a> = qlsq.<a class="code" href="a00511.html#ab8de83e0968609f28aa83854431e2593">get_b</a> ();
226
<a name="l00124"></a>00124 quadratics[segment].<a class="code" href="a00514.html#a7f5483bdd3103b7f4a732c8fbc0faa91">c</a> = qlsq.<a class="code" href="a00511.html#aeeedb9b168cef907167e4530dd4bf68a">get_c</a> ();
227
<a name="l00125"></a>00125 }
228
<a name="l00126"></a>00126 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(ptcounts);
229
<a name="l00127"></a>00127 }
230
<a name="l00128"></a>00128
231
<a name="l00129"></a>00129
232
<a name="l00130"></a>00130 <span class="comment">/**********************************************************************</span>
233
<a name="l00131"></a>00131 <span class="comment"> * QSPLINE::QSPLINE</span>
234
<a name="l00132"></a>00132 <span class="comment"> *</span>
235
<a name="l00133"></a>00133 <span class="comment"> * Constructor to build a QSPLINE from another.</span>
236
<a name="l00134"></a>00134 <span class="comment"> **********************************************************************/</span>
237
<a name="l00135"></a>00135
238
<a name="l00136"></a><a class="code" href="a00513.html#a388e310940b30950a7a5f32a9e4a41c2">00136</a> <a class="code" href="a00513.html#a550dbc1b50898c27be5237d72812d1c7">QSPLINE::QSPLINE</a>( <span class="comment">//constructor</span>
239
<a name="l00137"></a>00137 <span class="keyword">const</span> <a class="code" href="a00513.html">QSPLINE</a> &src) {
240
<a name="l00138"></a>00138 segments = 0;
241
<a name="l00139"></a>00139 xcoords = NULL;
242
<a name="l00140"></a>00140 quadratics = NULL;
243
<a name="l00141"></a>00141 *<span class="keyword">this</span> = src;
244
<a name="l00142"></a>00142 }
245
<a name="l00143"></a>00143
246
<a name="l00144"></a>00144
247
<a name="l00145"></a>00145 <span class="comment">/**********************************************************************</span>
248
<a name="l00146"></a>00146 <span class="comment"> * QSPLINE::~QSPLINE</span>
249
<a name="l00147"></a>00147 <span class="comment"> *</span>
250
<a name="l00148"></a>00148 <span class="comment"> * Destroy a QSPLINE.</span>
251
<a name="l00149"></a>00149 <span class="comment"> **********************************************************************/</span>
252
<a name="l00150"></a>00150
253
<a name="l00151"></a><a class="code" href="a00513.html#a7619e90d02d5127754c93073378a738a">00151</a> <a class="code" href="a00513.html#a7619e90d02d5127754c93073378a738a">QSPLINE::~QSPLINE</a> ( <span class="comment">//constructor</span>
254
<a name="l00152"></a>00152 ) {
255
<a name="l00153"></a>00153 <span class="keywordflow">if</span> (xcoords != NULL) {
256
<a name="l00154"></a>00154 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(xcoords);
257
<a name="l00155"></a>00155 xcoords = NULL;
258
<a name="l00156"></a>00156 }
259
<a name="l00157"></a>00157 <span class="keywordflow">if</span> (quadratics != NULL) {
260
<a name="l00158"></a>00158 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(quadratics);
261
<a name="l00159"></a>00159 quadratics = NULL;
262
<a name="l00160"></a>00160 }
263
<a name="l00161"></a>00161 }
264
<a name="l00162"></a>00162
265
<a name="l00163"></a>00163
266
<a name="l00164"></a>00164 <span class="comment">/**********************************************************************</span>
267
<a name="l00165"></a>00165 <span class="comment"> * QSPLINE::operator=</span>
268
<a name="l00166"></a>00166 <span class="comment"> *</span>
269
<a name="l00167"></a>00167 <span class="comment"> * Copy a QSPLINE</span>
270
<a name="l00168"></a>00168 <span class="comment"> **********************************************************************/</span>
271
<a name="l00169"></a>00169
272
<a name="l00170"></a><a class="code" href="a00513.html#a947d537ad3fd5ef782bc81055b291dcf">00170</a> <a class="code" href="a00513.html">QSPLINE</a> & <a class="code" href="a00513.html#a947d537ad3fd5ef782bc81055b291dcf">QSPLINE::operator= </a>( <span class="comment">//assignment</span>
273
<a name="l00171"></a>00171 <span class="keyword">const</span> <a class="code" href="a00513.html">QSPLINE</a> & source) {
274
<a name="l00172"></a>00172 <span class="keywordflow">if</span> (xcoords != NULL)
275
<a name="l00173"></a>00173 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(xcoords);
276
<a name="l00174"></a>00174 <span class="keywordflow">if</span> (quadratics != NULL)
277
<a name="l00175"></a>00175 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(quadratics);
278
<a name="l00176"></a>00176
279
<a name="l00177"></a>00177 segments = source.segments;
280
<a name="l00178"></a>00178 xcoords = (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((segments + 1) * <span class="keyword">sizeof</span> (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
281
<a name="l00179"></a>00179 quadratics = (<a class="code" href="a00514.html">QUAD_COEFFS</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> (segments * <span class="keyword">sizeof</span> (<a class="code" href="a00514.html">QUAD_COEFFS</a>));
282
<a name="l00180"></a>00180 memmove (xcoords, source.xcoords, (segments + 1) * sizeof (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a>));
283
<a name="l00181"></a>00181 memmove (quadratics, source.quadratics, segments * sizeof (<a class="code" href="a00514.html">QUAD_COEFFS</a>));
284
<a name="l00182"></a>00182 <span class="keywordflow">return</span> *<span class="keyword">this</span>;
285
<a name="l00183"></a>00183 }
286
<a name="l00184"></a>00184
287
<a name="l00185"></a>00185
288
<a name="l00186"></a>00186 <span class="comment">/**********************************************************************</span>
289
<a name="l00187"></a>00187 <span class="comment"> * QSPLINE::step</span>
290
<a name="l00188"></a>00188 <span class="comment"> *</span>
291
<a name="l00189"></a>00189 <span class="comment"> * Return the total of the step functions between the given coords.</span>
292
<a name="l00190"></a>00190 <span class="comment"> **********************************************************************/</span>
293
<a name="l00191"></a>00191
294
<a name="l00192"></a><a class="code" href="a00513.html#aaebd602e34d8b7f4e7b73884f1701a75">00192</a> <span class="keywordtype">double</span> <a class="code" href="a00513.html#aaebd602e34d8b7f4e7b73884f1701a75">QSPLINE::step</a>( <span class="comment">//find step functions</span>
295
<a name="l00193"></a>00193 <span class="keywordtype">double</span> x1, <span class="comment">//between coords</span>
296
<a name="l00194"></a>00194 <span class="keywordtype">double</span> x2) {
297
<a name="l00195"></a>00195 <span class="keywordtype">int</span> index1, index2; <span class="comment">//indices of coords</span>
298
<a name="l00196"></a>00196 <span class="keywordtype">double</span> total; <span class="comment">/*total steps */</span>
299
<a name="l00197"></a>00197
300
<a name="l00198"></a>00198 index1 = spline_index (x1);
301
<a name="l00199"></a>00199 index2 = spline_index (x2);
302
<a name="l00200"></a>00200 total = 0;
303
<a name="l00201"></a>00201 <span class="keywordflow">while</span> (index1 < index2) {
304
<a name="l00202"></a>00202 total +=
305
<a name="l00203"></a>00203 (double) quadratics[index1 + 1].<a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> ((<span class="keywordtype">float</span>) xcoords[index1 + 1]);
306
<a name="l00204"></a>00204 total -= (double) quadratics[index1].<a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> ((<span class="keywordtype">float</span>) xcoords[index1 + 1]);
307
<a name="l00205"></a>00205 index1++; <span class="comment">/*next segment */</span>
308
<a name="l00206"></a>00206 }
309
<a name="l00207"></a>00207 <span class="keywordflow">return</span> total; <span class="comment">/*total steps */</span>
310
<a name="l00208"></a>00208 }
311
<a name="l00209"></a>00209
312
<a name="l00210"></a>00210
313
<a name="l00211"></a>00211 <span class="comment">/**********************************************************************</span>
314
<a name="l00212"></a>00212 <span class="comment"> * QSPLINE::y</span>
315
<a name="l00213"></a>00213 <span class="comment"> *</span>
316
<a name="l00214"></a>00214 <span class="comment"> * Return the y value at the given x value.</span>
317
<a name="l00215"></a>00215 <span class="comment"> **********************************************************************/</span>
318
<a name="l00216"></a>00216
319
<a name="l00217"></a><a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">00217</a> <span class="keywordtype">double</span> <a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">QSPLINE::y</a>( <span class="comment">//evaluate</span>
320
<a name="l00218"></a>00218 <span class="keywordtype">double</span> x <span class="comment">//coord to evaluate at</span>
321
<a name="l00219"></a>00219 )<span class="keyword"> const </span>{
322
<a name="l00220"></a>00220 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> index; <span class="comment">//segment index</span>
323
<a name="l00221"></a>00221
324
<a name="l00222"></a>00222 index = spline_index (x);
325
<a name="l00223"></a>00223 <span class="keywordflow">return</span> quadratics[index].<a class="code" href="a00514.html#af7c4d23ebf42e6b74528204d1497d191">y</a> (x);<span class="comment">//in correct segment</span>
326
<a name="l00224"></a>00224 }
327
<a name="l00225"></a>00225
328
<a name="l00226"></a>00226
329
<a name="l00227"></a>00227 <span class="comment">/**********************************************************************</span>
330
<a name="l00228"></a>00228 <span class="comment"> * QSPLINE::spline_index</span>
331
<a name="l00229"></a>00229 <span class="comment"> *</span>
332
<a name="l00230"></a>00230 <span class="comment"> * Return the index to the largest xcoord not greater than x.</span>
333
<a name="l00231"></a>00231 <span class="comment"> **********************************************************************/</span>
334
<a name="l00232"></a>00232
335
<a name="l00233"></a>00233 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> QSPLINE::spline_index( <span class="comment">//evaluate</span>
336
<a name="l00234"></a>00234 <span class="keywordtype">double</span> x <span class="comment">//coord to evaluate at</span>
337
<a name="l00235"></a>00235 )<span class="keyword"> const </span>{
338
<a name="l00236"></a>00236 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> index; <span class="comment">//segment index</span>
339
<a name="l00237"></a>00237 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> bottom; <span class="comment">//bottom of range</span>
340
<a name="l00238"></a>00238 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> top; <span class="comment">//top of range</span>
341
<a name="l00239"></a>00239
342
<a name="l00240"></a>00240 bottom = 0;
343
<a name="l00241"></a>00241 top = segments;
344
<a name="l00242"></a>00242 <span class="keywordflow">while</span> (top - bottom > 1) {
345
<a name="l00243"></a>00243 index = (top + bottom) / 2; <span class="comment">//centre of range</span>
346
<a name="l00244"></a>00244 <span class="keywordflow">if</span> (x >= xcoords[index])
347
<a name="l00245"></a>00245 bottom = index; <span class="comment">//new min</span>
348
<a name="l00246"></a>00246 <span class="keywordflow">else</span>
349
<a name="l00247"></a>00247 top = index; <span class="comment">//new max</span>
350
<a name="l00248"></a>00248 }
351
<a name="l00249"></a>00249 <span class="keywordflow">return</span> bottom;
352
<a name="l00250"></a>00250 }
353
<a name="l00251"></a>00251
354
<a name="l00252"></a>00252
355
<a name="l00253"></a>00253 <span class="comment">/**********************************************************************</span>
356
<a name="l00254"></a>00254 <span class="comment"> * QSPLINE::move</span>
357
<a name="l00255"></a>00255 <span class="comment"> *</span>
358
<a name="l00256"></a>00256 <span class="comment"> * Reposition spline by vector</span>
359
<a name="l00257"></a>00257 <span class="comment"> **********************************************************************/</span>
360
<a name="l00258"></a>00258
361
<a name="l00259"></a><a class="code" href="a00513.html#a3d1808becf8a69d00747b444d9c76cf6">00259</a> <span class="keywordtype">void</span> <a class="code" href="a00513.html#a3d1808becf8a69d00747b444d9c76cf6">QSPLINE::move</a>( <span class="comment">// reposition spline</span>
362
<a name="l00260"></a>00260 <a class="code" href="a00409.html" title="integer coordinate">ICOORD</a> vec <span class="comment">// by vector</span>
363
<a name="l00261"></a>00261 ) {
364
<a name="l00262"></a>00262 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> segment; <span class="comment">//index of segment</span>
365
<a name="l00263"></a>00263 <a class="code" href="a00831.html#a8d41499d38c24d39b221ab0c158fe5a8">inT16</a> x_shift = vec.<a class="code" href="a00409.html#a59722a47c540007c58a539f0e35b3f33" title="access function">x</a> ();
366
<a name="l00264"></a>00264
367
<a name="l00265"></a>00265 <span class="keywordflow">for</span> (segment = 0; segment < segments; segment++) {
368
<a name="l00266"></a>00266 xcoords[segment] += x_shift;
369
<a name="l00267"></a>00267 quadratics[segment].<a class="code" href="a00514.html#a76a67699d3fe976ba0e6ea72f5c395a5">move</a> (vec);
370
<a name="l00268"></a>00268 }
371
<a name="l00269"></a>00269 xcoords[segment] += x_shift;
372
<a name="l00270"></a>00270 }
373
<a name="l00271"></a>00271
374
<a name="l00272"></a>00272
375
<a name="l00273"></a>00273 <span class="comment">/**********************************************************************</span>
376
<a name="l00274"></a>00274 <span class="comment"> * QSPLINE::overlap</span>
377
<a name="l00275"></a>00275 <span class="comment"> *</span>
378
<a name="l00276"></a>00276 <span class="comment"> * Return TRUE if spline2 overlaps this by no more than fraction less</span>
379
<a name="l00277"></a>00277 <span class="comment"> * than the bounds of this.</span>
380
<a name="l00278"></a>00278 <span class="comment"> **********************************************************************/</span>
381
<a name="l00279"></a>00279
382
<a name="l00280"></a><a class="code" href="a00513.html#af5a1e056b6b5c0348b640b5eaa987184">00280</a> <a class="code" href="a00831.html#a7712a7e28433d0ade59219a129549b6f">BOOL8</a> <a class="code" href="a00513.html#af5a1e056b6b5c0348b640b5eaa987184">QSPLINE::overlap</a>( <span class="comment">//test overlap</span>
383
<a name="l00281"></a>00281 <a class="code" href="a00513.html">QSPLINE</a> *spline2, <span class="comment">//2 cannot be smaller</span>
384
<a name="l00282"></a>00282 <span class="keywordtype">double</span> fraction <span class="comment">//by more than this</span>
385
<a name="l00283"></a>00283 ) {
386
<a name="l00284"></a>00284 <span class="keywordtype">int</span> leftlimit; <span class="comment">/*common left limit */</span>
387
<a name="l00285"></a>00285 <span class="keywordtype">int</span> rightlimit; <span class="comment">/*common right limit */</span>
388
<a name="l00286"></a>00286
389
<a name="l00287"></a>00287 leftlimit = xcoords[1];
390
<a name="l00288"></a>00288 rightlimit = xcoords[segments - 1];
391
<a name="l00289"></a>00289 <span class="comment">/*or too non-overlap */</span>
392
<a name="l00290"></a>00290 <span class="keywordflow">if</span> (spline2->segments < 3 || spline2->xcoords[1] > leftlimit + fraction * (rightlimit - leftlimit)
393
<a name="l00291"></a>00291 || spline2->xcoords[spline2->segments - 1] < rightlimit
394
<a name="l00292"></a>00292 - fraction * (rightlimit - leftlimit))
395
<a name="l00293"></a>00293 <span class="keywordflow">return</span> FALSE;
396
<a name="l00294"></a>00294 <span class="keywordflow">else</span>
397
<a name="l00295"></a>00295 <span class="keywordflow">return</span> TRUE;
398
<a name="l00296"></a>00296 }
399
<a name="l00297"></a>00297
400
<a name="l00298"></a>00298
401
<a name="l00299"></a>00299 <span class="comment">/**********************************************************************</span>
402
<a name="l00300"></a>00300 <span class="comment"> * extrapolate_spline</span>
403
<a name="l00301"></a>00301 <span class="comment"> *</span>
404
<a name="l00302"></a>00302 <span class="comment"> * Extrapolates the spline linearly using the same gradient as the</span>
405
<a name="l00303"></a>00303 <span class="comment"> * quadratic has at either end.</span>
406
<a name="l00304"></a>00304 <span class="comment"> **********************************************************************/</span>
407
<a name="l00305"></a>00305
408
<a name="l00306"></a><a class="code" href="a00513.html#a2459d3e247a16f9f01d6b77597193812">00306</a> <span class="keywordtype">void</span> <a class="code" href="a00513.html#a2459d3e247a16f9f01d6b77597193812">QSPLINE::extrapolate</a>( <span class="comment">//linear extrapolation</span>
409
<a name="l00307"></a>00307 <span class="keywordtype">double</span> gradient, <span class="comment">//gradient to use</span>
410
<a name="l00308"></a>00308 <span class="keywordtype">int</span> xmin, <span class="comment">//new left edge</span>
411
<a name="l00309"></a>00309 <span class="keywordtype">int</span> xmax <span class="comment">//new right edge</span>
412
<a name="l00310"></a>00310 ) {
413
<a name="l00311"></a>00311 <span class="keyword">register</span> <span class="keywordtype">int</span> segment; <span class="comment">/*current segment of spline */</span>
414
<a name="l00312"></a>00312 <span class="keywordtype">int</span> dest_segment; <span class="comment">//dest index</span>
415
<a name="l00313"></a>00313 <span class="keywordtype">int</span> *xstarts; <span class="comment">//new boundaries</span>
416
<a name="l00314"></a>00314 <a class="code" href="a00514.html">QUAD_COEFFS</a> *quads; <span class="comment">//new ones</span>
417
<a name="l00315"></a>00315 <span class="keywordtype">int</span> increment; <span class="comment">//in size</span>
418
<a name="l00316"></a>00316
419
<a name="l00317"></a>00317 increment = xmin < xcoords[0] ? 1 : 0;
420
<a name="l00318"></a>00318 <span class="keywordflow">if</span> (xmax > xcoords[segments])
421
<a name="l00319"></a>00319 increment++;
422
<a name="l00320"></a>00320 <span class="keywordflow">if</span> (increment == 0)
423
<a name="l00321"></a>00321 <span class="keywordflow">return</span>;
424
<a name="l00322"></a>00322 xstarts = (<span class="keywordtype">int</span> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((segments + 1 + increment) * <span class="keyword">sizeof</span> (int));
425
<a name="l00323"></a>00323 quads =
426
<a name="l00324"></a>00324 (<a class="code" href="a00514.html">QUAD_COEFFS</a> *) <a class="code" href="a00837.html#a4c14b091498fb7866a1dd0af51a591fe">alloc_mem</a> ((segments + increment) * <span class="keyword">sizeof</span> (<a class="code" href="a00514.html">QUAD_COEFFS</a>));
427
<a name="l00325"></a>00325 <span class="keywordflow">if</span> (xmin < xcoords[0]) {
428
<a name="l00326"></a>00326 xstarts[0] = xmin;
429
<a name="l00327"></a>00327 quads[0].<a class="code" href="a00514.html#aa9ab5e8947242fa1de50d9863165e032">a</a> = 0;
430
<a name="l00328"></a>00328 quads[0].<a class="code" href="a00514.html#a58eae9f9c9aa7b738c5493b8da560747">b</a> = gradient;
431
<a name="l00329"></a>00329 quads[0].<a class="code" href="a00514.html#a7f5483bdd3103b7f4a732c8fbc0faa91">c</a> = <a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> (xcoords[0]) - quads[0].<a class="code" href="a00514.html#a58eae9f9c9aa7b738c5493b8da560747">b</a> * xcoords[0];
432
<a name="l00330"></a>00330 dest_segment = 1;
433
<a name="l00331"></a>00331 }
434
<a name="l00332"></a>00332 <span class="keywordflow">else</span>
435
<a name="l00333"></a>00333 dest_segment = 0;
436
<a name="l00334"></a>00334 <span class="keywordflow">for</span> (segment = 0; segment < segments; segment++) {
437
<a name="l00335"></a>00335 xstarts[dest_segment] = xcoords[segment];
438
<a name="l00336"></a>00336 quads[dest_segment] = quadratics[segment];
439
<a name="l00337"></a>00337 dest_segment++;
440
<a name="l00338"></a>00338 }
441
<a name="l00339"></a>00339 xstarts[dest_segment] = xcoords[segment];
442
<a name="l00340"></a>00340 <span class="keywordflow">if</span> (xmax > xcoords[segments]) {
443
<a name="l00341"></a>00341 quads[dest_segment].<a class="code" href="a00514.html#aa9ab5e8947242fa1de50d9863165e032">a</a> = 0;
444
<a name="l00342"></a>00342 quads[dest_segment].<a class="code" href="a00514.html#a58eae9f9c9aa7b738c5493b8da560747">b</a> = gradient;
445
<a name="l00343"></a>00343 quads[dest_segment].<a class="code" href="a00514.html#a7f5483bdd3103b7f4a732c8fbc0faa91">c</a> = <a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> (xcoords[segments])
446
<a name="l00344"></a>00344 - quads[dest_segment].<a class="code" href="a00514.html#a58eae9f9c9aa7b738c5493b8da560747">b</a> * xcoords[segments];
447
<a name="l00345"></a>00345 dest_segment++;
448
<a name="l00346"></a>00346 xstarts[dest_segment] = xmax + 1;
449
<a name="l00347"></a>00347 }
450
<a name="l00348"></a>00348 segments = dest_segment;
451
<a name="l00349"></a>00349 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(xcoords);
452
<a name="l00350"></a>00350 <a class="code" href="a00837.html#a9ef1862b8e026e23eef238af032858e4">free_mem</a>(quadratics);
453
<a name="l00351"></a>00351 xcoords = (<a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> *) xstarts;
454
<a name="l00352"></a>00352 quadratics = quads;
455
<a name="l00353"></a>00353 }
456
<a name="l00354"></a>00354
457
<a name="l00355"></a>00355
458
<a name="l00356"></a>00356 <span class="comment">/**********************************************************************</span>
459
<a name="l00357"></a>00357 <span class="comment"> * QSPLINE::plot</span>
460
<a name="l00358"></a>00358 <span class="comment"> *</span>
461
<a name="l00359"></a>00359 <span class="comment"> * Draw the QSPLINE in the given colour.</span>
462
<a name="l00360"></a>00360 <span class="comment"> **********************************************************************/</span>
463
<a name="l00361"></a>00361
464
<a name="l00362"></a>00362 <span class="preprocessor">#ifndef GRAPHICS_DISABLED</span>
465
<a name="l00363"></a><a class="code" href="a00513.html#ae2fb284b33116e4d3a0f9e1eaa125dff">00363</a> <span class="preprocessor"></span><span class="keywordtype">void</span> <a class="code" href="a00513.html#ae2fb284b33116e4d3a0f9e1eaa125dff">QSPLINE::plot</a>( <span class="comment">//draw it</span>
466
<a name="l00364"></a>00364 <a class="code" href="a00532.html">ScrollView</a>* window, <span class="comment">//window to draw in</span>
467
<a name="l00365"></a>00365 <a class="code" href="a00532.html#a100504544a5423a94222149ee9ed0fe8">ScrollView::Color</a> colour <span class="comment">//colour to draw in</span>
468
<a name="l00366"></a>00366 )<span class="keyword"> const </span>{
469
<a name="l00367"></a>00367 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> segment; <span class="comment">//index of segment</span>
470
<a name="l00368"></a>00368 <a class="code" href="a00831.html#a8d41499d38c24d39b221ab0c158fe5a8">inT16</a> <a class="code" href="a00513.html#aaebd602e34d8b7f4e7b73884f1701a75">step</a>; <span class="comment">//index of poly piece</span>
471
<a name="l00369"></a>00369 <span class="keywordtype">double</span> increment; <span class="comment">//x increment</span>
472
<a name="l00370"></a>00370 <span class="keywordtype">double</span> x; <span class="comment">//x coord</span>
473
<a name="l00371"></a>00371
474
<a name="l00372"></a>00372 window-><a class="code" href="a00532.html#a79855c525ec660b452382e9813d2edb5">Pen</a>(colour);
475
<a name="l00373"></a>00373 <span class="keywordflow">for</span> (segment = 0; segment < segments; segment++) {
476
<a name="l00374"></a>00374 increment =
477
<a name="l00375"></a>00375 (double) (xcoords[segment + 1] -
478
<a name="l00376"></a>00376 xcoords[segment]) / <a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">QSPLINE_PRECISION</a>;
479
<a name="l00377"></a>00377 x = xcoords[segment];
480
<a name="l00378"></a>00378 <span class="keywordflow">for</span> (step = 0; step <= <a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">QSPLINE_PRECISION</a>; step++) {
481
<a name="l00379"></a>00379 <span class="keywordflow">if</span> (segment == 0 && step == 0)
482
<a name="l00380"></a>00380 window-><a class="code" href="a00532.html#a730b182e350dfe37b0986d7f0591fee5">SetCursor</a>(x, quadratics[segment].<a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> (x));
483
<a name="l00381"></a>00381 <span class="keywordflow">else</span>
484
<a name="l00382"></a>00382 window-><a class="code" href="a00532.html#a101de80a4722fe821969f61f6e7303a6">DrawTo</a>(x, quadratics[segment].<a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> (x));
485
<a name="l00383"></a>00383 x += increment;
486
<a name="l00384"></a>00384 }
487
<a name="l00385"></a>00385 }
488
<a name="l00386"></a>00386 }
489
<a name="l00387"></a>00387 <span class="preprocessor">#endif</span>
490
<a name="l00388"></a>00388 <span class="preprocessor"></span>
491
<a name="l00389"></a><a class="code" href="a00513.html#a390001e1e6af8b44213511435fec23bd">00389</a> <span class="keywordtype">void</span> <a class="code" href="a00513.html#ae2fb284b33116e4d3a0f9e1eaa125dff">QSPLINE::plot</a>(Pix *pix)<span class="keyword"> const </span>{
492
<a name="l00390"></a>00390 <span class="keywordflow">if</span> (pix == NULL) {
493
<a name="l00391"></a>00391 <span class="keywordflow">return</span>;
494
<a name="l00392"></a>00392 }
495
<a name="l00393"></a>00393
496
<a name="l00394"></a>00394 <a class="code" href="a00831.html#aba1f582fd0168f3ff9225d8c90fa9eb8">inT32</a> segment; <span class="comment">// Index of segment</span>
497
<a name="l00395"></a>00395 <a class="code" href="a00831.html#a8d41499d38c24d39b221ab0c158fe5a8">inT16</a> <a class="code" href="a00513.html#aaebd602e34d8b7f4e7b73884f1701a75">step</a>; <span class="comment">// Index of poly piece</span>
498
<a name="l00396"></a>00396 <span class="keywordtype">double</span> increment; <span class="comment">// x increment</span>
499
<a name="l00397"></a>00397 <span class="keywordtype">double</span> x; <span class="comment">// x coord</span>
500
<a name="l00398"></a>00398 <span class="keywordtype">double</span> height = <span class="keyword">static_cast<</span><span class="keywordtype">double</span><span class="keyword">></span>(pixGetHeight(pix));
501
<a name="l00399"></a>00399 Pta* points = ptaCreate(<a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">QSPLINE_PRECISION</a> * segments);
502
<a name="l00400"></a>00400 <span class="keyword">const</span> <span class="keywordtype">int</span> kLineWidth = 5;
503
<a name="l00401"></a>00401
504
<a name="l00402"></a>00402 <span class="keywordflow">for</span> (segment = 0; segment < segments; segment++) {
505
<a name="l00403"></a>00403 increment = <span class="keyword">static_cast<</span><span class="keywordtype">double</span><span class="keyword">></span>((xcoords[segment + 1] -
506
<a name="l00404"></a>00404 xcoords[segment])) / <a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">QSPLINE_PRECISION</a>;
507
<a name="l00405"></a>00405 x = xcoords[segment];
508
<a name="l00406"></a>00406 <span class="keywordflow">for</span> (step = 0; step <= <a class="code" href="a00785.html#a5c2f35a8400249d498f29982dfc63bf7">QSPLINE_PRECISION</a>; step++) {
509
<a name="l00407"></a>00407 <span class="keywordtype">double</span> <a class="code" href="a00513.html#ab9f77ca80057e3844e880a1bf003beca">y</a> = height - quadratics[segment].<a class="code" href="a00514.html#af7c4d23ebf42e6b74528204d1497d191">y</a>(x);
510
<a name="l00408"></a>00408 ptaAddPt(points, x, y);
511
<a name="l00409"></a>00409 x += increment;
512
<a name="l00410"></a>00410 }
513
<a name="l00411"></a>00411 }
514
<a name="l00412"></a>00412
515
<a name="l00413"></a>00413 <span class="keywordflow">switch</span> (pixGetDepth(pix)) {
516
<a name="l00414"></a>00414 <span class="keywordflow">case</span> 1:
517
<a name="l00415"></a>00415 pixRenderPolyline(pix, points, kLineWidth, L_SET_PIXELS, 1);
518
<a name="l00416"></a>00416 <span class="keywordflow">break</span>;
519
<a name="l00417"></a>00417 <span class="keywordflow">case</span> 32:
520
<a name="l00418"></a>00418 pixRenderPolylineArb(pix, points, kLineWidth, 255, 0, 0, 1);
521
<a name="l00419"></a>00419 <span class="keywordflow">break</span>;
522
<a name="l00420"></a>00420 <span class="keywordflow">default</span>:
523
<a name="l00421"></a>00421 pixRenderPolyline(pix, points, kLineWidth, L_CLEAR_PIXELS, 1);
524
<a name="l00422"></a>00422 <span class="keywordflow">break</span>;
525
<a name="l00423"></a>00423 }
526
<a name="l00424"></a>00424 ptaDestroy(&points);
527
<a name="l00425"></a>00425 }
528
</pre></div></div><!-- contents -->
530
<!-- window showing the filter options -->
531
<div id="MSearchSelectWindow"
532
onmouseover="return searchBox.OnSearchSelectShow()"
533
onmouseout="return searchBox.OnSearchSelectHide()"
534
onkeydown="return searchBox.OnSearchSelectKey(event)">
535
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark"> </span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark"> </span>Classes</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Namespaces</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark"> </span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark"> </span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark"> </span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark"> </span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(9)"><span class="SelectionMark"> </span>Friends</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(10)"><span class="SelectionMark"> </span>Defines</a></div>
537
<!-- iframe showing the search results (closed by default) -->
538
<div id="MSearchResultsWindow">
539
<iframe src="javascript:void(0)" frameborder="0"
540
name="MSearchResults" id="MSearchResults">
544
<div id="nav-path" class="navpath">
546
<li class="navelem"><a class="el" href="a00785.html">quspline.cpp</a> </li>
548
<li class="footer">Generated on Mon Feb 3 2014 10:59:07 for tesseract by
549
<a href="http://www.doxygen.org/index.html">
550
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.6.1 </li>