Performance improvement: use @mp_sqr@ and @mpmont_reduce@ instead of

[u/mdw/catacomb] / mpscan.h

/* -*-c-*-
 *
 * $Id: mpscan.h,v 1.2 1999/11/13 01:55:10 mdw Exp $
 *
 * Sequential bit scan of multiprecision integers
 *
 * (c) 1999 Straylight/Edgeware
 */

/*----- Licensing notice --------------------------------------------------* 
 *
 * This file is part of Catacomb.
 *
 * Catacomb is free software; you can redistribute it and/or modify
 * it under the terms of the GNU Library General Public License as
 * published by the Free Software Foundation; either version 2 of the
 * License, or (at your option) any later version.
 * 
 * Catacomb 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 Library General Public License for more details.
 * 
 * You should have received a copy of the GNU Library General Public
 * License along with Catacomb; if not, write to the Free
 * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 * MA 02111-1307, USA.
 */

/*----- Revision history --------------------------------------------------* 
 *
 * $Log: mpscan.h,v $
 * Revision 1.2  1999/11/13 01:55:10  mdw
 * Fixed so that they compile.  Minor interface changes.
 *
 * Revision 1.1  1999/09/03 08:41:12  mdw
 * Initial import.
 *
 */

#ifndef MPSCAN_H
#define MPSCAN_H

#ifdef __cplusplus
  extern "C" {
#endif

/*----- Header files ------------------------------------------------------*/

#ifndef MPW_H
#  include "mpw.h"
#endif

/*----- Data structures ---------------------------------------------------*/

typedef struct mpscan {
  const mpw *v, *vl;			/* Vector of words to scan */
  mpw w;				/* Current word to scan */
  int bits;				/* Number of bits left in @w@ */
} mpscan;

/*----- Functions provided ------------------------------------------------*/

/* --- @mpscan_initx@ --- *
 *
 * Arguments:	@mpscan *m@ = pointer to bitscanner structure
 *		@const mpw *v, *vl@ = vector of words to scan
 *
 * Returns:	---
 *
 * Use:		Initializes a bitscanner from a low-level base-and-limit
 *		representation of an integer.  Initially no bit is ready; you
 *		must call @mpscan_step@ before anything useful will come
 *		out.
 */

#define MPSCAN_INITX(m_, v_, vl_) do {					\
  mpscan *_m = (m_);							\
  _m->v = (v_);								\
  _m->vl = (vl_);							\
  _m->bits = 0;								\
} while (0)

extern void mpscan_initx(mpscan */*m*/, const mpw */*v*/, const mpw */*vl*/);

/* --- @mpscan_step@ --- *
 *
 * Arguments:	@mpscan *m@ = pointer to bitscanner
 *
 * Returns:	Nonzero if there is another bit to read.
 *
 * Use:		Steps on to the next bit in the integer.  The macro version
 *		evaluates its argument multiple times.
 */

#define MPSCAN_STEP(m)							\
  ((m)->bits ? ((m)->w >>= 1,						\
		(m)->bits--, 1) :					\
   (m)->v < (m)->vl ? ((m)->w = *(m)->v++,				\
		       (m)->bits = MPW_BITS - 1, 1) :			\
   0)

extern int mpscan_step(mpscan */*m*/);

/* --- @mpscan_bit@ --- *
 *
 * Arguments:	@const mpscan *m@ = pointer to bitscanner
 *
 * Returns:	The value of the current bit.
 *
 * Use:		Reads the value of the current bit looked at by a
 *		bitscanner.
 */

#define MPSCAN_BIT(m) ((m)->w & 1)

extern int mpscan_bit(const mpscan */*m*/);

/*----- That's all, folks -------------------------------------------------*/

#ifdef __cplusplus
  }
#endif

#endif