git.gag.com Git - fw/stlink/blob - exampleF4/CMSIS/Documentation/DSP_Lib/html/arm_matrix_example_f32_8c-example.html

   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">
   3 <head>
   4 <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
   5 <title>CMSIS DSP Software Library: arm_matrix_example_f32.c</title>
   6 <link href="tabs.css" rel="stylesheet" type="text/css"/>
   7 <link href="search/search.css" rel="stylesheet" type="text/css"/>
   8 <script type="text/javaScript" src="search/search.js"></script>
   9 <link href="doxygen.css" rel="stylesheet" type="text/css"/>
  10 </head>
  11 <body onload='searchBox.OnSelectItem(0);'>
  12 <!-- Generated by Doxygen 1.7.2 -->
  13 <script type="text/javascript"><!--
  14 var searchBox = new SearchBox("searchBox", "search",false,'Search');
  15 --></script>
  16 <div class="navigation" id="top">
  17   <div class="tabs">
  18     <ul class="tablist">
  19       <li><a href="index.html"><span>Main&#160;Page</span></a></li>
  20       <li><a href="modules.html"><span>Modules</span></a></li>
  21       <li><a href="annotated.html"><span>Data&#160;Structures</span></a></li>
  22       <li><a href="files.html"><span>Files</span></a></li>
  23       <li><a href="examples.html"><span>Examples</span></a></li>
  24       <li id="searchli">
  25         <div id="MSearchBox" class="MSearchBoxInactive">
  26         <span class="left">
  27           <img id="MSearchSelect" src="search/mag_sel.png"
  28                onmouseover="return searchBox.OnSearchSelectShow()"
  29                onmouseout="return searchBox.OnSearchSelectHide()"
  30                alt=""/>
  31           <input type="text" id="MSearchField" value="Search" accesskey="S"
  32                onfocus="searchBox.OnSearchFieldFocus(true)"
  33                onblur="searchBox.OnSearchFieldFocus(false)"
  34                onkeyup="searchBox.OnSearchFieldChange(event)"/>
  35           </span><span class="right">
  36             <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
  37           </span>
  38         </div>
  39       </li>
  40     </ul>
  41   </div>
  42 </div>
  43 <div class="header">
  44   <div class="headertitle">
  45 <h1>arm_matrix_example_f32.c</h1>  </div>
  46 </div>
  47 <div class="contents">
  48 <div class="fragment"><pre class="fragment"><span class="comment">/* ---------------------------------------------------------------------- </span>
  49 <span class="comment">* Copyright (C) 2010 ARM Limited. All rights reserved.   </span>
  50 <span class="comment">*  </span>
  51 <span class="comment">* $Date:        29. November 2010  </span>
  52 <span class="comment">* $Revision:    V1.0.3</span>
  53 <span class="comment">*  </span>
  54 <span class="comment">* Project:      CMSIS DSP Library  </span>
  55 <span class="comment">* Title:        arm_matrix_example_f32.c          </span>
  56 <span class="comment">*  </span>
  57 <span class="comment">* Description:  Example code demonstrating least square fit to data  </span>
  58 <span class="comment">*               using matrix functions  </span>
  59 <span class="comment">*                </span>
  60 <span class="comment">* Target Processor: Cortex-M4/Cortex-M3  </span>
  61 <span class="comment">*</span>
  62 <span class="comment">*</span>
  63 <span class="comment">* Version 1.0.3 2010/11/29 </span>
  64 <span class="comment">*    Re-organized the CMSIS folders and updated documentation. </span>
  65 <span class="comment">* </span>
  66 <span class="comment">* Version 1.0.1 2010/10/05 KK </span>
  67 <span class="comment">*    Production release and review comments incorporated.  </span>
  68 <span class="comment">*</span>
  69 <span class="comment">* Version 1.0.0 2010/09/20 KK</span>
  70 <span class="comment">*    Production release and review comments incorporated.</span>
  71 <span class="comment">* ------------------------------------------------------------------- */</span>
  72
  73 <span class="preprocessor">#include &quot;<a class="code" href="arm__math_8h.html">arm_math.h</a>&quot;</span>
  74 <span class="preprocessor">#include &quot;math_helper.h&quot;</span>
  75
  76 <span class="preprocessor">#define SNR_THRESHOLD   90 </span>
  77 <span class="preprocessor"></span>
  78 <span class="comment">/* -------------------------------------------------------------------------------- </span>
  79 <span class="comment">* Test input data(Cycles) taken from FIR Q15 module for differant cases of blockSize  </span>
  80 <span class="comment">* and tapSize </span>
  81 <span class="comment">* --------------------------------------------------------------------------------- */</span>
  82
  83 <span class="keyword">const</span> <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a0"></a><a class="code" href="arm__matrix__example__f32_8c.html#a974d5f0aace6a99e29ca767907fb3b9f">B_f32</a>[4] =
  84 {
  85     782.0, 7577.0, 470.0, 4505.0
  86 };
  87
  88 <span class="comment">/* -------------------------------------------------------------------------------- </span>
  89 <span class="comment">* Formula to fit is  C1 + C2 * numTaps + C3 * blockSize + C4 * numTaps * blockSize </span>
  90 <span class="comment">* -------------------------------------------------------------------------------- */</span>
  91
  92 <span class="keyword">const</span> <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a1"></a><a class="code" href="arm__matrix__example__f32_8c.html#aed27b92d9847194d9dcce40cecf2b48a">A_f32</a>[16] =
  93 {
  94     <span class="comment">/* Const,   numTaps,    blockSize,  numTaps*blockSize */</span>
  95     1.0,        32.0,       4.0,        128.0,
  96     1.0,        32.0,       64.0,       2048.0,
  97     1.0,        16.0,       4.0,        64.0,
  98     1.0,        16.0,       64.0,       1024.0,
  99 };
 100
 101
 102 <span class="comment">/* ---------------------------------------------------------------------- </span>
 103 <span class="comment">* Temporary buffers  for storing intermediate values </span>
 104 <span class="comment">* ------------------------------------------------------------------- */</span>
 105 <span class="comment">/* Transpose of A Buffer */</span>
 106 <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a2"></a><a class="code" href="arm__matrix__example__f32_8c.html#a46dc2aa6dfc692af7b4a1379d7329ccd">AT_f32</a>[16];
 107 <span class="comment">/* (Transpose of A * A) Buffer */</span>
 108 <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a3"></a><a class="code" href="arm__matrix__example__f32_8c.html#a867497c6bf86014513bf2ad3551aa896">ATMA_f32</a>[16];
 109 <span class="comment">/* Inverse(Transpose of A * A)  Buffer */</span>
 110 <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a4"></a><a class="code" href="arm__matrix__example__f32_8c.html#a44425c149c52b326a3b7a77676686f00">ATMAI_f32</a>[16];
 111 <span class="comment">/* Test Output Buffer */</span>
 112 <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a5"></a><a class="code" href="arm__matrix__example__f32_8c.html#a98c67c0fc0cb5f2df51b21482d31d21c">X_f32</a>[4];
 113
 114 <span class="comment">/* ---------------------------------------------------------------------- </span>
 115 <span class="comment">* Reference ouput buffer C1, C2, C3 and C4 taken from MATLAB  </span>
 116 <span class="comment">* ------------------------------------------------------------------- */</span>
 117 <span class="keyword">const</span> <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a6"></a><a class="code" href="arm__matrix__example__f32_8c.html#a6184758419722fa16bb883097c2f596b">xRef_f32</a>[4] = {73.0, 8.0, 21.25, 2.875};
 118
 119 <a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> <a name="a7"></a><a class="code" href="arm__convolution__example__f32_8c.html#af06013f588a7003278de222913c9d819">snr</a>;
 120
 121
 122 <span class="comment">/* ---------------------------------------------------------------------- </span>
 123 <span class="comment">* Max magnitude FFT Bin test </span>
 124 <span class="comment">* ------------------------------------------------------------------- */</span>
 125
 126 int32_t <a name="a8"></a><a class="code" href="arm__class__marks__example__f32_8c.html#a196718f834091385d38586a0ce4009dc">main</a>(<span class="keywordtype">void</span>)
 127 {
 128
 129     <a name="_a9"></a><a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> A;      <span class="comment">/* Matrix A Instance */</span>
 130     <a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> AT;     <span class="comment">/* Matrix AT(A transpose) instance */</span>
 131     <a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> ATMA;   <span class="comment">/* Matrix ATMA( AT multiply with A) instance */</span>
 132     <a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> ATMAI;  <span class="comment">/* Matrix ATMAI(Inverse of ATMA) instance */</span>
 133     <a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> B;      <span class="comment">/* Matrix B instance */</span>
 134     <a class="code" href="structarm__matrix__instance__f32.html" title="Instance structure for the floating-point matrix structure.">arm_matrix_instance_f32</a> X;      <span class="comment">/* Matrix X(Unknown Matrix) instance */</span>
 135
 136     uint32_t srcRows, srcColumns;   <span class="comment">/* Temporary variables */</span>
 137     <a class="code" href="arm__math_8h.html#a5e459c6409dfcd2927bb8a57491d7cf6" title="Error status returned by some functions in the library.">arm_status</a> <a name="a10"></a><a class="code" href="arm__dotproduct__example__f32_8c.html#a88ccb294236ab22b00310c47164c53c3">status</a>;
 138
 139     <span class="comment">/* Initialise A Matrix Instance with numRows, numCols and data array(A_f32) */</span>
 140     srcRows = 4;
 141     srcColumns = 4;
 142     <a name="a11"></a><a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;A, srcRows, srcColumns, (<a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> *)<a class="code" href="arm__matrix__example__f32_8c.html#aed27b92d9847194d9dcce40cecf2b48a">A_f32</a>);
 143
 144     <span class="comment">/* Initialise Matrix Instance AT with numRows, numCols and data array(AT_f32) */</span>
 145     srcRows = 4;
 146     srcColumns = 4;
 147     <a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;AT, srcRows, srcColumns, <a class="code" href="arm__matrix__example__f32_8c.html#a46dc2aa6dfc692af7b4a1379d7329ccd">AT_f32</a>);
 148
 149     <span class="comment">/* calculation of A transpose */</span>
 150     status = <a name="a12"></a><a class="code" href="group___matrix_trans.html#gad7dd9f108429da13d3864696ceeec789" title="Floating-point matrix transpose.">arm_mat_trans_f32</a>(&amp;A, &amp;AT);
 151
 152
 153     <span class="comment">/* Initialise ATMA Matrix Instance with numRows, numCols and data array(ATMA_f32) */</span>
 154     srcRows = 4;
 155     srcColumns = 4;
 156     <a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;ATMA, srcRows, srcColumns, <a class="code" href="arm__matrix__example__f32_8c.html#a867497c6bf86014513bf2ad3551aa896">ATMA_f32</a>);
 157
 158     <span class="comment">/* calculation of AT Multiply with A */</span>
 159     status = <a name="a13"></a><a class="code" href="group___matrix_mult.html#ga917bf0270310c1d3f0eda1fc7c0026a0" title="Floating-point matrix multiplication.">arm_mat_mult_f32</a>(&amp;AT, &amp;A, &amp;ATMA);
 160
 161     <span class="comment">/* Initialise ATMAI Matrix Instance with numRows, numCols and data array(ATMAI_f32) */</span>
 162     srcRows = 4;
 163     srcColumns = 4;
 164     <a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;ATMAI, srcRows, srcColumns, <a class="code" href="arm__matrix__example__f32_8c.html#a44425c149c52b326a3b7a77676686f00">ATMAI_f32</a>);
 165
 166     <span class="comment">/* calculation of Inverse((Transpose(A) * A) */</span>
 167     status = <a name="a14"></a><a class="code" href="group___matrix_inv.html#ga542be7aabbf7a2297a4b62cf212910e3" title="Floating-point matrix inverse.">arm_mat_inverse_f32</a>(&amp;ATMA, &amp;ATMAI);
 168
 169     <span class="comment">/* calculation of (Inverse((Transpose(A) * A)) *  Transpose(A)) */</span>
 170     status = <a class="code" href="group___matrix_mult.html#ga917bf0270310c1d3f0eda1fc7c0026a0" title="Floating-point matrix multiplication.">arm_mat_mult_f32</a>(&amp;ATMAI, &amp;AT, &amp;ATMA);
 171
 172     <span class="comment">/* Initialise B Matrix Instance with numRows, numCols and data array(B_f32) */</span>
 173     srcRows = 4;
 174     srcColumns = 1;
 175     <a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;B, srcRows, srcColumns, (<a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> *)<a class="code" href="arm__matrix__example__f32_8c.html#a974d5f0aace6a99e29ca767907fb3b9f">B_f32</a>);
 176
 177     <span class="comment">/* Initialise X Matrix Instance with numRows, numCols and data array(X_f32) */</span>
 178     srcRows = 4;
 179     srcColumns = 1;
 180     <a class="code" href="group___matrix_init.html#ga11e3dc41592a6401c13182fef9416a27" title="Floating-point matrix initialization.">arm_mat_init_f32</a>(&amp;X, srcRows, srcColumns, <a class="code" href="arm__matrix__example__f32_8c.html#a98c67c0fc0cb5f2df51b21482d31d21c">X_f32</a>);
 181
 182     <span class="comment">/* calculation ((Inverse((Transpose(A) * A)) *  Transpose(A)) * B) */</span>
 183     status = <a class="code" href="group___matrix_mult.html#ga917bf0270310c1d3f0eda1fc7c0026a0" title="Floating-point matrix multiplication.">arm_mat_mult_f32</a>(&amp;ATMA, &amp;B, &amp;X);
 184
 185     <span class="comment">/* Comparison of reference with test output */</span>
 186     <a class="code" href="arm__convolution__example__f32_8c.html#af06013f588a7003278de222913c9d819">snr</a> = arm_snr_f32((<a class="code" href="arm__math_8h.html#a4611b605e45ab401f02cab15c5e38715" title="32-bit floating-point type definition.">float32_t</a> *)<a class="code" href="arm__matrix__example__f32_8c.html#a6184758419722fa16bb883097c2f596b">xRef_f32</a>, <a class="code" href="arm__matrix__example__f32_8c.html#a98c67c0fc0cb5f2df51b21482d31d21c">X_f32</a>, 4);
 187
 188     <span class="comment">/*------------------------------------------------------------------------------ </span>
 189 <span class="comment">    *                   Initialise status depending on SNR calculations </span>
 190 <span class="comment">    *------------------------------------------------------------------------------*/</span>
 191     <span class="keywordflow">if</span>( <a class="code" href="arm__convolution__example__f32_8c.html#af06013f588a7003278de222913c9d819">snr</a> &gt; <a name="a15"></a><a class="code" href="arm__convolution__example__f32_8c.html#af08ec3fef897d77c6817638bf0e0c5c6">SNR_THRESHOLD</a>)
 192     {
 193         status = <a name="a16"></a><a class="code" href="arm__math_8h.html#a5e459c6409dfcd2927bb8a57491d7cf6a9f8b2a10bd827fb4600e77d455902eb0">ARM_MATH_SUCCESS</a>;
 194     }
 195     <span class="keywordflow">else</span>
 196     {
 197         status = <a name="a17"></a><a class="code" href="arm__math_8h.html#a5e459c6409dfcd2927bb8a57491d7cf6a09457f2be656b35015fd6d36202fa376">ARM_MATH_TEST_FAILURE</a>;
 198     }
 199
 200
 201     <span class="comment">/* ---------------------------------------------------------------------- </span>
 202 <span class="comment">    ** Loop here if the signals fail the PASS check. </span>
 203 <span class="comment">    ** This denotes a test failure </span>
 204 <span class="comment">    ** ------------------------------------------------------------------- */</span>
 205     <span class="keywordflow">if</span>( status != <a class="code" href="arm__math_8h.html#a5e459c6409dfcd2927bb8a57491d7cf6a9f8b2a10bd827fb4600e77d455902eb0">ARM_MATH_SUCCESS</a>)
 206     {
 207       <span class="keywordflow">while</span>(1);
 208     }
 209 }
 210
 211 </pre></div> </div>
 212 </div>
 213 <!--- window showing the filter options -->
 214 <div id="MSearchSelectWindow"
 215      onmouseover="return searchBox.OnSearchSelectShow()"
 216      onmouseout="return searchBox.OnSearchSelectHide()"
 217      onkeydown="return searchBox.OnSearchSelectKey(event)">
 218 <a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark">&#160;</span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark">&#160;</span>Data Structures</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark">&#160;</span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark">&#160;</span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark">&#160;</span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark">&#160;</span>Typedefs</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark">&#160;</span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark">&#160;</span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark">&#160;</span>Defines</a></div>
 219
 220 <!-- iframe showing the search results (closed by default) -->
 221 <div id="MSearchResultsWindow">
 222 <iframe src="" frameborder="0"
 223         name="MSearchResults" id="MSearchResults">
 224 </iframe>
 225 </div>
 226
 227 <hr class="footer"/><address class="footer"><small>Generated on Fri Jul 15 2011 13:16:15 for CMSIS DSP Software Library by&#160;
 228 <a href="http://www.doxygen.org/index.html">
 229 <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.2 </small></address>
 230 </body>
 231 </html>