00001 /* $Id$ */ 00002 00003 /* GenLib2 - General Library 2 00004 * Copyright (C) 2005 Dirk Stallmann 00005 * 00006 * This program is free software; you can redistribute it and/or 00007 * modify it under the terms of the GNU General Public License 00008 * as published by the Free Software Foundation; either version 2 00009 * of the License, or (at your option) any later version. 00010 * 00011 * This program is distributed in the hope that it will be useful, 00012 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00013 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00014 * GNU General Public License for more details. 00015 * 00016 * You should have received a copy of the GNU General Public License 00017 * along with this program; if not, write to the Free Software 00018 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 00019 */ 00020 00021 #ifndef GENLIB2_BLAS_INCLUDED 00022 #define GENLIB2_BLAS_INCLUDED 00023 00024 #include "vector.h" 00025 #include "matrix.h" 00026 00027 // external C functions 00028 extern "C" 00029 { 00030 // copy vector 00031 void dcopy(int&, double&, int&, double&, int&); 00032 00033 // dot product 00034 double ddot(int&, double&, int&, double&, int&); 00035 } 00036 00037 namespace GenLib2 00038 { 00039 00043 bool lapack_dpotrf(Matrix<double>& A) 00044 { 00045 bool rtn = true; 00046 char uplo = 'u'; 00047 int n = A.rows(); 00048 int lda = A.rows(); 00049 int info = 0; 00050 00051 F77_DPOTRF(uplo, n, A(0,0), lda, info); 00052 if (info != 0) rtn = false; 00053 00054 return rtn; 00055 } 00056 00057 } // namespace GenLib2 00058 00059 #endif /* GENLIB2_BLAS_INCLUDED */