}
}
/**************************************************************************************************/
+//build and swap functions based on pam.c by maechler from R cluster package
//sets Dp[0] does not set Dp[1]. chooses intial medoids.
int Pam::buildPhase() {
try {
for (int j = 0; j < numSamples; j++) {
if (m->control_pressed) { break; }
if (dists[i][j] == Dp[j][0]) {
- double small = 0.0;
- if (Dp[j][1] > dists[h][j]) { small = dists[h][j]; }
- else { small = Dp[j][1]; }
- dz += (- Dp[j][0]+ small);
+ double smallValue; smallValue = 0.0;
+ if (Dp[j][1] > dists[h][j]) { smallValue = dists[h][j]; }
+ else { smallValue = Dp[j][1]; }
+ dz += (- Dp[j][0]+ smallValue);
}else if (dists[h][j] < Dp[j][0]) {
dz += (- Dp[j][0] + dists[h][j]);
}