pnome-scal.c File Reference

#include <stdio.h>
#include <boolean.h>
#include <math.h>
#include "arithmetique.h"
#include "vecteur.h"
#include "polynome.h"

Include dependency graph for pnome-scal.c:

Go to the source code of this file.

Functions
void	polynome_scalar_mult (Ppolynome *ppp, float factor)
	void polynome_scalar_mult(Ppolynome* ppp, float factor) (*ppp) = factor * (*ppp) !usage: polynome_scalar_mult(&pp, factor); More...
Ppolynome	polynome_scalar_multiply (Ppolynome pp, float factor)
	Ppolynome polynome_scalar_multiply(Ppolynome pp, float factor) pp = factor * (pp) !usage: pp = polynome_scalar_mult(pp, factor);. More...
void	polynome_scalar_add (Ppolynome *ppp, float term)
	void polynome_scalar_add(Ppolynome* ppp, float term) (*ppp) = (*ppp) + term !usage: polynome_scalar_add(&pp, term); More...
Ppolynome	polynome_scalar_addition (Ppolynome pp, float term)
	Ppolynome polynome_scalar_addition(Ppolynome pp, float term) pp = pp + term !usage: pp = polynome_scalar_add(pp, term);. More...
Ppolynome	polynome_power_n (Ppolynome pp, int n)
	Ppolynome polynome_power_n(Ppolynome pp, int n) returns pp ^ n (n>=0) More...
Ppolynome	polynome_nth_root (Ppolynome p, int n)
	computes the n-root of polynomial if possible, that is if all exponents are multiple of n return POLYNOME_UNDEFINED if not possible symbolically More...
Ppolynome	number_replaces_var (Ppolynome pp, Variable var, float num)
	Ppolynome number_replaces_var(Ppolynome pp, Variable var, float num) returns a copy of polynomial pp where variable var is replaced by a floating-point number: num. More...
Ppolynome	polynome_incr (Ppolynome pp)
	Ppolynome polynome_incr(Ppolynome pp) returns pp + 1. More...
Ppolynome	polynome_decr (Ppolynome pp)
	Ppolynome polynome_decr(Ppolynome pp) returns pp - 1. More...

Function Documentation

number_replaces_var()

Ppolynome number_replaces_var	(	Ppolynome	pp,
		Variable	var,
		float	num
	)

Ppolynome number_replaces_var(Ppolynome pp, Variable var, float num) returns a copy of polynomial pp where variable var is replaced by a floating-point number: num.

removes also the monomial pmnum

Definition at line 201 of file pnome-scal.c.

{
    Pmonome pmnum;
    Ppolynome ppnum, ppnew;
 
    if (POLYNOME_UNDEFINED_P(pp)) 
        return (POLYNOME_UNDEFINED);
    else {
        pmnum = make_monome(num, TCST, VALUE_ONE);
        ppnum = monome_to_new_polynome(pmnum);
        ppnew = polynome_var_subst(pp, var, ppnum);
        polynome_rm(&ppnum);            /* removes also the monomial pmnum */
 
        return(ppnew);
    }
}

References make_monome(), monome_to_new_polynome(), num, polynome_rm(), POLYNOME_UNDEFINED, POLYNOME_UNDEFINED_P, polynome_var_subst(), TCST, and VALUE_ONE.

Here is the call graph for this function:

polynome_decr()

Ppolynome polynome_decr

(

Ppolynome

pp

)

Ppolynome polynome_decr(Ppolynome pp) returns pp - 1.

pp is NOT duplicated.

Parameters

pp

p

Definition at line 246 of file pnome-scal.c.

{
    if (POLYNOME_UNDEFINED_P(pp)) 
        return (POLYNOME_UNDEFINED);
    else {
        Pmonome minus_one = make_monome(-1.0, TCST, VALUE_ONE);
 
        polynome_monome_add(&pp, minus_one);
        monome_rm(&minus_one);
 
        return(pp);
    }
}

References make_monome(), monome_rm(), polynome_monome_add(), POLYNOME_UNDEFINED, POLYNOME_UNDEFINED_P, TCST, and VALUE_ONE.

Referenced by add_constraint_on_x(), and make_causal_external().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_incr()

Ppolynome polynome_incr

(

Ppolynome

pp

)

