helpers.h

Code for L2-optimal three view triangulation based on calculation of stationary points

#ifndef HELPERS
#define HELPERS

inline void multiply(double *eq, double *A, int *multiplicationMap, int mapSize,
	      int destRow, int m, int n) {
  for(int j = 0; j < mapSize; j++) {
    A[m * (multiplicationMap[j] - 1) + destRow] = // subtract by 1 to comply with matlab indices. 
      eq[n - mapSize + j];
  }
};

inline void divide(double *eq, double *A, int *multiplicationMap, int mapSize,
	    int destRow, int m, int n) {
  for(int j = 0; j < mapSize; j++) {
    A[m * (n - mapSize + j) + destRow] = eq[multiplicationMap[j] - 1]; // subtract by 1 to comply with matlab indices. 
  }
};

#endif