Boucle if qui ne fonctionne pas

**hbemar** · 27/03/2017, 11h16

Bonjour à tous!
Voilà je réalise un modèle en C que je compile après sous R. Mais voilà, j'ai deux de mes boucles if pour les paramètres gamma et beta qui ne
fonctionne pas car peut importe la valeur que j'attribue à ces deux paramètres, ils sont toujours de 0. Je crois que ça vient de la façon dont je les déclares mais je ne suis pas sûre et de toute façon je ne vois pas d'autre manière de les déclarés ( je suis novice en C).

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
 
 
static double parms[9];
#define K parms[0]
#define kapa parms[1]
#define t0 parms[2]
#define alpha parms[3]
#define sigma parms[4]
#define ts parms[5]
#define te parms[6]
#define gamma parms[7]
#define beta parms[8]
 
 
/* initializer  */
void initmod(void (* odeparms)(int *, double *))
{
    int N = 9;
    odeparms(&N, parms);
}
 
#define S y[0]
#define I y[1]
#define R y[2]
#define dSdt ydot[0]
#define dIdt ydot[1]
#define dRdt ydot[2]
 
 
/* Derivatives and 1 output variable */
void derivs (int *neq, double *t, double *y, double *ydot,
             double *yout, int *ip)
{
    if (ip[0] < 1) error("nout should be at least 1") ;
 
 
 
if (*t<=te) gamma=0* *t; 
if (*t>te) gamma=gamma;
 
 
if (*t<=te) beta=0* *t; 
if (*t>te) beta=beta;
 
 
double tau ;
tau =  alpha * exp(-kapa * (*t-t0));
 
 
double delta ;
if(*t<=ts) delta=0* *t;
if(*t>ts) delta=sigma *(1- exp(-kapa * (*t-ts)));
 
 
    dSdt = -beta*S*I + tau * S * (1 - S / K )- delta*S ; 
    dIdt = beta * S * I- (gamma + delta )*I;
    dRdt= delta*S +(gamma + delta)*I;
 
yout[0] = S + I + R ;
 
}

J'ai essayer d'autre façon de l'écrire avec des accolades pour les if ou encore en faisant:

Code :

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

1
2
3
4
5
6
7
double gamma1 ;
if(*t<=te) gamma1=0* *t;
if(*t>te) gamma1=gamma;
 
double beta1 ;
if(*t<=te) beta1=0* *t;
if(*t>te) beta1=beta;

Mais à chaque fois rien ne marche, gamma et beta reste à 0.

Avez vous une idée de ce qu'il conviendrait de faire?

Merci beaucoup pour votre aide.

**ternel** · 27/03/2017, 12h03

Question vocabulaire, if n'est pas une boucle, mais une "structure de controle de flux", qu'on appelle en langage normal une "conditionnelle" ou simplement "un if".

Il nous manque un peu de code.
Ta fonction d'initialisation (initmod) est-elle bien appelée avant ta fonction derivs?

**hbemar** · 27/03/2017, 13h11

Bonjour,
déjà merci de m'avoir répondu. Oui désolé pour mon langage, je ferais plus attention la prochaine fois.

Pour le script sous C, celui que je vous ai donné dans mon premier message est complet, mis à part:

Code :

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

1
2
3
4
5
6
7
8
/* file lorenz.c */
#include <R.h>
#include <Rinternals.h>
#include <Rdefines.h>
#include <R_ext/Rdynload.h>
#include <Rmath.h>
#include <stdio.h>
#include <math.h>

au début du script.

Après, je compile le sous-programmes C (appelé "epidemic") via la commande R CMD SHLIB sous linux.

Puis une fois ça fait, je l'ouvre sous R et je fais:

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
 
#sélection du modèle:
dyn.load("epidemic.so")
#vérifie si la fonction est bien chargée:
is.loaded("derivs")
 
#temps
time.out <- seq(0,1000,by =0.1)
#variables initiales:
init<-c(S=0.01, I=0.01, R=0.01)
#parametres
values<-c(K=1,kapa=1.029871e-03,t0=0,alpha=2.374821e-02,sigma=3.445653e-03,ts=4.504544e+02,beta=0.3,gamma=0.02,te=444.7)
 
