root/Minuit2/MatrixInverse.h

0001 // @(#)root/minuit2:$Id$
0002 // Authors: M. Winkler, F. James, L. Moneta, A. Zsenei   2003-2005
0003
0004 /**********************************************************************
0005  *                                                                    *
0006  * Copyright (c) 2005 LCG ROOT Math team,  CERN/PH-SFT                *
0007  *                                                                    *
0008  **********************************************************************/
0009
0010 #ifndef ROOT_Minuit2_MatrixInverse
0011 #define ROOT_Minuit2_MatrixInverse
0012
0013 #include "Minuit2/ABTypes.h"
0014 #include "Minuit2/ABObj.h"
0015
0016 namespace ROOT {
0017
0018 namespace Minuit2 {
0019
0020 template <class mtype, class M, class T>
0021 class MatrixInverse {
0022
0023 public:
0024    MatrixInverse(const M &obj) : fObject(obj) {}
0025
0026    ~MatrixInverse() {}
0027
0028    typedef mtype Type;
0029
0030    const M &Obj() const { return fObject; }
0031
0032 private:
0033    M fObject;
0034 };
0035
0036 template <class M, class T>
0037 class MatrixInverse<vec, M, T> {
0038
0039 private:
0040    MatrixInverse(const M &obj) : fObject(obj) {}
0041
0042 public:
0043    ~MatrixInverse() {}
0044
0045    typedef vec Type;
0046
0047    const M &Obj() const { return fObject; }
0048
0049 private:
0050    M fObject;
0051 };
0052
0053 template <class mt, class M, class T>
0054 inline ABObj<mt, MatrixInverse<mt, ABObj<mt, M, T>, T>, T> Inverse(const ABObj<mt, M, T> &obj)
0055 {
0056    return ABObj<mt, MatrixInverse<mt, ABObj<mt, M, T>, T>, T>(MatrixInverse<mt, ABObj<mt, M, T>, T>(obj));
0057 }
0058
0059 } // namespace Minuit2
0060
0061 } // namespace ROOT
0062
0063 #endif // ROOT_Minuit2_MatrixInverse