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