--- trunk/src/math/LU.hpp	2010/05/10 17:28:26	1442
+++ trunk/src/math/LU.hpp	2013/06/16 15:15:42	1879
@@ -35,8 +35,9 @@
  *                                                                      
  * [1]  Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).             
  * [2]  Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).          
- * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008).          
- * [4]  Vardeman & Gezelter, in progress (2009).                        
+ * [3]  Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).          
+ * [4]  Kuang & Gezelter,  J. Chem. Phys. 133, 164101 (2010).
+ * [5]  Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
  */
  
 /*=========================================================================
@@ -266,8 +267,8 @@ int LUFactorLinearSystem(MatrixType& A, int *index, in
  * lower triangular matrix and U is upper triangular matrix. 
  * @param A input square matrix
  * @param index pivot indices
+ * @param x vector
  * @param size size of the matrix and temporary arrays
- * @param tmpSize temporary array
  * @return true if inverse is computed, otherwise return false
  * @note A=LU and index[] are generated from method LUFactorLinearSystem). 
  * Also, solution vector is written directly over input load vector.