#include <eq_system_solver.h>
Definition at line 34 of file eq_system_solver.h.
◆ solve()
Method computes absolute shifts from relative shifts
- Parameters
-
bX | relative shifts in X dim |
bY | relative shifts in Y dim |
A | system matrix to be used |
shiftX | absolute shifts in X dim |
shiftY | absolute shifts in Y dim |
verbosity | level |
iterations | of the solver |
Definition at line 36 of file eq_system_solver.cpp.
43 std::vector<Matrix1D<double>> shifts(2);
67 std::cout <<
"Solving equation system ...\n";
76 ex = bX - helper.
A * shiftX;
77 ey = bY - helper.
A * shiftY;
84 double R2x = 1 - (stddeveX * stddeveX) / varbX;
85 double R2y = 1 - (stddeveY * stddeveY) / varbY;
87 std::cout <<
"Iteration " << it <<
" R2x=" << R2x <<
" R2y=" << R2y
91 double oldWeightSum = helper.
w.
sum();
97 double newWeightSum = helper.
w.
sum();
98 if ((newWeightSum == oldWeightSum) && (verbosity > 1)){
99 std::cout <<
"No outlier found\n\n";
101 }
else if (verbosity > 1)
102 std::cout <<
"Found " << (int) (oldWeightSum - newWeightSum)
106 }
while (it < iterations);
std::vector< Matrix1D< double > > bs
#define FOR_ALL_ELEMENTS_IN_MATRIX1D(v)
double sum(bool average=false) const
void typeCast(const Matrix1D< T1 > &v1, Matrix1D< T2 > &v2)
void resizeNoCopy(int Xdim)
void computeMeanAndStddev(double &mean, double &stddev) const
void weightedLeastSquares(WeightedLeastSquaresHelper &h, Matrix1D< double > &result)
The documentation for this class was generated from the following files: