factorizations.h

Go to the documentation of this file.

/* ------------------------------------------------------------
 * This is the file "factorizations.h" of the H2Lib package.
 * All rights reserved, Steffen Boerm 2010
 * ------------------------------------------------------------ */

#ifndef FACTORIZATIONS_H
#define FACTORIZATIONS_H

#include "blas.h"
#include "amatrix.h"
#include "truncation.h"

/* ------------------------------------------------------------
 * Diagonal matrices
 * ------------------------------------------------------------ */

HEADER_PREFIX void
diagsolve_amatrix_avector(bool atrans, pcamatrix a, pavector x);

HEADER_PREFIX void
diagsolve_amatrix(bool atrans, pcamatrix a, bool xtrans, pamatrix x);

HEADER_PREFIX void
diageval_amatrix_avector(bool atrans, pcamatrix a, pavector x);

HEADER_PREFIX void
diageval_amatrix(bool atrans, pcamatrix a, bool xtrans, pamatrix x);

HEADER_PREFIX void
diageval_realavector_amatrix(field alpha, bool atrans, pcrealavector a,
    bool xtrans, pamatrix x);

/* ------------------------------------------------------------
 * Triangular matrices
 * ------------------------------------------------------------ */

HEADER_PREFIX void
triangularsolve_amatrix_avector(bool alower, bool aunit, bool atrans,
    pcamatrix a, pavector x);

HEADER_PREFIX void
triangularsolve_amatrix(bool alower, bool aunit, bool atrans, pcamatrix a,
    bool xtrans, pamatrix x);

HEADER_PREFIX void
triangulareval_amatrix_avector(bool alower, bool aunit, bool atrans,
    pcamatrix a, pavector x);

HEADER_PREFIX void
triangulareval_amatrix(bool alower, bool aunit, bool atrans, pcamatrix a,
    bool xtrans, pamatrix x);

HEADER_PREFIX void
triangularaddmul_amatrix(field alpha, bool alower, bool atrans, pcamatrix a,
    bool blower, bool btrans, pcamatrix b, pamatrix c);

HEADER_PREFIX void
copy_lower_amatrix(pcamatrix a, bool aunit, pamatrix b);

HEADER_PREFIX void
copy_upper_amatrix(pcamatrix a, bool aunit, pamatrix b);

/* ------------------------------------------------------------
 * Triangular decompositions
 * ------------------------------------------------------------ */

HEADER_PREFIX uint
lrdecomp_amatrix(pamatrix a);

HEADER_PREFIX uint
lrdecomp_blocks_amatrix(pamatrix a, uint blocksize);

HEADER_PREFIX void
lrdecomp_tasks_amatrix(pamatrix a, uint blocksize);

HEADER_PREFIX void
lrsolve_n_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
lrsolve_t_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
lrsolve_amatrix_avector(bool atrans, pcamatrix a, pavector x);

HEADER_PREFIX void
lrsolve_amatrix(pcamatrix a, pamatrix x);

HEADER_PREFIX void
lreval_n_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
lreval_t_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
lreval_amatrix_avector(bool atrans, pcamatrix a, pavector x);

HEADER_PREFIX uint
choldecomp_amatrix(pamatrix a);

HEADER_PREFIX uint
choldecomp_blocks_amatrix(pamatrix a, uint blocksize);

HEADER_PREFIX void
choldecomp_tasks_amatrix(pamatrix a, uint blocksize);

HEADER_PREFIX void
cholsolve_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
cholsolve_amatrix(pcamatrix a, pamatrix x);

HEADER_PREFIX void
choleval_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX uint
ldltdecomp_amatrix(pamatrix a);

HEADER_PREFIX void
ldltsolve_amatrix_avector(pcamatrix a, pavector x);

HEADER_PREFIX void
ldltsolve_amatrix(pcamatrix a, pamatrix x);

/* ------------------------------------------------------------
 * Orthogonal decompositions
 * ------------------------------------------------------------ */

HEADER_PREFIX void
qrdecomp_amatrix(pamatrix a, pavector tau);

HEADER_PREFIX uint
qrdecomp_pivot_amatrix(pamatrix a, pavector tau, uint *colpiv);

HEADER_PREFIX uint
qrdecomp_rank_amatrix(pamatrix a, pavector tau, pctruncmode tm, real eps,
    uint *colpiv);

HEADER_PREFIX void
qreval_amatrix_avector(bool qtrans, pcamatrix a, pcavector tau, pavector x);

HEADER_PREFIX void
qreval_amatrix(bool qtrans, pcamatrix a, pcavector tau, pamatrix x);

HEADER_PREFIX void
qrsolve_amatrix_avector(pcamatrix a, pcavector tau, pavector x);

HEADER_PREFIX void
qrinvert_amatrix(pamatrix a);

HEADER_PREFIX void
qrexpand_amatrix(pcamatrix a, pcavector tau, pamatrix q);

#endif