[R-gregmisc-users] SF.net SVN: r-gregmisc: [1137] trunk/SASxport/src
Brought to you by:
warnes
Menu
▾
▴
[R-gregmisc-users] SF.net SVN: r-gregmisc: [1137] trunk/SASxport/src
|
From: <wa...@us...> - 2007-08-12 03:13:37
|
Revision: 1137
http://r-gregmisc.svn.sourceforge.net/r-gregmisc/?rev=1137&view=rev
Author: warnes
Date: 2007-08-11 20:13:36 -0700 (Sat, 11 Aug 2007)
Log Message:
-----------
Remove new stuff... use SAS's code instead
Removed Paths:
-------------
trunk/SASxport/src/B8.h
trunk/SASxport/src/IEEEtoIBM.c
trunk/SASxport/src/MASKS.h
trunk/SASxport/src/main.c
trunk/SASxport/src/reverse.c
Deleted: trunk/SASxport/src/B8.h
===================================================================
--- trunk/SASxport/src/B8.h 2007-08-12 03:12:21 UTC (rev 1136)
+++ trunk/SASxport/src/B8.h 2007-08-12 03:13:36 UTC (rev 1137)
@@ -1,60 +0,0 @@
-/* Binary constant generator macro
-By Tom Torfs - donated to the public domain
-*/
-
-#include <sys/types.h>
-
-
-/* All macro's evaluate to compile-time constants */
-
-/* *** helper macros *** /
-
-/* turn a numeric literal into a hex constant
-(avoids problems with leading zeroes)
-8-bit constants max value 0x11111111, always fits in unsigned long
-*/
-#define HEX__(n) 0x##n##LU
-
-/* 8-bit conversion function */
-#define B8__(x) ((x&0x0000000FLU)?1:0) \
-+((x&0x000000F0LU)?2:0) \
-+((x&0x00000F00LU)?4:0) \
-+((x&0x0000F000LU)?8:0) \
-+((x&0x000F0000LU)?16:0) \
-+((x&0x00F00000LU)?32:0) \
-+((x&0x0F000000LU)?64:0) \
-+((x&0xF0000000LU)?128:0)
-
-/* *** user macros *** /
-
-/* for upto 8-bit binary constants */
-#define B8(d) ((unsigned char)B8__(HEX__(d)))
-
-/* for upto 16-bit binary constants, MSB first */
-#define B16(dmsb,dlsb) (((unsigned short)B8(dmsb)<<8) \
-+ B8(dlsb))
-
-/* for upto 32-bit binary constants, MSB first */
-#define B32(dmsb,db2,db3,dlsb) ( \
-+ ((unsigned long)B8(dmsb)<<24) \
-+ ((unsigned long)B8(db2)<<16) \
-+ ((unsigned long)B8(db3)<<8) \
-+ B8(dlsb))
-
-/* for upto 32-bit binary constants, MSB first */
-#define B64(dmsb,db2,db3,db4,db5,db6,db7,dlsb) ( \
-+ ((unsigned int64_t)B8(dmsb<<56) \
-+ ((unsigned int64_t)B8(db2)<<48) \
-+ ((unsigned int64_t)B8(db3)<<40) \
-+ ((unsigned int64_t)B8(db4 <<32) \
-+ ((unsigned int64_t)B8(dm5)<<24) \
-+ ((unsigned int64_t)B8(db6)<<16) \
-+ ((unsigned int64_t)B8(db7)<<8 ) \
-+ B8(dlsb))
-
-
-/* Sample usage:
-B8(01010101) = 85
-B16(10101010,01010101) = 43605
-B32(10000000,11111111,10101010,01010101) = 2164238933
-*/
Deleted: trunk/SASxport/src/IEEEtoIBM.c
===================================================================
--- trunk/SASxport/src/IEEEtoIBM.c 2007-08-12 03:12:21 UTC (rev 1136)
+++ trunk/SASxport/src/IEEEtoIBM.c 2007-08-12 03:13:36 UTC (rev 1137)
@@ -1,109 +0,0 @@
-//#include "writeSAS.h"
-
-#include <stdio.h>
-#include <string.h>
-#include <stdlib.h>
-#import <assert.h>
-#import <sys/types.h>
-
-typedef struct { /* IBM floating point format */
- unsigned sign : 1; /* Sign bit */
- unsigned exponent: 7; /* Exponent */
- //unsigned fraction: 56; /* Fraction */
- unsigned fraction1: 24; /* Top half of fraction */
- unsigned fraction2: 32; /* Lower half of fraction */
-} IBM_fp_struct;
-
-
-typedef struct { /* IEEE floating point format */
- unsigned sign : 1; /* Sign bit */
- unsigned exponent: 11; /* Exponent */
- // unsigned fraction: 52; /* Fraction */
- unsigned fraction1: 20; /* Top half of fraction */
- unsigned fraction2: 32; /* Lower half of fraction */
-} IEEE_fp_struct;
-
-
-int IEEEtoIBM(double *ieee, double *ibm )
-{
- IEEE_fp_struct *ieee_fp = (IEEE_fp_struct*) ieee;
- IBM_fp_struct *ibm_fp = (IBM_fp_struct* ) ibm;
- int64_t sign;
- int64_t exponent;
- int64_t fraction;
- short low_exp_bits;
-
- /*** Extract Pieces ***/
- // Sign = Ieee & Ieee_Sign;
- // Exponent = *Ieee & Ieee_Exp;
- // Fraction = *Ieee & Ieee_Frac;
-
- printf("1\n");
-
- /*** copy IEEE sign to IBM sign ***/
- ibm_fp->sign = ieee_fp->sign;
-
- printf("2\n");
-
- /*** convert IEEE exponent to IBM exponent ***/
- exponent = ieee_fp->exponent;
-
- printf("3\n");
-
- /* store lowest 2 bits from exponent */
- low_exp_bits = exponent % 16;
-
-
- printf("4\n");
-
- /* divide exponent by 4 to get from IEEE pow(2) exponent to IBM pow(16) exponent */
- exponent = exponent >> 2;
-
- printf("5\n");
-
-
- /* put it into the return value */
- ibm_fp->exponent = exponent;
-
- printf("6\n");
-
-
- /*** convert IEEE fraction to IBM fraction */
- fraction = ((int64_t) ieee_fp->fraction1 ) << 32 | ((int64_t) ieee_fp->fraction2);
-
- printf("7\n");
-
-
- /* shift left by low order bits lost from exponent */
- fraction = fraction << (low_exp_bits && 0x03);
-
-
- printf("8\n");
-
- /* add leading 1 bit */
- fraction >> 1;
- fraction = fraction | ( ((int64_t) 1) <<55 );
-
- printf("9\n");
-
-
- ibm_fp->fraction1 = (int32_t) (fraction >> 32) | 0x0FFFFF;
- ibm_fp->fraction2 = (int32_t) fraction;
-
-
- printf("10\n");
-
-
- char *buf= (char*) ibm;
- char tmp;
- int i;
- for(i=0; i<8; i++)
- {
- tmp = buf[7-i];
- buf[7-i] = buf[i];
- buf[i] = tmp;
- }
-
- return 0;
-}
-
Deleted: trunk/SASxport/src/MASKS.h
===================================================================
--- trunk/SASxport/src/MASKS.h 2007-08-12 03:12:21 UTC (rev 1136)
+++ trunk/SASxport/src/MASKS.h 2007-08-12 03:13:36 UTC (rev 1137)
@@ -1,23 +0,0 @@
-/* Masks for IBM fields */
-#define IBM_SIGN B64(10000000, 0000000, 0000000, 00000000, \
- 00000000, 0000000, 0000000, 00000000)
-
-#define IBM_EXP B64(01111111, 0000000, 0000000, 00000000, \
- 00000000, 0000000, 0000000, 00000000)
-
-#define IBM_FRAC B64(00000000, 1111111, 1111111, 11111111, \
- 11111111, 1111111, 1111111, 11111111)
-
-#define IBM_FRAC_MSB B64(00000000, 1000000, 0000000, 00000000, \
- 00000000, 0000000, 0000000, 00000000)
-
-/* Masks for IEEE fields */
-#define IEEE_SIGN B64(10000000, 0000000, 0000000, 00000000, \
- 00000000, 0000000, 0000000, 00000000)
-
-#define IEEE_EXP B64(01111111, 1111000, 0000000, 00000000, \
- 00000000, 0000000, 0000000, 00000000)
-
-#define IEEE_FRAC B64(00000000, 0000111, 1111111, 11111111, \
- 11111111, 1111111, 1111111, 11111111)
-
Deleted: trunk/SASxport/src/main.c
===================================================================
--- trunk/SASxport/src/main.c 2007-08-12 03:12:21 UTC (rev 1136)
+++ trunk/SASxport/src/main.c 2007-08-12 03:13:36 UTC (rev 1137)
@@ -1,55 +0,0 @@
-#include <math.h>
-#include <stdio.h>
-#include <string.h>
-#include <stdlib.h>
-#import <assert.h>
-#import <sys/types.h>
-
-
-extern int cnxptiee(char *from, int fromtype, char *to, int totype);
-int IEEEtoIBM(double *ieee, double *ibm );
-
-int main(int argc, char *argv)
-{
- double ieee;
- double ibm_sas, ieee_sas;
- double ibm_our, ieee_our;
- char *ibm_sas_c = (char*) &ibm_sas;
- char *ibm_our_c = (char*) &ibm_our;
-
-
- for(int i=0; i<10; i++)
- {
- ieee = pow(2.0, ((double) i)) + 1.0 + 1.0 / ( ((double) i) + 1.0);
-
- /* from ieee to ibm, SAS code */
- cnxptiee( (char*) &ieee, 0, (char*) &ibm_sas, 1 );
-
- /* and back again, SAS code */
- cnxptiee( (char*) &ibm_sas, 1, (char*) &ieee_sas, 0 );
-
-
- /* from ieee to ibm, our code */
- IEEEtoIBM( &ieee, &ibm_our );
-
- /* and back again, SAS code */
- cnxptiee( (char*) &ibm_our, 1, (char*) &ieee_our, 0 );
-
-
- printf("i=%5d, ieee=%10f, ieee_sas=%10f ieee_our=%10f \n", i,
- ieee, ieee_sas, ieee_our);
-
- for(int index=0; index<8; index++)
- {
- printf(" ");
- printf("%2x ", ibm_sas_c[index]);
- printf("%2x ", ibm_our_c[index]);
- printf("\n");
- }
- printf("\n");
-
- }
-
- return 0;
-}
-
Deleted: trunk/SASxport/src/reverse.c
===================================================================
--- trunk/SASxport/src/reverse.c 2007-08-12 03:12:21 UTC (rev 1136)
+++ trunk/SASxport/src/reverse.c 2007-08-12 03:13:36 UTC (rev 1137)
@@ -1,83 +0,0 @@
-#include "writeSAS.h"
-
-#include <stdio.h>
-#include <string.h>
-#include <stdlib.h>
-#import <assert.h>
-#import <sys/types.h>
-
-/* reverse: convert current byte order to little endian */
-void reverse( u_char *intp, size_t size)
-{
- static u_char endianTest[2] = {0x01,0x00};
- size_t i;
- u_char tmp;
-
-#if !defined(BIG_ENDIAN) && !defined(LITTLE_ENDIAN)
- /* Test if we are on a big endian or little endian platform */
- if( (short) *endianTest == 1 )
- {
- /* Little Endian */
- /* Do nothing */
- return;
- }
-#endif
-
- /* Big Endian */
- /* Swap bytes */
- for(i=0; i < size/2; i++)
- {
- tmp = (u_char) intp[i];
- intp[i] = intp[size-i-1];
- intp[size-i] = tmp;
- }
-
- return;
-}
-
-
-/* test code */
-void test_reverse()
-{
- u_char byte_pattern[1] = { 0x00 };
- u_char byte_value = 0x00;
-
- u_char short_pattern[2] = { 0x00, 0x01 }; /* NB: little endian byte pattern */
- short short_value = 0x0100; /* NB: hex is written big endian */
-
- u_char int_pattern[4] = { 0x0, 0x01, 0x02, 0x03 };
- int int_value = 0x03020100;
-
- u_char long_pattern[4] = { 0x0, 0x01, 0x02, 0x03 };
- long long_value = 0x03020100;
-
- /* Do the reverse, then test */
-
- /* byte */
- REVERSE( &byte_value, sizeof(u_char) );
- assert( (u_char) *byte_pattern == byte_value );
-
- /* short */
- REVERSE( &short_value, sizeof(short) );
- assert( *((short *) short_pattern) == short_value );
-
- /* int */
- REVERSE( &int_value, sizeof(int) );
- assert( *((int *) int_pattern) == int_value );
-
- /* long */
- REVERSE( &long_value, sizeof(long) );
- assert( *((long*) long_pattern) == long_value );
-
-}
-
-
-#ifdef DO_TEST
-int main(int argc, u_char *argv)
-{
- test_reverse();
-}
-
-
-
-#endif
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
