truncation.h

Go to the documentation of this file.

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

#ifndef TRUNCATION_H
#define TRUNCATION_H

typedef struct _truncmode truncmode;

typedef truncmode *ptruncmode;

typedef const truncmode *pctruncmode;

#include "realavector.h"
#include "settings.h"

/* ------------------------------------------------------------
 Truncation strategy
 ------------------------------------------------------------ */

struct _truncmode {
  bool frobenius;
  bool absolute;
  bool blocks;

  real zeta_level;
  real zeta_age;
};

/* ------------------------------------------------------------
 Constructors for standard truncation strategies
 ------------------------------------------------------------ */

HEADER_PREFIX ptruncmode
new_truncmode();

HEADER_PREFIX void
del_truncmode(ptruncmode tm);

HEADER_PREFIX ptruncmode
new_releucl_truncmode();

HEADER_PREFIX ptruncmode
new_relfrob_truncmode();

HEADER_PREFIX ptruncmode
new_blockreleucl_truncmode();

HEADER_PREFIX ptruncmode
new_blockrelfrob_truncmode();

HEADER_PREFIX ptruncmode
new_abseucl_truncmode();

/* ------------------------------------------------------------
 Find minimal acceptable rank
 ------------------------------------------------------------ */

HEADER_PREFIX uint
findrank_truncmode(pctruncmode tm, real eps, pcrealavector sigma);

#endif