-rrange for byAdd

	modified:   src/Objects/DataManip/mrcImage/src/lmrcImageSymmetryFind.c

diff --git a/src/Objects/DataManip/mrcImage/src/lmrcImageSymmetryFind.c b/src/Objects/DataManip/mrcImage/src/lmrcImageSymmetryFind.c
index 843f80a..10a96f5 100755 (executable)
--- a/src/Objects/DataManip/mrcImage/src/lmrcImageSymmetryFind.c
+++ b/src/Objects/DataManip/mrcImage/src/lmrcImageSymmetryFind.c
@@ -356,6 +356,7 @@ lmrcImageSymmetryFindCheckByTTest(mrcImage* in, lmrcImageSymmetryFindInfo* info,
        mrcImage* pow1doutBacks;
        mrcImage  prob;
        mrcImage  ttest;
+       double ymin, ymax;
 
        lmrcImageSymmetryFindSubtraction(&subout, in, info, mode);
        subouts = (mrcImage*)memoryAllocate(sizeof(mrcImage)*info->flagImages, "in lmrcImageSymmetryFindTTest");
@@ -410,9 +411,17 @@ lmrcImageSymmetryFindCheckByTTest(mrcImage* in, lmrcImageSymmetryFindInfo* info,
                        pow1douts[j].HeaderN.z = 1;
                        pow1douts[j].HeaderMode = mrcFloatImage;
                        mrcInit(&(pow1douts[j]), NULL);
+                       
+                       if(info->rmin<=info->rmax) {
+                               ymin = info->rmin;
+                               ymax = info->rmax;
+                       } else {
+                               ymin = 0;
+                               ymax = fftouts[j].HeaderN.y-1;
+                       }
                        for(x=0; x<pow1douts[j].HeaderN.x; x++) { 
                                sum = 0.0;
-                               for(y=0; y<fftouts[j].HeaderN.y; y++) {
+                               for(y=ymin; y<=ymax; y++) {
                                        mrcPixelDataGet(&(fftouts[j]), x, y, 0.0, &data1, mrcPixelPow, mrcPixelHowNearest);  
                                        sum += data1*y;
                                }