Revision 795 branches/replacement/src/linad99/ludcmp_inv.cpp
ludcmp_inv.cpp (revision 795)  

3  3 
* 
4  4 
* Copyright (c) 2009 ADMB Foundation 
5  5 
*/ 
6 
/** 

7 
* \file 

8 
* This file contains the routines necessary to compute 

9 
* the inverse of a matrix 

10 
*/ 

6  11  
7  12 
#define HOME_VERSION 
8  13 
#include <fvar.hpp> 
...  ...  
185  190 
return vc; 
186  191 
} 
187  192  
188  
189  
190  
191  193 
/** Inverse of a varaiable matrix. 
192  194 
\param aa dvar_matrix conaining matrix to be inverted,\f$A\f$. 
193  195 
\return dvar_matrix containing \f$A^{1}\f$. 
...  ...  
332  334  
333  335 
dvector y(lb, ub); 
334  336 
dvector x(lb, ub); 
335 
//int lb=rowmin; 

336 
//int ub=rowmax; 

337  
337  338 
cltudecomp & b = bb; 
338  339 
ivector indxinv(lb, ub); 
339  340 
for (i = lb; i <= ub; i++) 
Also available in: Unified diff