tri2drt0.h

Go to the documentation of this file.

/* ------------------------------------------------------------
 * This is the file "tri2drt0.h" of the H2Lib package.
 * All rights reserved, Nadine Albrecht 2015
 * ------------------------------------------------------------ */

#ifndef TRI2DRT0_H
#define TRI2DRT0_H

typedef struct _tri2drt0 tri2drt0;

typedef tri2drt0 *ptri2drt0;

typedef const tri2drt0 *pctri2drt0;

#include "sparsematrix.h"
#include "tri2d.h"
#include "clustergeometry.h"

/* ------------------------------------------------------------
   Edge-elements for a triangular mesh
   ------------------------------------------------------------ */

struct _tri2drt0
{
  pctri2d t2;     
  
  uint ndof; 
  
  uint nfix;     
  
  uint *is_dof;   
  
  uint *idx2dof;  
};

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

HEADER_PREFIX ptri2drt0
new_tri2drt0(pctri2d rt);

HEADER_PREFIX void
update_tri2drt0(ptri2drt0 rt0);

HEADER_PREFIX void
del_tri2drt0(ptri2drt0 dc);

HEADER_PREFIX psparsematrix
build_tri2drt0_A_sparsematrix(pctri2drt0 dc);

HEADER_PREFIX psparsematrix
build_tri2drt0_B_sparsematrix(pctri2drt0 dc);

HEADER_PREFIX psparsematrix
build_tri2drt0_A_interaction_sparsematrix(pctri2drt0 dc);

HEADER_PREFIX psparsematrix
build_tri2drt0_B_interaction_sparsematrix(pctri2drt0 dc);


HEADER_PREFIX void
assemble_tri2drt0_darcy_A_sparsematrix(pctri2drt0 dc, psparsematrix A, psparsematrix Af, pavector K);

HEADER_PREFIX void
assemble_tri2drt0_darcy_B_sparsematrix(pctri2drt0 dc, psparsematrix A, psparsematrix Af);


HEADER_PREFIX void
assemble_tri2drt0_b_D_avector(pctri2drt0 dc,
      field (*f)(const real *e, void *fdata), void *fdata,
      pavector d);

HEADER_PREFIX void
assemble_tri2drt0_b_f_avector(ptri2drt0 dc,
                   field (*f)(const real *x, void *fdata), void *fdata,
                   pavector fv);

void assemble_tri2drt0_g_N_avector(ptri2drt0 dc, 
           field (*f)(const uint *e, void *data), void *data,
           pavector g);

HEADER_PREFIX real
norml2_pressure_centroid_tri2drt0(pctri2drt0 dc,
         field(*f) (const real * x, void *fdata), void *fdata,
         pcavector x2);

HEADER_PREFIX real
norml2_pressure_edgemidpoint_tri2drt0(pctri2drt0 dc,
         field(*f) (const real * x, void *fdata), void *fdata, pcavector x2);

HEADER_PREFIX real
norml2_flux_centroid_tri2drt0(pctri2drt0 dc,
         void (*f) (const real * x, void *fdata, pavector v), void *fdata,
         pcavector x1, pcavector g);

HEADER_PREFIX real
norml2_flux_edgemidpoint_tri2drt0(pctri2drt0 dc,
         void (*f) (const real * x, void *fdata, pavector v), void *fdata,
         pcavector x1, pcavector g);


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

HEADER_PREFIX pclustergeometry 
build_tri2drt0_A_clustergeometry(pctri2drt0 rt, uint *idx);

HEADER_PREFIX pclustergeometry 
build_tri2drt0_B_clustergeometry(pctri2drt0 rt, uint *idx);

#endif