tet3dp1.h

Go to the documentation of this file.

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

#ifndef TET3DP1_H
#define TET3DP1_H

typedef struct _tet3dp1 tet3dp1;

typedef tet3dp1 *ptet3dp1;

typedef const tet3dp1 *pctet3dp1;

#include "sparsematrix.h"
#include "tet3d.h"
#include "clustergeometry.h"

struct _tet3dp1 {
  pctet3d gr;

  uint ndof;

  uint nfix;

  bool *is_dof;

  uint *idx2dof;
};

HEADER_PREFIX ptet3dp1
new_tet3dp1(pctet3d gr);

HEADER_PREFIX void
del_tet3dp1(ptet3dp1 dc);

HEADER_PREFIX psparsematrix
build_tet3dp1_sparsematrix(pctet3dp1 dc);

HEADER_PREFIX psparsematrix
build_tet3dp1_interaction_sparsematrix(pctet3dp1 dc);

HEADER_PREFIX psparsematrix
build_tet3dp1_prolongation_sparsematrix(pctet3dp1 dfine, pctet3dp1 dcoarse,
         pctet3dref rf);

HEADER_PREFIX void
assemble_tet3dp1_laplace_sparsematrix(pctet3dp1 dc,
       psparsematrix A, psparsematrix Af);

HEADER_PREFIX void
assemble_tet3dp1_mass_sparsematrix(pctet3dp1 dc,
          psparsematrix M, psparsematrix Mf);

HEADER_PREFIX void
assemble_tet3dp1_dirichlet_avector(pctet3dp1 dc,
      field (*f)(const real *x, void *fdata), void *fdata,
      pavector d);

HEADER_PREFIX void
assemble_tet3dp1_functional_avector (pctet3dp1 dc,
       field (*f)(const real *x, void *fdata), void *fdata,
       pavector b);

HEADER_PREFIX real
normmax_tet3dp1(pctet3dp1 dc,
    field (*f)(const real *x, void *fdata), void *fdata,
    pavector xs);

HEADER_PREFIX real
norml2_tet3dp1(pctet3dp1 dc,
         field (*f)(const real *x, void *fdata), void *fdata,
         pcavector xs, pcavector xf);

/* ------------------------------------------------------------
   Clustergeometry
   ------------------------------------------------------------ */

HEADER_PREFIX pclustergeometry
build_tet3dp1_clustergeometry(pctet3dp1 p1, uint *idx);


#endif