1 |
|
#ifndef lint |
2 |
< |
static const char RCSid[] = "$Id$"; |
2 |
> |
static const char RCSid[] = "$Id$"; |
3 |
|
#endif |
4 |
|
/* |
5 |
|
* pf3.c - routines for gaussian and box filtering |
9 |
|
|
10 |
|
#include "standard.h" |
11 |
|
|
12 |
+ |
#include <string.h> |
13 |
+ |
|
14 |
|
#include "color.h" |
15 |
|
|
16 |
|
#define RSCA 1.13 /* square-radius multiplier: sqrt(4/PI) */ |
68 |
|
double d; |
69 |
|
register int x; |
70 |
|
|
71 |
< |
gtabsiz = 111*CHECKRAD*CHECKRAD/(rad*rad); |
71 |
> |
gtabsiz = 111*CHECKRAD*CHECKRAD; |
72 |
|
gausstable = (float *)malloc(gtabsiz*sizeof(float)); |
73 |
|
if (gausstable == NULL) |
74 |
|
goto memerr; |
192 |
|
register int c, x; |
193 |
|
register float *gscan; |
194 |
|
/* compute ring sums */ |
195 |
< |
bzero((char *)ringsum, (orad+1)*sizeof(float)); |
196 |
< |
bzero((char *)ringwt, (orad+1)*sizeof(short)); |
195 |
> |
memset((char *)ringsum, '\0', (orad+1)*sizeof(float)); |
196 |
> |
memset((char *)ringwt, '\0', (orad+1)*sizeof(short)); |
197 |
|
for (r = -orad; r <= orad; r++) { |
198 |
|
if (rcent+r < 0) continue; |
199 |
|
if (rcent+r >= nrows) break; |