tet3drt0.h

Go to the documentation of this file.

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

#ifndef TET3DRT0_H
#define TET3DRT0_H

typedef struct _tet3drt0 tet3drt0;

typedef tet3drt0 *ptet3drt0;

typedef const tet3drt0 *pctet3drt0;

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

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

struct _tet3drt0
{
  pctet3d t3;  
  
  uint ndof;    
  
  uint nfix;      
  
  uint *is_dof; 
  
  uint *idx2dof; 
};

HEADER_PREFIX ptet3drt0
new_tet3drt0(pctet3d gr);

HEADER_PREFIX void 
update_tet3drt0(ptet3drt0 rt0);

HEADER_PREFIX void
del_tet3drt0(ptet3drt0 rt);

HEADER_PREFIX psparsematrix
build_tet3drt0_A_sparsematrix(pctet3drt0 dc);

HEADER_PREFIX psparsematrix
build_tet3drt0_B_sparsematrix(pctet3drt0 dc);

HEADER_PREFIX psparsematrix
build_tet3drt0_A_interaction_sparsematrix(pctet3drt0 dc);


HEADER_PREFIX psparsematrix
build_tet3drt0_B_interaction_sparsematrix(pctet3drt0 dc);

HEADER_PREFIX void
assemble_tet3drt0_darcy_A_sparsematrix(pctet3drt0 dc, psparsematrix A, psparsematrix Af, pavector K);

HEADER_PREFIX void
assemble_tet3drt0_darcy_B_sparsematrix(pctet3drt0 dc, psparsematrix A, psparsematrix Af);

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


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

void assemble_tet3drt0_g_N_avector(ptet3drt0 dc, 
           field (*f)(const uint *e, void *data), void *data,
           pavector g);

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

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

HEADER_PREFIX real
norml2_flux_centroid_tet3drt0(pctet3drt0 dc,
         void (*q) (const real * x, void *fdata, pavector v), void *fdata,
         pcavector x1, pcavector g);

HEADER_PREFIX real
norml2_flux_facemidpoint_tet3drt0(pctet3drt0 dc,
         void (*q) (const real * x, void *fdata, pavector v), void *fdata,
         pcavector x1, pcavector g);


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

HEADER_PREFIX pclustergeometry 
build_tet3drt0_A_clustergeometry(pctet3drt0 rt, uint *idx);

HEADER_PREFIX pclustergeometry 
build_tet3drt0_B_clustergeometry(pctet3drt0 rt, uint *idx);

#endif