#calcul:
SIR.C <-as.data.frame(radau(init,time.out,"derivs",values,dllname="epidemic",initfunc="initmod",nout=1,outnames="N"))

Et c'est à ce moment que je remarque que peut importe les valeur de beta et gamma dans le vecteur "values", cela ne changent rien à la dynamique de S, I et R car en réalité beta et gamma ont toujours une valeur de 0 même lorsque t est supérieur à te=444.7.

Voilà, j'espère que cela permettra d'identifier le problème.

**hbemar** · 27/03/2017, 13h45

Huum je crois avoir trouvé le problème mais ça me parait bizarre,
j'ai écris dans le code c:

Code :

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

1
2
3
4
5
6
7
double gammabis ;
if(*t<=ts) gammabis=0* *t;
if(*t>ts) gammabis=gamma; 
 
double betabis ;
if(*t<=ts) betabis=0* *t;
if(*t>ts) betabis=beta;

Et sous R, le vecteur de paramètre "value" doit suivre exactement l'ordre dans lequel les paramètres sont déclarés sous C.
En effet, ça fonctionne si je fait:

Code :

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

values<-c(K=1,kapa=1.029871e-03,t0=0,alpha=2.374821e-02,sigma=3.445653e-03,ts=4.504544e+02,te=444.7,gamma=0.02, beta=0.1)

Mais ça ne marche plus si j'inverse beta et gamma:

Code :

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

values<-c(K=1,kapa=1.029871e-03,t0=0,alpha=2.374821e-02,sigma=3.445653e-03,ts=4.504544e+02,te=444.7,gamma=0.02, beta=0.1)

Étrange non? Peut être est ce normal.

**hbemar** · 27/03/2017, 14h37

Huum je crois avoir trouvé le problème mais ça me parait bizarre,
j'ai écris dans le code c:

Code :

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

1
2
3
4
5
6
7
double gammabis ;
if(*t<=ts) gammabis=0* *t;
if(*t>ts) gammabis=gamma; 
 
double betabis ;
if(*t<=ts) betabis=0* *t;
if(*t>ts) betabis=beta;

Et sous R, le vecteur de paramètre "value" doit suivre exactement l'ordre dans lequel les paramètres sont déclarés sous C.
En effet, ça fonctionne si je fait:

Code :

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

values<-c(K=1,kapa=1.029871e-03,t0=0,alpha=2.374821e-02,sigma=3.445653e-03,ts=4.504544e+02,te=444.7,gamma=0.02, beta=0.1)

Mais ça ne marche plus si j'inverse beta et gamma:

Code :

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

values<-c(K=1,kapa=1.029871e-03,t0=0,alpha=2.374821e-02,sigma=3.445653e-03,ts=4.504544e+02,te=444.7,gamma=0.02, beta=0.1)

Étrange non? Peut être est ce normal.

**gangsoleil** · 28/03/2017, 15h34

Code :

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

1
2
3
4
5
6
7
double gammabis ; 
if(*t<=ts) gammabis=0* *t;
if(*t>ts) gammabis=gamma; 
 
double betabis ;
if(*t<=ts) betabis=0* *t;
if(*t>ts) betabis=beta;

Si *t est inferieur a ts, tu affectes 0 a gammabis ou betabis, via la multiplication de 0 et de *t.... Est-ce bien ce que tu veux faire ?

**hbemar** · 30/03/2017, 11h51

Oui Oui, c'est bien ce que je voulais faire.

Mais en fait j'ai identifier le problème, c'est que je n'avais pas rentrer les paramètres sur R dans le vecteur "values" dans le même ordre que celui dans lequel j'ai déclaré les paramètres sur C. Je pensais que le fait de les nommés dans le vecteur sur R fonctionnerai mais en réalité non.

Boucle if qui ne fonctionne pas

C

Vue hybride

Discussions similaires

Partager

Partager