math_fce.h

Go to the documentation of this file.

/*!
 *
 * Copyright (C) 2015 Technical University of Liberec.  All rights reserved.
 * 
 * This program is free software; you can redistribute it and/or modify it under
 * the terms of the GNU General Public License version 3 as published by the
 * Free Software Foundation. (http://www.gnu.org/licenses/gpl-3.0.en.html)
 * 
 * This program is distributed in the hope that it will be useful, but WITHOUT
 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 * FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
 *
 * 
 * @file    math_fce.h
 * @brief   
 */
 
#ifndef MATH_H
#define MATH_H
 
// TODO: Only mh_fe_values depends on this file. This file can be removed.
 
//! small matrix types
typedef double SmallVec1_t[1];
typedef double SmallVec2_t[2];
typedef double SmallVec3_t[3];
typedef double SmallVec4_t[4];
 
typedef SmallVec1_t *SmallMtx1;
typedef SmallVec2_t *SmallMtx2;
typedef SmallVec3_t *SmallMtx3;
typedef SmallVec4_t *SmallMtx4;
 
 
//! Numerical helpers
#include <float.h>
// DBL_MIN approx.= 2.22507e-308, least nonzero double
 
// this should be used when we want to keep first digit nonzero
// such number have full double precision, smaller numbers loose precision
#define NUM_ZERO DBL_MIN/DBL_EPSILON
 
// each use of ESP(value) indicate that we use some small value
// this small value should be relative to the number we want to compare
#define EPS(value) (value)
 
// OBSOLETE
#define ZERO        EPS(1e-12)
 
#define DBL_EQ(i,j) (fabs((i)-(j))<NUM_ZERO)
#define DBL_GE(i,j) ((i)>(j)-NUM_ZERO)
#define DBL_LE(i,j) ((i)<(j)+NUM_ZERO)
#define DBL_GT(i,j) ((i)>(j)+NUM_ZERO)
#define DBL_LT(i,j) ((i)<(j)-NUM_ZERO)
 
//! Usefull math macros
#define SQUARE(x)   ((x) * (x))
#define SGN(x)      ( ((x)>ZERO)? (1) :  ( ((x)<(-ZERO))? (-1) : (0) ) )
#define SUBDET2(i,j,k,l) (a[(i)][(k)]*a[(j)][(l)]-a[(i)][(l)]*a[(j)][(k)])
#ifndef M_PI
    #define M_PI 3.14159265358979323846264338327950288f
#endif
 
// vector functions
double vector_length(double[]);
double scalar_product(double[],double[]);
void normalize_vector(double[]);
void scale_vector(double[],double);
 
void vector_product(double[],double[],double[]);
void vector_difference(double[],double[],double[]);
// small matrix operations
double Det3( SmallMtx3 a);
double MatrixInverse(double *a,double *b,int size);
void PrintSmallMatrix(double *mtx, int size );
#endif
//-----------------------------------------------------------------------------
// vim: set cindent: