Method Process

avida-core/source/actions/PrintActions.cc:1130–1169

  
  static const cString GetDescription() { return "Arguments: [string fname='lineage_counts.dat']\n  WARNING: This will only have the appropriate header if all lineages are present before this action is run for the first time."; }
  void Process(cAvidaContext&)
  {
    const int update = m_world->GetStats().GetUpdate();
    const double generation = m_world->GetStats().SumGeneration().Average();
    
    //only loop through living organisms
    const Apto::Array<cOrganism*, Apto::Smart>& living_orgs = m_world->GetPopulation().GetLiveOrgList();
    
    Apto::Map<int, int> lineage_label_counts;
    
    //build hash of lineage_label -> count
    for(int i = 0; i < living_orgs.GetSize(); i++) {
      const int cur_lineage_label = living_orgs[i]->GetLineageLabel();
      if (lineage_label_counts.Has(cur_lineage_label)) lineage_label_counts[cur_lineage_label]++;
      else lineage_label_counts[cur_lineage_label] = 1;
    }
    
    //setup lineage labels in the first pass
    if (first_run == false) {
      for (Apto::Map<int, int>::KeyIterator it = lineage_label_counts.Keys(); it.Next();) lineage_labels.Push(*it.Get());
      first_run = true;
    }
    
    Avida::Output::FilePtr df = Avida::Output::File::StaticWithPath(m_world->GetNewWorld(), (const char*)m_filename);
    df->Write(update, "Update");
    df->Write(generation, "Generation");
    
    //for each lineage label, output the counts
    //@LZ - handle dead lineages appropriately
    for (int i = 0; i < lineage_labels.GetSize(); i++) {
      //default to 0 in case this lineage is dead
      int count = 0;
      
      lineage_label_counts.Get(lineage_labels[i], count);
      
      df->Write(count, cStringUtil::Stringf("Lineage Label %d", lineage_labels[i]));
    }
    
    df->Endl();
  }
};