blm.h

Go to the documentation of this file.

 
#ifndef M4RIE_BLM_H
#define M4RIE_BLM_H
 
/******************************************************************************
*
*            M4RIE: Linear Algebra over GF(2^e)
*
*    Copyright (C) 2013 Martin Albrecht <martinralbrecht@googlemail.com>
*
*  Distributed under the terms of the GNU General Public License (GEL)
*  version 2 or higher.
*
*    This code is distributed in the hope that it will be useful,
*    but WITHOUT ANY WARRANTY; without even the implied warranty of
*    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
*    General Public License for more details.
*
*  The full text of the GPL is available at:
*
*                  http://www.gnu.org/licenses/
******************************************************************************/
 
#include <m4ri/m4ri.h>
#include "m4rie/gf2e.h"

 
#define M4RIE_CRT_LEN (M4RIE_MAX_DEGREE + 1)

 
typedef struct {
  mzd_t *H; 
  djb_t *h; 
 
  mzd_t *F; 
  djb_t *f; 
 
  mzd_t *G; 
  djb_t *g; 
} blm_t;

 
extern const int costs[17];

 
static inline int blm_cost_crt(const int p[M4RIE_CRT_LEN]) {
  int cost = costs[p[0]];
  for(deg_t d=1; d<M4RIE_CRT_LEN; d++)
    cost += costs[d] * p[d];
  return cost;
}

 
int *crt_init(const deg_t f_len, const deg_t g_len);

 
blm_t *blm_init_crt(const gf2e *ff, const deg_t f_ncols, const deg_t g_ncols, const int *p, int djb);

 
blm_t *_blm_finish_polymult(const gf2e *ff, blm_t *f);

 
void blm_free(blm_t *f);

 
blm_t *_blm_djb_compile(blm_t *f);

 
void _mzd_ptr_apply_blm_mzd(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f);

 
void _mzd_ptr_apply_blm_djb(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f);

 
static inline void _mzd_ptr_apply_blm(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f) {
  if (f->f!=NULL)
    _mzd_ptr_apply_blm_djb(X, A, B, f);
  else
    _mzd_ptr_apply_blm_mzd(X, A, B, f);
}
 
 
#endif //M4RIE_BLM_H