Changeset 12006

Timestamp:

05/19/08 23:22:11 (7 months ago)

Author:

ayu

Message:

• aded comment. etc.
  

Files:

• lang/cplusplus/misc/clustercpp/src/ward.cpp (modified) (3 diffs)

lang/cplusplus/misc/clustercpp/src/ward.cpp

@@ -126 +126 @@
       double mindist=-1.0;
       INDEX minj =0;
+      std::set<INDEX>::iterator kietaend = kieta.end();
+            
       for(INDEX j=0;j<nextnodenum; j++){
-        if(kieta.find(j) != kieta.end()){
+        if(kieta.find(j) != kietaend){
           continue; //すでに消えてるなら何もしない
         }
@@ -137 +139 @@
           (*mat)[j]->erase(p);  //今後使わなくなる成分の消去
         }else{
-          tmpdist = (*objp)[j] =((sp+sizes[j])*elm(j,p)+(sq+sizes[j])*elm(j,q)-sizes[j]*mqp)/(sqp+sizes[j]);
+          tmpdist = (*objp)[j] =((sp+sizes[j])*elm(j,p)+(sq+sizes[j])*elm(j,q)-sizes[j]*mqp)/(sqp+sizes[j]);      
           (*mat)[j]->erase(p); //今後使わなくなる成分の消去
           (*mat)[j]->erase(q); //今後使わなくなる成分の消去
@@ -145 +147 @@
           minj=j;
         }
-        if(hasrow(j)){
+        if(hasrow(j)){ //Q. この判定必要？ A. 必要。まだ合体させられてないのに、matの行自体は消えていることがあるよ。
           if((*mat)[j]->size()==0){
             caches->erase(j);
             mat->erase(j);
+            //一見ここでkieta.insert(j)したくなるが、そうしてはいけない！
+            //kietaにinsertするのは、真にその点が消えた時=他の点と合体させられた時のみである。
+            //そうでないのにここでinsertしてしまうと、まだ消滅していないのに、未来永劫計算対象からはずされてしまうのだっ！！
           }else{
-            Cache& c=*(getcache(j));
-            if(c.index==p || c.index==q){
+            CacheP c = getcache(j);
+            if(c->index==p || c->index==q){
               MATRIXROW::iterator it = std::min_element((*mat)[j]->begin(), (*mat)[j]->end(), compmatrow());
               savecache(j, it->second, it->first);