BlockSolver.h

Go to the documentation of this file.

#ifndef __CINT__
#ifndef BlockSolver_H
#define BlockSolver_H

#include <map>
#include <string>
#include <vector>

#include "CLHEP/Matrix/GenMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Random/RandGaussQ.h"

#include "Calibration/Tools/interface/InvMatrixCommonDefs.h"

struct BlockSolver {
  int operator()(const CLHEP::HepMatrix& matrix, const CLHEP::HepVector& vector, CLHEP::HepVector& result);

private:
  void shrink(const CLHEP::HepMatrix& matrix,
              CLHEP::HepMatrix& solution,
              const CLHEP::HepVector& result,
              CLHEP::HepVector& input,
              const std::vector<int>& where);
  void pour(CLHEP::HepVector& result, const CLHEP::HepVector& output, const std::vector<int>& where);
};

#endif
#endif