Index: doc/html/classldmat.html
===================================================================
--- doc/html/classldmat.html (revision 37)
+++ doc/html/classldmat.html (revision 79)
@@ -77,8 +77,8 @@
 <tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#e967b9425007f0cb6cd59b845f9756d8">mult_sym</a> (const mat &amp;C)</td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Inplace symmetric multiplication by a SQUARE matrix $C$, i.e. $V = C*V*C'$.  <a href="#e967b9425007f0cb6cd59b845f9756d8"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Inplace symmetric multiplication by a SQUARE matrix <img class="formulaInl" alt="$C$" src="form_24.png">, i.e. <img class="formulaInl" alt="$V = C*V*C'$" src="form_25.png">.  <a href="#e967b9425007f0cb6cd59b845f9756d8"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#4fd155f38eb6dd5af4bdf9c98a7999a9">mult_sym_t</a> (const mat &amp;C)</td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Inplace symmetric multiplication by a SQUARE transpose of matrix $C$, i.e. $V = C'*V*C$.  <a href="#4fd155f38eb6dd5af4bdf9c98a7999a9"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Inplace symmetric multiplication by a SQUARE transpose of matrix <img class="formulaInl" alt="$C$" src="form_24.png">, i.e. <img class="formulaInl" alt="$V = C'*V*C$" src="form_26.png">.  <a href="#4fd155f38eb6dd5af4bdf9c98a7999a9"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="a60f2c7e4f3c6a7738eaaaab81ffad20"></a><!-- doxytag: member="ldmat::add" ref="a60f2c7e4f3c6a7738eaaaab81ffad20" args="(const ldmat &amp;ld2, double w=1.0)" -->
 void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#a60f2c7e4f3c6a7738eaaaab81ffad20">add</a> (const <a class="el" href="classldmat.html">ldmat</a> &amp;ld2, double w=1.0)</td></tr>
@@ -92,5 +92,9 @@
 double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#d64f331b781903e913cb2ee836886f3f">qform</a> (const vec &amp;v) const </td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Evaluates quadratic form $x= v'*V*v$;. <br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Evaluates quadratic form <img class="formulaInl" alt="$x= v'*V*v$" src="form_29.png">;. <br></td></tr>
+<tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="d876c5f83e02b3e809b35c9de5068f14"></a><!-- doxytag: member="ldmat::invqform" ref="d876c5f83e02b3e809b35c9de5068f14" args="(const vec &amp;v) const " -->
+double&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#d876c5f83e02b3e809b35c9de5068f14">invqform</a> (const vec &amp;v) const </td></tr>
+
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Evaluates quadratic form <img class="formulaInl" alt="$x= v'*inv(V)*v$" src="form_30.png">;. <br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="4d6e401de9607332305c27e67972a07a"></a><!-- doxytag: member="ldmat::clear" ref="4d6e401de9607332305c27e67972a07a" args="()" -->
 void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#4d6e401de9607332305c27e67972a07a">clear</a> ()</td></tr>
@@ -107,5 +111,5 @@
 <tr><td class="memItemLeft" nowrap align="right" valign="top">vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#fc380626ced6f9244fb58c5f0231174d">sqrt_mult</a> (const vec &amp;v) const </td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Multiplies square root of $V$ by vector $x$.  <a href="#fc380626ced6f9244fb58c5f0231174d"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Multiplies square root of <img class="formulaInl" alt="$V$" src="form_27.png"> by vector <img class="formulaInl" alt="$x$" src="form_28.png">.  <a href="#fc380626ced6f9244fb58c5f0231174d"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top">virtual void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#2c160cb123c1102face7a50ec566a031">inv</a> (<a class="el" href="classldmat.html">ldmat</a> &amp;Inv) const </td></tr>
 
