Show
Ignore:
Timestamp:
06/02/09 10:24:26 (16 years ago)
Author:
smidl
Message:

doc - oprava

Files:
1 modified

Legend:

Unmodified
Added
Removed

  • doc/html/libDC_8h-source.html

    r353 r354  
    7575<a name="l00028"></a>00028  
    7676<a name="l00030"></a>00030 <span class="comment">//TODO can be done via: dtrtri.f from lapack</span> 
    77 <a name="l00031"></a>00031 <a class="code" href="classitpp_1_1Mat.html#6bba394f181c76fda12759568986c613" title="Default Matrix Type.">mat</a> <a class="code" href="group__math.html#g6715d039e6d5d97005cf9e2522dfa474" title="Auxiliary function ltuinv; inversion of a triangular matrix;.">ltuinv</a>( <span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html#6bba394f181c76fda12759568986c613" title="Default Matrix Type.">mat</a> &amp;L ); 
     77<a name="l00031"></a>00031 mat <a class="code" href="group__math.html#g6715d039e6d5d97005cf9e2522dfa474" title="Auxiliary function ltuinv; inversion of a triangular matrix;.">ltuinv</a>( <span class="keyword">const</span> mat &amp;L ); 
    7878<a name="l00032"></a>00032  
    7979<a name="l00037"></a><a class="code" href="classsqmat.html">00037</a> <span class="keyword">class </span><a class="code" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form...">sqmat</a> 
    … …  
    130130<a name="l00137"></a>00137 <span class="comment"></span>                <a class="code" href="classfsqmat.html#79e3f73e0ccd663c7f7e08083d272940" title="Default initialization.">fsqmat</a>(<span class="keyword">const</span> <span class="keywordtype">int</span> dim0); <span class="comment">// mat will be initialized OK</span> 
    131131<a name="l00139"></a>00139 <span class="comment"></span>                <a class="code" href="classfsqmat.html#79e3f73e0ccd663c7f7e08083d272940" title="Default initialization.">fsqmat</a> ( <span class="keyword">const</span> mat &amp;<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a> ); 
    132 <a name="l00141"></a><a class="code" href="classfsqmat.html#9e19a5e26083f4a493b6fad38cc40fab">00141</a>                 <a class="code" href="classfsqmat.html#9e19a5e26083f4a493b6fad38cc40fab" title="Constructor.">fsqmat</a> ( <span class="keyword">const</span> <a class="code" href="classfsqmat.html" title="Fake sqmat. This class maps sqmat operations to operations on full matrix.">fsqmat</a> &amp;<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>, <span class="keyword">const</span> ivec &amp;perm ):<a class="code" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form...">sqmat</a>(M.<a class="code" href="classsqmat.html#071e80ced9cc3b8cbb360fa7462eb646" title="Reimplementing common functions of mat: cols().">rows</a>()){<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a>(<span class="stringliteral">"not implemneted"</span>);}; 
    133 <a name="l00143"></a><a class="code" href="classfsqmat.html#c01f3e9bb590f2a2921369d672f3ce1e">00143</a>                 <a class="code" href="classfsqmat.html#c01f3e9bb590f2a2921369d672f3ce1e" title="Constructor.">fsqmat</a> ( <span class="keyword">const</span> vec &amp;d ):<a class="code" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form...">sqmat</a>(d.<a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>()){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>=<a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(d);}; 
     132<a name="l00141"></a><a class="code" href="classfsqmat.html#9e19a5e26083f4a493b6fad38cc40fab">00141</a>                 <a class="code" href="classfsqmat.html#9e19a5e26083f4a493b6fad38cc40fab" title="Constructor.">fsqmat</a> ( <span class="keyword">const</span> <a class="code" href="classfsqmat.html" title="Fake sqmat. This class maps sqmat operations to operations on full matrix.">fsqmat</a> &amp;<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>, <span class="keyword">const</span> ivec &amp;perm ):<a class="code" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form...">sqmat</a>(M.<a class="code" href="classsqmat.html#071e80ced9cc3b8cbb360fa7462eb646" title="Reimplementing common functions of mat: cols().">rows</a>()){it_error(<span class="stringliteral">"not implemneted"</span>);}; 
     133<a name="l00143"></a><a class="code" href="classfsqmat.html#c01f3e9bb590f2a2921369d672f3ce1e">00143</a>                 <a class="code" href="classfsqmat.html#c01f3e9bb590f2a2921369d672f3ce1e" title="Constructor.">fsqmat</a> ( <span class="keyword">const</span> vec &amp;d ):<a class="code" href="classsqmat.html" title="Virtual class for representation of double symmetric matrices in square-root form...">sqmat</a>(d.length()){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>=diag(d);}; 
    134134<a name="l00144"></a>00144  
    135135<a name="l00146"></a><a class="code" href="classfsqmat.html#2a8f104e4befbc2aa90d8b11edfedb2e">00146</a>                 <span class="keyword">virtual</span> <a class="code" href="classfsqmat.html#2a8f104e4befbc2aa90d8b11edfedb2e" title="Destructor for future use;.">~fsqmat</a>(){}; 
    … …  
    138138<a name="l00154"></a>00154                 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classfsqmat.html#9fa853e1ca28f2a1a1c43377e798ecb1" title="Matrix inversion preserving the chosen form.">inv</a> ( <a class="code" href="classfsqmat.html" title="Fake sqmat. This class maps sqmat operations to operations on full matrix.">fsqmat</a> &amp;Inv ); 
    139139<a name="l00155"></a>00155  
    140 <a name="l00156"></a><a class="code" href="classfsqmat.html#eb0d1358f536e4453b5f99d0418ca1e5">00156</a>                 <span class="keywordtype">double</span> <a class="code" href="classfsqmat.html#eb0d1358f536e4453b5f99d0418ca1e5" title="Logarithm of a determinant.">logdet</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="group__logexpfunc.html#g320dc9f6d27a6cfe4e149633fd1ea566" title="The natural logarithm of the elements.">log</a> ( <a class="code" href="group__determinant.html#gd9ded34e98aba6036329f53e8abf51a8" title="Determinant of complex square matrix.Calculate determinant of the complex matrix...">det</a> ( <a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a> ) );}; 
     140<a name="l00156"></a><a class="code" href="classfsqmat.html#eb0d1358f536e4453b5f99d0418ca1e5">00156</a>                 <span class="keywordtype">double</span> <a class="code" href="classfsqmat.html#eb0d1358f536e4453b5f99d0418ca1e5" title="Logarithm of a determinant.">logdet</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> log ( det ( <a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a> ) );}; 
    141141<a name="l00157"></a><a class="code" href="classfsqmat.html#a6c91b0389e73404324b2314b08d6e87">00157</a>                 <span class="keywordtype">double</span> <a class="code" href="classfsqmat.html#a6c91b0389e73404324b2314b08d6e87" title="Evaluates quadratic form ;.">qform</a> (<span class="keyword">const</span>  vec &amp;v )<span class="keyword"> const </span>{<span class="keywordflow">return</span> ( v* ( <a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>*v ) );}; 
    142 <a name="l00158"></a><a class="code" href="classfsqmat.html#58075da64ddadd4df40654c35b928c6f">00158</a>                 <span class="keywordtype">double</span> <a class="code" href="classfsqmat.html#58075da64ddadd4df40654c35b928c6f" title="Evaluates quadratic form ;.">invqform</a> (<span class="keyword">const</span>  vec &amp;v )<span class="keyword"> const </span>{<span class="keywordflow">return</span> ( v* ( <a class="code" href="group__inverse.html#g125b083397cc9450b8a67a485b545bc8" title="Inverse of real square matrix.Calculate the inverse of the real matrix .">itpp::inv</a>(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>)*v ) );}; 
    143 <a name="l00159"></a><a class="code" href="classfsqmat.html#842a774077ee34ac3c36d180ab33e103">00159</a>                 vec <a class="code" href="classfsqmat.html#842a774077ee34ac3c36d180ab33e103" title="Multiplies square root of  by vector .">sqrt_mult</a> (<span class="keyword">const</span> vec &amp;v )<span class="keyword"> const </span>{mat Ch=<a class="code" href="group__matrixdecomp.html#g8aa42136e884519422d4f2dd95e51cb5" title="Cholesky factorisation of real symmetric and positive definite matrix.">chol</a>(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>); <span class="keywordflow">return</span> Ch*v;}; 
     142<a name="l00158"></a><a class="code" href="classfsqmat.html#58075da64ddadd4df40654c35b928c6f">00158</a>                 <span class="keywordtype">double</span> <a class="code" href="classfsqmat.html#58075da64ddadd4df40654c35b928c6f" title="Evaluates quadratic form ;.">invqform</a> (<span class="keyword">const</span>  vec &amp;v )<span class="keyword"> const </span>{<span class="keywordflow">return</span> ( v* ( itpp::inv(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>)*v ) );}; 
     143<a name="l00159"></a><a class="code" href="classfsqmat.html#842a774077ee34ac3c36d180ab33e103">00159</a>                 vec <a class="code" href="classfsqmat.html#842a774077ee34ac3c36d180ab33e103" title="Multiplies square root of  by vector .">sqrt_mult</a> (<span class="keyword">const</span> vec &amp;v )<span class="keyword"> const </span>{mat Ch=chol(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>); <span class="keywordflow">return</span> Ch*v;}; 
    144144<a name="l00160"></a>00160  
    145145<a name="l00162"></a><a class="code" href="classfsqmat.html#a2e0bf7dbbbbe1d3358064c4ad455f1f">00162</a>                 <span class="keywordtype">void</span> <span class="keyword">add</span> ( <span class="keyword">const</span> <a class="code" href="classfsqmat.html" title="Fake sqmat. This class maps sqmat operations to operations on full matrix.">fsqmat</a> &amp;fsq2, <span class="keywordtype">double</span> w=1.0 ){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>+=fsq2.M;}; 
    146146<a name="l00163"></a>00163  
    147 <a name="l00165"></a><a class="code" href="classfsqmat.html#922f8190c13987cbcdb33ec2bf5cf105">00165</a>                 <span class="keywordtype">void</span> <a class="code" href="classfsqmat.html#922f8190c13987cbcdb33ec2bf5cf105" title="Access functions.">setD</a> (<span class="keyword">const</span> vec &amp;nD){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>=<a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(nD);} 
    148 <a name="l00167"></a><a class="code" href="classfsqmat.html#bcf837b2956745e8986044f5600dbd6e">00167</a>                 vec <a class="code" href="classfsqmat.html#bcf837b2956745e8986044f5600dbd6e" title="Access functions.">getD</a> (){<span class="keywordflow">return</span> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>);} 
     147<a name="l00165"></a><a class="code" href="classfsqmat.html#922f8190c13987cbcdb33ec2bf5cf105">00165</a>                 <span class="keywordtype">void</span> <a class="code" href="classfsqmat.html#922f8190c13987cbcdb33ec2bf5cf105" title="Access functions.">setD</a> (<span class="keyword">const</span> vec &amp;nD){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>=diag(nD);} 
     148<a name="l00167"></a><a class="code" href="classfsqmat.html#bcf837b2956745e8986044f5600dbd6e">00167</a>                 vec <a class="code" href="classfsqmat.html#bcf837b2956745e8986044f5600dbd6e" title="Access functions.">getD</a> (){<span class="keywordflow">return</span> diag(<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>);} 
    149149<a name="l00169"></a><a class="code" href="classfsqmat.html#03a8f49eb4d38a054ecc522be59cd2ad">00169</a>                 <span class="keywordtype">void</span> <a class="code" href="classfsqmat.html#03a8f49eb4d38a054ecc522be59cd2ad" title="Access functions.">setD</a> (<span class="keyword">const</span> vec &amp;nD, <span class="keywordtype">int</span> i){<span class="keywordflow">for</span>(<span class="keywordtype">int</span> j=i;j&lt;nD.length();j++){<a class="code" href="classfsqmat.html#a7a1fcb9aae19d1e4daddfc9c22ce453" title="Full matrix on which the operations are performed.">M</a>(j,j)=nD(j-i);}} <span class="comment">//Fixme can be more general</span> 
    150150<a name="l00170"></a>00170  
    … …  
    221221<a name="l00293"></a>00293 <span class="preprocessor">#endif // DC_H</span> 
    222222</pre></div></div> 
    223 <hr size="1"><address style="text-align: right;"><small>Generated on Tue Jun 2 10:02:11 2009 for mixpp by&nbsp; 
     223<hr size="1"><address style="text-align: right;"><small>Generated on Tue Jun 2 10:11:00 2009 for mixpp by&nbsp; 
    224224<a href="http://www.doxygen.org/index.html"> 
    225225<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.8 </small></address>