txlyre
/
avram


			
							123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128
							<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html401/loose.dtd">
<html>
<!-- Created on November 8, 2012 by texi2html 1.82
texi2html was written by: 
            Lionel Cons <[email protected]> (original author)
            Karl Berry  <[email protected]>
            Olaf Bachmann <[email protected]>
            and many others.
Maintained by: Many creative people.
Send bugs and suggestions to <[email protected]>
-->
<head>
<title>avram - a virtual machine code interpreter: D.10.3 Additional lapack notes</title>

<meta name="description" content="avram - a virtual machine code interpreter: D.10.3 Additional lapack notes">
<meta name="keywords" content="avram - a virtual machine code interpreter: D.10.3 Additional lapack notes">
<meta name="resource-type" content="document">
<meta name="distribution" content="global">
<meta name="Generator" content="texi2html 1.82">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<style type="text/css">
<!--
a.summary-letter {text-decoration: none}
blockquote.smallquotation {font-size: smaller}
pre.display {font-family: serif}
pre.format {font-family: serif}
pre.menu-comment {font-family: serif}
pre.menu-preformatted {font-family: serif}
pre.smalldisplay {font-family: serif; font-size: smaller}
pre.smallexample {font-size: smaller}
pre.smallformat {font-family: serif; font-size: smaller}
pre.smalllisp {font-size: smaller}
span.roman {font-family:serif; font-weight:normal;}
span.sansserif {font-family:sans-serif; font-weight:normal;}
ul.toc {list-style: none}
-->
</style>


</head>

<body lang="en" bgcolor="#FFFFFF" text="#000000" link="#0000FF" vlink="#800080" alink="#FF0000">

<a name="Additional-lapack-notes"></a>
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="lapack-exceptions.html#lapack-exceptions" title="Previous section in reading order"> &lt; </a>]</td>
<td valign="middle" align="left">[<a href="math.html#math" title="Next section in reading order"> &gt; </a>]</td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left">[<a href="External-Libraries.html#External-Libraries" title="Beginning of this chapter or previous chapter"> &lt;&lt; </a>]</td>
<td valign="middle" align="left">[<a href="lapack.html#lapack" title="Up section"> Up </a>]</td>
<td valign="middle" align="left">[<a href="Copying.html#Copying" title="Next chapter"> &gt;&gt; </a>]</td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left">[<a href="avram.html#Top" title="Cover (top) of document">Top</a>]</td>
<td valign="middle" align="left">[<a href="avram_toc.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Function-Index.html#Function-Index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="avram_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<hr size="1">
<a name="Additional-lapack-notes-1"></a>
<h3 class="subsection">D.10.3 Additional <code>lapack</code> notes</h3>

<p>The functions <code>dgesdd</code> and <code>zgesdd</code> are an effective
dimensionality reduction technique for a large database of time
<a name="index-dimensionality-reduction"></a>
series.  A set of basis vectors can be computed once for the database,
and then any time series in the database can be expressed as a linear
combination thereof. To the extent that the data embody any redundant
information, an approximate reconstruction of an individual series
from the database will require fewer coefficients (maybe far fewer) in
terms of the basis than original length of the series.
</p>
<p>The library functions <code>dgelsd</code> and <code>zgelsd</code> are good for
<a name="index-least-squares-2"></a>
finding least squares fits to empirical data. If the matrix parameter
<var>a</var> is interpreted as a list of inputs and the vector parameter
<var>b</var> as the list of corresponding output data from some unknown
linear function of <var>n</var> variables <var>f</var>, then <var>x</var> is the list
of coefficients whereby <var>f</var> achieves the optimum fit to the data
in the least squares sense.
</p>
<p>These functions solve a special case of the problem solved by
<a name="index-generalized-least-squares-1"></a>
<a name="index-least-squares-3"></a>
<code>dggglm</code> and <code>zggglm</code> where the parameter <var>B</var> is the
identity matrix. For the latter functions, the output vector <var>y</var>
can be interpreted as a measure of the error, and <var>B</var> can be
chosen to express unequal costs for errors at different points in
the fitted function.
</p>
<p>Cholesky decompositions obtained by <code>dpptrf</code> and <code>zpptrf</code>
<a name="index-Cholesky-decomposition-1"></a>
are useful for generating correlated random numbers. A population of
vectors of uncorrelated standard normally distributed random numbers
can be made to exhibit any correlations to order by multiplying all of
<a name="index-correlation"></a>
the vectors by the lower Cholesky factor of the desired covariance
<a name="index-covariance-matrix-1"></a>
matrix.
</p>
<hr size="1">
<table cellpadding="1" cellspacing="1" border="0">
<tr><td valign="middle" align="left">[<a href="lapack-exceptions.html#lapack-exceptions" title="Previous section in reading order"> &lt; </a>]</td>
<td valign="middle" align="left">[<a href="math.html#math" title="Next section in reading order"> &gt; </a>]</td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left">[<a href="External-Libraries.html#External-Libraries" title="Beginning of this chapter or previous chapter"> &lt;&lt; </a>]</td>
<td valign="middle" align="left">[<a href="lapack.html#lapack" title="Up section"> Up </a>]</td>
<td valign="middle" align="left">[<a href="Copying.html#Copying" title="Next chapter"> &gt;&gt; </a>]</td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left"> &nbsp; </td>
<td valign="middle" align="left">[<a href="avram.html#Top" title="Cover (top) of document">Top</a>]</td>
<td valign="middle" align="left">[<a href="avram_toc.html#SEC_Contents" title="Table of contents">Contents</a>]</td>
<td valign="middle" align="left">[<a href="Function-Index.html#Function-Index" title="Index">Index</a>]</td>
<td valign="middle" align="left">[<a href="avram_abt.html#SEC_About" title="About (help)"> ? </a>]</td>
</tr></table>
<p>
 <font size="-1">
  This document was generated on <i>November 8, 2012</i> using <a href="http://www.nongnu.org/texi2html/"><i>texi2html 1.82</i></a>.
 </font>
 <br>

</p>
</body>
</html>