Poly.h File Reference

#include "MultInd.h"
#include "Vector.h"
#include "Matrix.h"
#include <fstream>

Include dependency graph for Poly.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	Polynomial
struct	Polynomial::PolynomialDataTag

Functions
unsigned long	choose (unsigned n, unsigned k)
Polynomial	operator* (const double &dou, Polynomial &p)

Function Documentation

choose()

unsigned long choose	(	unsigned	n,
		unsigned	k
	)

Definition at line 38 of file tools.cpp.

{
    const unsigned long uupSize = 100;
    static unsigned long uup[uupSize];
    unsigned long *up;
    static unsigned long Nold = 0;
    static unsigned long Kold = 0;
    unsigned long l,m;
    unsigned i,j;

    if ( (n < k) || !n ) return 0;

    if ( (n == k) || !k )   // includes n == 1
        return 1;

    if ( k > (n >> 1) ) // Only lower half
        k = n-k;

    if ( (Nold == n) && (k < Kold) )    // We did it last time ...
        return *(uup + k - 1);

    if ( k > uupSize )
    {
        printf( "choose( unsigned, unsigned) :  overflow\n");
        getchar(); exit(-1);
    }

    Nold=n; Kold=k;

    *(up=uup)=2;
    for (i=2; i<n; i++)     // Pascal's triangle
    {
        // todo: remove next line:
        *(up+1)=1;
        l=1;
        m=*(up=uup);
        for (j=0; j<mmin(i,k); j++)
        {
            *up=m+l;
            l=m;
            m=*(++up);
        }
        // todo: remove next line:
        *up=1;
    }

    return *(uup + k - 1);
}

operator*()

Polynomial operator* ( const double &  dou, Polynomial &  p  )

inline

Definition at line 132 of file Poly.h.

{
    // we can use operator * defined on Polynomial because of commutativity
    return p * dou;     
}