Changeset 12090
- Timestamp:
- 09/06/07 12:59:41 (10 months ago)
- Files:
-
- spike/grib/degrib18/g2clib-1.0.4/compack.c (modified) (2 diffs)
- spike/grib/degrib18/g2clib-1.0.4/grib2.h (modified) (1 diff)
- spike/grib/degrib18/g2clib-1.0.4/jpcpack.c (modified) (4 diffs)
- spike/grib/degrib18/g2clib-1.0.4/misspack.c (modified) (2 diffs)
- spike/grib/degrib18/g2clib-1.0.4/mkieee.c (modified) (1 diff)
- spike/grib/degrib18/g2clib-1.0.4/pngpack.c (modified) (4 diffs)
- spike/grib/degrib18/g2clib-1.0.4/simpack.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
spike/grib/degrib18/g2clib-1.0.4/compack.c
r10748 r12090 109 109 // 110 110 if (idrstmpl[1] == 0) { // No binary scaling 111 imin=(g2int) rint(rmin*dscale);111 imin=(g2int)RINT(rmin*dscale); 112 112 //imax=(g2int)rint(rmax*dscale); 113 113 rmin=(g2float)imin; 114 114 for (j=0;j<ndpts;j++) 115 ifld[j]=(g2int) rint(fld[j]*dscale)-imin;115 ifld[j]=(g2int)RINT(fld[j]*dscale)-imin; 116 116 } 117 117 else { // Use binary scaling factor … … 119 119 //rmax=rmax*dscale; 120 120 for (j=0;j<ndpts;j++) 121 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);121 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 122 122 } 123 123 // spike/grib/degrib18/g2clib-1.0.4/grib2.h
r10754 r12090 202 202 typedef struct gribfield gribfield; 203 203 204 #define rint(d) (floor(d+0.5))204 #define RINT(d) (floor(d+0.5)) 205 205 206 206 /* Prototypes for unpacking API */ spike/grib/degrib18/g2clib-1.0.4/jpcpack.c
r12077 r12090 93 93 } 94 94 if (idrstmpl[1] == 0) 95 maxdif = (g2int) ( rint(rmax*dscale) - rint(rmin*dscale));95 maxdif = (g2int) (RINT(rmax*dscale) - RINT(rmin*dscale)); 96 96 else 97 maxdif = (g2int) rint( (rmax-rmin)*dscale*bscale );97 maxdif = (g2int)RINT( (rmax-rmin)*dscale*bscale ); 98 98 // 99 99 // If max and min values are not equal, pack up field. … … 113 113 // bits in which the data will fit. 114 114 // 115 imin=(g2int) rint(rmin*dscale);116 imax=(g2int) rint(rmax*dscale);115 imin=(g2int)RINT(rmin*dscale); 116 imax=(g2int)RINT(rmax*dscale); 117 117 maxdif=imax-imin; 118 118 temp=log((double)(maxdif+1))/alog2; … … 121 121 // scale data 122 122 for(j=0;j<ndpts;j++) 123 ifld[j]=(g2int) rint(fld[j]*dscale)-imin;123 ifld[j]=(g2int)RINT(fld[j]*dscale)-imin; 124 124 } 125 125 else { … … 130 130 rmin=rmin*dscale; 131 131 rmax=rmax*dscale; 132 maxdif=(g2int) rint((rmax-rmin)*bscale);132 maxdif=(g2int)RINT((rmax-rmin)*bscale); 133 133 temp=log((double)(maxdif+1))/alog2; 134 134 nbits=(g2int)ceil(temp); 135 135 // scale data 136 136 for (j=0;j<ndpts;j++) 137 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);137 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 138 138 } 139 139 // spike/grib/degrib18/g2clib-1.0.4/misspack.c
r10748 r12090 145 145 nonmiss=0; 146 146 if (idrstmpl[1] == 0) { // No binary scaling 147 imin=(g2int) rint(rmin*dscale);147 imin=(g2int)RINT(rmin*dscale); 148 148 //imax=(g2int)rint(rmax*dscale); 149 149 rmin=(g2float)imin; 150 150 for ( j=0; j<ndpts; j++) { 151 151 if (ifldmiss[j] == 0) { 152 jfld[nonmiss]=(g2int) rint(fld[j]*dscale)-imin;152 jfld[nonmiss]=(g2int)RINT(fld[j]*dscale)-imin; 153 153 nonmiss++; 154 154 } … … 160 160 for ( j=0; j<ndpts; j++ ) { 161 161 if (ifldmiss[j] == 0) { 162 jfld[nonmiss]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);162 jfld[nonmiss]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 163 163 nonmiss++; 164 164 } spike/grib/degrib18/g2clib-1.0.4/mkieee.c
r10748 r12090 104 104 else 105 105 atemp=atemp*two126; 106 imant=(g2int) rint(atemp*two23);106 imant=(g2int)RINT(atemp*two23); 107 107 } 108 108 else { spike/grib/degrib18/g2clib-1.0.4/pngpack.c
r12075 r12090 76 76 if (fld[j] < rmin) rmin=fld[j]; 77 77 } 78 maxdif = (g2int) rint( (rmax-rmin)*dscale*bscale );78 maxdif = (g2int)RINT( (rmax-rmin)*dscale*bscale ); 79 79 // 80 80 // If max and min values are not equal, pack up field. … … 94 94 // bits in which the data will fit. 95 95 // 96 imin=(g2int) rint(rmin*dscale);97 imax=(g2int) rint(rmax*dscale);96 imin=(g2int)RINT(rmin*dscale); 97 imax=(g2int)RINT(rmax*dscale); 98 98 maxdif=imax-imin; 99 99 temp=log((double)(maxdif+1))/alog2; … … 102 102 // scale data 103 103 for(j=0;j<ndpts;j++) 104 ifld[j]=(g2int) rint(fld[j]*dscale)-imin;104 ifld[j]=(g2int)RINT(fld[j]*dscale)-imin; 105 105 } 106 106 else { … … 111 111 rmin=rmin*dscale; 112 112 rmax=rmax*dscale; 113 maxdif=(g2int) rint((rmax-rmin)*bscale);113 maxdif=(g2int)RINT((rmax-rmin)*bscale); 114 114 temp=log((double)(maxdif+1))/alog2; 115 115 nbits=(g2int)ceil(temp); 116 116 // scale data 117 117 for (j=0;j<ndpts;j++) 118 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);118 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 119 119 } 120 120 // spike/grib/degrib18/g2clib-1.0.4/simpack.c
r10748 r12090 100 100 // bits in which the data will fit. 101 101 // 102 imin=(g2int) rint(rmin*dscale);103 imax=(g2int) rint(rmax*dscale);102 imin=(g2int)RINT(rmin*dscale); 103 imax=(g2int)RINT(rmax*dscale); 104 104 maxdif=imax-imin; 105 105 temp=log((double)(maxdif+1))/alog2; … … 108 108 // scale data 109 109 for(j=0;j<ndpts;j++) 110 ifld[j]=(g2int) rint(fld[j]*dscale)-imin;110 ifld[j]=(g2int)RINT(fld[j]*dscale)-imin; 111 111 } 112 112 else if (nbits!=0 && idrstmpl[1]==0) { … … 123 123 // scale data 124 124 for (j=0;j<ndpts;j++) 125 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);125 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 126 126 } 127 127 else if (nbits==0 && idrstmpl[1]!=0) { … … 132 132 rmin=rmin*dscale; 133 133 rmax=rmax*dscale; 134 maxdif=(g2int) rint((rmax-rmin)*bscale);134 maxdif=(g2int)RINT((rmax-rmin)*bscale); 135 135 temp=log((double)(maxdif+1))/alog2; 136 136 nbits=(g2int)ceil(temp); 137 137 // scale data 138 138 for (j=0;j<ndpts;j++) 139 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);139 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 140 140 } 141 141 else if (nbits!=0 && idrstmpl[1]!=0) { … … 149 149 // scale data 150 150 for (j=0;j<ndpts;j++) 151 ifld[j]=(g2int) rint(((fld[j]*dscale)-rmin)*bscale);151 ifld[j]=(g2int)RINT(((fld[j]*dscale)-rmin)*bscale); 152 152 } 153 153 //
