MultinomialLogisticRegression.hpp

#ifndef _LDAPLUSPLUS_OPTIMIZATION_MULTINOMIAL_LOGISTIC_REGRESSION
#define _LDAPLUSPLUS_OPTIMIZATION_MULTINOMIAL_LOGISTIC_REGRESSION

#include <cmath>

#include <Eigen/Core>

namespace ldaplusplus {
namespace optimization {


template <typename Scalar>
class MultinomialLogisticRegression
{
    typedef Eigen::Matrix<Scalar, Eigen::Dynamic, Eigen::Dynamic> MatrixX;
    typedef Eigen::Matrix<Scalar, Eigen::Dynamic, 1> VectorX;

    public:
        MultinomialLogisticRegression(const MatrixX &X, const Eigen::VectorXi &y, VectorX Cy, Scalar L);
        MultinomialLogisticRegression(const MatrixX &X, const Eigen::VectorXi &y, Scalar L);

        Scalar value(const MatrixX &eta) const;

        void gradient(const MatrixX &eta, Eigen::Ref<MatrixX> grad) const;

    private:
        const MatrixX &X_;
        const Eigen::VectorXi &y_;
        Scalar L_;
        VectorX Cy_;
};


}  // namespace optimization
}  // namespace ldaplusplus
#endif // _LDAPLUSPLUS_OPTIMIZATION_MULTINOMIAL_LOGISTIC_REGRESSION