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