Changeset 82


Ignore:
Timestamp:
May 26, 2004, 4:16:03 PM (19 years ago)
Author:
cecilia
Message:

annot-map

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/Merge.cc

    r81 r82  
    1717    : data_(data), error_(error), nof_assays_(nof_assays)
    1818  {
    19     // skapa annot-vektor (string?)
    20     while(annotations){
     19    // skapa annot-vektor
     20    while(!annotations.eof()){
    2121      string str;
    2222      getline(annotations,str,'\n');
     
    4040    gslapi::matrix weights;
    4141    weights = calculate_weights();
     42
    4243    // för varje annot: kolla vilka rader den finns på
    4344    map<string, vector<u_int> > name;
    44     cout << annots_.size() << " size" << endl;
    45     for (u_int i=0; i<annots_.size(); i++) {
    46       map<string,vector<u_int> >::iterator p;
    47       p = name.find(annots_[i]);
    48       vector<u_int> index;
    49       if(p != name.end()){
    50         index = (*p).second;
    51       }
    52       index.push_back(i);
    53       for (vector<u_int>::const_iterator j=index.begin(); j!=index.end(); j++){
    54         cout << *j << "* ";
    55       }
    56       cout << endl;
    57       name.insert(pair<string, vector<u_int> >(annots_[i], index));
    58     }
     45    for (u_int i=0; i<annots_.size(); i++)
     46      name[annots_[i]].push_back(i);
    5947   
    6048    // för varje namn:
    6149   
    62     for (map<string,vector<u_int> >::const_iterator k=name.begin(); k!=name.end(); k++) {
    63 //      merging((*k).second, weights);
    64       vector<u_int> t = (*k).second;
    65       for (vector<u_int>::const_iterator i=t.begin(); i!=t.end(); i++){
    66         cout << *i << " ";
    67       }
    68       cout << endl;
    69     }
     50    for (map<string,vector<u_int> >::const_iterator k=name.begin(); k!=name.end(); k++)
     51      merging((*k).second, weights);
    7052  }
    7153
Note: See TracChangeset for help on using the changeset viewer.