matrixutil.hpp

Go to the documentation of this file.

#ifndef MATRIXUTIL_HPP
#define MATRIXUTIL_HPP


template <typename T>
void reserve_2D( T ** &arr, int xsize, int ysize) {
  // zusammenhängendes 2D-Array reservieren
  arr=new T* [ysize];
  arr[0]=new T [xsize*ysize];
  for (int y=1; y<ysize; y++) 
    arr[y]=arr[0]+y*xsize;
}

template <typename T>
void delete_2D( T ** &arr) {
  delete[] arr[0];
  delete[] arr;
}  

bool QR (int zeilen, int spalten, double **A, double *d, double *b, const double epsilon);
void rdsolve (int spalten, double** R, double *d, double *c, double *x);
void Invert_RtR(double **R, double **I, int n);
void compute_errors(int num_param, int num_data, double **R, double *d, double *b, double *errors);

inline double signum(double a) { 
 if (a>0) return 1.0;
 else return -1.0; 
}



#endif

---