Discussion :

[yanndublanche]

Bonjour,

Je possede une classe SimEvent construite comme suit:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
class SimEvent{
           public:
	// Enumeration of simulation event type
	enum Type {...};
 
	// Constructor & destructor
	SimEvent(){}
	SimEvent(const SimEvent& s){operator = (s);}
	~SimEvent(){}
 
	// Operator
	void operator = (const SimEvent& s){}
	bool operator < (const SimEvent& b) const{}
 
	// Set functions
	void setTau(double a_tau){}
	void setEndTau(double a_endTau){}
	void setDT(double a_dt){}
	void setType(Type a_type){}
	void setId(wstring a_id){}
	void setIdx(unsigned int a_idx){}
 
	void addEntity(Entity* a_entity){}
 
	void setAbsAmount(bool a_absAmount){}
	void setDimensionality(unsigned int a_dimensionality){}
	void setEntityName(wstring a_entity){}
	void setCompartmentName(wstring a_compartmentName){}
	void setTimeseriesID(unsigned int a_timeseriesID){}
	void setTimeseriesEnd(bool a_timeseriesEnd){}
	void setConstraint(Expression a_expr){}
	void setIsDynamicConstraint(bool a_isDynamicConstraint){}
 
	// Get functions
	wstring getId(){}
	unsigned int getDimensionality(){}
 
         //private:
	// Execution time
	double m_tau;
	double m_endTau;
	double m_dt;
 
	// Simulation event type
	Type m_type;
 
	wstring m_id;
 
	// Vector of entity index
	vector<Entity*> m_entityVector;
 
	// Constraint element
	Expression m_constraint;
	bool m_isDynamicConstraint;
 
	// Request for absolute amount 
	bool m_absAmount;
 
	// Dimension of the compartment
	unsigned int m_dimensionality;
 
	// Parameter / constraint / expression idx
	unsigned int m_idx;
 
	// EntityName
	wstring m_entity;
 
	// Compartment
	wstring m_compartmentName;
 
	// Time serie element
	unsigned int m_timeseriesID;
	bool m_timeseriesEnd;		
};

Ensuite dans une autre classe, je cree un tableau de 5000 elements de SimEvent:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

SimEvent m_heap[5000];

Lorsque je compile le programme, il ne semble pas y avoir de probleme.
Par contre, lorsque je l'execute, j'ai une fenetre d'erreur windows qui apparait (j'utilise CodeWarrior pour programmer et Windows XP).

Si je fais le debug pas a pas de Code Warrior, il me trouve une erreur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Stack Overflow at address 0x453832

qui semble etre lie a la class C++ stackall.c

J'ai ensuite essaye de tester la meme chose avec un tableau de taille 10 et dans ce cas la, cela fonctionne correctement.

A quoi cela peut etre du? Est ce un probleme de capacite memoire maximum pour un tableau? Que puis je faire pour resoudre ce probleme.

Merci d'avance pour votre aide,

Yann

[Aurelien.Regat-Barrel]

Ton tableau est tros gros pour être stocké sur la pile. Utilise un tableau dynamique, et mieux encore, std::vector.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
std::vector<SimEvent> m_heap( 5000 );

[yanndublanche]

Merci pour ces renseignements.
Le but de mon tableau est de pouvoir trier les simevents par ordre croissant defini par l'operator <.
Est ce qu'avec un vecteur je pourrais aussi intervertir les elements comme dans un tableau?

[Aurelien.Regat-Barrel]

Haha... beaucoup mieux, tu peux trier ton tableau avec une ligne de plus.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
std::sort( m_heap.begin(), m_heap.end() );

Et voilà. Vive la STL!

[yanndublanche]

Merci pour ton aide

Le probleme du sort est que ce serait trop long a faire.
En fait, on cherche a simuler a l'aide de ce tableau un arbre binaire, on a donc sur le tableau une fonction PercolateUp:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
void SimEventPriorityQueue::percolateUp(const SimEvent& a_element, unsigned int a_pos){
	unsigned int pos = a_pos;
 
	// Percolate up
	while(pos > 0 && (a_element < m_heap[(pos-1)/2])){
		// The element is smaller than its parent, so push down the parent and move the hole
		m_heap[pos] = m_heap[(pos-1)/2];
		pos = (pos-1)/2;
	}
	// Found location: insert element
	m_heap[pos] = a_element;
}

Je pense que ca va fonctionner aussi un vecteur. Non?

[Aurelien.Regat-Barrel]

Oui un vecteur ça s'utilise comme un tableau classique. Quel est le but de ton arbre binaire ?

[yanndublanche]

Le but de mon arbre binaire est d'avoir a chaque fois en premier position l'element le plus petit sans avoir necessairement toute la liste triee, car a chque cycle d'execution, on execute le premier element, puis on le met a jour et on doit ensuite le repositionner dans la liste donc ce serait trop long si on faisait un sort de la liste a chaque fois.
Merci pour tes infos et tes conseils

[Luc Hermitte]

Avec un set, tu retires l'élément et tu le réinjectes après modification de son critère (pourle tri).

[Aurelien.Regat-Barrel]

Et pourquoi pas carrement faire uen simple recherche du plus petit (min_element) ? Car d'après ton algo, si j'ai bien compris, tu as un tableau de taille 2^n pout n éléments à stocker...

[yanndublanche]

J'ai un tableau de taille 2^n et 2^n elements a classer.
L'algorithme optimal est celui utilise avec arbre binaire car pour placer un element on a au plus n permutations a effectuer.