怎么将ibm浮点数转换成float浮点，请高手指点，最好用c#

dorothyorsusie

最近想把segy中的ibm浮点数转换成float浮点，很多高手已经做过这个了吧，求指点

Franklin

我也想做

dorothyorsusie

Franklin 发表于 2012-9-18 20:22
我也想做

把低位的字节和高位的调换，然后转换成float

Franklin

dorothyorsusie 发表于 2012-9-19 09:05
把低位的字节和高位的调换，然后转换成float

我从su里边找了个ibm_to_float解决了

dorothyorsusie

Franklin 发表于 2012-9-19 21:44
我从su里边找了个ibm_to_float解决了

su是什么？请讲中文，朋友

Franklin

su是cwp（克罗拉多矿院）的一款开源软件

Franklin

dorothyorsusie 发表于 2012-9-19 22:14
su是什么？请讲中文，朋友

su是cwp（克罗拉多矿院）的一款开源软件

dorothyorsusie

Franklin 发表于 2012-9-21 19:25
su是cwp（克罗拉多矿院）的一款开源软件

能不能贴出来看看，怎么写的

Franklin

dorothyorsusie 发表于 2012-9-24 10:34
能不能贴出来看看，怎么写的

static void ibm_to_float(int from[], int to[], int n, int endian)
/***********************************************************************
ibm_to_float - convert between 32 bit IBM and IEEE floating numbers
************************************************************************
Input::
from                input vector
to                output vector, can be same as input vector
endian                byte order =0 little endian (DEC, PC's)
                            =1 other systems
*************************************************************************
Notes:
Up to 3 bits lost on IEEE -> IBM

Assumes sizeof(int) == 4

IBM -> IEEE may overflow or underflow, taken care of by
substituting large number or zero

Only integer shifting and masking are used.
*************************************************************************
Credits: CWP: Brian Sumner,  c.1985
*************************************************************************/
{
    register int fconv, fmant, i, t;

    for (i = 0;i < n; ++i) {

        fconv = from;

        /* if little endian, i.e. endian=0 do this */
        if (endian == 0) fconv = (fconv << 24) | ((fconv >> 24) & 0xff) |
                ((fconv & 0xff00) << 8) | ((fconv & 0xff0000) >> 8);

        if (fconv) {
            fmant = 0x00ffffff & fconv;
            /* The next two lines were added by Toralf Foerster */
            /* to trap non-IBM format data i.e. conv=0 data  */
            if (fmant == 0)
                warn("mantissa is zero data may not be in IBM FLOAT Format !");
            t = (int) ((0x7f000000 & fconv) >> 22) - 130;
            while (!(fmant & 0x00800000)) { --t; fmant <<= 1; }
            if (t > 254) fconv = (0x80000000 & fconv) | 0x7f7fffff;
            else if (t <= 0) fconv = 0;
            else fconv =   (0x80000000 & fconv) | (t << 23)
                         | (0x007fffff & fmant);
        }
        to = fconv;
    }
    return;
}

dorothyorsusie

Franklin 发表于 2012-9-24 16:34
static void ibm_to_float(int from[], int to[], int n, int endian)
/****************************** ...

谢谢了