Ppolynome polynome_incr(Ppolynome pp) returns pp + 1.

pp is NOT duplicated.

Parameters

pp

p

Definition at line 226 of file pnome-scal.c.

{
    if (POLYNOME_UNDEFINED_P(pp)) 
        return (POLYNOME_UNDEFINED);
    else {
        Pmonome one = make_monome(1.0, TCST, VALUE_ONE);
 
        polynome_monome_add(&pp, one);
        monome_rm(&one);
 
        return(pp);
    }
}

References make_monome(), monome_rm(), polynome_monome_add(), POLYNOME_UNDEFINED, POLYNOME_UNDEFINED_P, TCST, and VALUE_ONE.

Here is the call graph for this function:

polynome_nth_root()

Ppolynome polynome_nth_root	(	Ppolynome	p,
		int	n
	)

computes the n-root of polynomial if possible, that is if all exponents are multiple of n return POLYNOME_UNDEFINED if not possible symbolically

Definition at line 177 of file pnome-scal.c.

                                                {
    Ppolynome pp = polynome_dup(p);
    for(p=pp;!POLYNOME_NUL_P(p);p=polynome_succ(p)) {
        Pmonome m = polynome_monome(p);
        Pvecteur v ;
        monome_coeff(m)=powf(monome_coeff(m),1.f/n);
        for(v = monome_term(m); !VECTEUR_NUL_P(v); v=vecteur_succ(v)) {
            if(vecteur_val(v)%n == 0) {
                vecteur_val(v)/=n; 
            }
            else if(vecteur_var(v)!=(Variable)TCST){
                polynome_rm(&pp);
                return POLYNOME_UNDEFINED;
            }
        }
    }
    return pp;
}

References monome_coeff, monome_term, polynome_dup(), polynome_monome, POLYNOME_NUL_P, polynome_rm(), polynome_succ, POLYNOME_UNDEFINED, TCST, VECTEUR_NUL_P, vecteur_succ, vecteur_val, and vecteur_var.

Referenced by polynome_roots().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_power_n()

Ppolynome polynome_power_n	(	Ppolynome	pp,
		int	n
	)

Ppolynome polynome_power_n(Ppolynome pp, int n) returns pp ^ n (n>=0)

Modification:

• treat n < 0 if pp is a monomial. LZ 6 Nov. 92

FI: a unique return would be welcome! No enough time for cleaning

Parameters

pp

p

Definition at line 121 of file pnome-scal.c.

{
    if (POLYNOME_UNDEFINED_P(pp)) 
        return (POLYNOME_UNDEFINED);
    else if (POLYNOME_NUL_P(pp)) {
        if(n>0)
            return POLYNOME_NUL;
        else if (n == 0)
            return make_polynome(1.0, TCST, VALUE_ONE);
        else
            return POLYNOME_UNDEFINED;
    }
    else if (n < 0) {
        if ( is_polynome_a_monome(pp) ) {
            int i,m=-n;
            Ppolynome pptemp, ppresult = polynome_dup(pp);
 
            for(i = 1; i < m; i++)      {
                pptemp = polynome_mult(ppresult, pp);
                polynome_rm(&ppresult);
                ppresult = pptemp;
            }
            return(polynome_div(make_polynome(1.0, TCST, VALUE_ONE), 
                                ppresult));
        }
        else
            polynome_error("polynome_power_n",
                           "negative power n=%d"
                           " and polynome is not a monome\n",n);
    }
    else if (n == 0) 
        return(make_polynome(1.0, TCST, VALUE_ONE));
    else if (n == 1) 
        return(polynome_dup(pp));
    else if (n > 1) {
        int i;
        Ppolynome pptemp, ppresult = polynome_dup(pp);
 
        for(i = 1; i < n; i++)  {
            pptemp = polynome_mult(ppresult, pp);
            polynome_rm(&ppresult);
            ppresult = pptemp;
        }
        return(ppresult);
    }
    /* FI: a unique return would be welcome! No enough time for cleaning */
    polynome_error("polynome_power_n", "Cannot happen!\n");
    return POLYNOME_UNDEFINED;
}

