In DSCSOLV changed if (ALFA0[I-1] >= 0.0) to if (ALFA0[I-1} > 0.0).
authorilb@NIH.GOV <ilb@NIH.GOV@ba61647d-9d00-f842-95cd-605cb4296b96>
Thu, 30 Nov 2017 21:18:17 +0000 (21:18 +0000)
committerilb@NIH.GOV <ilb@NIH.GOV@ba61647d-9d00-f842-95cd-605cb4296b96>
Thu, 30 Nov 2017 21:18:17 +0000 (21:18 +0000)
git-svn-id: https://citdcbmipav.cit.nih.gov/repos-pub/mipav/trunk@15287 ba61647d-9d00-f842-95cd-605cb4296b96

mipav/src/gov/nih/mipav/model/algorithms/DoublyConnectedSC.java

index 0492ff1..934ad45 100644 (file)
@@ -1388,7 +1388,7 @@ double neweps;
                if (ISHAPE != 0) {\r
                        NSHAPE = 1;\r
                        for (I = 2; I <= M-1; I++) {\r
-                           if (ALFA0[I-1] >= 0.0) {\r
+                           if (ALFA0[I-1] > 0.0) {\r
                                continue;\r
                            }\r
                            NSHAPE = NSHAPE + 1;\r
@@ -3864,14 +3864,14 @@ double neweps;
                          ratio = X3MAX/XABS;\r
                          S3 = ONE + S3* ratio * ratio;\r
                          X3MAX = XABS;\r
-                 }\r
-                 else {\r
-                     if (XABS != ZERO) {\r
-                         ratio = XABS/X3MAX;\r
-                         S3 = S3 + ratio * ratio;\r
-                     }\r
-                 }\r
-                 continue;\r
+                         }\r
+                         else {\r
+                             if (XABS != ZERO) {\r
+                                 ratio = XABS/X3MAX;\r
+                                 S3 = S3 + ratio * ratio;\r
+                             }\r
+                         }\r
+                     continue;\r
 \r
                  } // if (XABS <= RDWARF || XABS >= AGIANT)\r
        \r
@@ -4357,7 +4357,7 @@ double neweps;
                } // if (AJNORM != 0)\r
                RDIAG[J-1] = -AJNORM;\r
            } // for (J = 1; J <= MINMN; J++)\r
-             return;\r
+           return;\r
        }\r
 \r
 \r