fast_frommer.h

Go to the documentation of this file.

#ifndef FAST_FROMMER_2627D46B
#define FAST_FROMMER_2627D46B

#if _MSC_VER > 1000
#pragma once
#endif // _MSC_VER > 1000


#include "FrommerTimer.h"

//Problem soll ich Frommer durch "generic_Ring" oder durch zZahl parametrisieren?? Eher durch ein Ring,
// Der Ring sollte ein TypeDef für den ElementTyp haben!
//
//Koerper ist momentan durch ein #define austauschbar...


namespace nCenterFocus
{


template < class TFrommerDefs, int variant>
class fast_Frommer //has to implement the  IFrommer Interface
{

        public:
                typedef  typename TFrommerDefs::RingType        TRing;

                typedef  typename TFrommerDefs::RingType        RingType;

                typedef typename TFrommerDefs::PolynomXYType    TPolynomXY;

                typedef typename TFrommerDefs::TMatrixXYType    TMatrixXYType;
                typedef typename TFrommerDefs::TMatrixXY_0Type  TMatrixXY_0Type;

                typedef typename TFrommerDefs::TMatrixXYPairType TMatrixXYPairType;

                typedef typename TFrommerDefs::PolynomXYPairType TPolynomXYPair;


        private:

        std::string name_m;

        #ifdef FROMMERSTAGE_TIMER
                Timer tim_m;
                Timer tim2_m;
                FrommerTimer frommerTimer_m;
                MtcpManager_g.connectTimer(tim_m);
                MtcpManager_g.connectTimer(tim2_m);
        #endif

        #ifdef SAFE
                bool initialized;
        #endif


        const TRing * const     ring1; 

        const typename RingType::FieldType & field_ref_m;

        const short             maxFocalValuesToCompute;

        const short             maxDegree; 




        TPolynomXY const *      q; 

        TPolynomXY const *      minusp; 

         
        TPolynomXYPair          pq;


        TMatrixXYType * const   A;
        
        

        TMatrixXYPairType* const        dA; 


        TMatrixXYType    * const        B;

 
        TMatrixXY_0Type  * const B0;

        const typename TRing::ScalarType        * const         a_table; 
        
        short int       currentDegree;

        short int       lastDegree;

        short           currFocalValuePos;

        

protected:
        
                void wellDefined();
                inline void init();
                static const typename TRing::ScalarType * const         createLowATable(int maxDegree, const TRing * _ring );
                template <class TemplateMatrix>
                inline void     outputMatrix(TemplateMatrix *mat) const;

        
                inline void     outputMatrices() const;
                inline void     compute_dyA_and_dxA();
        
                 void   perform_the_A_Step();
                inline void     perform_the_even_A_Step();
                inline void     perform_the_odd_A_Step();
                
                inline void     perform_generic_C_Step();
                inline void     perform_C_Step();
                void                    perform_fast_C_Step();
                 void           perform_first_C_Step();

                //void   perform_second_C_Step();


                inline void     perform_C_StepDeg3();
                
                inline void     perform_inner_C_Step(   register typename TRing::ElementType * const b0, 
                                                        register const typename TRing::ElementType * const p0,
                                                        register const typename TRing::ElementType * const q0,
                                                        register const typename TRing::ElementType * const da0,
                                                        short _n
                                                );

                inline void     perform_fast_inner_C_Step(      register typename TMatrixXY_0Type::CoefficientType &  b0,
                                                        register const typename TRing::ElementType &  p0,
                                                        register const typename TRing::ElementType &  q0,
                                                        register const typename TRing::ElementType &  da0,
                                                        short _n
                                                );

public:
        void setName(std::string name) {name_m=name;    }
        
        fast_Frommer(short m, const TRing *rr);
        virtual ~fast_Frommer();


                  void  setPolynoms(const TPolynomXY * const minuspp, 
                                    const TPolynomXY * const qq);


                inline void     doit(short bComputeAllFocalValues, int jacobian=false);

                void    doitNew(short bComputeAllFocalValues, int jacobian);
        
                void    doitOld(short bComputeAllFocalValues, int jacobian);

                

        

        template<int >
        inline    void  doitx(short howManyFocalValues);

        template<int >
                  void  doitxNew(short howManyFocalValues);

        
        template<int >
                void    doitxOld(short howManyFocalValues);
        

                inline const TRing*     getRing() const;

                inline const TRing&     getRingRef() const;
                inline typename TRing::ElementType &    getFocalValue( int pos) const;

                void    getComputedFocalValues(vector<typename TRing::ElementType> &computedFV ) const;

                inline short    getSuccessiveVanishedFocalValuesCount() const;

                short   getComputedFocalValuesNum() const {return  currFocalValuePos-1 ; }

                short   getMaxFocalValuesToCompute() const {return  maxFocalValuesToCompute ; }
         typename TRing::ElementType*   computedFocalValues;

        void printStageTimings(std::ostream & os) const;

        #ifdef CF_TEST
                TMatrixXYType getA()    { return *A; }
                TMatrixXYType getB()    { return *B; }
        #endif
};



        #include "fast_frommer.cpp"



} // end Namespace FrommerAlgorithm
#endif // !defined(FAST_FROMMER_2627D46B)