regress_linear.h

Go to the documentation of this file.

 
#ifndef LIBA_REGRESS_LINEAR_H
#define LIBA_REGRESS_LINEAR_H
 
#include "a.h"
 
typedef struct a_regress_linear a_regress_linear;
 
#if defined(__cplusplus)
extern "C" {
#endif /* __cplusplus */
 
A_EXTERN void a_regress_linear_init(a_regress_linear *ctx, a_float *coef_p, a_size coef_n, a_float bias);
 
A_EXTERN a_float a_regress_linear_eval(a_regress_linear const *ctx, a_float const *val);
 
A_EXTERN void a_regress_linear_err1(a_regress_linear const *ctx, a_size n, a_float const *x, a_float const *y, a_float *err);
 
A_EXTERN void a_regress_linear_err2(a_regress_linear const *ctx, a_size n, a_float const *const *x, a_float const *y, a_float *err);
 
A_EXTERN void a_regress_linear_pdm1(a_regress_linear const *ctx, a_size n, a_float const *x, a_float *pdm, a_float y_mean);
 
A_EXTERN void a_regress_linear_pdm2(a_regress_linear const *ctx, a_size n, a_float const *const *x, a_float *pdm, a_float y_mean);
 
A_EXTERN void a_regress_linear_gd(a_regress_linear *ctx, a_float const *input, a_float error, a_float alpha);
 
A_EXTERN void a_regress_linear_sgd(a_regress_linear *ctx, a_float const *x, a_float y, a_float alpha);
 
A_EXTERN void a_regress_linear_sgd1(a_regress_linear *ctx, a_size n, a_float const *x, a_float const *y, a_float alpha);
 
A_EXTERN void a_regress_linear_sgd2(a_regress_linear *ctx, a_size n, a_float const *const *x, a_float const *y, a_float alpha);
 
A_EXTERN void a_regress_linear_bgd1(a_regress_linear *ctx, a_size n, a_float const *x, a_float const *err, a_float alpha);
 
A_EXTERN void a_regress_linear_bgd2(a_regress_linear *ctx, a_size n, a_float const *const *x, a_float const *err, a_float alpha);
 
A_EXTERN a_float a_regress_linear_mgd1(a_regress_linear *ctx, a_size n, a_float const *x, a_float const *y, a_float *err,
                                       a_float delta, a_float lrmax, a_float lrmin, a_size lrtim, a_size epoch, a_size batch);
 
A_EXTERN a_float a_regress_linear_mgd2(a_regress_linear *ctx, a_size n, a_float const *const *x, a_float const *y, a_float *err,
                                       a_float delta, a_float lrmax, a_float lrmin, a_size lrtim, a_size epoch, a_size batch);
 
A_EXTERN void a_regress_linear_zero(a_regress_linear *ctx);
 
#if defined(__cplusplus)
} /* extern "C" */
namespace a
{
typedef struct a_regress_linear regress_linear;
} /* namespace a */
#endif /* __cplusplus */
 
struct a_regress_linear
{
    a_float *coef_p; // points to regression coefficients
    a_size coef_n; // number of regression coefficients
    a_float bias; // intercept
#if defined(__cplusplus)
    A_INLINE void init(a_float *p, a_size n, a_float b = 0)
    {
        a_regress_linear_init(this, p, n, b);
    }
    A_INLINE a_float eval(a_float const *val) const
    {
        return a_regress_linear_eval(this, val);
    }
    A_INLINE void err1(a_size n, a_float const *x, a_float const *y, a_float *err) const
    {
        a_regress_linear_err1(this, n, x, y, err);
    }
    A_INLINE void err2(a_size n, a_float const *const *x, a_float const *y, a_float *err) const
    {
        a_regress_linear_err2(this, n, x, y, err);
    }
    A_INLINE void pdm1(a_size n, a_float const *x, a_float *pdm, a_float y_mean) const
    {
        a_regress_linear_pdm1(this, n, x, pdm, y_mean);
    }
    A_INLINE void pdm2(a_size n, a_float const *const *x, a_float *pdm, a_float y_mean) const
    {
        a_regress_linear_pdm2(this, n, x, pdm, y_mean);
    }
    A_INLINE void gd(a_float const *input, a_float error, a_float alpha)
    {
        a_regress_linear_gd(this, input, error, alpha);
    }
    A_INLINE void sgd(a_float const *x, a_float y, a_float alpha)
    {
        a_regress_linear_sgd(this, x, y, alpha);
    }
    A_INLINE void sgd1(a_size n, a_float const *x, a_float const *y, a_float alpha)
    {
        a_regress_linear_sgd1(this, n, x, y, alpha);
    }
    A_INLINE void sgd2(a_size n, a_float const *const *x, a_float const *y, a_float alpha)
    {
        a_regress_linear_sgd2(this, n, x, y, alpha);
    }
    A_INLINE void bgd1(a_size n, a_float const *x, a_float const *err, a_float alpha)
    {
        a_regress_linear_bgd1(this, n, x, err, alpha);
    }
    A_INLINE void bgd2(a_size n, a_float const *const *x, a_float const *err, a_float alpha)
    {
        a_regress_linear_bgd2(this, n, x, err, alpha);
    }
    A_INLINE a_float mgd1(a_size n, a_float const *x, a_float const *y, a_float *err, a_float delta,
                          a_float lrmax, a_float lrmin, a_size lrtim, a_size epoch, a_size batch)
    {
        return a_regress_linear_mgd1(this, n, x, y, err, delta, lrmax, lrmin, lrtim, epoch, batch);
    }
    A_INLINE a_float mgd2(a_size n, a_float const *const *x, a_float const *y, a_float *err, a_float delta,
                          a_float lrmax, a_float lrmin, a_size lrtim, a_size epoch, a_size batch)
    {
        return a_regress_linear_mgd2(this, n, x, y, err, delta, lrmax, lrmin, lrtim, epoch, batch);
    }
    A_INLINE void zero() { a_regress_linear_zero(this); }
#endif /* __cplusplus */
};
 
#endif /* a/regress_linear.h */