You cannot select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
CMSIS-DSP/v1.12.0/group__MatrixHouseholder.html

344 lines
17 KiB
HTML

<!-- HTML header for doxygen 1.9.2-->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://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=11"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>CMSIS-DSP: Householder transform of a vector</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>
<script type="text/javascript" src="footer.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(document).ready(function() { init_search(); });
/* @license-end */
</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 type="text/javascript" async="async" src="https://cdn.jsdelivr.net/npm/mathjax@2/MathJax.js"></script>
<link href="$extra_stylesheet" rel="stylesheet" type="text/css" />
<link href="extra_stylesheet.css" rel="stylesheet" type="text/css"/>
<link href="version.css" rel="stylesheet" type="text/css" />
<script type="text/javascript" src="../version.js"></script>
</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: 55px;">
<td id="projectlogo" style="padding: 1.5em;"><img alt="Logo" src="cmsis_logo_white_small.png"/></td>
<td style="padding-left: 1em; padding-bottom: 1em;padding-top: 1em;">
<div id="projectname">CMSIS-DSP
&#160;<span id="projectnumber"><script type="text/javascript">
<!--
writeHeader.call(this);
writeVersionDropdown(this);
//-->
</script>
</span>
</div>
<div id="projectbrief">CMSIS DSP Software Library</div>
</td>
<!--END !PROJECT_NAME-->
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.2 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
var searchBox = new SearchBox("searchBox", "search",'Search','.html');
/* @license-end */
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main&#160;Page</span></a></li>
<li><a href="pages.html"><span>Usage&#160;and&#160;Description</span></a></li>
<li class="current"><a href="modules.html"><span>&#160;API&#160;Reference</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.svg"
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.svg" 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">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(document).ready(function(){initNavTree('group__MatrixHouseholder.html',''); initResizable(); });
/* @license-end */
</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)">
</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="summary">
<a href="#func-members">Functions</a> </div>
<div class="headertitle"><div class="title">Householder transform of a vector<div class="ingroups"><a class="el" href="group__groupMatrix.html">Matrix Functions</a></div></div></div>
</div><!--header-->
<div class="contents">
<table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a id="func-members" name="func-members"></a>
Functions</h2></td></tr>
<tr class="memitem:ga2ac1a8e5e186599ae8ca5a67206c19c4"><td class="memItemLeft" align="right" valign="top">float16_t&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="group__MatrixHouseholder.html#ga2ac1a8e5e186599ae8ca5a67206c19c4">arm_householder_f16</a> (const float16_t *pSrc, const float16_t threshold, uint32_t <a class="el" href="arm__variance__example__f32_8c.html#ab6558f40a619c2502fbc24c880fd4fb0">blockSize</a>, float16_t *pOut)</td></tr>
<tr class="memdesc:ga2ac1a8e5e186599ae8ca5a67206c19c4"><td class="mdescLeft">&#160;</td><td class="mdescRight">Householder transform of a half floating point vector. <a href="group__MatrixHouseholder.html#ga2ac1a8e5e186599ae8ca5a67206c19c4">More...</a><br /></td></tr>
<tr class="separator:ga2ac1a8e5e186599ae8ca5a67206c19c4"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ga4dcb8af2a65875bf84da6ed35c436537"><td class="memItemLeft" align="right" valign="top"><a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="group__MatrixHouseholder.html#ga4dcb8af2a65875bf84da6ed35c436537">arm_householder_f32</a> (const <a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> *pSrc, const <a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> threshold, uint32_t <a class="el" href="arm__variance__example__f32_8c.html#ab6558f40a619c2502fbc24c880fd4fb0">blockSize</a>, <a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> *pOut)</td></tr>
<tr class="memdesc:ga4dcb8af2a65875bf84da6ed35c436537"><td class="mdescLeft">&#160;</td><td class="mdescRight">Householder transform of a floating point vector. <a href="group__MatrixHouseholder.html#ga4dcb8af2a65875bf84da6ed35c436537">More...</a><br /></td></tr>
<tr class="separator:ga4dcb8af2a65875bf84da6ed35c436537"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ga62826e5b9e3c7787147e4e4de21015fc"><td class="memItemLeft" align="right" valign="top"><a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="group__MatrixHouseholder.html#ga62826e5b9e3c7787147e4e4de21015fc">arm_householder_f64</a> (const <a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> *pSrc, const <a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> threshold, uint32_t <a class="el" href="arm__variance__example__f32_8c.html#ab6558f40a619c2502fbc24c880fd4fb0">blockSize</a>, <a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> *pOut)</td></tr>
<tr class="memdesc:ga62826e5b9e3c7787147e4e4de21015fc"><td class="mdescLeft">&#160;</td><td class="mdescRight">Householder transform of a double floating point vector. <a href="group__MatrixHouseholder.html#ga62826e5b9e3c7787147e4e4de21015fc">More...</a><br /></td></tr>
<tr class="separator:ga62826e5b9e3c7787147e4e4de21015fc"><td class="memSeparator" colspan="2">&#160;</td></tr>
</table>
<a name="details" id="details"></a><h2 class="groupheader">Description</h2>
<p >Computes the Householder transform of a vector x.</p>
<p >The Householder transform of x is a vector v with</p>
<p class="formulaDsp">
\[ v_0 = 1 \]
</p>
<p >and a scalar \(\beta\) such that:</p>
<p class="formulaDsp">
\[ P = I - \beta v v^T \]
</p>
<p >is an orthogonal matrix and</p>
<p class="formulaDsp">
\[ P x = ||x||_2 e_1 \]
</p>
<p >So P is an hyperplane reflection such that the image of x is proportional to \(e_1\).</p>
<p >\(e_1\) is the vector of coordinates:</p>
<p class="formulaDsp">
\[ \begin{pmatrix} 1 \\ 0 \\ \vdots \\ \end{pmatrix} \]
</p>
<p >If x is already proportional to \(e_1\) then the matrix P should be the identity.</p>
<p >Thus, \(\beta\) should be 0 and in this case the vector v can also be null.</p>
<p >But how do we detect that x is already proportional to \(e_1\).</p>
<p >If x </p><p class="formulaDsp">
\[ x = \begin{pmatrix} x_0 \\ xr \\ \end{pmatrix} \]
</p>
<p >where \(xr\) is a vector.</p>
<p >The algorithm is computing the norm squared of this vector:</p>
<p class="formulaDsp">
\[ ||xr||^2 \]
</p>
<p >and this value is compared to a <code>threshold</code>. If the value is smaller than the <code>threshold</code>, the algorithm is returning 0 for \(\beta\) and the householder vector.</p>
<p >This <code>threshold</code> is an argument of the function.</p>
<p >Default values are provided in the header <code><a class="el" href="matrix__functions_8h.html">dsp/matrix_functions.h</a></code> like for instance <code>DEFAULT_HOUSEHOLDER_THRESHOLD_F32</code> </p>
<h2 class="groupheader">Function Documentation</h2>
<a id="ga2ac1a8e5e186599ae8ca5a67206c19c4" name="ga2ac1a8e5e186599ae8ca5a67206c19c4"></a>
<h2 class="memtitle"><span class="permalink"><a href="#ga2ac1a8e5e186599ae8ca5a67206c19c4">&#9670;&nbsp;</a></span>arm_householder_f16()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname">float16_t arm_householder_f16 </td>
<td>(</td>
<td class="paramtype">const float16_t *&#160;</td>
<td class="paramname"><em>pSrc</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const float16_t&#160;</td>
<td class="paramname"><em>threshold</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">uint32_t&#160;</td>
<td class="paramname"><em>blockSize</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">float16_t *&#160;</td>
<td class="paramname"><em>pOut</em>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
</tr>
</table>
</div><div class="memdoc">
<dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">pSrc</td><td>points to the input vector. </td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">threshold</td><td>norm2 threshold. <br />
</td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">blockSize</td><td>dimension of the vector space. </td></tr>
<tr><td class="paramdir">[out]</td><td class="paramname">pOut</td><td>points to the output vector. </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>beta return the scaling factor beta </dd></dl>
</div>
</div>
<a id="ga4dcb8af2a65875bf84da6ed35c436537" name="ga4dcb8af2a65875bf84da6ed35c436537"></a>
<h2 class="memtitle"><span class="permalink"><a href="#ga4dcb8af2a65875bf84da6ed35c436537">&#9670;&nbsp;</a></span>arm_householder_f32()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> arm_householder_f32 </td>
<td>(</td>
<td class="paramtype">const <a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> *&#160;</td>
<td class="paramname"><em>pSrc</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a>&#160;</td>
<td class="paramname"><em>threshold</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">uint32_t&#160;</td>
<td class="paramname"><em>blockSize</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype"><a class="el" href="arm__math__types_8h.html#a4611b605e45ab401f02cab15c5e38715">float32_t</a> *&#160;</td>
<td class="paramname"><em>pOut</em>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
</tr>
</table>
</div><div class="memdoc">
<dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">pSrc</td><td>points to the input vector. </td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">threshold</td><td>norm2 threshold. <br />
</td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">blockSize</td><td>dimension of the vector space. </td></tr>
<tr><td class="paramdir">[out]</td><td class="paramname">pOut</td><td>points to the output vector. </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>beta return the scaling factor beta </dd></dl>
</div>
</div>
<a id="ga62826e5b9e3c7787147e4e4de21015fc" name="ga62826e5b9e3c7787147e4e4de21015fc"></a>
<h2 class="memtitle"><span class="permalink"><a href="#ga62826e5b9e3c7787147e4e4de21015fc">&#9670;&nbsp;</a></span>arm_householder_f64()</h2>
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> arm_householder_f64 </td>
<td>(</td>
<td class="paramtype">const <a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> *&#160;</td>
<td class="paramname"><em>pSrc</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a>&#160;</td>
<td class="paramname"><em>threshold</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">uint32_t&#160;</td>
<td class="paramname"><em>blockSize</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype"><a class="el" href="arm__math__types_8h.html#ac55f3ae81b5bc9053760baacf57e47f4">float64_t</a> *&#160;</td>
<td class="paramname"><em>pOut</em>&#160;</td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td></td>
</tr>
</table>
</div><div class="memdoc">
<dl class="params"><dt>Parameters</dt><dd>
<table class="params">
<tr><td class="paramdir">[in]</td><td class="paramname">pSrc</td><td>points to the input vector. </td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">threshold</td><td>norm2 threshold. <br />
</td></tr>
<tr><td class="paramdir">[in]</td><td class="paramname">blockSize</td><td>dimension of the vector space. </td></tr>
<tr><td class="paramdir">[out]</td><td class="paramname">pOut</td><td>points to the output vector. </td></tr>
</table>
</dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>beta return the scaling factor beta </dd></dl>
</div>
</div>
</div><!-- contents -->
</div><!-- doc-content -->
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
<ul>
<li class="footer">
<script type="text/javascript">
<!--
writeFooter.call(this);
//-->
</script>
</li>
</ul>
</div>
</body>
</html>