singquad2d.h

Go to the documentation of this file.

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

#ifndef SINGQUAD2D_H_
#define SINGQUAD2D_H_

/* C STD LIBRARY */
#include <assert.h>
/* CORE 0 */
#ifdef USE_SIMD
#include "simd.h"
#endif
#include "settings.h"
/* CORE 1 */
#include "gaussquad.h"
/* CORE 2 */
/* CORE 3 */
/* SIMPLE */
/* PARTICLES */
/* BEM */
#include "surface3d.h"

struct _singquad2d {
  real *x_id;
  real *y_id;
  real *w_id;
  real base_id;
  uint n_id;
  real *x_edge;
  real *y_edge;
  real *w_edge;
  real base_edge;
  uint n_edge;
  real *x_vert;
  real *y_vert;
  real *w_vert;
  real base_vert;
  uint n_vert;
  real *x_dist;
  real *y_dist;
  real *w_dist;
  real base_dist;
  uint n_dist;
  real *x_single;
  real *y_single;
  real *w_single;
  real base_single;
  uint n_single;

#ifdef USE_TRIQUADPOINTS

  real *tri_x;

  real *tri_y;

  real *tri_z;
#endif

  uint q;
  uint q2;
  uint nmax;
};

typedef struct _singquad2d singquad2d;

typedef singquad2d* psingquad2d;

typedef const singquad2d *pcsingquad2d;

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

HEADER_PREFIX psingquad2d
build_singquad2d(pcsurface3d gr, uint q, uint q2);

HEADER_PREFIX void
del_singquad2d(psingquad2d sq);

/* ------------------------------------------------------------
 Weighting quadrature rules
 ------------------------------------------------------------ */

HEADER_PREFIX void
weight_basisfunc_ll_singquad2d(real *x, real *y, real *w, uint nq);

HEADER_PREFIX void
weight_basisfunc_cl_singquad2d(real *x, real *y, real *w, uint nq);

HEADER_PREFIX void
weight_basisfunc_lc_singquad2d(real *x, real *y, real *w, uint nq);

HEADER_PREFIX void
weight_basisfunc_l_singquad2d(real *x, real *y, real*w, uint nq);

/* ------------------------------------------------------------
 Select quadrature rule
 ------------------------------------------------------------ */

HEADER_PREFIX uint
fast_select_quadrature(uint (*geo_t)[3], uint t, uint s);


HEADER_PREFIX uint
select_quadrature_singquad2d(pcsingquad2d sq, const uint *tv, const uint *sv,
    uint *tp, uint *sp, real **x, real **y, real **w, uint *n, real *base);

#endif /* SINGQUAD2D_H_ */