h2arith.h

Go to the documentation of this file.

/*@todo extend the function prototypes to general cases*/
#ifndef H2ARITH_H
#define H2ARITH_H

#include "clusteroperator.h"
#include "h2matrix.h"
#include "h2compression.h"
#include "h2update.h"

/* ------------------------------------------------------------
   Auxiliary functions
   ------------------------------------------------------------ */

HEADER_PREFIX ph2matrix
build_from_block_lower_h2matrix(pcblock b, pclusterbasis rb, pclusterbasis cb);

HEADER_PREFIX ph2matrix
build_from_block_upper_h2matrix(pcblock b, pclusterbasis rb, pclusterbasis cb);


HEADER_PREFIX void
tolower_h2matrix(ph2matrix a);

/* ------------------------------------------------------------
   Sum
   ------------------------------------------------------------ */

HEADER_PREFIX void
add_hmatrix_h2matrix(pchmatrix a,
                     ph2matrix B, pclusteroperator rwf, pclusteroperator cwf,
                     ptruncmode tm, real tol);

/* ------------------------------------------------------------
   Convert functions
   ------------------------------------------------------------ */

HEADER_PREFIX pamatrix
convert_h2matrix_amatrix(bool atrans, pch2matrix a);

/* static or in clusterbasis ??? */
HEADER_PREFIX void
fastaddmul_clusterbasis_amatrix(pcclusterbasis cb, pamatrix Yt, pamatrix Y);

HEADER_PREFIX prkmatrix
convert_uniform_rkmatrix(bool atrans, pcuniform a);

HEADER_PREFIX phmatrix
convert_h2matrix_hmatrix(pch2matrix a);

/* ------------------------------------------------------------
   Multiplication
   ------------------------------------------------------------ */

HEADER_PREFIX prkmatrix
mul_h2matrix_rkmatrix(pch2matrix a, bool btrans, pch2matrix B, real tol);

HEADER_PREFIX void
addmul_h2matrix(field alpha, pch2matrix a, bool btrans, pch2matrix B,
                              ph2matrix C, pclusteroperator rwf, pclusteroperator cwf,
                              ptruncmode tm, real tol);

HEADER_PREFIX ph2matrix
mul_h2matrix(field alpha, ph2matrix a, ph2matrix B,
             ph2matrix h2, ptruncmode tm, real tol);

/* ------------------------------------------------------------
   Inversion
   ------------------------------------------------------------ */

HEADER_PREFIX void
invert_h2matrix(ph2matrix a, pclusteroperator arwf, pclusteroperator acwf,
                ph2matrix B, pclusteroperator brwf, pclusteroperator bcwf,
                ptruncmode tm, real tol);

/* ------------------------------------------------------------
   Forward / backward substitution
   ------------------------------------------------------------ */

HEADER_PREFIX void
lowersolve_h2matrix_avector(bool unit, bool atrans, pch2matrix a, pavector x);

HEADER_PREFIX void
uppersolve_h2matrix_avector(bool unit, bool atrans, pch2matrix a, pavector x);



HEADER_PREFIX void
lowersolve_h2matrix_amatrix(bool unit, bool atrans, pch2matrix a,
                            bool xtrans, pamatrix X);


HEADER_PREFIX void
uppersolve_h2matrix_amatrix(bool unit, bool atrans, pch2matrix a,
                            bool xtrans, pamatrix X);

HEADER_PREFIX void
lowersolve_amatrix_h2matrix(bool aunit, bool atrans, pcamatrix a,
                            bool xytrans, pch2matrix X,
                            ph2matrix Y, pclusteroperator rwf, pclusteroperator cwf,
                            ptruncmode tm, real tol);

HEADER_PREFIX void
uppersolve_amatrix_h2matrix(bool aunit, bool atrans, pcamatrix a,
                            bool ytrans, ph2matrix Y, pclusteroperator rwf, pclusteroperator cwf,
                            ptruncmode tm, real tol);

HEADER_PREFIX void
lowersolve_h2matrix(bool aunit, bool atrans, pch2matrix a,
  bool xytrans, ph2matrix X, pclusteroperator xrwf, pclusteroperator xcwf,
  ph2matrix Y, pclusteroperator yrwf, pclusteroperator ycwf,
  ptruncmode tm, real tol);

HEADER_PREFIX void
uppersolve_h2matrix(bool aunit, bool atrans, pch2matrix a,
  bool ytrans, ph2matrix Y, pclusteroperator yrwf, pclusteroperator ycwf,
  ptruncmode tm, real tol);

/* ------------------------------------------------------------
   LR decomposition
   ------------------------------------------------------------ */

HEADER_PREFIX void
lrdecomp_h2matrix(ph2matrix a, pclusteroperator arwf, pclusteroperator acwf,
                  ph2matrix L, pclusteroperator lrwf, pclusteroperator lcwf,
                  ph2matrix R, pclusteroperator rrwf, pclusteroperator rcwf,
                  ptruncmode tm, real tol);

HEADER_PREFIX void
lrsolve_h2matrix_avector(pch2matrix L, pch2matrix R, pavector x);


/* ------------------------------------------------------------
   Cholesky decomposition
   ------------------------------------------------------------ */

HEADER_PREFIX void
init_cholesky_h2matrix(ph2matrix a, pclusteroperator *prwf, pclusteroperator *pcwf,
                  ptruncmode tm);

HEADER_PREFIX void
choldecomp_h2matrix(ph2matrix a, pclusteroperator arwf, pclusteroperator acwf,
                  ph2matrix L, pclusteroperator lrwf, pclusteroperator lcwf,
                  ptruncmode tm, real tol);

HEADER_PREFIX void
cholsolve_h2matrix_avector(pch2matrix a, pavector x);

#endif