Function maximum

src/isomodel/SimModel.cpp:198–206

    return vd;
  }
  double maximum(const Vector& v1) {
    double max = -std::numeric_limits<double>::max();
    for (size_t i = 0; i < v1.size(); i++) {
      if (v1[i] > max) {
        max = v1[i];
      }
    }
    return max;
  }
  Vector maximum(const Vector& v1, const Vector& v2) {
    assert(v1.size() == v2.size());
    Vector vx = Vector(v1.size());