Linear.h

Go to the documentation of this file.

#ifndef LINEAR
#define LINEAR

#include "Shape.h"

template <size_t SPD>
class Linear: public Shape {
public:

  Linear (const size_t * iMap = 0);

  inline virtual ~Linear () {}

  Linear (const Linear<SPD> &);

  virtual inline Linear<SPD> * clone () const {
    return new Linear<SPD> (*this);
  }

  // Accessors/Mutators
  inline size_t getNumFunctions () const {
    return SPD + 1;
  }
  inline size_t getNumVariables () const {
    return SPD;
  }

  // Functionality

  double getVal (size_t a, const double *x) const;
  double getDVal (size_t a, const double *x, size_t i) const;

private:
  size_t bctMap[SPD + 1];
};

template<size_t SPD>
Linear<SPD>::Linear (const size_t * iMap) {
  for (size_t a = 0; a < SPD + 1; a++)
    bctMap[a] = a;

  if (iMap != 0) {
    for (size_t a = 0; a < SPD + 1; a++) {
      bctMap[a] = iMap[a];
    }
  }

  return;
}

template<size_t SPD>
Linear<SPD>::Linear (const Linear<SPD> &Lin) {
  for (size_t a = 0; a < SPD + 1; a++) {
    bctMap[a] = Lin.bctMap[a];
  }
}

template<size_t SPD>
double Linear<SPD>::getVal (size_t a, const double *x) const {
  if (bctMap[a] != SPD) {
    return x[bctMap[a]];
  }
  else {
    double va = 0;

    for (size_t k = 0; k < SPD; k++) {
      va += x[k];
    }

    return 1 - va;
  }
}

template<size_t SPD>
double Linear<SPD>::getDVal (size_t a, const double *x, size_t i) const {
  if (bctMap[a] != SPD) {
    return bctMap[a] == i ? 1 : 0;
  } else {
    return -1;
  }
}

#endif

---