(file) Return to update_sharp_keys.c CVS log (file) (dir) Up to [Development] / JSOC / proj / sharp / apps

Diff for /JSOC/proj/sharp/apps/update_sharp_keys.c between version 1.6 and 1.7

version 1.6, 2014/05/16 21:56:33 version 1.7, 2014/06/02 19:47:26
Line 268  int DoIt(void)
Line 268  int DoIt(void)
         int scale = round(2.0/cdelt1);         int scale = round(2.0/cdelt1);
         int nx1 = nx/scale;         int nx1 = nx/scale;
         int ny1 = ny/scale;         int ny1 = ny/scale;
           int nxp = nx1+40;
           int nyp = ny1+40;
  
         if (nx1 > floor((nx-1)/scale + 1) )         if (nx1 > floor((nx-1)/scale + 1) )
                 DIE("X-dimension of output array in fsample() is too large.");                 DIE("X-dimension of output array in fsample() is too large.");
Line 281  int DoIt(void)
Line 283  int DoIt(void)
         float *p1p = (float *)malloc(nx1*ny1*sizeof(float));         float *p1p = (float *)malloc(nx1*ny1*sizeof(float));
         float *p1n = (float *)malloc(nx1*ny1*sizeof(float));         float *p1n = (float *)malloc(nx1*ny1*sizeof(float));
         float *p1 = (float *)malloc(nx1*ny1*sizeof(float));         float *p1 = (float *)malloc(nx1*ny1*sizeof(float));
         float *pmap = (float *)malloc(nx1*ny1*sizeof(float));          float *pmap    = (float *)malloc(nxp*nyp*sizeof(float));
           float *p1pad   = (float *)malloc(nxp*nyp*sizeof(float));
           float *pmapn   = (float *)malloc(nx1*ny1*sizeof(float));
  
         for (irec=0;irec<nrecs;irec++)         for (irec=0;irec<nrecs;irec++)
         {         {
Line 385  int DoIt(void)
Line 389  int DoIt(void)
            if (r_valueflag)            if (r_valueflag)
            {            {
               if (computeR(bz_err, los , dims, &Rparam, cdelt1, rim, p1p0, p1n0,               if (computeR(bz_err, los , dims, &Rparam, cdelt1, rim, p1p0, p1n0,
                      p1p, p1n, p1, pmap, nx1, ny1))                             p1p, p1n, p1, pmap, nx1, ny1, scale, p1pad, nxp, nyp, pmapn))
   
               {               {
                  Rparam            = DRMS_MISSING_FLOAT;                  Rparam            = DRMS_MISSING_FLOAT;
               }               }
Line 697  int DoIt(void)
Line 702  int DoIt(void)
         free(bt_err); free(bh_err);  free(jz_err);         free(bt_err); free(bh_err);  free(jz_err);
         free(jz_err_squared); free(jz_rms_err);         free(jz_err_squared); free(jz_rms_err);
         free(cvsinfoall);         free(cvsinfoall);
   
         free(rim);         free(rim);
         free(p1p0);         free(p1p0);
         free(p1n0);         free(p1n0);
Line 704  int DoIt(void)
Line 710  int DoIt(void)
         free(p1n);         free(p1n);
         free(p1);         free(p1);
         free(pmap);         free(pmap);
           free(p1pad);
           free(pmapn);
  
 return 0; return 0;
 }       // DoIt }       // DoIt

Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
Karen Tian
 Powered by
ViewCVS 0.9.4