@@ -113,11 +117,11 @@
 <tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#e7207748909325bb0f99b43f090a2b7e">mult_sym</a> (const mat &amp;C, <a class="el" href="classldmat.html">ldmat</a> &amp;U) const </td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Symmetric multiplication of $U$ by a general matrix $C$, result of which is stored in the current class.  <a href="#e7207748909325bb0f99b43f090a2b7e"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Symmetric multiplication of <img class="formulaInl" alt="$U$" src="form_31.png"> by a general matrix <img class="formulaInl" alt="$C$" src="form_24.png">, result of which is stored in the current class.  <a href="#e7207748909325bb0f99b43f090a2b7e"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#f94dc3a233f3d40fc853d8d4ac3b8eab">mult_sym_t</a> (const mat &amp;C, <a class="el" href="classldmat.html">ldmat</a> &amp;U) const </td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Symmetric multiplication of $U$ by a transpose of a general matrix $C$, result of which is stored in the current class.  <a href="#f94dc3a233f3d40fc853d8d4ac3b8eab"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Symmetric multiplication of <img class="formulaInl" alt="$U$" src="form_31.png"> by a transpose of a general matrix <img class="formulaInl" alt="$C$" src="form_24.png">, result of which is stored in the current class.  <a href="#f94dc3a233f3d40fc853d8d4ac3b8eab"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top">void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#f291faa073e7bc8dfafc7ae93daa2506">ldform</a> (const mat &amp;A, const vec &amp;D0)</td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Transforms general $A'D0 A$ into pure $L'DL$.  <a href="#f291faa073e7bc8dfafc7ae93daa2506"></a><br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Transforms general <img class="formulaInl" alt="$A'D0 A$" src="form_32.png"> into pure <img class="formulaInl" alt="$L'DL$" src="form_33.png">.  <a href="#f291faa073e7bc8dfafc7ae93daa2506"></a><br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="0884a613b94fde61bfc84288e73ce57f"></a><!-- doxytag: member="ldmat::setD" ref="0884a613b94fde61bfc84288e73ce57f" args="(const vec &amp;nD)" -->
 void&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#0884a613b94fde61bfc84288e73ce57f">setD</a> (const vec &amp;nD)</td></tr>
@@ -146,9 +150,9 @@
 vec&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#4cce04824539c4a8d062d9a36d6e014e">D</a></td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Positive vector $D$. <br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Positive vector <img class="formulaInl" alt="$D$" src="form_22.png">. <br></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="f74a64b99fe58a75ebd37bb679e121ea"></a><!-- doxytag: member="ldmat::L" ref="f74a64b99fe58a75ebd37bb679e121ea" args="" -->
 mat&nbsp;</td><td class="memItemRight" valign="bottom"><a class="el" href="classldmat.html#f74a64b99fe58a75ebd37bb679e121ea">L</a></td></tr>
 
-<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Lower-triangular matrix $L$. <br></td></tr>
+<tr><td class="mdescLeft">&nbsp;</td><td class="mdescRight">Lower-triangular matrix <img class="formulaInl" alt="$L$" src="form_21.png">. <br></td></tr>
 <tr><td colspan="2"><br><h2>Friends</h2></td></tr>
 <tr><td class="memItemLeft" nowrap align="right" valign="top"><a class="anchor" name="eaaa0baa6026b84cfcbced41c84599d1"></a><!-- doxytag: member="ldmat::operator&lt;&lt;" ref="eaaa0baa6026b84cfcbced41c84599d1" args="(std::ostream &amp;os, const ldmat &amp;sq)" -->
@@ -163,5 +167,5 @@
 <img class="formulaDsp" alt="\[M = L'DL\]" src="form_5.png">
 <p>
- where only $L$ and $D$ matrices are stored. All inplace operations modifies only these and the need to compose and decompose the matrix is avoided. <hr><h2>Member Function Documentation</h2>
+ where only <img class="formulaInl" alt="$L$" src="form_21.png"> and <img class="formulaInl" alt="$D$" src="form_22.png"> matrices are stored. All inplace operations modifies only these and the need to compose and decompose the matrix is avoided. <hr><h2>Member Function Documentation</h2>
 <a class="anchor" name="0f0f6e083e6d947cf58097ffce3ccd1a"></a><!-- doxytag: member="ldmat::opupdt" ref="0f0f6e083e6d947cf58097ffce3ccd1a" args="(const vec &amp;v, double w)" -->
 <div class="memitem">
@@ -190,5 +194,5 @@
 
 <p>
