Can't Copy and Paste this?

Click here for a copy-and-paste friendly version of this code!

//**************************************

//     

//INCLUDE files for :Convert MSBIN forma

//     t to/from IEEE float

//**************************************

//     

/* +++Date last modified: 05-Jul-1997 */

/*

** SNIPMATH.H - Header file for SNIPPETS math functions and macros

*/

#ifndef SNIPMATH__H

#define SNIPMATH__H

#include 

#include "sniptype.h"

#include "round.h"

/*

** Callable library functions begin here

*/

voidSetBCDLen(int n); /* Bcdl.C */

longBCDtoLong(char *BCDNum); /* Bcdl.C */

voidLongtoBCD(long num, char BCDNum[]);/* Bcdl.C */

double bcd_to_double(void *buf, size_t len, /* Bcdd.C */

int digits);

int double_to_bcd(double arg, char *buf, /* Bcdd.C */

size_t length, size_t digits );

DWORDncomb1 (int n, int m);/* Combin.C*/

DWORDncomb2 (int n, int m);/* Combin.C*/

voidSolveCubic(double a, double b, double c, /* Cubic.C*/

double d, int *solutions,

double *x);

DWORDdbl2ulong(double t); /* Dbl2Long.C */

longdbl2long(double t);/* Dbl2Long.C */

double dround(double x); /* Dblround.C */

/* Use #defines for Permutations and Combinations -- Factoryl.C */

#define log10P(n,r) (log10factorial(n)-log10factorial((n)-(r)))

#define log10C(n,r) (log10P((n),(r))-log10factorial(r))

double log10factorial(double N); /* Factoryl.C */

double fibo(unsigned short term);/* Fibo.C */

double frandom(int n);/* Frand.C*/

double ipow(double x, int n);/* Ipow.C */

int ispow2(int x);/* Ispow2.C*/

longdouble ldfloor(long double a);/* Ldfloor.C */

int initlogscale(long dmax, long rmax);/* Logscale.C */

longlogscale(long d); /* Logscale.C */

floatMSBINToIEEE(float f); /* Msb2Ieee.C */

floatIEEEToMSBIN(float f); /* Msb2Ieee.C */

int perm_index (char pit[], int size);/* Perm_Idx.C */

int round_div(int n, int d); /* Rnd_Div.C */

longround_ldiv(long n, long d);/* Rnd_Div.C */

double rad2deg(double rad); /* Rad2Deg.C */

double deg2rad(double deg); /* Rad2Deg.C */

#include "pi.h"

#ifndef PHI

#define PHI ((1.0+sqrt(5.0))/2.0) /* the golden number*/

#define INV_PHI (1.0/PHI) /* the golden ratio */

#endif

/*

** File: ISQRT.C

*/

    struct int_sqrt {

    unsigned sqrt,

    frac;

};

void usqrt(unsigned long x, struct int_sqrt *q);

#endif /* SNIPMATH__H */

code:

Can't Copy and Paste this?

Click here for a copy-and-paste friendly version of this code!

Terms of Agreement:   

By using this code, you agree to the following terms...   

1) You may use this code in your own programs (and may compile it into a program and distribute it in compiled format for langauges that allow it) freely and with no charge.   

2) You MAY NOT redistribute this code (for example to a web site) without written permission from the original author. Failure to do so is a violation of copyright laws.   

3) You may link to this code from another website, but ONLY if it is not wrapped in a frame. 

4) You will abide by any additional copyright restrictions which the author may have placed in the code or code's description.  

//     

/* +++Date last modified: 05-Jul-1997 */

/*** MSBIN conversion routines***/

/*** public domain by Jeffery Foy ***/

#include "snipmath.h"

    union Converter {

    unsigned char uc[10];

    unsigned int ui[5];

    unsigned long ul[2];

    float f[2];

    double d[1];

};

/* MSBINToIEEE - Converts an MSBIN floating point number */

/*to IEEE floating point format*/

/**/

/* Input: f - floating point number in MSBIN format */

/* Output: Same number in IEEE format*/

float MSBINToIEEE(float f)

    {

    union Converter t;

    int sign, exp;/* sign and exponent */

    t.f[0] = f;

    /* extract the sign & move exponent bias from 0x81 to 0x7f */

    sign = t.uc[2] / 0x80;

    exp = (t.uc[3] - 0x81 + 0x7f) & 0xff;

    /* reassemble them in IEEE 4 byte real number format */

    t.ui[1] = (t.ui[1] & 0x7f) | (exp << 7) | (sign << 15);

    return t.f[0];

} /* End of MSBINToIEEE */

/* IEEEToMSBIN - Converts an IEEE floating point number */

/*to MSBIN floating point format */

/**/

/* Input: f - floating point number in IEEE format */

/* Output: Same number in MSBIN format*/

float IEEEToMSBIN(float f)

    {

    union Converter t;

    int sign, exp;/* sign and exponent */

    t.f[0] = f;

    /* extract sign & change exponent bias from 0x7f to 0x81 */

    sign = t.uc[3] / 0x80;

    exp = ((t.ui[1] >> 7) - 0x7f + 0x81) & 0xff;

    /* reassemble them in MSBIN format */

    t.ui[1] = (t.ui[1] & 0x7f) | (sign << 7) | (exp << 8);

    return t.f[0];

} /* End of IEEEToMSBIN */