auxiliary.h

// copyright by benedikt oswald, all rights reserved, 2002-2006.
// project - aqhdes2003
// file name - auxiliary.c
// file type - C include file file
// objective - povide Whitney base functionality
// creation - 2004 jan 31 ~ 00:00:00 by holger gerhards
// modified - 2004 feb 12 ~ 15:15:00 by benedikt oswald, adapted to overall solver structure
// modified - 2004 feb 23 ~ 14:00:00 by holger gerhards, include new functions (efld2dslice, efld3dvoxel,GlobalVectorOut)
// modified - 2004 feb 23 ~ 14:00:00 by holger gerhards, include modi for some functions
// modified - 2004 mar 04 ~ 13:00:00 by holger gerhards, include new functions (IntegerOut, DoubleOut)
// modified - 2004 mar 08 ~ 10:30:00 by holger gerhards, added Integer "timestep" in GlobalMatrixOut, 
//                                                             GlobalVectorOut, IntegerOut, DoubleOut  
// required software - UG 3.1 or higher
// file id -
// feature - implements function declarations from whitney.h

#ifndef __AUX__
#define __AUX__

#include <stddef.h>
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>

#include "gm.h"
#include "parallel.h"
#include "evm.h"
#include "shapes.h"
#include "misc.h"
#include "ugdevices.h"
#include "np.h"
#include "assemble.h"
#include "disctools.h"
#include "general.h"
#include "ugstruct.h"
#include "quadrature.h"
#include "cmdline.h"
#include "udm.h"
#include "commands.h"
#include "cmdint.h"

#include "aqhdPhysicoMath.h"
#include "aqhdSysParams.h"
#include "aqhdMaterials.h"
#include "aqhdBaseTypes.h"
#include "disconst.h"
#include "whitney.h"

/* PUBLIC CONSTANT DECLARATION */

#define NORM_OK   0
#define NORM_BAD  1
#define DIST_OK   0

#define MATHEMATICA      1
#define ONLY_NUM         2
#define ABSOLUTE         3
#define X_COMPONENT      4
#define Y_COMPONENT      5
#define Z_COMPONENT      6
#define ALL_COMPONENTS   7
#define BRACKETS         8

#define CHECK_LOCAL      9
#define EFIELD          10
#define OUTPUT_ONLY     11

#define MAT_CLEAR       12
#define MAT_SCALE       13
#define MAT_ADD         14

/* function recuperates the number of the corner that is opposite a given side in a tetrahedron */
INT CornerOppositeToSideOnTetrahedron(ELEMENT *t,INT s);

/* function recuperates the number of the side that is opposite a given corner in a tetrahedron */
INT SideOppositeToCorner(ELEMENT *t,INT i);

/* proprietary implementation of global-to-local coordinate conversion function */
INT GlobalToLocalTetrahedron(INT ncoe,DOUBLE **cornercoord, DOUBLE *globalcoord,DOUBLE *localcoord);

/* proprietary implementation of FindElementFromPosition */
ELEMENT *FindTetrahedronFromPosition (GRID *theGrid,DOUBLE *pos);

/* define a normal vector on a boundary segment */
INT NormalVecOnBounSeg(INT bn , DOUBLE_VECTOR norm);

/* define a distance to a boundary segment    */
INT DistanceToBounSegments( DOUBLE *cartcoords , DOUBLE *distance );

/* give output of a local matrix */
void LocalMatrixOut(ELEMENT *t, MATDATA_DESC *M, INT mode, char *name, char *zusatz);

/* give output of a local vector */
void LocalVectorOut(ELEMENT *t, VECDATA_DESC *V, INT mode, char *name, char *zusatz);

/* gives global matrix out */
INT GlobalMatrixOut(MATDATA_DESC *A, INT mode, INT timestep, char *name, char *zusatz);

/* gives global vector out */
INT GlobalVectorOut(VECDATA_DESC *A, INT mode, INT timestep, char *name, char *zusatz); 

/* gives global matrix out - for ARPACK */
INT GlobalMatrixARPACK(MATDATA_DESC *A, char *name, INT mode, INT anzahl); 

/* modify matrix according to mode */
INT ModifyMatrix(MATDATA_DESC *A, INT mode, DOUBLE value);

/* global matrix clear */
INT aqmatclear(MATDATA_DESC *A);

/* gives integer value out */
void IntegerOut(INT a, INT mode, INT timestep, char *name, char *zusatz);

/* gives double value out */
void DoubleOut(DOUBLE a, INT mode, INT timestep, char *name, char *zusatz);

#endif

---