duniform.h

Go to the documentation of this file.

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

#ifndef DUNIFORM_H
#define DUNIFORM_H

typedef struct _duniform duniform;

typedef duniform *pduniform;

typedef const duniform *pcduniform;

#include "settings.h"
#include "amatrix.h"
#include "dclusterbasis.h"

struct _duniform {
  pdclusterbasis rb;
  pdclusterbasis cb;

  uint rd;
  uint cd;

  amatrix S;
};

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

HEADER_PREFIX pduniform
new_duniform(pdclusterbasis rb, uint rd, pdclusterbasis cb, uint cd);

HEADER_PREFIX void
del_duniform(pduniform u);

/* ------------------------------------------------------------
 * Statistics
 * ------------------------------------------------------------ */

HEADER_PREFIX size_t
getsize_duniform(pcduniform u);

/* ------------------------------------------------------------
 * Simple utility functions
 * ------------------------------------------------------------ */

HEADER_PREFIX void
clear_duniform(pduniform u);

HEADER_PREFIX void
copy_duniform(pcduniform src, pduniform trg);

/* ------------------------------------------------------------
 * Matrix-vector multiplication
 * ------------------------------------------------------------ */

HEADER_PREFIX void
fastaddeval_duniform_avector(field alpha, pcduniform u,
           pcavector xt, pavector yt);

HEADER_PREFIX void
fastaddevaltrans_duniform_avector(field alpha, pcduniform u,
          pcavector xt, pavector yt);

HEADER_PREFIX void
slowaddeval_duniform_avector(field alpha, pcduniform u,
           pcavector x, pavector y);

HEADER_PREFIX void
slowaddevaltrans_duniform_avector(field alpha, pcduniform u,
          pcavector x, pavector y);

/* ------------------------------------------------------------
 * Conversion to a full matrix
 * ------------------------------------------------------------ */

HEADER_PREFIX void
expand_duniform(field alpha, pcduniform u, pamatrix G);

/* ------------------------------------------------------------
 * Matrix norm 
 * ------------------------------------------------------------ */

HEADER_PREFIX real
norm2_fast_duniform(pcduniform u, pcdclusteroperator ro, pcdclusteroperator co);

HEADER_PREFIX real
normfrob_fast_duniform(pcduniform u, pcdclusteroperator ro, pcdclusteroperator co);


/* ---------------------------------------------------------- 
 * Projection
 * ---------------------------------------------------------- */

HEADER_PREFIX void
add_projected_duniform(pcduniform u, pcdclusteroperator ro, pcdclusteroperator co, pduniform unew);

#endif