blas.h Source File

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 */