# IEEE 754 and float 範例程式

功能：互轉浮點數與二進制表示法<br>

```c

/* ieee754.c  
 * Author: Aaron Liao (http://applezulab.netdpi.net)
 * Last updated: 2010/12/29 
 * IEEE754 sample code
 * LICENSE: GPL
 * 32 bits:  [1][8][23]
 * 64 bits:  [1][11][52]
 */

#include <strings.h>
#include <string.h>
#include <stdio.h>

int ieee754(char *data, double real)
{
    int exp=0, i=0, exp_tmp;
    double cal, frag_tmp;
    char buf[32];

    bzero(buf, sizeof(buf));

    cal = real;
    /* get sign */
    if( cal < 0 ){

        buf[0] = 1;
        cal *= (double) -1;

    } else {

        buf[0] = 0;
    }

    /* get exp */
    while( cal >= 1.0 ) {
        cal /= 2;
        exp++;
    }

    for( i=8, exp_tmp=1; i>0; i--){

        if( (exp & exp_tmp) ) {

            buf[i] = 1;

        } else {

            buf[i] = 0;

        }

        exp_tmp *=2 ;
    }

    /* get fraction */
    for(i=9, frag_tmp = (double) 1/2; i< 32; i++) {

        if( (cal - frag_tmp) >= 0 ){
            buf[i] = 1;
            cal = (double) cal - (double) frag_tmp;
        } else {
            buf[i] = 0;
        }
        frag_tmp /= (double) 2;

    }

    /* To ascii */
    for(i=0; i<32; i++) {

        buf[i] += '0';
    }

    strncpy(data, buf, 32);

    return 0;
}

double rev_ieee754(char *buf)
{
    int exp=0, i=0, exp_tmp, sign;
    double real, frag_tmp;
    int len = 32;

    /* get sign */
    if( buf[0] == '1' ){

        sign = -1;

    } else {

        sign = 1;
    }

    /* get exp */
    for( exp=0, i=8, exp_tmp=1; i>0; i--){

        if( buf[i] == '1') {
            exp += exp_tmp;
        }

        exp_tmp*=2;
    }

    /* get fraction */
    for(i=9, real=(double)0, frag_tmp = (double) 1/2; i< len; i++) {

        if( buf[i] == '1') {
            real += frag_tmp;
        }
        
        frag_tmp /= (double) 2;

    }

    while(exp>0) {

        real *= 2;
        exp--;
    }

    real *= (double) sign;

    return real;
}

int ieee754_64(char *data, double real)
{
    int exp=0, i=0, exp_tmp;
    double cal, frag_tmp;
    static char buf[64];

    bzero(buf, sizeof(buf));

    cal = real;
    /* get sign */
    if( cal < 0 ){

        cal *= (double) -1;
        buf[0] = 1;

    } else {

        buf[0] = 0;
    }

    /* get exp */
    while( cal >= 1.0 ) {
        cal /= 2;
        exp++;
    }

    for( i=11, exp_tmp=1; i>0; i--){

        if( (exp & exp_tmp) ) {

            buf[i] = 1;

        } else {

            buf[i] = 0;

        }

        exp_tmp *=2 ;
    }

    /* get fraction */
    for(i=12, frag_tmp = (double) 1/2; i< 64; i++) {

        if( (cal - frag_tmp) >= 0 ){
            buf[i] = 1;
            cal = (double) cal - (double) frag_tmp;
        } else {
            buf[i] = 0;
        }
        frag_tmp /= (double) 2;

    }

    /* To ascii */
    for(i=0; i<64; i++) {

        buf[i] += '0';
    }

    strncpy(data, buf, 64);

    return 0;
}

double rev_ieee754_64(char *buf)
{
    int exp=0, i=0, exp_tmp, sign;
    double real, frag_tmp;
    int len=64;

    /* get sign */
    if( buf[0] == '1' ){

        sign = -1;

    } else {

        sign = 1;
    }

    /* get exp */
    for( exp=0, i=11, exp_tmp=1; i>0; i--){

        if( buf[i] == '1') {
            exp += exp_tmp;
        }

        exp_tmp*=2;
    }

    /* get fraction */
    for(i=12, real=(double)0, frag_tmp = (double) 1/2; i< len; i++) {

        if( buf[i] == '1') {
            real += frag_tmp;
        }
        
        frag_tmp /= (double) 2;

    }

    while(exp>0) {

        real *= 2;
        exp--;
    }

    real *= (double) sign;

    return real;
}

/* test */
int main(void)
{
    int i, j;
    char buf[32], buf2[64];

    for(i=1; i<=1000000; i*=2) {
        
        bzero(buf, sizeof(buf));
        bzero(buf2, sizeof(buf2));

        ieee754( (char*)buf, (double) i);
        printf("[32bit] %20f ->\t", rev_ieee754( (char*)buf));;

        for(j=0; j<sizeof(buf); j++) {
            printf("%c ", buf[j]);
        }
        printf("\n");

        ieee754_64( buf2, (double) i);
        printf("[64bit] %20f ->\t", rev_ieee754_64( (char*)buf2));;

        for(j=0; j<sizeof(buf2);j++) {
            printf("%c ", buf2[j]);
        }
        printf("\n");
    }

    return 0;
}
```

<br>


---

# Agent Instructions: Querying This Documentation

If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question.

Perform an HTTP GET request on the current page URL with the `ask` query parameter:

```
GET https://applezulab.netdpi.net/linux-prog/ieee-754-and-float-sample-code.md?ask=<question>
```

The question should be specific, self-contained, and written in natural language.
The response will contain a direct answer to the question and relevant excerpts and sources from the documentation.

Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.