IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

C Discussion :

Boucle if qui ne fonctionne pas


Sujet :

C

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre averti
    Femme Profil pro
    Étudiant
    Inscrit en
    Janvier 2017
    Messages
    12
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 33
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Janvier 2017
    Messages : 12
    Par défaut Boucle if qui ne fonctionne pas
    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.

  2. #2
    Expert éminent

    Femme Profil pro
    Ingénieur développement logiciels
    Inscrit en
    Juin 2007
    Messages
    5 202
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Essonne (Île de France)

    Informations professionnelles :
    Activité : Ingénieur développement logiciels

    Informations forums :
    Inscription : Juin 2007
    Messages : 5 202
    Par défaut
    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?

  3. #3
    Membre averti
    Femme Profil pro
    Étudiant
    Inscrit en
    Janvier 2017
    Messages
    12
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 33
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Janvier 2017
    Messages : 12
    Par défaut
    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.

  4. #4
    Membre averti
    Femme Profil pro
    Étudiant
    Inscrit en
    Janvier 2017
    Messages
    12
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 33
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Janvier 2017
    Messages : 12
    Par défaut
    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.

  5. #5
    Membre averti
    Femme Profil pro
    Étudiant
    Inscrit en
    Janvier 2017
    Messages
    12
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 33
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Janvier 2017
    Messages : 12
    Par défaut
    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.

  6. #6
    Modérateur
    Avatar de gangsoleil
    Homme Profil pro
    Manager / Cyber Sécurité
    Inscrit en
    Mai 2004
    Messages
    10 150
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Haute Savoie (Rhône Alpes)

    Informations professionnelles :
    Activité : Manager / Cyber Sécurité

    Informations forums :
    Inscription : Mai 2004
    Messages : 10 150
    Par défaut
    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 ?
    "La route est longue, mais le chemin est libre" -- https://framasoft.org/
    Les règles du forum

  7. #7
    Membre averti
    Femme Profil pro
    Étudiant
    Inscrit en
    Janvier 2017
    Messages
    12
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Âge : 33
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Janvier 2017
    Messages : 12
    Par défaut
    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.

Discussions similaires

  1. Une boucle While qui ne fonctionne pas
    Par seremei dans le forum Langage
    Réponses: 6
    Dernier message: 25/10/2012, 11h33
  2. [XL-2010] Boucle For qui ne fonctionne pas pour toutes les valeurs
    Par pleskyd dans le forum Macros et VBA Excel
    Réponses: 15
    Dernier message: 13/02/2012, 15h16
  3. une boucle for qui ne fonctionne pas
    Par piffeo dans le forum Général JavaScript
    Réponses: 3
    Dernier message: 06/11/2008, 21h52
  4. boucle if-then qui ne fonctionne pas
    Par les4c77 dans le forum Macros et VBA Excel
    Réponses: 3
    Dernier message: 17/09/2007, 16h24
  5. [MySQL] Boucle qui ne fonctionne pas bien que la 1ère fois
    Par R.L. dans le forum PHP & Base de données
    Réponses: 8
    Dernier message: 27/05/2006, 21h59

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo