1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283
|
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.6"/>
<title>ViennaCL - The Vienna Computing Library: tests/src/bisect.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
$(document).ready(initResizable);
$(window).load(resizeHeight);
</script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
$(document).ready(function() { searchBox.OnSelectItem(0); });
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td style="padding-left: 0.5em;">
<div id="projectname">ViennaCL - The Vienna Computing Library
 <span id="projectnumber">1.7.1</span>
</div>
<div id="projectbrief">Free open-source GPU-accelerated linear algebra and solver library.</div>
</td>
<td> <div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.6 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
<div id="nav-tree">
<div id="nav-tree-contents">
<div id="nav-sync" class="sync"></div>
</div>
</div>
<div id="splitbar" style="-moz-user-select:none;"
class="ui-resizable-handle">
</div>
</div>
<script type="text/javascript">
$(document).ready(function(){initNavTree('tests_2src_2bisect_8cpp_source.html','');});
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
<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>Macros</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(11)"><span class="SelectionMark"> </span>Pages</a></div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div class="header">
<div class="headertitle">
<div class="title">bisect.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="tests_2src_2bisect_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/* =========================================================================</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> Copyright (c) 2010-2016, Institute for Microelectronics,</span></div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment"> Institute for Analysis and Scientific Computing,</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> TU Wien.</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> Portions of this software are copyright by UChicago Argonne, LLC.</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"></span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> -----------------</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> ViennaCL - The Vienna Computing Library</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> -----------------</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment"></span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> Project Head: Karl Rupp rupp@iue.tuwien.ac.at</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment"></span></div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"> (A list of authors and contributors can be found in the PDF manual)</span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"></span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> License: MIT (X11), see file LICENSE in the base directory</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment">============================================================================= */</span></div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> </div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment">// includes, system</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="preprocessor">#include <stdlib.h></span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="preprocessor">#include <stdio.h></span></div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#include <string.h></span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> </div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> </div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment">// includes, project</span></div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> </div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor">#include "<a class="code" href="scalar_8hpp.html">viennacl/scalar.hpp</a>"</span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="preprocessor">#include "<a class="code" href="vector_8hpp.html">viennacl/vector.hpp</a>"</span></div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> </div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span> </div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="preprocessor">#include "<a class="code" href="bisect_8hpp.html">viennacl/linalg/bisect.hpp</a>"</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="preprocessor">#include "<a class="code" href="bisect__gpu_8hpp.html">viennacl/linalg/bisect_gpu.hpp</a>"</span></div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="preprocessor">#include "<a class="code" href="tql2_8hpp.html">viennacl/linalg/tql2.hpp</a>"</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> </div>
<div class="line"><a name="l00038"></a><span class="lineno"><a class="line" href="tests_2src_2bisect_8cpp.html#a6ebf6899d6c1c8b7b9d09be872c05aae"> 38</a></span> <span class="preprocessor">#define EPS 10.0e-4</span></div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> </div>
<div class="line"><a name="l00040"></a><span class="lineno"><a class="line" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c"> 40</a></span> <span class="keyword">typedef</span> <span class="keywordtype">float</span> <a class="code" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a>;</div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="comment">// declaration, forward</span></div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="keywordtype">bool</span> <a class="code" href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab">runTest</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> mat_size);</div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> </div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="keyword">template</span><<span class="keyword">typename</span> NumericT></div>
<div class="line"><a name="l00053"></a><span class="lineno"><a class="line" href="tests_2src_2bisect_8cpp.html#a4ff9461d32b6b78b32e2615bd0fb5363"> 53</a></span> <span class="keywordtype">void</span> <a class="code" href="tests_2src_2bisect_8cpp.html#a4ff9461d32b6b78b32e2615bd0fb5363">initInputData</a>(std::vector<NumericT> &diagonal, std::vector<NumericT> &superdiagonal, <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> mat_size)</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> {</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  </div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  srand(278217421);</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span> </div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span> <span class="preprocessor">#define RANDOM_VALUES false</span></div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span> </div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  <span class="keywordflow">if</span> (<a class="code" href="tests_2src_2bisect_8cpp.html#a429ac4051f4f4f60df33112d161f8ca6">RANDOM_VALUES</a> == <span class="keyword">true</span>)</div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  {</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  <span class="comment">// Initialize diagonal and superdiagonal elements with random values</span></div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i < mat_size; ++i)</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  {</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  diagonal[i] = <span class="keyword">static_cast<</span><a class="code" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a><span class="keyword">></span>(2.0 * (((double)rand()</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  / (double) RAND_MAX) - 0.5));</div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  superdiagonal[i] = <span class="keyword">static_cast<</span><a class="code" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a><span class="keyword">></span>(2.0 * (((double)rand()</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  / (double) RAND_MAX) - 0.5));</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  }</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  }</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  { </div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="comment">// Initialize diagonal and superdiagonal elements with modulo values</span></div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  <span class="comment">// This will cause in many multiple eigenvalues.</span></div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i < mat_size; ++i)</div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  {</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  diagonal[i] = ((<a class="code" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a>)(i % 3)) - 4.5f;</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  superdiagonal[i] = ((<a class="code" href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a>)(i % 3)) - 5.5f;</div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  }</div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  }</div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <span class="comment">// the first element of s is used as padding on the device (thus the</span></div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="comment">// whole vector is copied to the device but the kernels are launched</span></div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="comment">// with (s+1) as start address</span></div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  superdiagonal[0] = 0.0f; </div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span> }</div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span> </div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span> </div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> <span class="comment">// Program main</span></div>
<div class="line"><a name="l00091"></a><span class="lineno"><a class="line" href="tests_2src_2bisect_8cpp.html#ae66f6b31b5ad750f1fe042a706a4e3d4"> 91</a></span> <span class="comment"></span><span class="keywordtype">int</span> <a class="code" href="tests_2src_2bisect_8cpp.html#ae66f6b31b5ad750f1fe042a706a4e3d4">main</a>()</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span> {</div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="keywordtype">bool</span> test_result = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span> </div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <span class="comment">// run test for large matrix</span></div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  test_result = <a class="code" href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab">runTest</a>(550);</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="keywordflow">if</span>(test_result == <span class="keyword">true</span>)</div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  {</div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  std::cout << <span class="stringliteral">"First Test Succeeded!"</span> << std::endl << std::endl;</div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  }</div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  {</div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  std::cout << <span class="stringliteral">"---TEST FAILED---"</span> << std::endl;</div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  <span class="keywordflow">return</span> EXIT_FAILURE;</div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  }</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span> </div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  <span class="comment">// run test for small matrix</span></div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  test_result = <a class="code" href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab">runTest</a>(96);</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> </div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  <span class="keywordflow">if</span>(test_result == <span class="keyword">true</span>)</div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  {</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  std::cout << std::endl << <span class="stringliteral">"---TEST SUCCESSFULLY COMPLETED---"</span> << std::endl;</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="keywordflow">return</span> EXIT_SUCCESS;</div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  }</div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  {</div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  std::cout << <span class="stringliteral">"---TEST FAILED---"</span> << std::endl;</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="keywordflow">return</span> EXIT_FAILURE;</div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  }</div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span> }</div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span> </div>
<div class="line"><a name="l00125"></a><span class="lineno"><a class="line" href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab"> 125</a></span> <span class="keywordtype">bool</span> <a class="code" href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab">runTest</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> mat_size)</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span> {</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <span class="keywordtype">bool</span> bResult = <span class="keyword">false</span>;</div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span> </div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  std::vector<NumericT> diagonal(mat_size);</div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  std::vector<NumericT> superdiagonal(mat_size);</div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  std::vector<NumericT> eigenvalues_bisect(mat_size);</div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span> </div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  <span class="comment">// -------------Initialize data-------------------</span></div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="comment">// Fill the diagonal and superdiagonal elements of the vector</span></div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <a class="code" href="tests_2src_2bisect_8cpp.html#a4ff9461d32b6b78b32e2615bd0fb5363">initInputData</a>(diagonal, superdiagonal, mat_size);</div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span> </div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  <span class="comment">// -------Start the bisection algorithm------------</span></div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  std::cout << <span class="stringliteral">"Start the bisection algorithm"</span> << std::endl;</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  std::cout << <span class="stringliteral">"Matrix size: "</span> << mat_size << std::endl;</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  bResult = <a class="code" href="namespaceviennacl_1_1linalg.html#a54d70c731aed90556e228b7f14ac3a52">viennacl::linalg::bisect</a>(diagonal, superdiagonal, eigenvalues_bisect);</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  <span class="comment">// Exit if an error occured during the execution of the algorithm</span></div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  <span class="keywordflow">if</span> (bResult == <span class="keyword">false</span>)</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  <span class="keywordflow">return</span> <span class="keyword">false</span>;</div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span> </div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <span class="comment">// ---------------Check the results---------------</span></div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  <span class="comment">// The results of the bisection algorithm will be checked with the tql algorithm</span></div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  <span class="comment">// Initialize Data for tql1 algorithm</span></div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  std::vector<NumericT> diagonal_tql(mat_size);</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  std::vector<NumericT> superdiagonal_tql(mat_size);</div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  diagonal_tql = diagonal;</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  superdiagonal_tql = superdiagonal;</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span> </div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  <span class="comment">// Start the tql algorithm</span></div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  std::cout << <span class="stringliteral">"Start the tql algorithm..."</span> << std::endl;</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  viennacl::linalg::tql1<NumericT>(mat_size, diagonal_tql, superdiagonal_tql);</div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> </div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="comment">// Ensure that eigenvalues from tql1 algorithm are sorted in ascending order</span></div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  std::sort(diagonal_tql.begin(), diagonal_tql.end());</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span> </div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  <span class="comment">// Compare the results from the bisection algorithm with the results</span></div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  <span class="comment">// from the tql algorithm.</span></div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  std::cout << <span class="stringliteral">"Start comparison..."</span> << std::endl;</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="keywordflow">for</span> (<span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> i = 0; i < mat_size; i++)</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  {</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="keywordflow">if</span> (std::abs(diagonal_tql[i] - eigenvalues_bisect[i]) > <a class="code" href="tests_2src_2bisect_8cpp.html#a6ebf6899d6c1c8b7b9d09be872c05aae">EPS</a>)</div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  {</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  std::cout << std::setprecision(12) << diagonal_tql[i] << <span class="stringliteral">" != "</span> << eigenvalues_bisect[i] << <span class="stringliteral">"\n"</span>;</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  <span class="keywordflow">return</span> <span class="keyword">false</span>;</div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  }</div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  }</div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span> <span class="comment"> // ------------Print the results---------------</span></div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span> <span class="comment"> std::cout << "mat_size = " << mat_size << std::endl;</span></div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span> <span class="comment"> for (unsigned int i = 0; i < mat_size; ++i)</span></div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span> <span class="comment"> {</span></div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span> <span class="comment"> std::cout << "Eigenvalue " << i << ": \tbisect: " << std::setprecision(14) << eigenvalues_bisect[i] << "\ttql: " << diagonal_tql[i] << std::endl;</span></div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span> <span class="comment"> }</span></div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span> <span class="comment">*/</span></div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordflow">return</span> bResult;</div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  </div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span> }</div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_a6ebf6899d6c1c8b7b9d09be872c05aae"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#a6ebf6899d6c1c8b7b9d09be872c05aae">EPS</a></div><div class="ttdeci">#define EPS</div><div class="ttdef"><b>Definition:</b> <a href="tests_2src_2bisect_8cpp_source.html#l00038">bisect.cpp:38</a></div></div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_ae66f6b31b5ad750f1fe042a706a4e3d4"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#ae66f6b31b5ad750f1fe042a706a4e3d4">main</a></div><div class="ttdeci">int main()</div><div class="ttdef"><b>Definition:</b> <a href="tests_2src_2bisect_8cpp_source.html#l00091">bisect.cpp:91</a></div></div>
<div class="ttc" id="tql2_8hpp_html"><div class="ttname"><a href="tql2_8hpp.html">tql2.hpp</a></div><div class="ttdoc">Implementation of the tql2-algorithm for eigenvalue computations. </div></div>
<div class="ttc" id="bisect__gpu_8hpp_html"><div class="ttname"><a href="bisect__gpu_8hpp.html">bisect_gpu.hpp</a></div><div class="ttdoc">Implementation of an bisection algorithm for eigenvalues. </div></div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_a52b5d30a2d7b064678644a3bf49b7f6c"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#a52b5d30a2d7b064678644a3bf49b7f6c">NumericT</a></div><div class="ttdeci">float NumericT</div><div class="ttdef"><b>Definition:</b> <a href="tests_2src_2bisect_8cpp_source.html#l00040">bisect.cpp:40</a></div></div>
<div class="ttc" id="namespaceviennacl_1_1linalg_html_a54d70c731aed90556e228b7f14ac3a52"><div class="ttname"><a href="namespaceviennacl_1_1linalg.html#a54d70c731aed90556e228b7f14ac3a52">viennacl::linalg::bisect</a></div><div class="ttdeci">std::vector< typename viennacl::result_of::cpu_value_type< typename VectorT::value_type >::type > bisect(VectorT const &alphas, VectorT const &betas)</div><div class="ttdoc">Implementation of the bisect-algorithm for the calculation of the eigenvalues of a tridiagonal matrix...</div><div class="ttdef"><b>Definition:</b> <a href="bisect_8hpp_source.html#l00078">bisect.hpp:78</a></div></div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_adc4051332fbb2e2ed482d04b70e3a9ab"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#adc4051332fbb2e2ed482d04b70e3a9ab">runTest</a></div><div class="ttdeci">bool runTest(unsigned int mat_size)</div><div class="ttdoc">Run a simple test. </div><div class="ttdef"><b>Definition:</b> <a href="tests_2src_2bisect_8cpp_source.html#l00125">bisect.cpp:125</a></div></div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_a4ff9461d32b6b78b32e2615bd0fb5363"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#a4ff9461d32b6b78b32e2615bd0fb5363">initInputData</a></div><div class="ttdeci">void initInputData(std::vector< NumericT > &diagonal, std::vector< NumericT > &superdiagonal, unsigned int mat_size)</div><div class="ttdoc">initInputData Initialize the diagonal and superdiagonal elements of the matrix </div><div class="ttdef"><b>Definition:</b> <a href="tests_2src_2bisect_8cpp_source.html#l00053">bisect.cpp:53</a></div></div>
<div class="ttc" id="vector_8hpp_html"><div class="ttname"><a href="vector_8hpp.html">vector.hpp</a></div><div class="ttdoc">The vector type with operator-overloads and proxy classes is defined here. Linear algebra operations ...</div></div>
<div class="ttc" id="tests_2src_2bisect_8cpp_html_a429ac4051f4f4f60df33112d161f8ca6"><div class="ttname"><a href="tests_2src_2bisect_8cpp.html#a429ac4051f4f4f60df33112d161f8ca6">RANDOM_VALUES</a></div><div class="ttdeci">#define RANDOM_VALUES</div></div>
<div class="ttc" id="bisect_8hpp_html"><div class="ttname"><a href="bisect_8hpp.html">bisect.hpp</a></div><div class="ttdoc">Implementation of the algorithm for finding eigenvalues of a tridiagonal matrix. </div></div>
<div class="ttc" id="scalar_8hpp_html"><div class="ttname"><a href="scalar_8hpp.html">scalar.hpp</a></div><div class="ttdoc">Implementation of the ViennaCL scalar class. </div></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
<ul>
<li class="navelem"><a class="el" href="dir_4cce297acab7995b556697564a5c781a.html">tests</a></li><li class="navelem"><a class="el" href="dir_95ddc301de7dff4151ad233d69616cc0.html">src</a></li><li class="navelem"><a class="el" href="tests_2src_2bisect_8cpp.html">bisect.cpp</a></li>
<li class="footer">Generated on Wed Jan 20 2016 22:32:38 for ViennaCL - The Vienna Computing Library by
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.6 </li>
</ul>
</div>
</body>
</html>
|