Exemplo n.º 1
0
  @Override
  public void trigger(double time) {
    Measure m = new Measure();
    m.time = time;
    int cnt = 0;
    double dsCnt = 0;
    double usCnt = 0;
    double dsRho = 0;
    double usRho = 0;
    for (int lane = 0; lane < this.lanes; lane++) {
      CAMoveableEntity[] p = this.l.getParticles(lane);
      for (int i = this.from; i <= this.to; i++) {
        if (p[i] != null) {
          cnt++;
          m.in.add(p[i].getId());
          AgentInfo ai = this.ais.get(p[i].getId());
          // updateTimeRho(i, ai.dir, ai, time, p);
          if (p[i].getDir() == 1) {

            dsCnt++;
            // dsRho += ai.timeRho.get(time);
          } else {
            usCnt++;
            // usRho += ai.timeRho.get(time);
          }
        }
      }
    }
    double rho = cnt / this.a;

    m.rho = rho;
    m.dsRho = dsCnt / this.a;
    m.usRho = usCnt / this.a;
    // m.dsRho = dsRho / dsCnt;
    // m.usRho = usRho / usCnt;

    for (int lane = 0; lane < this.lanes; lane++) {
      CAMoveableEntity[] p = this.l.getParticles(lane);
      for (int i = 0; i < this.l.getNumOfCells(); i++) {
        CAMoveableEntity part = p[i];
        if (part != null) {
          Id id = part.getId();
          AgentInfo ai = this.ais.get(id);
          if (ai == null && (i == 0 || i == this.l.getNumOfCells() - 1)) {
            ai = new AgentInfo(part.getDir());
            // ai.lastPos = i;
            this.ais.put(id, ai);
          }
          if (ai.lastPos != i) {
            ai.timePos.put(time, i);
            ai.lastPos = i;
            if (i == 0) {
              if (part.getDir() == 1) {
                this.dsIn++;
              } else {
                this.usOut++;
              }
            } else if (i == this.l.getNumOfCells() - 1) {
              if (part.getDir() == 1) {
                this.dsOut++;
              } else {
                this.usIn++;
              }
            }
          }
        }
      }
    }
    if (m.rho > 0) {
      this.ms.add(m);
    }
  }
Exemplo n.º 2
0
  @Override
  public void report(BufferedWriter bw) throws IOException {
    int in = this.dsIn + this.usIn - this.dsOut - this.usOut;
    System.out.println(
        in
            + " "
            + " dsIn:"
            + this.dsIn
            + " dsOut:"
            + this.dsOut
            + " usIn:"
            + this.usIn
            + " usOut:"
            + this.usOut);

    List<Tuple<Integer, Integer>> ranges = new ArrayList<>();
    Variance vd = new Variance();
    Variance vs = new Variance();
    Variance vf = new Variance();
    int cnt = 0;
    double tm = 15;
    int from = 0;
    int to = 0;
    for (Measure m : this.ms) {
      cpdSpd(m);
      bw.append(
          m.time + " " + m.dsRho + " " + m.dsSpd + " " + m.usRho + " " + m.usSpd + " " + m.rho + " "
              + " " + m.spd + "\n");

      if (true) continue;
      if (m.time < 15 || m.time > 100) {
        to++;
        from++;
        continue;
      }
      cpdSpd(m);
      vd.addVar(m.rho);
      vs.addVar(m.spd);
      vf.addVar(m.rho * m.spd);
      cnt++;
      if (cnt < 5) {
        to++;
        continue;
      }
      double md = vd.getMean();
      double std = Math.sqrt(vd.getVar());

      double ms = vs.getMean();
      double sts = Math.sqrt(vs.getVar());

      double mf = vf.getMean();
      double stf = Math.sqrt(vf.getVar());

      double bal = m.dsRho / m.usRho;
      bal = bal > 1 ? 1 / bal : bal;
      // bal = 1;
      double timespan = m.time - tm;
      if (md > 1 && std > 0.1 * md || sts > 0.1 * ms || bal < 0.5 || timespan > 20) {
        if (timespan > 5) {
          Tuple<Integer, Integer> t = new Tuple<>(from, to);
          ranges.add(t);
        }
        tm = m.time;
        from = to;
        cnt = 0;
        vd = new Variance();
        vs = new Variance();
      }
      to++;
    }

    int rng = ranges.size();
    int cntt = 0;
    for (Tuple<Integer, Integer> r : ranges) {
      cntt++;
      // if (cntt < rng / 2 + rng / 4 || cntt > rng / 2 + rng / 4 + rng /
      // 8) {
      // continue;
      // }
      from = r.getFirst();
      to = r.getSecond();
      Measure m = new Measure();

      double range = to - from;
      double ccnt = 0;
      double lastTime = 0;
      for (int idx = from; idx < to; idx++) {
        Measure mm = this.ms.get(idx);
        if (lastTime + 1 > mm.time) {
          continue;
        }
        ccnt++;
        m.rho += mm.rho;
        m.spd += mm.spd;
        m.time += mm.time;
        m.dsRho += mm.dsRho;
        m.usRho += mm.usRho;
        m.dsSpd += mm.dsSpd;
        m.usSpd += mm.usSpd;
        bw.append(
            mm.time + " " + mm.dsRho + " " + mm.dsSpd + " " + mm.usRho + " " + mm.usSpd + " "
                + mm.rho + " " + " " + mm.spd + "\n");
        lastTime = mm.time;
      }
      m.rho /= ccnt;
      m.spd /= ccnt;
      m.dsRho /= ccnt;
      m.dsSpd /= ccnt;
      m.usRho /= ccnt;
      m.usSpd /= ccnt;
      m.time /= ccnt;
      // bw.append(m.time + " " + m.dsRho + " " + m.dsSpd + " " + m.usRho
      // + " " + m.usSpd + " " + m.rho + " " + " " + m.spd + "\n");

    }
    // bw.append(m.time + " " + m.rho + " " + m.spd + "\n");

  }