plvar-ecart.c File Reference

#include <stdio.h>
#include <stdlib.h>
#include "boolean.h"
#include "arithmetique.h"
#include "vecteur.h"
#include "contrainte.h"
#include "sc.h"
#include "sommet.h"
#include "ray_dte.h"
#include "polyedre.h"
#include "matrix.h"
#include "plint.h"

Include dependency graph for plvar-ecart.c:

Go to the source code of this file.

Macros
#define	TRACE
	package plint More...
#define	MALLOC(s, t, f)   malloc(s)
	pour recuperer les declarations des fonctions de conversion de sc en liste de sommets et reciproquement, bien que ca casse le DAG des types de donnees More...

Functions
Psommet	eq_in_ineq (Psommet *sys, int *nb_som, Pvecteur *lvbase)
	Psommet eq_in_ineq(Psommet * sys, int * nb_som, Pvecteur * lvbase): Transformation des egalites du systeme en inegalites. More...
Psommet	var_ecart_sup (Psommet sys, int nb_som, Pvecteur *lvbase, int *nbvars, Pbase *b)
	Psommet var_ecart_sup(Psommet sys, int nb_som, Pvecteur * lvbase, int * nbvars, Pbase *b): ajout des variables d'ecart necessaires a la transformation du systeme A . More...
Psommet	add_var_sup (Psommet sys, int nb_som, Pvecteur *lvbase, Pvecteur *lvsup, int *nbvars, Pbase *b, Psommet fonct)

Macro Definition Documentation

MALLOC

#define MALLOC	(		s,
			t,
			f
	)		malloc(s)

pour recuperer les declarations des fonctions de conversion de sc en liste de sommets et reciproquement, bien que ca casse le DAG des types de donnees

Definition at line 52 of file plvar-ecart.c.

TRACE

#define TRACE

package plint

Definition at line 41 of file plvar-ecart.c.

Function Documentation

add_var_sup()

Psommet add_var_sup	(	Psommet	sys,
		int	nb_som,
		Pvecteur *	lvbase,
		Pvecteur *	lvsup,
		int *	nbvars,
		Pbase *	b,
		Psommet	fonct
	)

??? 999 ???

Definition at line 210 of file plvar-ecart.c.

{
    Psysteme ps1=NULL;
    Psommet ineg;
    Psommet sys2 = sys;
    Variable nv;
 
#ifdef TRACE
    Pvecteur pvec;
 
    printf(" ** ajout des variables suupplementaire \n");
#endif
    if (sys) {
        ps1 = som_sys_conv(sys);
        if (ps1 != NULL) {
            ps1->dimension = *nbvars;
            ps1->base = base_dup(*b);
                
            for (ineg = sys2;ineg!= NULL;ineg=ineg->succ) {
                if (*(ineg->eq_sat) == -1) {
                    nv = creat_new_var(ps1);
                    vect_add_elem(&(ineg->vecteur),nv,VALUE_ONE);
                    *b = vect_add_variable(*b,nv);
                    (*nbvars) ++;
                    if (value_pos_p(vect_coeff(TCST,ineg->vecteur))) {
                        nv = creat_new_var(ps1);
                        vect_add_elem(&(ineg->vecteur),nv,VALUE_MONE);
                        vect_add_elem(&(fonct->vecteur),nv,
                                      int_to_value(999)); /* ??? 999 ??? */
                        lvbase_add(nv,nb_som,lvsup);
                    }
                    lvbase_add(nv,nb_som,lvbase);
 
                    *(ineg->eq_sat) = 1;
                    *b = vect_add_variable(*b,nv);
                    (*nbvars) ++;
                }
                nb_som --;
            }
 
#ifdef TRACE
            printf (" -- variables de base :");
            for (pvec = *lvbase;pvec!= NULL;pvec=pvec->succ)
            {
                printf(" (0x%x,",(unsigned int)pvec->var);
                print_Value(pvec->val);
                printf("), ");
            }
            printf("\n");
            ps1 = som_sys_conv(sys2);
            sc_fprint(stdout,ps1,*variable_default_name);
#endif
 
            sc_rm(ps1);
        }
    }
    return(sys2);
}

References Ssysteme::base, base_dup(), creat_new_var(), Ssysteme::dimension, typ_som::eq_sat, int_to_value, lvbase_add(), print_Value(), printf(), sc_fprint(), sc_rm(), typ_som::succ, Svecteur::succ, TCST, Svecteur::val, VALUE_MONE, VALUE_ONE, value_pos_p, Svecteur::var, variable_default_name(), vect_add_elem(), vect_add_variable(), vect_coeff(), and typ_som::vecteur.

Referenced by primal(), and primal_positive().

Here is the call graph for this function:

Here is the caller graph for this function:

eq_in_ineq()

Psommet eq_in_ineq	(	Psommet *	sys,
		int *	nb_som,
		Pvecteur *	lvbase
	)

Psommet eq_in_ineq(Psommet * sys, int * nb_som, Pvecteur * lvbase): Transformation des egalites du systeme en inegalites.

N egalites du systeme seront transformees en N+1 inegalites.

resultat retourne par la fonction :

Psommet : systeme lineaire initial modifie.

Les parametres de la fonction :

Psommet sys : systeme lineaire Pvecteur lvbase: liste des variables de base du systeme int nb_som : nombre de contraintes du systeme

Definition at line 68 of file plvar-ecart.c.

