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