Marsyas
0.6.0-alpha
|
00001 /****************************************************** 00002 * LU decomposition routines used by test_lu.cpp * 00003 * with dynamic allocations * 00004 * * 00005 * C++ version by J-P Moreau, Paris * 00006 * --------------------------------------------------- * 00007 * Reference: * 00008 * * 00009 * "Numerical Recipes by W.H. Press, B. P. Flannery, * 00010 * S.A. Teukolsky and W.T. Vetterling, Cambridge * 00011 * University Press, 1986". * 00012 * --------------------------------------------------- * 00013 * Uses: basis_r.cpp and vmblock.cpp * 00014 ******************************************************/ 00015 00016 #if !defined(_LU_CPP_) 00017 #define _LU_CPP_ 00018 00019 #include <marsyas/basis.h> 00020 #include "vmblock.h" 00021 00022 #define NMAX 100 00023 #define TINY 1.5e-16 00024 00025 int LUDCMP(REAL **A, int n, int *INDX, int *d); 00026 void LUBKSB(REAL **A, int n, int *INDX, REAL *B); 00027 00028 00029 #endif //_LU_CPP_ 00030 00031