<!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.8"/>
<title>SphinxBase: src/libsphinxbase/feat/lda.c 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="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">SphinxBase
 <span id="projectnumber">0.6</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.8 -->
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>Globals</span></a></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('lda_8c_source.html','');});
</script>
<div id="doc-content">
<div class="header">
<div class="headertitle">
<div class="title">lda.c</div> </div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/* -*- c-basic-offset: 4; indent-tabs-mode: nil -*- */</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"> * Copyright (c) 2006 Carnegie Mellon University. All rights</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment"> * reserved.</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> * Redistribution and use in source and binary forms, with or without</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> * modification, are permitted provided that the following conditions</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> * are met:</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"> * 1. Redistributions of source code must retain the above copyright</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> * notice, this list of conditions and the following disclaimer. </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"> * 2. Redistributions in binary form must reproduce the above copyright</span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"> * notice, this list of conditions and the following disclaimer in</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> * the documentation and/or other materials provided with the</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment"> * distribution.</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"> * This work was supported in part by funding from the Defense Advanced </span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment"> * Research Projects Agency and the National Science Foundation of the </span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment"> * United States of America, and the CMU Sphinx Speech Consortium.</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment"> * THIS SOFTWARE IS PROVIDED BY CARNEGIE MELLON UNIVERSITY ``AS IS'' AND </span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment"> * ANY EXPRESSED OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, </span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="comment"> * THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR</span></div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="comment"> * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL CARNEGIE MELLON UNIVERSITY</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="comment"> * NOR ITS EMPLOYEES BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,</span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment"> * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT </span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment"> * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, </span></div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment"> * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY </span></div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment"> * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT </span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment"> * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE </span></div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="comment"> * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.</span></div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="comment"> *</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"> *</span></div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment"> * lda.c -- Read and apply LDA matrices to features.</span></div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment"> * Author: David Huggins-Daines <dhuggins@cs.cmu.edu></span></div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> </div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="preprocessor">#include <assert.h></span></div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="preprocessor">#include <string.h></span></div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="preprocessor">#ifdef HAVE_CONFIG_H</span></div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="preprocessor">#include <config.h></span></div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span> </div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="preprocessor">#ifdef _MSC_VER</span></div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="preprocessor">#pragma warning (disable: 4018)</span></div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> </div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="preprocessor">#include "<a class="code" href="feat_8h.html">sphinxbase/feat.h</a>"</span></div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="preprocessor">#include "<a class="code" href="ckd__alloc_8h.html">sphinxbase/ckd_alloc.h</a>"</span></div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span> <span class="preprocessor">#include "<a class="code" href="bio_8h.html">sphinxbase/bio.h</a>"</span></div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="preprocessor">#include "<a class="code" href="err_8h.html">sphinxbase/err.h</a>"</span></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 MATRIX_FILE_VERSION "0.1"</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> int32</div>
<div class="line"><a name="l00061"></a><span class="lineno"><a class="line" href="feat_8h.html#a7639fbc34e3328d7e3259e61372c64b8"> 61</a></span> <a class="code" href="feat_8h.html#a7639fbc34e3328d7e3259e61372c64b8">feat_read_lda</a>(<a class="code" href="structfeat__t.html">feat_t</a> *feat, <span class="keyword">const</span> <span class="keywordtype">char</span> *ldafile, int32 dim)</div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span> {</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  FILE *fh;</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  int32 byteswap;</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  uint32 chksum, i, m, n;</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <span class="keywordtype">char</span> **argname, **argval;</div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span> </div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  assert(feat);</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  <span class="keywordflow">if</span> (feat->n_stream != 1) {</div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  <a class="code" href="err_8h.html#a5f7b2f58f5a663a6bdd51f197ae21993">E_ERROR</a>(<span class="stringliteral">"LDA incompatible with multi-stream features (n_stream = %d)\n"</span>,</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  feat->n_stream);</div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="keywordflow">return</span> -1;</div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  }</div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span> </div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  <span class="keywordflow">if</span> ((fh = fopen(ldafile, <span class="stringliteral">"rb"</span>)) == NULL) {</div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  <a class="code" href="err_8h.html#a54ffbfe898d74595c586a1f48f32ef03">E_ERROR_SYSTEM</a>(<span class="stringliteral">"Failed to open transform file '%s' for reading"</span>, ldafile);</div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span>  <span class="keywordflow">return</span> -1;</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  }</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>  <span class="keywordflow">if</span> (<a class="code" href="bio_8h.html#a7a155ff51740f1d9a31f7bb7b6e4bc3a">bio_readhdr</a>(fh, &argname, &argval, &byteswap) < 0) {</div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <a class="code" href="err_8h.html#a5f7b2f58f5a663a6bdd51f197ae21993">E_ERROR</a>(<span class="stringliteral">"Failed to read header from transform file '%s'\n"</span>, ldafile);</div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  fclose(fh);</div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="keywordflow">return</span> -1;</div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  }</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>  <span class="keywordflow">for</span> (i = 0; argname[i]; i++) {</div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <span class="keywordflow">if</span> (strcmp(argname[i], <span class="stringliteral">"version"</span>) == 0) {</div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="keywordflow">if</span> (strcmp(argval[i], MATRIX_FILE_VERSION) != 0)</div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  <a class="code" href="err_8h.html#a6a794bec721b555ac1f2167f9e12f662">E_WARN</a>(<span class="stringliteral">"%s: Version mismatch: %s, expecting %s\n"</span>,</div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  ldafile, argval[i], MATRIX_FILE_VERSION);</div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  }</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> </div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <a class="code" href="bio_8h.html#adb3c858a50381aeff022ccd4271af2d9">bio_hdrarg_free</a>(argname, argval);</div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  argname = argval = NULL;</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>  chksum = 0;</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>  <span class="keywordflow">if</span> (feat->lda)</div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <a class="code" href="ckd__alloc_8h.html#a443accde4663a4a07c6d147088614d57">ckd_free_3d</a>((<span class="keywordtype">void</span> ***)feat->lda);</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>  <span class="comment">/* Use a temporary variable to avoid strict-aliasing problems. */</span></div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  <span class="keywordtype">void</span> ***outlda;</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>  <span class="keywordflow">if</span> (<a class="code" href="bio_8h.html#aca6a6b71d02b9e5ca6b136a28640e0c1">bio_fread_3d</a>(&outlda, <span class="keyword">sizeof</span>(float32),</div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  &feat->n_lda, &m, &n,</div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  fh, byteswap, &chksum) < 0) {</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  <a class="code" href="err_8h.html#a54ffbfe898d74595c586a1f48f32ef03">E_ERROR_SYSTEM</a>(<span class="stringliteral">"%s: bio_fread_3d(lda) failed\n"</span>, ldafile);</div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  fclose(fh);</div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  <span class="keywordflow">return</span> -1;</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  }</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  feat->lda = (<span class="keywordtype">void</span> *)outlda;</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>  fclose(fh);</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> <span class="preprocessor">#ifdef FIXED_POINT</span></div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="comment">/* FIXME: This is a fragile hack that depends on mfcc_t and</span></div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span> <span class="comment"> * float32 being the same size (which they are, but...) */</span></div>
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <span class="keywordflow">for</span> (i = 0; i < feat->n_lda * m * n; ++i) {</div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  feat->lda[0][0][i] = FLOAT2MFCC(((<span class="keywordtype">float</span> *)feat->lda[0][0])[i]);</div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  }</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  <span class="comment">/* Note that SphinxTrain stores the eigenvectors as row vectors. */</span></div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  <span class="keywordflow">if</span> (n != feat->stream_len[0])</div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  <a class="code" href="err_8h.html#a1a4495946ab2449d61108fe829a94613">E_FATAL</a>(<span class="stringliteral">"LDA matrix dimension %d doesn't match feature stream size %d\n"</span>, n, feat->stream_len[0]);</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>  <span class="comment">/* Override dim from file if it is 0 or greater than m. */</span></div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <span class="keywordflow">if</span> (dim > m || dim <= 0) {</div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  dim = m;</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>  feat->out_dim = dim;</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span> </div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="keywordflow">return</span> 0;</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> </div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span> <span class="keywordtype">void</span></div>
<div class="line"><a name="l00139"></a><span class="lineno"><a class="line" href="feat_8h.html#a7452702cb435ea9244c7e1f61bbf5da1"> 139</a></span> <a class="code" href="feat_8h.html#a7452702cb435ea9244c7e1f61bbf5da1">feat_lda_transform</a>(<a class="code" href="structfeat__t.html">feat_t</a> *fcb, mfcc_t ***inout_feat, uint32 nfr)</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span> {</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  mfcc_t *tmp;</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  uint32 i, j, k;</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span> </div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  tmp = <a class="code" href="ckd__alloc_8h.html#aa00ef21903bc4f8a972488417adc8d2e">ckd_calloc</a>(fcb->stream_len[0], <span class="keyword">sizeof</span>(mfcc_t));</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <span class="keywordflow">for</span> (i = 0; i < nfr; ++i) {</div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  <span class="comment">/* Do the matrix multiplication inline here since fcb->lda</span></div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span> <span class="comment"> * is transposed (eigenvectors in rows not columns). */</span></div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  <span class="comment">/* FIXME: In the future we ought to use the BLAS. */</span></div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  memset(tmp, 0, <span class="keyword">sizeof</span>(mfcc_t) * fcb->stream_len[0]);</div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="keywordflow">for</span> (j = 0; j < <a class="code" href="feat_8h.html#a66dc16aeea4307bf36ec4e08e05d0c7c">feat_dimension</a>(fcb); ++j) {</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  <span class="keywordflow">for</span> (k = 0; k < fcb->stream_len[0]; ++k) {</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  tmp[j] += MFCCMUL(inout_feat[i][0][k], fcb->lda[0][j][k]);</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  }</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  }</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  memcpy(inout_feat[i][0], tmp, fcb->stream_len[0] * <span class="keyword">sizeof</span>(mfcc_t));</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>  <a class="code" href="ckd__alloc_8h.html#a31c6b405558620ac37599737b5722fbf">ckd_free</a>(tmp);</div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span> }</div>
<div class="ttc" id="ckd__alloc_8h_html_aa00ef21903bc4f8a972488417adc8d2e"><div class="ttname"><a href="ckd__alloc_8h.html#aa00ef21903bc4f8a972488417adc8d2e">ckd_calloc</a></div><div class="ttdeci">#define ckd_calloc(n, sz)</div><div class="ttdoc">Macros to simplify the use of above functions. </div><div class="ttdef"><b>Definition:</b> <a href="ckd__alloc_8h_source.html#l00248">ckd_alloc.h:248</a></div></div>
<div class="ttc" id="feat_8h_html_a7639fbc34e3328d7e3259e61372c64b8"><div class="ttname"><a href="feat_8h.html#a7639fbc34e3328d7e3259e61372c64b8">feat_read_lda</a></div><div class="ttdeci">SPHINXBASE_EXPORT int32 feat_read_lda(feat_t *feat, const char *ldafile, int32 dim)</div><div class="ttdoc">Add an LDA transformation to the feature module from a file. </div><div class="ttdef"><b>Definition:</b> <a href="lda_8c_source.html#l00061">lda.c:61</a></div></div>
<div class="ttc" id="feat_8h_html_a66dc16aeea4307bf36ec4e08e05d0c7c"><div class="ttname"><a href="feat_8h.html#a66dc16aeea4307bf36ec4e08e05d0c7c">feat_dimension</a></div><div class="ttdeci">#define feat_dimension(f)</div><div class="ttdoc">Total dimensionality of feature output. </div><div class="ttdef"><b>Definition:</b> <a href="feat_8h_source.html#l00237">feat.h:237</a></div></div>
<div class="ttc" id="ckd__alloc_8h_html"><div class="ttname"><a href="ckd__alloc_8h.html">ckd_alloc.h</a></div><div class="ttdoc">Sphinx's memory allocation/deallocation routines. </div></div>
<div class="ttc" id="bio_8h_html"><div class="ttname"><a href="bio_8h.html">bio.h</a></div><div class="ttdoc">Cross platform binary IO to process files in sphinx3 format. </div></div>
<div class="ttc" id="bio_8h_html_aca6a6b71d02b9e5ca6b136a28640e0c1"><div class="ttname"><a href="bio_8h.html#aca6a6b71d02b9e5ca6b136a28640e0c1">bio_fread_3d</a></div><div class="ttdeci">SPHINXBASE_EXPORT int32 bio_fread_3d(void ****arr, size_t e_sz, uint32 *d1, uint32 *d2, uint32 *d3, FILE *fp, uint32 swap, uint32 *chksum)</div><div class="ttdoc">Read a 3-d array (set of matrices) </div><div class="ttdef"><b>Definition:</b> <a href="bio_8c_source.html#l00430">bio.c:430</a></div></div>
<div class="ttc" id="structfeat__t_html"><div class="ttname"><a href="structfeat__t.html">feat_t</a></div><div class="ttdoc">Structure for describing a speech feature type Structure for describing a speech feature type (no...</div></div>
<div class="ttc" id="ckd__alloc_8h_html_a31c6b405558620ac37599737b5722fbf"><div class="ttname"><a href="ckd__alloc_8h.html#a31c6b405558620ac37599737b5722fbf">ckd_free</a></div><div class="ttdeci">SPHINXBASE_EXPORT void ckd_free(void *ptr)</div><div class="ttdoc">Test and free a 1-D array. </div><div class="ttdef"><b>Definition:</b> <a href="ckd__alloc_8c_source.html#l00241">ckd_alloc.c:241</a></div></div>
<div class="ttc" id="err_8h_html_a6a794bec721b555ac1f2167f9e12f662"><div class="ttname"><a href="err_8h.html#a6a794bec721b555ac1f2167f9e12f662">E_WARN</a></div><div class="ttdeci">#define E_WARN</div><div class="ttdoc">Print warning information to standard error stream. </div><div class="ttdef"><b>Definition:</b> <a href="err_8h_source.html#l00164">err.h:164</a></div></div>
<div class="ttc" id="feat_8h_html"><div class="ttname"><a href="feat_8h.html">feat.h</a></div><div class="ttdoc">compute the dynamic coefficients from the cepstral vector. </div></div>
<div class="ttc" id="err_8h_html"><div class="ttname"><a href="err_8h.html">err.h</a></div><div class="ttdoc">Implementation of logging routines. </div></div>
<div class="ttc" id="ckd__alloc_8h_html_a443accde4663a4a07c6d147088614d57"><div class="ttname"><a href="ckd__alloc_8h.html#a443accde4663a4a07c6d147088614d57">ckd_free_3d</a></div><div class="ttdeci">SPHINXBASE_EXPORT void ckd_free_3d(void *ptr)</div><div class="ttdoc">Free a 3-D array (ptr) previously allocated by ckd_calloc_3d. </div><div class="ttdef"><b>Definition:</b> <a href="ckd__alloc_8c_source.html#l00294">ckd_alloc.c:294</a></div></div>
<div class="ttc" id="err_8h_html_a1a4495946ab2449d61108fe829a94613"><div class="ttname"><a href="err_8h.html#a1a4495946ab2449d61108fe829a94613">E_FATAL</a></div><div class="ttdeci">#define E_FATAL</div><div class="ttdoc">Exit with non-zero status after error message. </div><div class="ttdef"><b>Definition:</b> <a href="err_8h_source.html#l00127">err.h:127</a></div></div>
<div class="ttc" id="feat_8h_html_a7452702cb435ea9244c7e1f61bbf5da1"><div class="ttname"><a href="feat_8h.html#a7452702cb435ea9244c7e1f61bbf5da1">feat_lda_transform</a></div><div class="ttdeci">SPHINXBASE_EXPORT void feat_lda_transform(feat_t *fcb, mfcc_t ***inout_feat, uint32 nfr)</div><div class="ttdoc">Transform a block of features using the feature module's LDA transform. </div><div class="ttdef"><b>Definition:</b> <a href="lda_8c_source.html#l00139">lda.c:139</a></div></div>
<div class="ttc" id="err_8h_html_a5f7b2f58f5a663a6bdd51f197ae21993"><div class="ttname"><a href="err_8h.html#a5f7b2f58f5a663a6bdd51f197ae21993">E_ERROR</a></div><div class="ttdeci">#define E_ERROR</div><div class="ttdoc">Print error message to standard error stream. </div><div class="ttdef"><b>Definition:</b> <a href="err_8h_source.html#l00169">err.h:169</a></div></div>
<div class="ttc" id="err_8h_html_a54ffbfe898d74595c586a1f48f32ef03"><div class="ttname"><a href="err_8h.html#a54ffbfe898d74595c586a1f48f32ef03">E_ERROR_SYSTEM</a></div><div class="ttdeci">#define E_ERROR_SYSTEM</div><div class="ttdoc">Print error text; Call perror("");. </div><div class="ttdef"><b>Definition:</b> <a href="err_8h_source.html#l00142">err.h:142</a></div></div>
<div class="ttc" id="bio_8h_html_adb3c858a50381aeff022ccd4271af2d9"><div class="ttname"><a href="bio_8h.html#adb3c858a50381aeff022ccd4271af2d9">bio_hdrarg_free</a></div><div class="ttdeci">SPHINXBASE_EXPORT void bio_hdrarg_free(char **name, char **val)</div><div class="ttdoc">Free name and value strings previously allocated and returned by bio_readhdr. </div><div class="ttdef"><b>Definition:</b> <a href="bio_8c_source.html#l00121">bio.c:121</a></div></div>
<div class="ttc" id="bio_8h_html_a7a155ff51740f1d9a31f7bb7b6e4bc3a"><div class="ttname"><a href="bio_8h.html#a7a155ff51740f1d9a31f7bb7b6e4bc3a">bio_readhdr</a></div><div class="ttdeci">SPHINXBASE_EXPORT int32 bio_readhdr(FILE *fp, char ***name, char ***val, int32 *swap)</div><div class="ttdoc">Read binary file format header: has the following format. </div><div class="ttdef"><b>Definition:</b> <a href="bio_8c_source.html#l00187">bio.c:187</a></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_68267d1309a1af8e8297ef4c3efbcdba.html">src</a></li><li class="navelem"><a class="el" href="dir_3b5e6589d8d7d3e6424fe0fcd4bf9c2a.html">libsphinxbase</a></li><li class="navelem"><a class="el" href="dir_4d78eab29b32367e0c8277af3fb65de0.html">feat</a></li><li class="navelem"><b>lda.c</b></li>
<li class="footer">Generated on Sat Oct 18 2014 15:21:16 for SphinxBase by
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.8 </li>
</ul>
</div>
</body>
</html>
distrib
>
Mageia
>
5
>
x86_64
>
media
>
core-release
>
by-pkgid
>
b56c9eed3fa0f00740b6f09a2460a97f
>
files
>
323