-Perfroms a rank-1 update by outer product of vectors: $V = V + w v v'$. <dl compact><dt><b>Parameters:</b></dt><dd>
+Perfroms a rank-1 update by outer product of vectors: <img class="formulaInl" alt="$V = V + w v v'$" src="form_23.png">. <dl compact><dt><b>Parameters:</b></dt><dd>
   <table border="0" cellspacing="2" cellpadding="0">
     <tr><td valign="top"></td><td valign="top"><em>v</em>&nbsp;</td><td>Vector forming the outer product to be added </td></tr>
@@ -218,5 +222,5 @@
 
 <p>
-Inplace symmetric multiplication by a SQUARE matrix $C$, i.e. $V = C*V*C'$. 
+Inplace symmetric multiplication by a SQUARE matrix <img class="formulaInl" alt="$C$" src="form_24.png">, i.e. <img class="formulaInl" alt="$V = C*V*C'$" src="form_25.png">. 
 <p>
 <dl compact><dt><b>Parameters:</b></dt><dd>
@@ -247,5 +251,5 @@
 
 <p>
-Inplace symmetric multiplication by a SQUARE transpose of matrix $C$, i.e. $V = C'*V*C$. 
+Inplace symmetric multiplication by a SQUARE transpose of matrix <img class="formulaInl" alt="$C$" src="form_24.png">, i.e. <img class="formulaInl" alt="$V = C'*V*C$" src="form_26.png">. 
 <p>
 <dl compact><dt><b>Parameters:</b></dt><dd>
@@ -276,5 +280,5 @@
 
 <p>
-Multiplies square root of $V$ by vector $x$. 
+Multiplies square root of <img class="formulaInl" alt="$V$" src="form_27.png"> by vector <img class="formulaInl" alt="$x$" src="form_28.png">. 
 <p>
 Used e.g. in generating normal samples. 
@@ -336,5 +340,5 @@
 
 <p>
-Symmetric multiplication of $U$ by a general matrix $C$, result of which is stored in the current class. 
+Symmetric multiplication of <img class="formulaInl" alt="$U$" src="form_31.png"> by a general matrix <img class="formulaInl" alt="$C$" src="form_24.png">, result of which is stored in the current class. 
 <p>
 <dl compact><dt><b>Parameters:</b></dt><dd>
@@ -373,5 +377,5 @@
 
 <p>
-Symmetric multiplication of $U$ by a transpose of a general matrix $C$, result of which is stored in the current class. 
+Symmetric multiplication of <img class="formulaInl" alt="$U$" src="form_31.png"> by a transpose of a general matrix <img class="formulaInl" alt="$C$" src="form_24.png">, result of which is stored in the current class. 
 <p>
 <dl compact><dt><b>Parameters:</b></dt><dd>
@@ -410,7 +414,7 @@
 
 <p>
-Transforms general $A'D0 A$ into pure $L'DL$. 
-<p>
-The new decomposition fullfills: $A'*diag(D)*A = self.L'*diag(self.D)*self.L$ <dl compact><dt><b>Parameters:</b></dt><dd>
+Transforms general <img class="formulaInl" alt="$A'D0 A$" src="form_32.png"> into pure <img class="formulaInl" alt="$L'DL$" src="form_33.png">. 
+<p>
+The new decomposition fullfills: <img class="formulaInl" alt="$A'*diag(D)*A = self.L'*diag(self.D)*self.L$" src="form_34.png"> <dl compact><dt><b>Parameters:</b></dt><dd>
   <table border="0" cellspacing="2" cellpadding="0">
     <tr><td valign="top"></td><td valign="top"><em>A</em>&nbsp;</td><td>general matrix </td></tr>
@@ -467,5 +471,5 @@
 <hr>The documentation for this class was generated from the following files:<ul>
 <li>work/mixpp/bdm/math/<a class="el" href="libDC_8h-source.html">libDC.h</a><li>work/mixpp/bdm/math/libDC.cpp</ul>
-<hr size="1"><address style="text-align: right;"><small>Generated on Wed Mar 12 16:15:52 2008 for mixpp by&nbsp;
+<hr size="1"><address style="text-align: right;"><small>Generated on Fri Apr 18 11:15:21 2008 for mixpp by&nbsp;
 <a href="http://www.doxygen.org/index.html">
 <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.3 </small></address>
