harith2.h

Go to the documentation of this file.

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

#ifndef HARITH2_H
#define HARITH2_H

typedef struct _haccum haccum;

typedef haccum *phaccum;

typedef const haccum *pchaccum;

typedef struct _hprodentry hprodentry;

typedef hprodentry *phprodentry;

#include "harith.h"

struct _haccum {
  phmatrix z;

  pctruncmode tm;

  real eps;
  
  prkmatrix r;

  phprodentry xy;
};

HEADER_PREFIX phaccum
new_haccum(phmatrix z, pctruncmode tm, real eps);

HEADER_PREFIX void
del_haccum(phaccum ha);

HEADER_PREFIX void
addproduct_haccum(field alpha,
      bool xtrans, pchmatrix x, bool ytrans, pchmatrix y,
      phaccum za);

HEADER_PREFIX phaccum *
split_haccum(phmatrix z, phaccum za);

HEADER_PREFIX void
flush_haccum(phaccum ha);

HEADER_PREFIX void
add_amatrix_destructive_rkmatrix(field alpha, bool atrans, pamatrix a,
         pctruncmode tm, real eps, prkmatrix b);

HEADER_PREFIX void
add_amatrix_destructive_hmatrix(field alpha, bool atrans, pamatrix a,
        pctruncmode tm, real eps, phmatrix b);

HEADER_PREFIX void
addmul2_hmatrix(field alpha,
    bool xtrans, pchmatrix x, bool ytrans, pchmatrix y,
    pctruncmode tm, real eps, phmatrix z);

HEADER_PREFIX void
lowersolve_haccum(bool aunit, bool atrans, pchmatrix a,
      bool xtrans, phaccum xa);

HEADER_PREFIX void
lowersolve2_hmatrix_hmatrix(bool aunit, bool atrans, pchmatrix a,
          pctruncmode tm, real eps,
          bool xtrans, phmatrix x);

HEADER_PREFIX void
uppersolve_haccum(bool aunit, bool atrans, pchmatrix a,
      bool xtrans, phaccum xa);

HEADER_PREFIX void
uppersolve2_hmatrix_hmatrix(bool aunit, bool atrans, pchmatrix a,
          pctruncmode tm, real eps,
          bool xtrans, phmatrix x);

HEADER_PREFIX void
lrdecomp_haccum(phaccum aa);

HEADER_PREFIX void
lrdecomp2_hmatrix(phmatrix a, pctruncmode tm, real eps);

HEADER_PREFIX void
choldecomp_haccum(phaccum aa);

HEADER_PREFIX void
choldecomp2_hmatrix(phmatrix a, pctruncmode tm, real eps);

#endif