whitney.h

// copyright by benedikt oswald, all rights reserved, 2002-2006.
// project - aqhdes2003
// file name - whitney.h
// file type - C include file
// objective - povide Whitney base functionality
// creation - 2003 feb 03 ~ 18:12:00 by benedikt oswald
// modified - 2004 feb 12 ~ 15:15:00 by benedikt oswald, transferred functions into auxiliary.h
// required software - UG 3.1 or higher
// file id -
// feature - provides types, macros and functions definitions required
// feature - to implement finite element base functions of the Whitney 1-forms
// feature - and also Whitney 2-forms

#ifndef AQ_WHITNEY_ELEMENT_BASE_FUNCTIONS
#define AQ_WHITNEY_ELEMENT_BASE_FUNCTIONS

#ifndef __COMPILER__
#include "compiler.h" /* compiler specific data type definitions */
#endif

/* include standars C header files */
#include <stddef.h>
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>

/* include UG related header files */
#include "gm.h"          /* inserted as is from osterz          */
#include "parallel.h"
#include "evm.h"         /* inserted as is from osterz          */
#include "shapes.h"      /* inserted as is from osterz          */
#include "misc.h"        /* inserted as is from osterz          */
#include "ugdevices.h"   /* inserted as is from osterz          */
#include "np.h"          /* inserted as is from osterz          */            
#include "assemble.h"    /* inserted as is from osterz          */
#include "disctools.h"   /* inserted as is from osterz          */
#include "general.h"     /* inserted as is from osterz          */
#include "ugstruct.h"    /* inserted as is from osterz          */
#include "quadrature.h"  /* include Gaussian quadrature routine */
#include "udm.h"         /* contains VECDATA_DESC data structure */
#include "compiler.h"
#include "namespace.h"


/* include our own header files */
#include "aqhdPhysicoMath.h"
#include "disconst.h"


USING_UG_NAMESPACE
USING_UGDIM_NAMESPACE


/* PUBLIC CONSTANT DECLARATION */

/* define numerical constants */

#define WHITNEY1ST 0
#define WHITNEY2ND 1

#define NUM_WHITNEY_ELEM_NODES (DIM+1)
#define NUM_WHITNEY_ELEM_EDGES (3*(DIM-1))
#define NUM_WHITNEY_ELEM_FACETS (DIM+1)

#define NUM_WHITNEY_ELEM_EDGE_COEFFS NUM_WHITNEY_ELEM_EDGES
#define NUM_WHITNEY_ELEM_FACET_COEFFS NUM_WHITNEY_ELEM_FACETS

#define WHITNEY_FE_OK 0

#define THREE 3

/* define function constants */

#define KP3D(k,x,y)    k[0]=(x[1]*y[2]) - (x[2]*y[1]); k[1]=-((x[0]*y[2]) - (x[2]*y[0])); k[2]=(x[0]*y[1]) - (x[1]*y[0]) /* 3D cross product */
//
//#define DOTPRODUCT2D(s,a,b) s=a[0]*b[0] + a[1]*b[1]
#define DOTPRODUCT3D(s,a,b) s=a[0]*b[0] + a[1]*b[1] + a[2]*b[2]
//
//#define VECLEN2D(l,x)       l=sqrt(x[0]*x[0] + x[1]*x[1]); /* 2D vector length */
#define VECLEN3D(l,x)       l=sqrt(x[0]*x[0] + x[1]*x[1] + x[2]*x[2]); /* 3D vector length */
//
//#define VECSCALE2D(v,x,s)   v[0]=x[0]*s;v[1]=x[1]*s /* scale 2D vector by factor s */ 
#define VECSCALE3D(v,x,s)   v[0]=x[0]*s;v[1]=x[1]*s; v[2]=x[2]*s /* scale 3D vector by factor s */ 
//
#define ELEMID(t)           (t)->ge.id
//#define NODEID(t)           (t)->id

/* PUBLIC VARIABLE DECLARATION */


/* PUBLIC FUNCTION DECLARATIONS */

/* integrate Whitney 1st kind of edge e along edge ae using quadrature formula provided by UG      */
DOUBLE ConsistencyW1(const ELEMENT *t, /* for which element                                        */
                     INT e,            /* Whitney 1st kind of edge e                               */
                     INT ae            /* along edge which we want to integrate Whitney 1st kind   */
    );

/* integrate Whitney function of the 2nd kind on facet (l,m,n) of element t along edge with corners i and j */
DOUBLE IntegrateAlongEdge(
    const ELEMENT *t,
    INT eos,      /* Whitney function of edge e or facet s        */
    INT e,        /* edge along which element shall be integrated */
    INT whichkind /* which Whitney function: 1st or 2nd kind      */
    );

/* integrate Whitney function of the 2nd kind on facet (l,m,n) of element t along edge with corners i and j */
DOUBLE IntegrateOverFacet(
    const ELEMENT *t,
    INT eos,      /* Whitney function of edge e or facet s        */
    INT f,        /* facet over which base function shall be integrated */
    INT whichkind /* which Whitney function: 1st or 2nd kind      */
    );

/* function returns barycentric coordinates of a point specified in 3D cartesian coordinates */
void BarycentricCoordinates(
    const ELEMENT *t,
    const DOUBLE *cartcoords,
    DOUBLE *barycoords
    );

/*
  function calculates gradient that belongs to barycentric coordinate attached
  to corner c === side s with respect to the numbering scheme within the element t
*/

void BarycentricGradient(
    const ELEMENT *t,
    const INT s,
    DOUBLE *gradient
    );

/* function computes the value of the whitney function of the 1st kind along the edge with corners c0 and c1 */

void Whitney1(
    ELEMENT *t,
    INT e,
    const DOUBLE *localcoords,
    DOUBLE *w1
    );

/* function computes the curl of the Whitney 1st aling edge e, expressed in terms of the respective gradients */

void CurlWhitney1(ELEMENT *t,
                  const INT e,
                  DOUBLE *curl
    );

/* function computes the value of the whitney function of the 2nd kind */

/* void Whitney2( */
/*     ELEMENT *t, */
/*     const INT s, */
/*     const DOUBLE *cartcoords, */
/*     DOUBLE *whitney2 */
/*     ); */

#endif

---