| <!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.4"/> |
| <title>MADlib: conjugate_gradient.sql_in 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> |
| <script type="text/x-mathjax-config"> |
| MathJax.Hub.Config({ |
| extensions: ["tex2jax.js", "TeX/AMSmath.js", "TeX/AMSsymbols.js"], |
| jax: ["input/TeX","output/HTML-CSS"], |
| }); |
| </script><script src="../mathjax/MathJax.js"></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">MADlib |
|  <span id="projectnumber">0.5</span> |
| </div> |
| <div id="projectbrief">User Documentation</div> |
| </td> |
| </tr> |
| </tbody> |
| </table> |
| </div> |
| <!-- end header part --> |
| <!-- Generated by Doxygen 1.8.4 --> |
| <script type="text/javascript"> |
| var searchBox = new SearchBox("searchBox", "search",false,'Search'); |
| </script> |
| <div id="navrow1" class="tabs"> |
| <ul class="tablist"> |
| <li><a href="index.html"><span>Main Page</span></a></li> |
| <li><a href="modules.html"><span>Modules</span></a></li> |
| <li> |
| <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> |
| </li> |
| </ul> |
| </div> |
| </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('conjugate__gradient_8sql__in_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>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Groups</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">conjugate_gradient.sql_in</div> </div> |
| </div><!--header--> |
| <div class="contents"> |
| <a href="conjugate__gradient_8sql__in.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><span class="comment">/** </span></div> |
| <div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment"> *</span></div> |
| <div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment"> * @file conjugate_gradient.sql_in</span></div> |
| <div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> *</span></div> |
| <div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> * @brief SQL function computing Conjugate Gradient</span></div> |
| <div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> * @date March 2011</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"> *</span></div> |
| <div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> */</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">/**</span></div> |
| <div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">@addtogroup grp_cg</span></div> |
| <div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment">@about</span></div> |
| <div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment">This function uses the iterative conjugate gradient method [1] to find a solution to the function: \f[ \boldsymbol Ax = \boldsymbol b \f]</span></div> |
| <div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment">where \f$ \boldsymbol A \f$ is a symmetric, positive definite matrix and \f$x\f$ and \f$ \boldsymbol b \f$ are vectors. </span></div> |
| <div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment">@input</span></div> |
| <div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment">Matrix \f$ \boldsymbol A \f$ is assumed to be stored in a table where each row consists of at least two columns: array containing values of a given row, row number:</span></div> |
| <div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment"><pre>{TABLE|VIEW} <em>matrix_A</em> (</span></div> |
| <div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> <em>row_number</em> FLOAT,</span></div> |
| <div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment"> <em>row_values</em> FLOAT[],</span></div> |
| <div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment">)</pre></span></div> |
| <div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="comment">The number of elements in each row should be the same.</span></div> |
| <div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="comment">\f$ \boldsymbol b \f$ is passed as a FLOAT[] to the function.</span></div> |
| <div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment">@usage</span></div> |
| <div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment">Conjugate gradient can be called as follows:</span></div> |
| <div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment"><pre>SELECT \ref conjugate_gradient('<em>table_name</em>', </span></div> |
| <div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment"> '<em>name_of_row_values_col</em>', '<em>name_of_row_number_col</em>', '<em>aray_of_b_values</em>', </span></div> |
| <div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="comment"> '<em>desired_precision</em>');</pre></span></div> |
| <div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="comment">Function returns x as an array. </span></div> |
| <div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment"> </span></div> |
| <div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment">@examp</span></div> |
| <div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment">-# Construct matrix A according to structure:</span></div> |
| <div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment">\code</span></div> |
| <div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment">sql> SELECT * FROM data;</span></div> |
| <div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment"> row_num | row_val </span></div> |
| <div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment">---------+---------</span></div> |
| <div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"> 1 | {2,1}</span></div> |
| <div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="comment"> 2 | {1,4}</span></div> |
| <div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="comment">(2 rows)</span></div> |
| <div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment">\endcode</span></div> |
| <div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="comment">-# Call conjugate gradient function:</span></div> |
| <div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="comment">\code</span></div> |
| <div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="comment">sql> SELECT conjugate_gradient('data','row_val','row_num','{2,1}',1E-6,1);</span></div> |
| <div class="line"><a name="l00048"></a><span class="lineno"> 48</span> <span class="comment">INFO: COMPUTE RESIDUAL ERROR 14.5655661859659</span></div> |
| <div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="comment">INFO: ERROR 0.144934004246004</span></div> |
| <div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="comment">INFO: ERROR 3.12963615962926e-31</span></div> |
| <div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="comment">INFO: TEST FINAL ERROR 2.90029642185163e-29</span></div> |
| <div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="comment"> conjugate_gradient </span></div> |
| <div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="comment">---------------------------</span></div> |
| <div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="comment"> {1,-1.31838984174237e-15}</span></div> |
| <div class="line"><a name="l00055"></a><span class="lineno"> 55</span> <span class="comment">(1 row)</span></div> |
| <div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="comment">\endcode</span></div> |
| <div class="line"><a name="l00057"></a><span class="lineno"> 57</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00058"></a><span class="lineno"> 58</span> <span class="comment">@literature</span></div> |
| <div class="line"><a name="l00059"></a><span class="lineno"> 59</span> <span class="comment">[1] "Conjugate gradient method" Wikipedia - http://en.wikipedia.org/wiki/Conjugate_gradient_method</span></div> |
| <div class="line"><a name="l00060"></a><span class="lineno"> 60</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="comment">@sa File conjugate_gradient.sql_in documenting the SQL function.</span></div> |
| <div class="line"><a name="l00062"></a><span class="lineno"> 62</span> <span class="comment">*/</span></div> |
| <div class="line"><a name="l00063"></a><span class="lineno"> 63</span> <span class="comment"></span></div> |
| <div class="line"><a name="l00064"></a><span class="lineno"> 64</span> <span class="comment">/**</span></div> |
| <div class="line"><a name="l00065"></a><span class="lineno"> 65</span> <span class="comment"> * @brief Compute conjugate gradient</span></div> |
| <div class="line"><a name="l00066"></a><span class="lineno"> 66</span> <span class="comment"> * </span></div> |
| <div class="line"><a name="l00067"></a><span class="lineno"> 67</span> <span class="comment"> * @param matrix Name of the table containing argument matrix A</span></div> |
| <div class="line"><a name="l00068"></a><span class="lineno"> 68</span> <span class="comment"> * @param val_id Name of the column contains row values</span></div> |
| <div class="line"><a name="l00069"></a><span class="lineno"> 69</span> <span class="comment"> * @param row_id Name of the column contains row number</span></div> |
| <div class="line"><a name="l00070"></a><span class="lineno"> 70</span> <span class="comment"> * @param b Array containing values of b</span></div> |
| <div class="line"><a name="l00071"></a><span class="lineno"> 71</span> <span class="comment"> * @param precision_limit Precision threshold after which process will terminate</span></div> |
| <div class="line"><a name="l00072"></a><span class="lineno"> 72</span> <span class="comment"> * @param verbosity Verbose flag (0 = false, 1 = true)</span></div> |
| <div class="line"><a name="l00073"></a><span class="lineno"> 73</span> <span class="comment"> * @returns Array containing values of x</span></div> |
| <div class="line"><a name="l00074"></a><span class="lineno"> 74</span> <span class="comment"> *</span></div> |
| <div class="line"><a name="l00075"></a><span class="lineno"> 75</span> <span class="comment"> */</span></div> |
| <div class="line"><a name="l00076"></a><span class="lineno"><a class="code" href="conjugate__gradient_8sql__in.html#ab94e8a12cc9c2d0dac35168721ec8585"> 76</a></span> CREATE OR REPLACE FUNCTION MADLIB_SCHEMA.conjugate_gradient(Matrix TEXT, val_id TEXT, row_id TEXT, b FLOAT[], precision_limit FLOAT, verbosity INT) RETURNS FLOAT[] AS $$</div> |
| <div class="line"><a name="l00077"></a><span class="lineno"> 77</span> declare</div> |
| <div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  r FLOAT[];</div> |
| <div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  p FLOAT[];</div> |
| <div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  x FLOAT[];</div> |
| <div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  k INT;</div> |
| <div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  iter INT = 0;</div> |
| <div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  recidual_refresh INT := 30;</div> |
| <div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  alpha FLOAT;</div> |
| <div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  r_size FLOAT;</div> |
| <div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  r_new_size FLOAT;</div> |
| <div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  Ap FLOAT[];</div> |
| <div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  Ax FLOAT[];</div> |
| <div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  pAp_size FLOAT;</div> |
| <div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  beta FLOAT;</div> |
| <div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  exit_if_no_progess_in INT := 15;</div> |
| <div class="line"><a name="l00092"></a><span class="lineno"> 92</span> begin </div> |
| <div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  DROP TABLE IF EXISTS X_val;</div> |
| <div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  CREATE TEMP TABLE X_val(</div> |
| <div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  value FLOAT[]</div> |
| <div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  );</div> |
| <div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  </div> |
| <div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  DROP TABLE IF EXISTS P_val;</div> |
| <div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  CREATE TEMP TABLE P_val(</div> |
| <div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  value FLOAT[]</div> |
| <div class="line"><a name="l00101"></a><span class="lineno"> 101</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>  SELECT INTO k array_upper(b,1);</div> |
| <div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  --INSERT INTO X_val SELECT ARRAY(SELECT random() FROM generate_series(1, k));</div> |
| <div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  EXECUTE <span class="stringliteral">'INSERT INTO X_val SELECT ARRAY(SELECT random() FROM generate_series(1,'</span> ||k|| <span class="stringliteral">'))'</span>;</div> |
| <div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  LOOP</div> |
| <div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  IF(iter%recidual_refresh = 0)THEN </div> |
| <div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  EXECUTE 'SELECT ARRAY(SELECT MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>('||val_id||', j.x) FROM (SELECT value AS x FROM X_val) AS j, '|| Matrix ||' ORDER BY '||row_id||' LIMIT '|| k ||')' INTO Ax;</div> |
| <div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  SELECT INTO r MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a2875a161a01c7dcdea9a4997b074eefc" title="Subtracts two arrays. It requires that all the values are NON-NULL. Return type is the same as the in...">array_sub</a>(b, Ax);</div> |
| <div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  SELECT INTO r_size MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>(r, r);</div> |
| <div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  IF(verbosity > 0) THEN</div> |
| <div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  RAISE INFO 'COMPUTE RESIDUAL ERROR %', r_size;</div> |
| <div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  END IF;</div> |
| <div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  SELECT INTO p r; </div> |
| <div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  END IF;</div> |
| <div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  iter = iter + 1;</div> |
| <div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  TRUNCATE TABLE P_val;</div> |
| <div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  --INSERT INTO P_val VALUES(p);</div> |
| <div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  EXECUTE 'INSERT INTO P_val(value) VALUES(array['|| array_to_string(p,',') ||'])';</div> |
| <div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  EXECUTE 'SELECT ARRAY(SELECT MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>('||val_id||', j.p) FROM (SELECT value AS p FROM P_val) AS j,'|| Matrix ||' ORDER BY '||row_id||' LIMIT '|| k ||')' INTO Ap;</div> |
| <div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  SELECT INTO pAp_size MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>(p, Ap);</div> |
| <div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  alpha = r_size/pAp_size;</div> |
| <div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  </div> |
| <div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  --SELECT INTO x MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a91c8d3715142927b3967f05a4fbf1575" title="Adds two arrays. It requires that all the values are NON-NULL. Return type is the same as the input t...">array_add</a>(value, MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#ae6881cc5c86941b6ffca35d7f3cd5c12" title="This function takes an array as the input and executes element with multiplication by the scalar prov...">array_scalar_mult</a>(p,alpha)) FROM X_val;</div> |
| <div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  EXECUTE 'SELECT MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a91c8d3715142927b3967f05a4fbf1575" title="Adds two arrays. It requires that all the values are NON-NULL. Return type is the same as the input t...">array_add</a>(value, MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#ae6881cc5c86941b6ffca35d7f3cd5c12" title="This function takes an array as the input and executes element with multiplication by the scalar prov...">array_scalar_mult</a>(array['|| array_to_string(p,',') ||']::<span class="keywordtype">float</span>[],'||alpha||'::<span class="keywordtype">float</span>)) FROM X_val' INTO x;</div> |
| <div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  TRUNCATE TABLE X_val;</div> |
| <div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  --INSERT INTO X_val VALUES(x);</div> |
| <div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  EXECUTE 'INSERT INTO X_val VALUES(array['|| array_to_string(x,',') ||'])';</div> |
| <div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  </div> |
| <div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  SELECT INTO r MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a91c8d3715142927b3967f05a4fbf1575" title="Adds two arrays. It requires that all the values are NON-NULL. Return type is the same as the input t...">array_add</a>(r,MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#ae6881cc5c86941b6ffca35d7f3cd5c12" title="This function takes an array as the input and executes element with multiplication by the scalar prov...">array_scalar_mult</a>(Ap, -alpha));</div> |
| <div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  SELECT INTO r_new_size MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>(r,r);</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>  IF(verbosity > 0) THEN</div> |
| <div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  RAISE INFO 'ERROR %',r_new_size; </div> |
| <div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  END IF;</div> |
| <div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  IF (r_new_size < precision_limit) THEN</div> |
| <div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  EXECUTE 'SELECT ARRAY(SELECT MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>('||val_id||', j.x) FROM (SELECT value AS x FROM X_val) AS j, '|| Matrix ||' ORDER BY '||row_id||' LIMIT '|| k ||')' INTO Ax;</div> |
| <div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  SELECT INTO r MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a2875a161a01c7dcdea9a4997b074eefc" title="Subtracts two arrays. It requires that all the values are NON-NULL. Return type is the same as the in...">array_sub</a>(b, Ax);</div> |
| <div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  SELECT INTO r_new_size MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#acde10964ed23b7c8da515fb84cb8d5e0" title="Dot-product of two arrays. It requires that all the values are NON-NULL. Return type is the same as t...">array_dot</a>(r, r);</div> |
| <div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  IF(verbosity > 0) THEN</div> |
| <div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  RAISE INFO 'TEST FINAL ERROR %', r_new_size;</div> |
| <div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  END IF;</div> |
| <div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  IF (r_new_size < precision_limit) THEN</div> |
| <div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  EXIT;</div> |
| <div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  END IF;</div> |
| <div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  END IF;</div> |
| <div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  SELECT INTO p MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#a91c8d3715142927b3967f05a4fbf1575" title="Adds two arrays. It requires that all the values are NON-NULL. Return type is the same as the input t...">array_add</a>(r, MADLIB_SCHEMA.<a class="code" href="array__ops_8sql__in.html#ae6881cc5c86941b6ffca35d7f3cd5c12" title="This function takes an array as the input and executes element with multiplication by the scalar prov...">array_scalar_mult</a>(p, r_new_size/r_size));</div> |
| <div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  IF(r_size < r_new_size) THEN</div> |
| <div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  exit_if_no_progess_in = exit_if_no_progess_in-1;</div> |
| <div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  RAISE INFO 'No progress! count = %',exit_if_no_progess_in;</div> |
| <div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  </div> |
| <div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  IF(exit_if_no_progess_in <= 0) THEN</div> |
| <div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  RAISE EXCEPTION 'Algorithm failed to converge. Check is input is positive definate.';</div> |
| <div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  END IF;</div> |
| <div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  ELSE</div> |
| <div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  exit_if_no_progess_in = 15;</div> |
| <div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  END IF;</div> |
| <div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  r_size = r_new_size;</div> |
| <div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  END LOOP; </div> |
| <div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  IF(verbosity > 1) THEN</div> |
| <div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  RETURN ARRAY[r_new_size];</div> |
| <div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  END IF;</div> |
| <div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  --SELECT INTO x value FROM X_val;</div> |
| <div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  EXECUTE 'SELECT value FROM X_val' INTO x;</div> |
| <div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  RETURN x;</div> |
| <div class="line"><a name="l00166"></a><span class="lineno"> 166</span> end</div> |
| <div class="line"><a name="l00167"></a><span class="lineno"> 167</span> $$ LANGUAGE plpgsql;</div> |
| <div class="line"><a name="l00168"></a><span class="lineno"> 168</span> </div> |
| <div class="line"><a name="l00169"></a><span class="lineno"> 169</span> CREATE OR REPLACE FUNCTION MADLIB_SCHEMA.<a class="code" href="conjugate__gradient_8sql__in.html#ab94e8a12cc9c2d0dac35168721ec8585" title="Compute conjugate gradient. ">conjugate_gradient</a>(Matrix TEXT, val_id TEXT, row_id TEXT, b FLOAT[], precision_limit FLOAT) RETURNS FLOAT[] AS $$</div> |
| <div class="line"><a name="l00170"></a><span class="lineno"> 170</span> declare</div> |
| <div class="line"><a name="l00171"></a><span class="lineno"> 171</span> begin</div> |
| <div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  RETURN MADLIB_SCHEMA.<a class="code" href="conjugate__gradient_8sql__in.html#ab94e8a12cc9c2d0dac35168721ec8585" title="Compute conjugate gradient. ">conjugate_gradient</a>(Matrix, val_id, row_id, b, precision_limit,0);</div> |
| <div class="line"><a name="l00173"></a><span class="lineno"> 173</span> end</div> |
| <div class="line"><a name="l00174"></a><span class="lineno"> 174</span> $$ LANGUAGE plpgsql;</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_68267d1309a1af8e8297ef4c3efbcdba.html">src</a></li><li class="navelem"><a class="el" href="dir_efbcf68973d247bbf15f9eecae7f24e3.html">ports</a></li><li class="navelem"><a class="el" href="dir_a4a48839224ef8488facbffa8a397967.html">postgres</a></li><li class="navelem"><a class="el" href="dir_dc596537ad427a4d866006d1a3e1fe29.html">modules</a></li><li class="navelem"><a class="el" href="dir_42e4eb27424bda9fbbfa95509de09bad.html">conjugate_gradient</a></li><li class="navelem"><a class="el" href="conjugate__gradient_8sql__in.html">conjugate_gradient.sql_in</a></li> |
| <li class="footer">Generated on Tue Sep 10 2013 13:36:13 for MADlib by |
| <a href="http://www.doxygen.org/index.html"> |
| <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.4 </li> |
| </ul> |
| </div> |
| </body> |
| </html> |
