bem2d.h

Go to the documentation of this file.

/* ------------------------------------------------------------
 * This is the file "bem2d.h" of the H2Lib package.
 * All rights reserved, Sven Christophersen 2011
 * ------------------------------------------------------------ */

#ifndef BEM2D_H_
#define BEM2D_H_

/* C STD LIBRARY */
#include <stdio.h>
/* CORE 0 */
#include "settings.h"
/* CORE 1 */
#include "amatrix.h"
#include "gaussquad.h"
#include "factorizations.h"
/* CORE 2 */
#include "cluster.h"
#include "clustergeometry.h"
#include "hmatrix.h"
#include "clusterbasis.h"
#include "h2matrix.h"
/* CORE 3 */
#include "hcoarsen.h"
#include "h2update.h"
#include "aca.h"
/* SIMPLE */
/* PARTICLES */
/* BEM */
#include "curve2d.h"
#include "singquad1d.h"

/* ------------------------------------------------------------ *
 * Type definitions                                             *
 * ------------------------------------------------------------ */

typedef struct _bem2d bem2d;
typedef bem2d *pbem2d;
typedef const bem2d* pcbem2d;

typedef struct _aprxbem2d aprxbem2d;
typedef aprxbem2d *paprxbem2d;
typedef const aprxbem2d *pcaprxbem2d;

typedef struct _parbem2d parbem2d;
typedef parbem2d *pparbem2d;
typedef const parbem2d *pcparbem2d;

typedef struct _kernelbem2d kernelbem2d;

typedef kernelbem2d *pkernelbem2d;

typedef const kernelbem2d *pckernelbem2d;

typedef void (*quadpoints2d)(pcbem2d bem, const real bmin[2],
    const real bmax[2], const real delta, real (**Z)[2], real (**N)[2]);

enum _basisfunctionbem2d {
  BASIS_NONE_BEM2D = 0,
  BASIS_CONSTANT_BEM2D = 'c',
  BASIS_LINEAR_BEM2D = 'l'
};

typedef enum _basisfunctionbem2d basisfunctionbem2d;

typedef field (*boundary_func2d)(const real *x, const real *n);

struct _bem2d {
pccurve2d gr;

psingquad1d sq;

uint N_neumann;

basisfunctionbem2d basis_neumann;

uint N_dirichlet;

basisfunctionbem2d basis_dirichlet;

real *mass;

field alpha;

uint *v2t;

uint *v2ts;

void (*nearfield)(const uint *ridx, const uint *cidx, pcbem2d bem, bool ntrans,
    pamatrix N);

void (*farfield_rk)(pccluster rc, uint rname, pccluster cc, uint cname,
    pcbem2d bem, prkmatrix R);

void (*farfield_u)(uint rname, uint cname, pcbem2d bem, puniform U);

void (*leaf_row)(pcbem2d bem, pclusterbasis rb, uint rname);

void (*leaf_col)(pcbem2d bem, pclusterbasis cb, uint cname);

void (*transfer_row)(pcbem2d bem, pclusterbasis rb, uint rname);

void (*transfer_col)(pcbem2d bem, pclusterbasis cb, uint cname);

paprxbem2d aprx;

pparbem2d par;

pkernelbem2d kernels;
};

struct _kernelbem2d {
void (*fundamental)(const real (*X)[2], const real (*Y)[2], pamatrix V);

void (*dny_fundamental)(const real (*X)[2], const real (*Y)[2],
    const real (*NY)[2], pamatrix V);

void (*dnx_dny_fundamental)(const real (*X)[2], const real (*NX)[2],
    const real (*Y)[2], const real (*NY)[2], pamatrix V);

void (*kernel_row)(const uint *idx, const real (*Z)[2], pcbem2d bem, pamatrix A);

void (*kernel_col)(const uint *idx, const real (*Z)[2], pcbem2d bem, pamatrix B);

void (*dnz_kernel_row)(const uint *idx, const real (*Z)[2], const real (*NZ)[2],
    pcbem2d bem, pamatrix A);

void (*dnz_kernel_col)(const uint *idx, const real (*Z)[2], const real (*NZ)[2],
    pcbem2d bem, pamatrix B);

void (*fundamental_row)(const uint *idx, const real (*Z)[2], pcbem2d bem,
    pamatrix A);

void (*fundamental_col)(const uint *idx, const real (*Z)[2], pcbem2d bem,
    pamatrix B);

void (*dnz_fundamental_row)(const uint *idx, const real (*Z)[2],
    const real (*NZ)[2], pcbem2d bem, pamatrix A);

void (*dnz_fundamental_col)(const uint *idx, const real (*Z)[2],
    const real (*NZ)[2], pcbem2d bem, pamatrix B);

void (*lagrange_row)(const uint *idx, pcavector px, pcavector py, pcbem2d bem,
    pamatrix V);

void (*lagrange_col)(const uint *idx, pcavector px, pcavector py, pcbem2d bem,
    pamatrix W);
};

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

HEADER_PREFIX pbem2d new_bem2d(pccurve2d gr);

HEADER_PREFIX void del_bem2d(pbem2d bem);

/* ------------------------------------------------------------
 * Methods to build cluster trees
 * ------------------------------------------------------------ */

HEADER_PREFIX pclustergeometry build_bem2d_const_clustergeometry(pcbem2d bem,
  uint **idx);

HEADER_PREFIX pclustergeometry build_bem2d_linear_clustergeometry(pcbem2d bem,
  uint **idx);

HEADER_PREFIX pclustergeometry build_bem2d_clustergeometry(pcbem2d bem,
  uint **idx, basisfunctionbem2d basis);

HEADER_PREFIX pcluster build_bem2d_cluster(pcbem2d bem, uint clf,
  basisfunctionbem2d basis);