References is_polynome_a_monome(), make_polynome(), polynome_div(), polynome_dup(), polynome_error(), polynome_mult(), POLYNOME_NUL, POLYNOME_NUL_P, polynome_rm(), POLYNOME_UNDEFINED, POLYNOME_UNDEFINED_P, TCST, and VALUE_ONE.

Referenced by polynome_sigma(), polynome_sum_of_power(), polynome_var_subst(), and power_op_handler().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_scalar_add()

void polynome_scalar_add	(	Ppolynome *	ppp,
		float	term
	)

void polynome_scalar_add(Ppolynome* ppp, float term) (*ppp) = (*ppp) + term !usage: polynome_scalar_add(&pp, term);

Definition at line 86 of file pnome-scal.c.

{
    if ((term != 0) && (!POLYNOME_UNDEFINED_P(*ppp))) {
        Pmonome pmtoadd = make_monome((float) term, TCST, VALUE_ONE);
        polynome_monome_add(ppp, pmtoadd);
        monome_rm(&pmtoadd);
    }
}

References make_monome(), monome_rm(), polynome_monome_add(), POLYNOME_UNDEFINED_P, TCST, term(), and VALUE_ONE.

Referenced by do_computation_intensity().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_scalar_addition()

Ppolynome polynome_scalar_addition	(	Ppolynome	pp,
		float	term
	)

Ppolynome polynome_scalar_addition(Ppolynome pp, float term) pp = pp + term !usage: pp = polynome_scalar_add(pp, term);.

Definition at line 101 of file pnome-scal.c.

{
    if ((term != 0) && (!POLYNOME_UNDEFINED_P(pp))) {
        Pmonome pmtoadd = make_monome((float) term, TCST, VALUE_ONE);
        pp = polynome_monome_addition(pp, pmtoadd);
        monome_rm(&pmtoadd);
    }
    return pp;
}

References make_monome(), monome_rm(), polynome_monome_addition(), POLYNOME_UNDEFINED_P, TCST, term(), and VALUE_ONE.

Referenced by complexity_float_add().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_scalar_mult()

void polynome_scalar_mult	(	Ppolynome *	ppp,
		float	factor
	)

void polynome_scalar_mult(Ppolynome* ppp, float factor) (*ppp) = factor * (*ppp) !usage: polynome_scalar_mult(&pp, factor);

returns *ppp pointing to POLYNOME_NUL

Definition at line 46 of file pnome-scal.c.

{
    Ppolynome curpp;
 
    if ((!POLYNOME_NUL_P(*ppp)) && (!POLYNOME_UNDEFINED_P(*ppp))) {
        if (factor == 0)
            polynome_rm(ppp);   /* returns *ppp pointing to POLYNOME_NUL */
        else
            for(curpp = *ppp; curpp != POLYNOME_NUL; curpp = polynome_succ(curpp))
                monome_coeff(polynome_monome(curpp)) *= factor;
    }
}

References monome_coeff, polynome_monome, POLYNOME_NUL, POLYNOME_NUL_P, polynome_rm(), polynome_succ, and POLYNOME_UNDEFINED_P.

Referenced by do_computation_intensity(), make_causal_external(), make_causal_internal(), polynome_roots(), and polynome_sum_of_power().

Here is the call graph for this function:

Here is the caller graph for this function:

polynome_scalar_multiply()

Ppolynome polynome_scalar_multiply	(	Ppolynome	pp,
		float	factor
	)

Ppolynome polynome_scalar_multiply(Ppolynome pp, float factor) pp = factor * (pp) !usage: pp = polynome_scalar_mult(pp, factor);.

Definition at line 65 of file pnome-scal.c.

{
    Ppolynome curpp;
 
    if ((!POLYNOME_NUL_P(pp)) && (!POLYNOME_UNDEFINED_P(pp))) {
        if (factor == 0)
            pp = polynome_free(pp);
        else
            for(curpp = pp; curpp != POLYNOME_NUL; curpp = polynome_succ(curpp))
                monome_coeff(polynome_monome(curpp)) *= factor;
    }
    return pp;
}

References monome_coeff, polynome_free(), polynome_monome, POLYNOME_NUL, POLYNOME_NUL_P, polynome_succ, and POLYNOME_UNDEFINED_P.

Referenced by complexity_scalar_mult(), and polynome_roots().

Here is the call graph for this function:

Here is the caller graph for this function: