dclusterbasis.h

Go to the documentation of this file.

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

#ifndef DCLUSTERBASIS_H
#define DCLUSTERBASIS_H

typedef struct _dclusterbasis dclusterbasis;

typedef dclusterbasis *pdclusterbasis;

typedef const dclusterbasis *pcdclusterbasis;

#include "dcluster.h"
#include "dblock.h"
#include "amatrix.h"
#include "dclusteroperator.h"

struct _dclusterbasis
{
  pcdcluster t;

  uint directions;

  uint *k;

  uint *koff;

  uint ktree;

  uint kbranch;
        
  pamatrix V;

  pamatrix *E;

  uint sons;
  pdclusterbasis *son;
  uint **dirson;
};

/* ------------------------------------------------------------
 * Constructors and destructors
 * ------------------------------------------------------------ */

HEADER_PREFIX pdclusterbasis
init_dclusterbasis(pdclusterbasis cb, pcdcluster t);

HEADER_PREFIX void
uninit_dclusterbasis(pdclusterbasis cb);

HEADER_PREFIX pdclusterbasis
new_dclusterbasis(pcdcluster t);

HEADER_PREFIX void
del_dclusterbasis(pdclusterbasis cb);

/* ------------------------------------------------------------
 * Low-level management
 * ------------------------------------------------------------ */

HEADER_PREFIX void
update_dclusterbasis(pdclusterbasis cb);

HEADER_PREFIX void
setrank_dclusterbasis(pdclusterbasis cb, uint iota, uint k);

HEADER_PREFIX void
initmatrices_dclusterbasis(pdclusterbasis cb);

HEADER_PREFIX void 
  findranks_dclusterbasis(uint rank, pcdblock b, pdclusterbasis rb,
                          pdclusterbasis cb); 


/* ------------------------------------------------------------
 * Statistics
 * ------------------------------------------------------------ */
HEADER_PREFIX uint
getactives_dclusterbasis();

HEADER_PREFIX size_t
getsize_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX size_t
getsize_nonrecursive_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX uint
getmaxrank_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX uint
getactivedirections_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX void
print_tree_dclusterbasis(pcdclusterbasis cb);

/* ------------------------------------------------------------
 * Build cluster basis based on directional cluster tree
 * ------------------------------------------------------------ */

HEADER_PREFIX pdclusterbasis
buildfromdcluster_dclusterbasis(pcdcluster t);

/* ------------------------------------------------------------
 * Forward and backward transformation
 * ------------------------------------------------------------ */
HEADER_PREFIX pavector
newcoeffs_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX void
forward_dclusterbasis(pcdclusterbasis cb, pcavector x, pavector xt);

HEADER_PREFIX void
backward_dclusterbasis(pcdclusterbasis cb, pavector yt, pavector y);

HEADER_PREFIX void
slowforward_dclusterbasis(pcdclusterbasis cb, pcavector x, pavector xt);

HEADER_PREFIX void
slowbackward_dclusterbasis(pcdclusterbasis cb, pcavector yt, pavector y);

/* ------------------------------------------------------------
 * Forward and backward transformation for the root only
 * ------------------------------------------------------------ */

HEADER_PREFIX void
compress_dclusterbasis(pcdclusterbasis cb,
           pcavector xp, pavector xt);

HEADER_PREFIX void
expand_dclusterbasis(field alpha, pcdclusterbasis cb,
         pavector yt, pavector yp);

HEADER_PREFIX void
blockcompress_dclusterbasis(pcdclusterbasis cb,
          pcamatrix Xp, pamatrix Xt);

HEADER_PREFIX void
blockexpand_dclusterbasis(field alpha, pcdclusterbasis cb,
        pamatrix Xt, pamatrix Xp);

/* ------------------------------------------------------------
 * Enumeration by cluster number
 * ------------------------------------------------------------ */

HEADER_PREFIX pdclusterbasis *
enumerate_dclusterbasis(pcdcluster t, pdclusterbasis cb);

/* ------------------------------------------------------------
 * Hierarchical iterator
 * ------------------------------------------------------------ */

HEADER_PREFIX void
iterate_dclusterbasis(pdclusterbasis cb, uint tname,
          uint pardepth,
          void (*pre)(pdclusterbasis cb, uint tname,
          uint pardepth, void *data),
          void (*post)(pdclusterbasis cb, uint tname,
           uint pardepth, void *data),
          void *data);

/* ------------------------------------------------------------
 * Orthogonalization
 * ------------------------------------------------------------ */


HEADER_PREFIX void
ortho_dclusterbasis(pdclusterbasis cb, pdclusteroperator co);

HEADER_PREFIX void
weight_dclusterbasis_dclusteroperator(pdclusterbasis cb, pdclusteroperator co);

HEADER_PREFIX real
check_ortho_dclusterbasis(pcdclusterbasis cb);


HEADER_PREFIX pdclusterbasis
clone_structure_dclusterbasis(pcdclusterbasis cb);

HEADER_PREFIX pdclusterbasis
duplicate_dclusterbasis(pcdclusterbasis cb);

#endif