subst_mul_series_rep< C, V, UV > Class Template Reference

#include <series_implicit.hpp>

List of all members.

Public Member Functions

subst_mul_series_rep (const series< unknown< C, V > > &f, const series< unknown< C, V > > &g)
unknown< C, V > next ()

Public Attributes

nat f_sh
nat g_sh
series< C, V > f_kn
series< C, V > g_kn
series< C, V > inner

Detailed Description

template<typename C, typename V, typename UV>
class mmx::subst_mul_series_rep< C, V, UV >

Definition at line 374 of file series_implicit.hpp.

Constructor & Destructor Documentation

subst_mul_series_rep	(	const series< unknown< C, V > > &	f,
		const series< unknown< C, V > > &	g
	)			`[inline]`

Definition at line 380 of file series_implicit.hpp.

00380                                                                   :
00381     implementation<series_abstractions,UV>
00382      ::template binary_series_rep<mul_op,UC,UV> (f, g),
00383     f_sh (1), g_sh (1), f_kn (known (f)), g_kn (known (g)),
00384     inner (rshiftz (f_kn, (int) f_sh) * rshiftz (g_kn, (int) g_sh)) {}
  UC next () {

Member Function Documentation

unknown<C,V> next ( ) [inline]

Definition at line 385 of file series_implicit.hpp.

References subst_mul_series_rep< C, V, UV >::f_kn, subst_mul_series_rep< C, V, UV >::f_sh, subst_mul_series_rep< C, V, UV >::g_kn, subst_mul_series_rep< C, V, UV >::g_sh, subst_mul_series_rep< C, V, UV >::inner, mmx::is_known(), mmx::rshiftz(), mmx::substitute(), and UC.

00386 00387 00388 00389 00390     } 00391 00392 00393 00394 00395     } 00396 00397 00398 00399 00400 00401     } 00402 00403 00404 00405 00406 00407 00408 00409     } 00410   }

class="fragment">00385 { if (this->n >= g_sh && !is_known (substitute (this->f [this->n - g_sh]))) { g_sh = g_sh << 1; inner= rshiftz (f_kn, (int) f_sh) * rshiftz (g_kn, (int) g_sh); if (this->n >= f_sh && !is_known (substitute (this->g [this->n - f_sh]))) { f_sh = f_sh << 1; inner= rshiftz (f_kn, (int) f_sh) * rshiftz (g_kn, (int) g_sh); if (this->n < f_sh + g_sh) { UC acc= this->f[0] * this->g[this->n]; for (nat i=1; i<=this->n; i++) acc += this->f[i] * this->g[this->n-i]; return acc; else { UC acc= f_kn[0] * this->g[this->n] + this->f[this->n] * g_kn[0]; for (nat i=1; i<f_sh; i++) acc += f_kn[i] * this->g[this->n-i]; for (nat i=1; i<g_sh; i++) acc += this->f[this->n-i] * g_kn[i]; return acc + inner[this->n - f_sh - g_sh];