c:/Users/efeyhl/devs/math2/c_basic_math.h

Go to the documentation of this file.

#ifndef __C_BASIC_MATH_H__
#define __C_BASIC_MATH_H__

#ifdef  __cplusplus
extern "C" {
#endif

#include <stdio.h> //FILE

#ifndef __INT64_DEFINED__

#ifdef WIN32
typedef __int64 int64; 
typedef unsigned __int64 uint64; 
#else
typedef long long int64; 
typedef unsigned long long uint64; 
#endif //WIN32

#define __INT64_DEFINED__

#endif //__INT64_DEFINED__


#define EXA (1.E18)

#define PETA (1.E15)

#define TERA (1.E12)

#define GIGA (1.E9)

#define MEGA (1.E6)

#define KILO (1.E3)

#define MILLI (1.E-3)

#define MICRO (1.E-6)

#define NANO (1.E-9)

#define PICO (1.E-12)

#define FEMTO (1.E-15)

#define ATTO (1.E-18)


#define ONE_DIV_3 (1.0 / 3.0)

#define ONE_DIV_4 (1.0 / 4.0)


#define PI (3.1415926535897932384626433832795)

#define PI1 (PI / 180.0)

#define PI18 (PI / 10.0)

#define PI20 (PI / 9.0)

#define PI30 (PI / 6.0)

#define PI36 (PI / 5.0)

#define PI40 (2.0 * PI / 9.0)

#define PI45 (PI / 4.0)

#define PI54 (3.0 * PI / 10.0)

#define PI60 (PI / 3.0)

#define PI72 (2.0 * PI / 5.0)

#define PI90 (PI / 2.0)

#define PI108 (3.0 * PI / 5.0)

#define PI120 (2.0 * PI / 3.0)

#define PI135 (3.0 * PI / 4.0)

#define PI180 (PI)

#define PI240 (2.0 * PI120)

#define PI270 (3.0 * PI / 2.0)

#define PI360 (2.0 * PI)

#define DEG1 (180.0 / PI)

#define SQRT2 (1.4142135623730950488016887242097)

#define SQRT3 (1.7320508075688772935274463415059)


typedef union
{

  int status;
  //switch off compiler warning for union without name
#pragma warning( disable : 4201 )
  struct // status as bit-field
  {

    int undef : 1;


    int overflow : 1;


    int divideZero : 1;


    int zeroDivZero : 1;


    int baseExpZero : 1;


    int zeroExpNeg : 1;


    int negRadicand : 1;


    int trigonArcErr : 1;


    int trigonCotErr : 1;


    int rootExpZero : 1;


    int logBaseZero : 1;


    int logBaseOne : 1;


    int int64overflow : 1;


    int strErr : 1;


    int intPowNegExp : 1;


    int unused : 15;


    int underflow : 1;


    int negative : 1;
  }; //end of struct for the status-bits
  //take the default for compiler warning "union without name"
#pragma warning( default : 4201 )
}BM_STATUS; //end of unnamed union for status


void printBM_STATUS( const BM_STATUS bms, FILE* file );


int isDouble( const double* value, BM_STATUS* bms );


int isI64( const int64* value, BM_STATUS* bms );


double doubleAddDouble(
  const double augend,
  const double addend,
  BM_STATUS* bms );


int64 i64AddI64(
  const int64 augend,
  const int64 addend,
  BM_STATUS* bms );


double doubleSubDouble(
  const double minuend,
  const double subtrahend,
  BM_STATUS* bms );


int64 i64SubI64(
  const int64 minuend,
  const int64 subtrahend,
  BM_STATUS* bms );


double doubleMulDouble(
  const double multiplier,
  const double multiplicand,
  BM_STATUS* bms );



int64 i64MulI64(
  const int64 multiplier,
  const int64 multiplicand,
  BM_STATUS* bms );


double doubleDivDouble(
  const double numerator,
  const double denominator,
  BM_STATUS* bms );


int64 i64DivI64(
  const int64 numerator,
  const int64 denominator,
  BM_STATUS* bms );


double doubleModDouble(
  const double numerator,
  const double denominator,
  BM_STATUS* bms );


int64 i64ModI64(
  const int64 numerator,
  const int64 denominator,
  BM_STATUS* bms);


double reciprocDouble(
  const double value,
  BM_STATUS* bms );


double ceilDouble( const double value, BM_STATUS* bms );


double floorDouble( const double value, BM_STATUS* bms );


int64 doubleToInt64( const double value, double* intDbl, BM_STATUS* bms );


double i64ToDouble( const int64 value );


double roundDouble( const double value, const int decPlaces, BM_STATUS* bms );


int doubleCompareDouble( const double value1, const double value2, const double diff, BM_STATUS* bms );


int i64CompareI64( const int64 value1, const int64 value2, const int64 diff );


int isEvenI64( const int64 value );


double positiveDouble( const double value, BM_STATUS* bms );


int64 positiveI64( const int64 value, BM_STATUS* bms );


double negativeDouble( const double value, BM_STATUS* bms );


int64 negativeI64( const int64 value );


double invertSignDouble( const double value, BM_STATUS* bms );


int64 invertSignI64( const int64 value, BM_STATUS* bms );


double doublePow_2( const double base, BM_STATUS* bms );


int64 i64Pow_2( const int64 base, BM_STATUS* bms );


double doublePow_3( const double base, BM_STATUS* bms );


int64 i64Pow_3( const int64 base, BM_STATUS* bms );


double doublePow_4( const double base, BM_STATUS* bms );


int64 i64Pow_4( const int64 base, BM_STATUS* bms );


double doublePowDouble(
  const double base,
  const double exponent,
  BM_STATUS* bms ); //status is optional (can be null)


int64 i64PowI64(
  const int64 base,
  const int64 exponent,
  BM_STATUS* bms ); //status is optional (can be null)


double doubleRoot_2(
  const double radicand,
  BM_STATUS* bms );


int64 i64Root_2(const int64 radicand, BM_STATUS* bms);


double doubleRoot_3(
  const double radicant,
  BM_STATUS* bms );


double doubleRoot_4(
  const double radicant,
  BM_STATUS* bms );


double doubleRootDouble(
  const double radicant,
  const double rootExponent,
  BM_STATUS* bms );


int64 i64RootI64(
  const int64 radicant,
  const int64 rootExponent,
  BM_STATUS* bms );


double doubleLogDouble( 
  const double base,
  const double value,
  BM_STATUS* bms ); //status is optional



int64 i64LogI64( 
  const int64 base,
  const int64 value,
  BM_STATUS* bms ); //status is optional


void sortAscendantDouble(
  double* values,
  const unsigned long count,
  BM_STATUS* bms );


void sortAscendantI64(
  int64* values,
  const unsigned long count );


long doubleSolveP1norm(
  double* xr,
  const double* pf );


long doubleSolveP1(
  double* xr,
  const double* pf );


long doubleSolveP2norm(
  double* xr,
  const double* pf );


long doubleSolveP2(
  double* xr,
  const double* pf );


long doubleSolveP3norm(
  double* xr,
  const double* pf );


long doubleSolveP3(
  double* xr,
  const double* pf );


long doubleSolveP4norm(
  double* xr,
  const double* pf );


long doubleSolveP4(
  double* xr,
  const double* pf );


long doubleGetPolyFactors1(
  double* pf,
  const double* xz);


BM_STATUS doubleGetPolyFactors2(
  double* pf,
  const double* xz);


BM_STATUS doubleGetPolyFactors3(
  double* pf,
  const double* xz);


BM_STATUS doubleGetPolyFactors4(
  double* pf,
  const double* xz);


double deg2RadDouble( const double degree, BM_STATUS* bms );


double rad2DegDouble( const double radian, BM_STATUS* bms );


double stringToDouble( const char* stringValue, BM_STATUS* bms ); //status is optional


int64 stringToI64( const char* stringValue, char** stopChar, BM_STATUS* bms );


double strToDouble(
  const char* str,
  BM_STATUS* bms );


BM_STATUS doubleSolveL1(
  double* xr,
  const double* lf );


BM_STATUS doubleSolveL2(
  double* xr,
  double* lf[2] );


void printDoubleBitpattern( double value, FILE* file );


void printDoubleBitpatternFormatted( double value, FILE* file );


double sinDouble( const double rad, BM_STATUS* bms );


double asinDouble( const double value, BM_STATUS* bms );


double cosDouble( const double rad, BM_STATUS* bms );


double acosDouble( const double value, BM_STATUS* bms );


double tanDouble( const double rad, BM_STATUS* bms );


double atanDouble( const double value, BM_STATUS* bms );


double cotDouble( const double rad, BM_STATUS* bms );


double acotDouble( const double value, BM_STATUS* bms );



int findDouble( const double* element1, const unsigned int count, const double elem );


int findI64( const __int64* element1, const unsigned int count, const __int64 elem );

#ifdef  __cplusplus
} //end of extern "C"
#endif

#endif //__C_BASIC_MATH_H__

---