sco_common.hpp

#pragma once
#include <vector>
#include <cmath>
#include <algorithm>

namespace sco {

using std::vector;
typedef vector<double> DblVec;
typedef vector<unsigned char> BoolVec;

inline double vecSum(const DblVec& v) {
  double out = 0;
  for (int i=0; i < v.size(); ++i) out += v[i];
  return out;
}
inline double vecAbsSum(const DblVec& v) {
  double out = 0;
  for (int i=0; i < v.size(); ++i) out += fabs(v[i]);
  return out;
}
inline double pospart(double x) {
  return (x > 0) ? x : 0;
}
inline double sq(double x) {
  return x*x;
}
inline double vecHingeSum(const DblVec& v) {
  double out = 0;
  for (int i=0; i < v.size(); ++i) out += pospart(v[i]);
  return out;
}
inline double vecMax(const DblVec& v) {
  return *std::max_element(v.begin(), v.end());
}
inline double vecDot(const DblVec& a, const DblVec& b) {
  assert(a.size() == b.size());
  double out=0;
  for (int i=0; i < a.size(); ++i) out += a[i]*b[i];
  return out;
}



}