jacobian_analytic.h

Go to the documentation of this file.

/*
    This software is Copyright by the Board of Trustees of Michigan
    State University (c) Copyright 2017.
 
    You may use this software under the terms of the GNU public license
    (GPL).  The terms of this license are described at:
 
     http://www.gnu.org/licenses/gpl.txt
 
     Authors:
             Ron Fox
             Girodano Cerizza
         Aaron Chester
         FRIB
         Michigan State University
         East Lansing, MI 48824-1321
*/
 
#ifndef JACOBIAN_ANALYTIC_H
#define JACOBIAN_ANALYTIC_H
 
#include "CFitEngine.h"
 
namespace ddastoys {
    namespace analyticfit {
 
    // Concrete classes
 
    class SerialFitEngine1 : public CFitEngine {
    public:
        SerialFitEngine1(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        ~SerialFitEngine1() {}
        virtual void jacobian(const gsl_vector* p,  gsl_matrix *J);
        virtual void residuals(const gsl_vector*p, gsl_vector* r);
    };
 
    class CudaFitEngine1 : public CFitEngine {
    private:
        void* m_dXtrace;    
        void* m_dYtrace;    
        void* m_dResiduals; 
        void* m_dJacobian;  
        unsigned m_npts;    
    
    public:
        CudaFitEngine1(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        ~CudaFitEngine1();
        virtual void jacobian(const gsl_vector* p, gsl_matrix *J);
        virtual void residuals(const gsl_vector*p, gsl_vector* r);
    
    private:
        void throwCudaError(const char* msg);
    };
 
    class SerialFitEngine2 : public CFitEngine {
    public:
        SerialFitEngine2(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        virtual void jacobian(const gsl_vector* p,  gsl_matrix *J);
        virtual void residuals(const gsl_vector*p, gsl_vector* r);
    };
 
    class CudaFitEngine2 : public CFitEngine {
    private:
        void* m_dXtrace;    
        void* m_dYtrace;    
        void* m_dResiduals; 
        void* m_dJacobian;  
        unsigned m_npts;    
    
    public:
        CudaFitEngine2(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        ~CudaFitEngine2();
        virtual void jacobian(const gsl_vector* p,  gsl_matrix *J);
        virtual void residuals(const gsl_vector*p, gsl_vector* r);
    
    private:
        void throwCudaError(const char* msg);
    };
 
    class FitEngineFactory {
    public:
        CFitEngine* createSerialFitEngine1(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        CFitEngine* createCudaFitEngine1(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        CFitEngine* createSerialFitEngine2(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
        CFitEngine* createCudaFitEngine2(
        std::vector<std::pair<uint16_t, uint16_t>>&  data
        );
    };
 
    
    }
};
 
#endif