{
    Psommet eq=NULL;
    Psommet sys2;
    Psommet som2;
    Pvecteur pv3=NULL;
    Pvecteur pvsom=NULL;
    Pvecteur pvsom2 = NULL;
    bool egalite = false;
    Value den = VALUE_ONE;
 
#ifdef TRACE
 
    Psysteme ps1=NULL;
    Pvecteur pvec;
 
    printf(" ** transformation des equations en inequations \n");
#endif
 
 
    sys2 = *sys;
 
    for (eq = sys2;eq!= NULL;eq=eq->succ) {
        if (*(eq->eq_sat) == 1) {
            egalite = true;
            *(eq->eq_sat) = -1;
            pv3 = vect_dup(eq->vecteur);
            vect_chg_sgn(pv3);
            (void) vect_multiply(pv3,den);
            value_product(den,eq->denominateur);
            (void) vect_multiply(pvsom,eq->denominateur);
            pvsom2 = vect_add(pvsom,pv3);
            vect_rm(pv3);
            vect_rm(pvsom);
            pvsom = pvsom2;
        }
    }
    if (egalite) {
        som2 = (Psommet)MALLOC(sizeof(Ssommet),SOMMET,"eq_in_ineq");
        som2->eq_sat =(int *)MALLOC(sizeof(int),INTEGER,"eq_in_ineq");
        *(som2->eq_sat) = -1;
        som2->denominateur =den ;
        som2->vecteur = pvsom;
        som2->succ = NULL;
        sommet_add(&sys2,som2,nb_som);
    }
 
#ifdef TRACE
    ps1 = som_sys_conv(sys2);
 
    sc_fprint(stdout,ps1,*variable_default_name);
    printf (" -- variables de base :");
    for (pvec = *lvbase;pvec!= NULL;pvec=pvec->succ)
    {
        printf(" (0x%x,", (unsigned int) pvec->var);
        print_Value(pvec->val);
        printf("), ");
    }
    printf("\n");
    sc_rm(ps1);
 
#endif
 
    sommets_normalize(sys2);
    return(sys2);
 
}

References typ_som::denominateur, eq, Scontrainte::eq_sat, typ_som::eq_sat, MALLOC, print_Value(), printf(), sc_fprint(), sc_rm(), SOMMET, sommet_add(), sommets_normalize(), Scontrainte::succ, typ_som::succ, Svecteur::succ, Svecteur::val, VALUE_ONE, value_product, Svecteur::var, variable_default_name(), vect_add(), vect_chg_sgn(), vect_dup(), vect_multiply(), vect_rm(), Scontrainte::vecteur, and typ_som::vecteur.

Referenced by dual_pivot(), plint(), plreal(), primal(), primal_positive(), and sc_supress_same_constraints().

Here is the call graph for this function:

Here is the caller graph for this function:

var_ecart_sup()

Psommet var_ecart_sup	(	Psommet	sys,
		int	nb_som,
		Pvecteur *	lvbase,
		int *	nbvars,
		Pbase *	b
	)

Psommet var_ecart_sup(Psommet sys, int nb_som, Pvecteur * lvbase, int * nbvars, Pbase *b): ajout des variables d'ecart necessaires a la transformation du systeme A .

X <= B en A'. X' = B .

resultat retourne par la fonction :

Psommet : systeme lineaire initial modifie

Les parametres de la fonction :

Psommet sys : systeme lineaire Pvecteur lvbase: liste des variables de base du systeme int nb_som : nombre de contraintes du systeme int nbvars : nombre de variables du systeme Pbase b : liste des variables du systeme

Definition at line 156 of file plvar-ecart.c.

{
    Psysteme ps1=NULL;
    Psommet ineg;
    Psommet sys2 = sys;
    Variable nv;
 
#ifdef TRACE
    Pvecteur pvec;
 
    printf(" ** ajout des variables d'ecart \n");
#endif
    if (sys) {
        ps1 = som_sys_conv(sys);
        if (ps1 != NULL) {
            ps1->dimension = *nbvars;
            ps1->base = base_dup(*b);
                
            for (ineg = sys2;ineg!= NULL;ineg=ineg->succ) {
                if (*(ineg->eq_sat) == -1) {
                    nv = creat_new_var(ps1);
                    vect_add_elem(&(ineg->vecteur),nv,VALUE_ONE);
                    lvbase_add(nv,nb_som,lvbase);
                    *(ineg->eq_sat) = 1;
                    *b = vect_add_variable(*b,nv);
                    (*nbvars)++; 
                }
                nb_som --;
            }
 
#ifdef TRACE
            printf (" -- variables de base :");
            for (pvec = *lvbase;pvec!= NULL;pvec=pvec->succ)
            {
                printf(" (0x%x,",(unsigned int)pvec->var);
                print_Value(pvec->val);
                printf("), ");
            }
            printf("\n");
            ps1 = som_sys_conv(sys2);
            sc_fprint(stdout,ps1,*variable_default_name);
#endif
            sc_rm(ps1);
        }
    }
    return(sys2);
}

References Ssysteme::base, base_dup(), creat_new_var(), Ssysteme::dimension, typ_som::eq_sat, lvbase_add(), print_Value(), printf(), sc_fprint(), sc_rm(), typ_som::succ, Svecteur::succ, Svecteur::val, VALUE_ONE, Svecteur::var, variable_default_name(), vect_add_elem(), vect_add_variable(), and typ_som::vecteur.

Referenced by dual_pivot(), dual_pivot_pas(), plint(), plint_degen(), plint_pas(), plreal(), and syst_smith().

Here is the call graph for this function:

Here is the caller graph for this function: