OurLibLapack.h

Go to the documentation of this file.

/* Copyright (C) INRIA 1999-2005
** 
** This program is free software; you can redistribute it and/or modify it
** under the terms of the GNU General Public License version 2 as published
** by the Free Software Foundation.
** 
** This program is distributed in the hope that it will be useful, but
** WITHOUT ANY WARRANTY; without even the implied warranty of
** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
** Public License for more details.
** 
** You should have received a copy of the GNU General Public License along
** with this program; if not, write to the Free Software Foundation, Inc.,
** 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
**
*/
#ifndef __ourliblapack_h
#define __ourliblapack_h

#ifdef WINDOWS
#define extern __declspec (dllexport)
#endif


// Fonctions BLAS et LAPACK

// Decomposition (P)LU de la matrice A, stockee dans A et ipiv

extern void
dgetrf_(const int *m,
        const int *n,
        double *A,
        const int *lda,
        int *ipiv,
        int *info);


// Inversion de la matrice A apres avoir fait sa decomposition (P)LU,
// stockee dans A et ipiv
extern void
dgetri_(const int *n,
        double *A,
        const int *lda,
        int *ipiv,
        double *work,
        int *lwork,
        int *info);

// Multiplication par inv(A) apres avoir fait sa decomposition (P)LU,
// stockee dans A et ipiv
extern void
dgetrs_(const char *trans,
        const int *n,
        const int *nrhs,
        const double *A,
        const int *lda,
        const int *ipiv,
        double *B,
        const int *ldb,
        int *info);

// B = inv(A)*B, strictement equivalent a dgetrf_ suivi de dgetrs_
extern int
dgesv_(const int *n,
       const int *nrhs,
       double *A,
       const int *lda,
       int *ipiv,
       double *B,
       const int *ldb,
       int *info);

// Multiplication de matrices: C = alpha*A*B+beta*C
extern int
dgemm_(const char *transa,
       const char *transb,
       const int *m,
       const int *n,
       const int *k,
       const double *alpha,
       const double *A,
       const int *lda,
       const double *B,
       const int *ldb,
       const double *beta,
       double *C,
       const int *ldc);

// Multiplication matrice/vecteur: y = alpha*A*x+beta*y
extern int
dgemv_(const char *trans,
       const int *m,
       const int *n,
       const double *alpha,
       const double *A,
       const int *lda,
       const double *x,
       const int *incx,
       const double *beta,
       double *y,
       const int *incy);

// Addition de vecteurs: y = alpha*x+y
extern int
daxpy_(const int *n,
       const double *alpha,
       const double *x,
       const int *incx,
       double *y,
       const int *incy);

// Multiplication scalaire/vecteur: x = alpha*x
extern int
dscal_(const int *n,
       const double *alpha,
       double *x,
       const int *incx);

// Copie d'un vecteur: y = x
extern int
dcopy_(const int *n,
       const double *x,
       const int *incx,
       double *y,
       const int *incy);


// Algo de resolution de QP, qld
extern int 
ql0001_(const int *m,
        const int *me,
        const int *mmax,
        const int *n,
        const int *nmax,
        const int *mnn,
        const double *c,
        const double *d,
        const double *a,
        const double *b,
        const double *xl,
        const double *xu,
        double *x,
        double *u,
        const int *iout,
        int *ifail,
        const int *iprint,
        double *war,
        int *lwar,
        int *iwar,
        int *liwar,
        const double *eps1);



#endif

HuMAnS

humans-users@inrialpes.fr

---