mul_mod_helper< C, D, m > Struct Template Reference

#include <modular_int.hpp>

List of all members.

Static Public Member Functions

static void op (C &dest, const C &src, const C &p, const C &q, nat r, nat s, nat t)
static void op (C &dest, const C &src, const C &p, C &carry, const C &q, nat r, nat s, nat t)

Static Public Attributes

static const nat n = 8 * sizeof(C)

Detailed Description

template<typename V>
template<typename C, typename D, nat m>
struct mmx::modulus_mul_int_preinverse< V >::mul_mod_helper< C, D, m >

Definition at line 700 of file modular_int.hpp.

Member Function Documentation

static void op	(	C &	dest,
		const C &	src,
		const C &	p,
		C &	carry,
		const C &	q,
		nat	r,
		nat	s,
		nat	t
	)			`[inline, static]`

Definition at line 740 of file modular_int.hpp.

References mmx::C, and mul_mod_helper< C, D, m >::n.

00742 00743 00744 00745 00746 00747 00748 00749 00750 00751       } 00752 00753 00754 00755 00756 00757 00758 00759 00760         } 00761 00762 00763 00764 00765       } 00766 00767 00768 00769 00770 00771 00772         } 00773       } 00774 00775 00776     }

class="fragment">00741 { D a (dest); a *= src; if (m + 2 <= n) { C hi_a = a >> s; C b = (((D) hi_a) * q) >> t, c = 0;; dest = ((C) a) - b * p; if (dest >= p) { dest -= p; b++; } V::add_mod_core (dest, carry, p, c); carry = b + c; return; D hi_a = a >> s; C b = (hi_a * q) >> t, z= 0; D c = a - ((D) b) * p; if (m + 1 <= n) { if (c >= p) { dest = c - p; b++; if (dest >= p) { dest -= p; b++; } else dest = c; V::add_mod_core (dest, carry, p, z); carry = b + z; return; if (c >= p) { c -= p; b++; if (c >= p) { c -= p; b++; if (c >= p) { c -= p; b++; } dest = c; V::add_mod_core (dest, carry, p, z); carry = b + z;

static void op	(	C &	dest,
		const C &	src,
		const C &	p,
		const C &	q,
		nat	r,
		nat	s,
		nat	t
	)			`[inline, static]`

Definition at line 704 of file modular_int.hpp.

References mmx::C, and mul_mod_helper< C, D, m >::n.

Referenced by modulus_mul_int_preinverse< modulus_add_int_naive< modulus_reduction_int_naive< modulus_normalization_int_naive< modulus_maximum_size_int< m > > > > >::mul_mod(), and mul_mod_helper< C, void, m >::op().

00705                                          {
00706       D a (dest);
00707       a *= src;
00708       if (m + 2 <= n) {
00709         C hi_a = a >> s;
00710         C b = (((D) hi_a) * q) >> t;
00711         dest = ((C) a) - b * p;
00712         if (dest >= p) dest -=p;
00713         return;
00714       }
00715       D hi_a = a >> s;
00716       C b = (hi_a * q) >> t;
00717       D c = a - ((D) b) * p;
00718       if (m + 1 <= n) {
00719         if (c >= p) {
00720           dest = c - p;
00721           if (dest >= p)
00722             dest -= p;
00723         }
00724         else
00725           dest = c;
00726         return;
00727       }
00728       if (c >= p) {
00729         c -= p;
00730         if (c >= p) {
00731           c -= p;
00732           if (c >= p) 
00733             c -= p;
00734         }
00735       }
00736       dest = c;
00737     }

Member Data Documentation

const nat n = 8 * sizeof(C) [static]

Definition at line 701 of file modular_int.hpp.

Referenced by mul_mod_helper< C, void, m >::op(), and mul_mod_helper< C, D, m >::op().

The documentation for this struct was generated from the following file:

include/numerix/modular_int.hpp