00001 // This is oxl/mvl/PMatrixDecompAa.cxx 00002 #ifdef VCL_NEEDS_PRAGMA_INTERFACE 00003 #pragma implementation 00004 #endif 00005 00006 #include "PMatrixDecompAa.h" 00007 00008 void 00009 PMatrixDecompAa::set(const vnl_matrix<double>& P) 00010 { 00011 for (int r = 0; r < 3; ++r) { 00012 for (int c = 0; c < 3; ++c) 00013 A(r,c) = P(r,c); 00014 a[r] = P(r,3); 00015 } 00016 }