version 1.31, 2014/06/16 04:45:18
|
version 1.32, 2015/01/23 19:21:09
|
Line 1973 void computeSWIndex(struct swIndex *swKe |
|
Line 1973 void computeSWIndex(struct swIndex *swKe |
|
float *jz_err_squared = (float *) (malloc(nxny * sizeof(float))); | float *jz_err_squared = (float *) (malloc(nxny * sizeof(float))); |
float *jz_err_squared_smooth = (float *) (malloc(nxny * sizeof(float))); | float *jz_err_squared_smooth = (float *) (malloc(nxny * sizeof(float))); |
float *jz_rms_err = (float *) (malloc(nxny * sizeof(float))); | float *jz_rms_err = (float *) (malloc(nxny * sizeof(float))); |
|
float *err_term1 = (float *) (malloc(nxny * sizeof(float))); |
|
float *err_term2 = (float *) (malloc(nxny * sizeof(float))); |
| |
// define some values for the R calculation | // define some values for the R calculation |
int scale = round(2.0/cdelt1); | int scale = round(2.0/cdelt1); |
Line 2041 void computeSWIndex(struct swIndex *swKe |
|
Line 2043 void computeSWIndex(struct swIndex *swKe |
|
} | } |
| |
computeJz(bx_err, by_err, bx, by, dims, jz, jz_err, jz_err_squared, mask, bitmask, cdelt1, rsun_ref, rsun_obs, | computeJz(bx_err, by_err, bx, by, dims, jz, jz_err, jz_err_squared, mask, bitmask, cdelt1, rsun_ref, rsun_obs, |
derx, dery); |
derx, dery, err_term1, err_term2); |
| |
| |
if(computeJzsmooth(bx, by, dims, jz, jz_smooth, jz_err, jz_rms_err, jz_err_squared_smooth, &(swKeys_ptr->mean_jz), | if(computeJzsmooth(bx, by, dims, jz, jz_smooth, jz_err, jz_rms_err, jz_err_squared_smooth, &(swKeys_ptr->mean_jz), |