macros.h

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           macros.h
// author(s):      Jan Novak, Jan Zeman
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 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 Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************
#ifndef MUMECH_SHARE_H
#define MUMECH_SHARE_H

// There is conflict name XMLDocument in msxml and tinyxml.
// Solution 1:
// This define have to be before first win32 include.
#define WIN32_LEAN_AND_MEAN


/* Default header files */
#include <stdio.h>
#include <stdlib.h>
#include <iostream>  // Standard header defining std::cout
#include <float.h>   //real number limits
#include <limits.h>  //integer number limits
#include "gelib.h"

using namespace gelibspace;

namespace mumech {

#define _NOCASE

#ifdef _NOCASE
#define CASE false
#define _STRCMP _STRCASECMP
#define _STRNCMP strncmp
#else
#define CASE true
#define _STRCMP strcmp
#define _STRNCMP strncmp
#endif



/* Symbolic constants */
#define MACHINE_EPS 1E-16 //alternatively FLT_EPSILON = 1e-5
#define MACHINE_EPS_PRN 1E-8 //alternatively FLT_EPSILON = 1e-5

#define PI 3.14159265358979323846

#define       SQRT_2 1.41421356237309504880168872
#define _1_BY_SQRT_2 0.70710678118654752440084436

/* File manipulating macros */
#define FileError(fmode,name) fprintf(stderr,"\033[31mCannot %s file %s in file %s line %d\033[0m\n", fmode, name, __FILE__,__LINE__); exit( 1 );
#define FOPEN(file, name, mode) if( ( file = fopen(name,mode) ) == NULL) {FileError("open", name)}
#define FCLOSE(file) if (( fclose(file) ) == EOF) {FileError("close","")}

#define FREAD( var, number, file ) fread ( &var, sizeof ( var ), number, file )
#define FWRITE( var, number, file ) fwrite ( &var, sizeof ( var ), number, file )

/* Error functions */
//#define Error(string,value) fprintf( stderr, "\033[31mError : " );fprintf(stderr,string,value); fprintf(stderr,"in file %s line %d.\033[0m\n", __FILE__, __LINE__ );exit(FALSE)                   
#define NotImplementedFunction( name ) _errorr2( "Function %s : not implemented yet\n", name );
#define Dprintf( format, var ) printf (" %s = ", #var ); printf ( format, var )
#define Dfprintf( file, format, var ) fprintf (file," %s = ", #var ); fprintf ( file, format, var )
//#define Warning(string,value) fprintf( stderr, "\033[31mWarning : " ); fprintf(stderr,string,value); fprintf(stderr,"in file %s line %d.\033[0m\n", __FILE__, __LINE__ )
#define NfoPrint(string,value) fprintf( stderr, "\033[31mNotice : \033[0m" ); fprintf(stderr,string,value);

/* Other macros */
#ifndef SWAP
#define SWAP(a, b, dummy) dummy=a;a=b;b=dummy
#endif

#define SQR( a ) ( ( a ) * ( a ) )
#define CUB( a ) ( ( a ) * ( a ) * ( a ) )

#define MIN(a,b) ( ( ( a ) < ( b ) ) ? ( a ) : ( b ) )
#define MAX(a,b) ( ( ( a ) > ( b ) ) ? ( a ) : ( b ) )

#define MID(a,b,c) MIN( MIN(MAX(a,b), MAX(b,c)), MIN(MAX(b,c), MAX(c,a)) )

#ifndef CDMAX //collation of complex numbers' distances, requires <cmath> to be included
#define CDMAX(a,b) (((SQR(real(a))+SQR(imag(a)))>(SQR(real(b))+SQR(imag(b)))) ? (a) : (b))
#endif

#ifndef CReMAX //collation of complex numbers' real parts, requires <cmath> to be included
#define CReMAX(a,b) ( ( real( a ) > real( b ) ) ? ( a ) : ( b ) )
#endif

#ifndef CImMAX //collation of complex numbers' imaginary parts, requires <cmath> to be included
#define CImMAX(a,b) ( ( imag( a ) > imag( b ) ) ? ( a ) : ( b ) )
#endif

#ifndef ABS
#define ABS( a ) ( ( ( a ) > 0 ) ? ( a ) : ( -1 * ( a ) ) )
#endif

//#ifndef sgn
//#define sgn( a ) ( ( ( a ) < 0 ) ? -1 : 1 )
//#endif

#ifndef odd_or_even
#define odd_or_even( a ) ( ( ( ( a ) % 2 ) > 0 ) ? ODD : EVEN )
#endif

/* angle unit transformation*/

#ifndef DEG2RAD
#define DEG2RAD( a ) ( ( a ) * PI / 180. )
#endif

#ifndef RAD2DEG
#define RAD2DEG( a ) ( ( a ) * 180. / PI )
#endif

/* infinity number constant */

#ifndef INFTY
//#define INFTY FLT_MAX //alternatively DBL_MAX for douple precision
#define INFTY 1e+36 //alternatively DBL_MAX for douple precision
#endif


/* memory alocation macros */
/*-------------------------*/
//#define EXIT_SEQUENCE printf("\033[31mExiting\033[0m\n" ); printf( "\033[31mUnexpected exit!\033[0m\n" ); exit( 1 );
//extern void exit_throw( void );

//#define AllocError(var) printf("\033[31mCannot allocate %s in file %s line %d, stack's capacity exceeded\033[0m\n", #var, __FILE__, __LINE__ ); EXIT_SEQUENCE;
//#define NotNULLError(var) printf("\033[31mYou tried to allocate not NULL variable %s in file %s line %d\033[0m\n",#var,__FILE__,__LINE__);

/* Memory manipulating macros */
//#define MALLOC(var,number,type) if (var != NULL){NotNULLError(var)} else if( (var = (type*) calloc( (number), sizeof(type) ) ) == NULL ){AllocError(var)}
//#define FREE(var) if((var) != NULL) {free((void *) var); var = NULL;}
//#define REALLOC( var, number, type ) if ( ( var = ( type * ) realloc( var, (number) * sizeof( type ) ) ) == NULL ){AllocError(var)}
//#define MALLOC_IF_NULL(var,number,type) if (var == NULL) if( (var = (type*) calloc( (number), sizeof(type) ) ) == NULL ){AllocError(var)}

/* "Safe" allocations */
//#define MALLOC_IFNOT_EMPTY(var, number, type) if ( (number) != 0 ) { MALLOC( var, number, type ); }
//#define FREE_IFNOT_EMPTY(var, number)  if ( (number) != 0 ) { FREE ( var ); number = 0;}

/* "Vector" memory allocation */
//#define _CreateVector(type,var,number) if (var != NULL) {NotNULLError(var)} else if ( (var = (type *) calloc((number),sizeof(type))) == NULL ){AllocError(var)}
//#define _DeleteVector(var) if ((var) != NULL) {free( (void*) var); var=NULL;}

/* "Matrix" memory allocation */
//#define _CreateArrayStatic(type,name,rows,columns) {int _ii_; for(_ii_=0; _ii_<rows;_ii_++){ _CreateVector(type,name[_ii_],columns)}}
//#define _DeleteArrayStatic(name,rows) {int _ii_; for( _ii_=0; _ii_<rows;_ii_++) _DeleteVector(name[_ii_]);}

/* "Matrix" memory allocation - pure dynamic */
//#define _CreateArray(type,name,rows,columns) {if (name) NotNULLError(name) else {MALLOC(name,rows,type*); for (int _ii_=0; _ii_<rows;_ii_++) {name[_ii_]=NULL;_CreateVector(type,name[_ii_],columns)}}}
//#define _DeleteArray(name,rows) if ( name != NULL ) {for(int _ii_=0; _ii_<rows;_ii_++) _DeleteVector(name[_ii_]); _DeleteVector( name );}

/* "pointer" alocation */
//#define _CreatePointer(type,var) _CreateVector(type,var,1)

/* delete pointer - the same as _DeleteVector */
//#define _DeletePointer(var) _DeleteVector(var)



/* kronecker delta */
#ifndef delta
#define delta( i, j )  ( ( i == j ) ? 1. : 0. )
#endif

/* continuous reading from string (as from the file) */
#ifndef jump_empties
#define jump_empties( p_str )  while( ( * p_str ) <= ' ' && ( * p_str ) != '\0' ){ p_str++; }
#endif

#ifndef jump_word
#define jump_word( p_str )  while( ( * p_str ) > ' ' && ( * p_str ) != '\0' ){ p_str++; }
#endif

#ifndef sscanff
#define sscanff( p_str, format, p_var ) jump_empties( p_str ); sscanf( p_str, format, p_var ); jump_word( p_str )
#endif

//vector length
#define VEC_SQR_NORM_2D( x, y )    ( SQR( x ) + SQR( y ) )

#define VEC_SQR_NORM_3D( x, y, z ) ( SQR( x ) + SQR( y ) + SQR( z ) )
#define VEC_NORM( x, y, z ) sqrt( VEC_SQR_NORM_3D( x, y, z ) )

#define DIST_POINTS_SQR_2D( P1, P2 ) ( SQR( P1[0] - P2[0] ) + SQR( P1[1] - P2[1] ) )
#define DIST_POINTS_SQR_3D( P1, P2 ) ( SQR( P1[0] - P2[0] ) + SQR( P1[1] - P2[1] ) + SQR( P1[2] - P2[2] ) )


} // end of namespace mumech

#endif

/*end of file*/