version 1.27, 2014/03/05 19:51:19
|
version 1.28, 2014/03/13 18:40:43
|
Line 1057 int computeR(float *bz_err, float *los, |
|
Line 1057 int computeR(float *bz_err, float *los, |
|
// set up convolution kernel | // set up convolution kernel |
init_fresize_gaussian(&fresgauss,sigma,20,1); | init_fresize_gaussian(&fresgauss,sigma,20,1); |
| |
// make sure convolution kernel is smaller than or equal to array size |
|
if ( (nx < 41.) || (ny < 41.) ) return -1; |
|
|
|
fsample(los, rim, nx, ny, nx, nx1, ny1, nx1, scale, 0, 0, 0.0); | fsample(los, rim, nx, ny, nx, nx1, ny1, nx1, scale, 0, 0, 0.0); |
for (i = 0; i < nx1; i++) | for (i = 0; i < nx1; i++) |
{ | { |