Tree/transit.cpp at master · chelseah/Tree · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
#include "transit.h"
Item::Item(){
//  std::string s1;
//  line_= new std::string();
  line_=NULL;
  p_=NULL;
  dsp_=NULL;
  harm_flag_=1;
}
Item::Item(std::string * line, double p,double dsp,int harm_flag)
  : p_(p),dsp_(dsp),harm_flag_(harm_flag){
  //swap(*line_,*line);
  line_ = new string(*line);
//  line_->replace(line_->begin(),line_->end(),*line);
}

Item::~Item(){
  delete line_;
}
void Item::copy(const Item *nb){
  p_=nb->p_;
  if(!line_){
    delete line_;
  }
  line_ = new string(*(nb->line_));
  //  line_= nb->line_;
  dsp_=nb->dsp_;
  harm_flag_= nb->harm_flag_;
}

int Item::Ge(const Item *nb){
  if(harm_flag_){
    if(fabs(p_-nb->p_)/p_<1e-3){
      return 0;
    } else {
      if(p_<nb->p_){
        return 1;
      } else {
        return -1;
      }
    }
  } else {
    if(Isharm(nb->p_)){
      return 0;
    } else {
       if(p_<nb->p_){
        return 1;
      } else {
        return -1;
      }

    }
  }
}

bool Item::Isharm(const double p){
  int order;
  double err;
  if(p_>p){
    order = floor(p_/p+0.5);
    err = fabs((p_-order*p)/p_);
  } else {
    order = floor(p/p_+0.5);
    err = fabs((p/order-p_)/p_);
  }
  //cout << order << " " << err << endl;
  if (err<1.e-3) {
    return true;
  } else {
    return false;
  }
}

int Item::Better(const Item *nb) {
  if(dsp_<nb->dsp_){
    return 1;
  } else {
      return 0;
  }
}
void Item::PrintState(vector <string> &outstr){
  if(outstr.size()==0){
    std::cout << *line_ << '\n';
  } else {
//    std::cout << outstr.size() << endl;
    outstr.push_back(*line_);
  //  std::cout << "inprint"<<outstr.size() << endl;
  }
  //  std::cout << *line_ << '\n';
}