/* ------------------------------------------------------------
 * Initializer functions for h-matrix approximations
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_recomp_bem2d(pbem2d bem, bool recomp,
  real accur_recomp, bool coarsen, real accur_coarsen);

/* ------------------------------------------------------------
 * Interpolation
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_aprx_inter_row_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, uint m);

HEADER_PREFIX void setup_hmatrix_aprx_inter_col_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, uint m);

HEADER_PREFIX void setup_hmatrix_aprx_inter_mixed_bem2d(pbem2d bem,
  pccluster rc, pccluster cc, pcblock tree, uint m);

/* ------------------------------------------------------------
 * Green
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_aprx_green_row_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, uint m, uint l, real delta,
  quadpoints2d quadpoints);

HEADER_PREFIX void setup_hmatrix_aprx_green_col_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, uint m, uint l, real delta,
  quadpoints2d quadpoints);

HEADER_PREFIX void setup_hmatrix_aprx_green_mixed_bem2d(pbem2d bem,
  pccluster rc, pccluster cc, pcblock tree, uint m, uint l, real delta,
  quadpoints2d quadpoints);

/* ------------------------------------------------------------
 * Greenhybrid
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_aprx_greenhybrid_row_bem2d(pbem2d bem,
  pccluster rc, pccluster cc, pcblock tree, uint m, uint l, real delta,
  real accur, quadpoints2d quadpoints);

HEADER_PREFIX void setup_hmatrix_aprx_greenhybrid_col_bem2d(pbem2d bem,
  pccluster rc, pccluster cc, pcblock tree, uint m, uint l, real delta,
  real accur, quadpoints2d quadpoints);

HEADER_PREFIX void setup_hmatrix_aprx_greenhybrid_mixed_bem2d(pbem2d bem,
  pccluster rc, pccluster cc, pcblock tree, uint m, uint l, real delta,
  real accur, quadpoints2d quadpoints);

/* ------------------------------------------------------------
 * ACA
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_aprx_aca_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, real accur);

HEADER_PREFIX void setup_hmatrix_aprx_paca_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, real accur);

/* ------------------------------------------------------------
 * HCA
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_hmatrix_aprx_hca_bem2d(pbem2d bem, pccluster rc,
  pccluster cc, pcblock tree, uint m, real accur);

/* ------------------------------------------------------------
 * Initializer functions for h2-matrix approximations
 * ------------------------------------------------------------ */

HEADER_PREFIX void setup_h2matrix_recomp_bem2d(pbem2d bem, bool hiercomp,
  real accur_hiercomp);

HEADER_PREFIX void setup_h2matrix_aprx_inter_bem2d(pbem2d bem,
  pcclusterbasis rb, pcclusterbasis cb, pcblock tree, uint m);

HEADER_PREFIX void setup_h2matrix_aprx_greenhybrid_bem2d(pbem2d bem,
  pcclusterbasis rb, pcclusterbasis cb, pcblock tree, uint m, uint l,
  real delta, real accur, quadpoints2d quadpoints);

HEADER_PREFIX void setup_h2matrix_aprx_greenhybrid_ortho_bem2d(pbem2d bem,
    pcclusterbasis rb, pcclusterbasis cb, pcblock tree, uint m, uint l,
    real delta, real accur, quadpoints2d quadpoints);

/* ------------------------------------------------------------
 * Fill hmatrix
 * ------------------------------------------------------------ */

HEADER_PREFIX void assemble_bem2d_hmatrix(pbem2d bem, pblock b, phmatrix G);

HEADER_PREFIX void assemblecoarsen_bem2d_hmatrix(pbem2d bem, pblock b,
  phmatrix G);

/* ------------------------------------------------------------
 * Fill h2-matrix
 * ------------------------------------------------------------ */

HEADER_PREFIX void assemble_bem2d_h2matrix_row_clusterbasis(pcbem2d bem,
  pclusterbasis rb);

HEADER_PREFIX void assemble_bem2d_h2matrix_col_clusterbasis(pcbem2d bem,
  pclusterbasis cb);

HEADER_PREFIX void assemble_bem2d_h2matrix(pbem2d bem, pblock b, ph2matrix G);

HEADER_PREFIX void assemblehiercomp_bem2d_h2matrix(pbem2d bem, pblock b,
  ph2matrix G);

/* ------------------------------------------------------------
 * Lagrange polynomials
 * ------------------------------------------------------------ */

HEADER_PREFIX void assemble_bem2d_lagrange_const_amatrix(const uint *idx,
  pcavector px, pcavector py, pcbem2d bem, pamatrix V);

HEADER_PREFIX void assemble_bem2d_dn_lagrange_const_amatrix(const uint *idx,
  pcavector px, pcavector py, pcbem2d bem, pamatrix V);

HEADER_PREFIX void assemble_bem2d_lagrange_amatrix(const real (*X)[2],
  pcavector px, pcavector py, pamatrix V);

/* ------------------------------------------------------------
 * some useful functions
 * ------------------------------------------------------------ */

HEADER_PREFIX void projectl2_bem2d_const_avector(pbem2d bem,
  field (*rhs)(const real *x, const real *n), pavector f);

HEADER_PREFIX void build_bem2d_rect_quadpoints(pcbem2d bem, const real a[2],
  const real b[2], const real delta, real (**Z)[2], real (**N)[2]);

HEADER_PREFIX void build_bem2d_circle_quadpoints(pcbem2d bem, const real a[2],
  const real b[2], const real delta, real (**Z)[2], real (**N)[2]);

HEADER_PREFIX prkmatrix build_bem2d_rkmatrix(pccluster row, pccluster col,
  void* data);

HEADER_PREFIX pamatrix build_bem2d_amatrix(pccluster row, pccluster col,
  void* data);

#endif /* BEM2D_H_ */