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

Python Discussion :

Boucle while infinie après passage de C++ à python [Python 3.X]


Sujet :

Python

Vue hybride

psycopatate35 Boucle while infinie après... 02/01/2022, 21h08
wiztricks Pour trouver la solution, il... 02/01/2022, 23h23
psycopatate35 C'est effectivement tout le... 03/01/2022, 08h58
wiztricks C'est peut être un défaut... 03/01/2022, 09h27
Sve@r Bonjour Comme le dit... 03/01/2022, 11h13
lg_53 Dans ce cas pourquoi cherchez... 04/01/2022, 09h49
psycopatate35 Je cherche à modifier des... 04/01/2022, 10h01
jurassic pork hello, le souci c'est... 03/01/2022, 13h30
marco056 A priori, cela ne fonctionne... 03/01/2022, 13h58
jurassic pork heu c'est bizarre chez moi... 03/01/2022, 14h00
marco056 Le programme se termine chez... 03/01/2022, 14h22
marco056 Je n'y connais rien en C++... 03/01/2022, 14h21
Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre à l'essai
    Homme Profil pro
    Étudiant
    Inscrit en
    Janvier 2022
    Messages
    4
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 30
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant
    Secteur : Enseignement

    Informations forums :
    Inscription : Janvier 2022
    Messages : 4
    Par défaut Boucle while infinie après passage de C++ à python
    Bonjour et bonne année à tous !

    J'ai récemment tenté de convertir en python un programme créé par un collègue en C++ (ci-après le début) afin de le modifier un peu et de le rendre plus automatique.
    Code C++ : 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
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127
    128
    129
    130
    131
    132
    133
    134
    135
    136
    137
    138
    139
    140
    141
    142
    143
    144
    145
    146
    147
    148
    149
    150
    151
    #include <stdlib.h>
    #include <stdio.h>
    #include <math.h>
     
    void pause();
     
     
    int main(int argc, char *argv[])
     
    {
       int Nmax = 4000;
       int I = 0;
       int J = 0;
       int N = 0;
       long double Psif[8100];
       long double PsiTotal[8100];
       long double Rqd[8100];
       long double Pota[8100];
       long double Dtau;
       long double Csi;
       long double Tau = 0;
       long double Tauf = 0;
       long double u;
       long double Pi = 3.14159265358979323846264338327950288419716939937510;
       long double Sum;
       long double Ct1;
       long double Ct2;
       long double Ct3;
       long double Ct4;
       long double Ct5;
       long double Lambda;
       long double Alpha = 0.5;
       long double PsiProv = 0;
       long double CsiPrime;
       long double Beta = 0; // Zero pour lest tests
       long double Gamma = 0; // Zero pour les tests
       long double Rau;
       long double CsiP;
       long double PsiP;
       long double CsiM;
       long double PsiM;
       double Ei=-1.0;
       double Ef=1.0;
       double E0=0.39;
       double S=0.0706858;
       S=S/10000;
       double C0=0.001;
       C0=C0*1000;
       double T=294.65;
       double D=0.0001;
       D=D/10000;
       double Ru=165;
       double Cd=0.000017;
       double ks=20;
       ks=ks/100;
       double v=0.05;
       long double FRT;
       long double FSC0;
       long double DFRT;
       long double F = 96485.3329;
       long double R = 8.3144621;
       long double Ep1;
       long double Ep2;
       long double Ip1;
       long double Ip2;
        FRT = F / R / T;
        FSC0 = F * S * C0;
        DFRT = D * FRT;
     
        //parameters are defined for oxidation
        u = -(Ei - E0) * FRT;
        CsiM = (Ef - E0) * FRT;
        Lambda = ks / sqrt (v * DFRT);
        Rau = FRT * Ru * FSC0 * sqrt (v *DFRT);
        Gamma = Cd * sqrt (v) / FSC0 /sqrt (DFRT);
     
        printf (" Ei = %f ; Ef = %f ; S = %e ;  C0 = %e  \n", Ei, Ef, S, C0);
        printf (" T = %f ; D = %e ; Ru = %e ;  Cd = %e  \n", T, D, Ru, Cd);
        printf (" ks = %f ; v = %f \n", ks, v);
     
     
        Psif[0] = 0;
        Psif[1] = 0;
        Dtau = (CsiM + u) / Nmax; // calculation of the time step.
        Tau = Dtau;
        Csi = -u;
        Beta = Rau * Gamma;
        Ct3 = Beta/Dtau;
     
        printf (" Lambda = %LF ; Rau = %Lf ; Gamma = %Lf ; Dtau = %Lf \n", Lambda, Rau, Gamma, Dtau);
     
        for (I = 1 ; I <= 2.0 * Nmax + 1; I++)
     
            {
                Rqd[I] = (powl ((I + 1.0) , 1.5)) - (2.0 * powl (I, 1.5)) + (powl ((I - 1) , 1.5));
            }
     
     // Forward Scan
     
     
        for (I = 1 ; I <= Nmax ; I++)
            {
                Ct1 = expl (-Alpha * Csi) / Lambda;
                Ct2 = 4.0 / 3.0 * sqrtl (Dtau / Pi) * (1.0 + expl (-Csi));
                Pota[I] = Csi;
     
                Sum = 0;
     
                for (N = 1 ; N < I ; N++)
     
                    {
                        Sum = Sum + (Psif[N] * (1.0 + Ct3) - Psif[N - 1] * Ct3) * Rqd[I - N];
                    }
     
                N = 0;
     
                Psif[I] = Psif[I - 1];
     
                do
     
                    {
                        Psif[I] = (Psif[I] + 100.0 * PsiProv) /101.0; // Il faut calmer la convergence
                        PsiProv = Psif[I];
                        CsiPrime = Csi - Rau * Psif[I] - Beta * (1.0 - expl (-Tau/Beta));
                        Ct1 = expl (-Alpha * CsiPrime) / Lambda;
                        Ct2 = 4.0 / 3.0 * sqrtl (Dtau / Pi) * (1.0 + expl (-CsiPrime));
     
                        Psif[I] = (1.0 - Sum * Ct2 + (Ct1 * Ct3 + Ct2 * Ct3) * Psif[I - 1] ) / (Ct1 + Ct2 + Ct1 * Ct3 + Ct2 * Ct3);
     
                        N++;
     
                    }
     
                while ( fabs ( Psif[I] - PsiProv ) > 0.0001 );
     
                PsiTotal[I] = Psif[I] + Gamma * (1.0 - expl (-Tau / Beta));
     
                if ((I % 200) == 0)
     
                {
     
                    printf (" I= %d Csi= %Lf Psif= %Lf PsiTotal= %Lf N= %d \n", I, Csi, Psif[I], PsiTotal[I], N);
     
                }
     
                Csi = Csi + Dtau;
                Tau = Tau + Dtau;
            }
     
        Tauf = Tau - Dtau;
    }

    Cependant, je me heurte à un problème de boucle infinie dans mon programme python et qui n'existe pas en C++.
    Code Python : 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
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    import math
    import numpy as np
    import statistics
     
    Nmax = int(8000)
    I = int(0)
    J = int(0)
    N = int(0)
     
    Psif = [0,0]
    PsiTotal = [0,0]
    Rqd = [0,0]
    Pota = [0,0]
    for i in range(3*Nmax):
        Psif.append(np.longdouble(0e0))
        PsiTotal.append(np.longdouble(0e0))
        Rqd.append(np.longdouble(0e0))
        Pota.append(np.longdouble(0e0))
     
    Tau = float(0)
    Tauf = 0
    PsiProv = np.longdouble(0e0)
    F = float(96485.3329)
    R = float(8.3144621)
     
    Ei=-1
    Ef=1
    E0=0.39
    S=0.0706858
    S=S/10000
    C0=0.001
    C0=C0*1000
    T=294.65
    D=0.0001
    D=D/10000
    Ru=165
    Cd=0.000017
    ks=20
    ks=ks/100
    Alpha = 0.5
    v=0.05
     
    FRT=F/R/T
    FSC0=F*S*C0
    DFRT=D*FRT
    u=-(Ei-E0)*FRT
    CsiM=(Ef-E0)*FRT
    Lambda=ks/math.sqrt(v*DFRT)
    Rau=FRT*Ru*FSC0*math.sqrt(v*DFRT)
    Gamma=Cd*math.sqrt(v)/FSC0/math.sqrt(DFRT)
     
    "calculation of step time"
    Dtau = (CsiM+u)/Nmax
    Tau=Dtau
    Csi=-u
    Beta=Rau*Gamma
    Ct3=Beta/Dtau
     
    print(f"Lambda = {Lambda} ; Rau = {Rau} ; Gamma = {Gamma} ; Dtau = {Dtau}")
     
    for i in range(2*Nmax+1):
        Rqd[i]=math.pow((i+2),1.5)-(2*math.pow(i+1,1.5))+(math.pow((i),1.5))
     
    "Forward Scan"
     
    for i in range(Nmax+1):
        if i==0:
            pass
        else:
            Ct1=math.exp(-Alpha*Csi/Lambda)
            Ct2=np.longdouble(4/3*math.sqrt(Dtau/math.pi)*(1+math.exp(-Csi)))
            Pota[i]=Csi
            Sum=0
            for N in range(i):
                if N==0:
                    pass
                else :
                    Sum=Sum+(Psif[N]*(1+Ct3)-Psif[N-1]*Ct3)*Rqd[i-N]
            N=int(0)
            Psif[i]=Psif[i-1]
            while True:
                Psif
                Psif[i] = np.longdouble((Psif[i] + 100*PsiProv)/101)
                PsifProv = np.longdouble(Psif[i])
                CsiPrime = Csi - Rau*Psif[i] -Beta*(1+math.exp(-Tau/Beta))
                Ct1 = math.exp(-Alpha*CsiPrime)/Lambda
                Ct2 = np.longdouble(4/3*math.sqrt(Dtau/math.pi)*(1+math.exp(-CsiPrime)))
                Psif[i] = np.longdouble((1-Sum*Ct2+(Ct1*Ct3+Ct2*Ct3)*Psif[i-1])/(Ct1+Ct2+Ct1*Ct3+Ct2*Ct3))
                N=N+1
                if abs(Psif[i]-PsiProv)<=0.0001:
                    break
            PsiTotal[i]=Psif[i]+Gamma*(1-math.exp(-Tau/Beta))
            if i%200==0 and i!=0:
                print(f"I= {i}, Csi= {Csi}, Psif= {Psif[i]}, PsiTotal= {PsiTotal[i]}, N= {N}")
            Csi=Csi+Dtau
            Tau=Tau+Dtau
    Tauf=Tau-Dtau

    Quelqu'un aurait la solution svp ?

  2. #2
    Expert éminent
    Homme Profil pro
    Architecte technique retraité
    Inscrit en
    Juin 2008
    Messages
    21 723
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Manche (Basse Normandie)

    Informations professionnelles :
    Activité : Architecte technique retraité
    Secteur : Industrie

    Informations forums :
    Inscription : Juin 2008
    Messages : 21 723
    Par défaut
    Citation Envoyé par psycopatate35 Voir le message
    Quelqu'un aurait la solution svp ?
    Pour trouver la solution, il faut essayer de comprendre le problème.

    Au pif, il y a un while True dont on ne sort qu'à la condition abs(Psif[i]-PsiProv)<=0.0001. Déjà vérifier si cette condition est réalisée sinon voir ce qui a été loupé dans la traduction du C++...

    C'est juste le boulot qu'on fait pour mettre au point un programme.... et je ne vais pas le faire pour vous.

    - W
    Architectures post-modernes.
    Python sur DVP c'est aussi des FAQs, des cours et tutoriels

  3. #3
    Membre à l'essai
    Homme Profil pro
    Étudiant
    Inscrit en
    Janvier 2022
    Messages
    4
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 30
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant
    Secteur : Enseignement

    Informations forums :
    Inscription : Janvier 2022
    Messages : 4
    Par défaut
    C'est effectivement tout le problème !

    La boucle Do While de C++ fonctionne parfaitement et arrive à sortir pour l'ensemble de la boucle for définie par I = [1;4000].

    En python, ça fonctionne également mais seulement jusqu'à I = 2310. Jusque la, on obtient sensiblement les mêmes réponses qu'en C++. Par contre, dès que I = 2311, la boucle While n'arrive plus à atteindre la condition de sortie abs(Psif[i]-PsiProv)<=0.0001, qui pourtant ce fait bien en C++.

    La programmation n'étant pas mon domaine principal d'expertise, j'ai bien tenté de comprendre avec des proches ce qui coince mais sans succès.

  4. #4
    Expert éminent
    Homme Profil pro
    Architecte technique retraité
    Inscrit en
    Juin 2008
    Messages
    21 723
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Manche (Basse Normandie)

    Informations professionnelles :
    Activité : Architecte technique retraité
    Secteur : Industrie

    Informations forums :
    Inscription : Juin 2008
    Messages : 21 723
    Par défaut
    C'est peut être un défaut côté précision.
    Là ou il y a du boulot c'est dans la vérification de la cohérence entre types de données C++, leur traduction en types numpy et s'assurer que les opérations dans le code ne le dégradent pas.

    - W
    Architectures post-modernes.
    Python sur DVP c'est aussi des FAQs, des cours et tutoriels

  5. #5
    Membre prolifique
    Avatar de Sve@r
    Homme Profil pro
    Ingénieur développement logiciels
    Inscrit en
    Février 2006
    Messages
    12 830
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Oise (Picardie)

    Informations professionnelles :
    Activité : Ingénieur développement logiciels
    Secteur : Aéronautique - Marine - Espace - Armement

    Informations forums :
    Inscription : Février 2006
    Messages : 12 830
    Billets dans le blog
    1
    Par défaut
    Bonjour
    Citation Envoyé par psycopatate35 Voir le message
    En python, ça fonctionne également mais seulement jusqu'à I = 2310. Jusque la, on obtient sensiblement les mêmes réponses qu'en C++. Par contre, dès que I = 2311, la boucle While n'arrive plus à atteindre la condition de sortie abs(Psif[i]-PsiProv)<=0.0001, qui pourtant ce fait bien en C++.
    Comme le dit wiztricks, il peut y avoir un souci de précision. Souci qui est d'ailleurs parfaitement connu et qui n'est pas spécifique à Python mais dû à la façon dont un ordinateur convertit un décimal en binaire (tape par exemple 0.6 / 3 puis 0.2 * 3 dans une console Python et tu comprendras).
    Bon ce qui m'étonne c'est que ce souci étant aussi présent donc en C++, ça n'aurait pas dû fonctionner en C++ mais peut-être qu'en C++ la profondeur de la précision n'est pas la même et que cette différence fait justement la différence.

    Ce que je peux te proposer en retour c'est d'utiliser le module decimal, lequel travaille chiffre par chiffre et non pas en binaire. C'est plus long mais c'est garanti sans erreur. Le même exemple que tout à l'heure en utilisant ce module: decimal.Decimal("0.6") / 3 ou decimal.Decimal("0.2") * 3.
    Mon Tutoriel sur la programmation «Python»
    Mon Tutoriel sur la programmation «Shell»
    Sinon il y en a pleins d'autres. N'oubliez pas non plus les différentes faq disponibles sur ce site
    Et on poste ses codes entre balises [code] et [/code]

  6. #6
    Membre Expert

    Homme Profil pro
    Ingénieur calcul scientifique
    Inscrit en
    Mars 2013
    Messages
    1 229
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France, Alpes Maritimes (Provence Alpes Côte d'Azur)

    Informations professionnelles :
    Activité : Ingénieur calcul scientifique

    Informations forums :
    Inscription : Mars 2013
    Messages : 1 229
    Par défaut
    Citation Envoyé par psycopatate35 Voir le message
    La programmation n'étant pas mon domaine principal d'expertise
    Dans ce cas pourquoi cherchez vous à réécrire en python ce qui fonctionne très bien en C++ ?

  7. #7
    Membre à l'essai
    Homme Profil pro
    Étudiant
    Inscrit en
    Janvier 2022
    Messages
    4
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 30
    Localisation : France, Ille et Vilaine (Bretagne)

    Informations professionnelles :
    Activité : Étudiant
    Secteur : Enseignement

    Informations forums :
    Inscription : Janvier 2022
    Messages : 4
    Par défaut
    Citation Envoyé par lg_53 Voir le message
    Dans ce cas pourquoi cherchez vous à réécrire en python ce qui fonctionne très bien en C++ ?
    Je cherche à modifier des choses dans le programme, notamment l'automatiser, et créer une petite interface graphique, et python est l'un des langages les plus simple pour cela apparemment.

  8. #8
    Expert confirmé
    Avatar de jurassic pork
    Homme Profil pro
    Bidouilleur
    Inscrit en
    Décembre 2008
    Messages
    4 210
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Bidouilleur
    Secteur : Industrie

    Informations forums :
    Inscription : Décembre 2008
    Messages : 4 210
    Par défaut
    hello,
    Citation Envoyé par psycopatate35 Voir le message
    Je cherche à modifier des choses dans le programme, notamment l'automatiser, et créer une petite interface graphique, et python est l'un des langages les plus simple pour cela apparemment.
    pour les interfaces graphiques en C++ il y a Qt :
    Qt, bien connu pour son framework C++, avec son environnement de développement intégré (IDE), appelé Qt Creator. Cette application intègre un éditeur de code C++ avancé avec détection d'erreur syntaxique et intégration automatique de fonction ou de variables, un éditeur de formulaire et d'interface graphique, un outil de gestion de projet, une aide contextuelle intégrée, un débugueur visuel, etc. Idéal pour les développement cross-plateforme.
    Ami calmant, J.P

  9. #9
    Expert confirmé
    Avatar de jurassic pork
    Homme Profil pro
    Bidouilleur
    Inscrit en
    Décembre 2008
    Messages
    4 210
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Bidouilleur
    Secteur : Industrie

    Informations forums :
    Inscription : Décembre 2008
    Messages : 4 210
    Par défaut
    hello,
    Citation Envoyé par psycopatate35 Voir le message
    Bonjour et bonne année à tous !

    J'ai récemment tenté de convertir en python un programme créé par un collègue en C++ (ci-après le début) afin de le modifier un peu et de le rendre plus automatique.
    Code C++ : 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
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
    100
    101
    102
    103
    104
    105
    106
    107
    108
    109
    110
    111
    112
    113
    114
    115
    116
    117
    118
    119
    120
    121
    122
    123
    124
    125
    126
    127
    128
    129
    130
    131
    132
    133
    134
    135
    136
    137
    138
    139
    140
    141
    142
    143
    144
    145
    146
    147
    148
    149
    150
    151
    #include <stdlib.h>
    #include <stdio.h>
    #include <math.h>
     
    void pause();
     
     
    int main(int argc, char *argv[])
     
    {
       int Nmax = 4000;
       int I = 0;
       int J = 0;
       int N = 0;
       long double Psif[8100];
       long double PsiTotal[8100];
       long double Rqd[8100];
       long double Pota[8100];
       long double Dtau;
       long double Csi;
       long double Tau = 0;
       long double Tauf = 0;
       long double u;
       long double Pi = 3.14159265358979323846264338327950288419716939937510;
       long double Sum;
       long double Ct1;
       long double Ct2;
       long double Ct3;
       long double Ct4;
       long double Ct5;
       long double Lambda;
       long double Alpha = 0.5;
       long double PsiProv = 0;
       long double CsiPrime;
       long double Beta = 0; // Zero pour lest tests
       long double Gamma = 0; // Zero pour les tests
       long double Rau;
       long double CsiP;
       long double PsiP;
       long double CsiM;
       long double PsiM;
       double Ei=-1.0;
       double Ef=1.0;
       double E0=0.39;
       double S=0.0706858;
       S=S/10000;
       double C0=0.001;
       C0=C0*1000;
       double T=294.65;
       double D=0.0001;
       D=D/10000;
       double Ru=165;
       double Cd=0.000017;
       double ks=20;
       ks=ks/100;
       double v=0.05;
       long double FRT;
       long double FSC0;
       long double DFRT;
       long double F = 96485.3329;
       long double R = 8.3144621;
       long double Ep1;
       long double Ep2;
       long double Ip1;
       long double Ip2;
        FRT = F / R / T;
        FSC0 = F * S * C0;
        DFRT = D * FRT;
     
        //parameters are defined for oxidation
        u = -(Ei - E0) * FRT;
        CsiM = (Ef - E0) * FRT;
        Lambda = ks / sqrt (v * DFRT);
        Rau = FRT * Ru * FSC0 * sqrt (v *DFRT);
        Gamma = Cd * sqrt (v) / FSC0 /sqrt (DFRT);
     
        printf (" Ei = %f ; Ef = %f ; S = %e ;  C0 = %e  \n", Ei, Ef, S, C0);
        printf (" T = %f ; D = %e ; Ru = %e ;  Cd = %e  \n", T, D, Ru, Cd);
        printf (" ks = %f ; v = %f \n", ks, v);
     
     
        Psif[0] = 0;
        Psif[1] = 0;
        Dtau = (CsiM + u) / Nmax; // calculation of the time step.
        Tau = Dtau;
        Csi = -u;
        Beta = Rau * Gamma;
        Ct3 = Beta/Dtau;
     
        printf (" Lambda = %LF ; Rau = %Lf ; Gamma = %Lf ; Dtau = %Lf \n", Lambda, Rau, Gamma, Dtau);
     
        for (I = 1 ; I <= 2.0 * Nmax + 1; I++)
     
            {
                Rqd[I] = (powl ((I + 1.0) , 1.5)) - (2.0 * powl (I, 1.5)) + (powl ((I - 1) , 1.5));
            }
     
     // Forward Scan
     
     
        for (I = 1 ; I <= Nmax ; I++)
            {
                Ct1 = expl (-Alpha * Csi) / Lambda;
                Ct2 = 4.0 / 3.0 * sqrtl (Dtau / Pi) * (1.0 + expl (-Csi));
                Pota[I] = Csi;
     
                Sum = 0;
     
                for (N = 1 ; N < I ; N++)
     
                    {
                        Sum = Sum + (Psif[N] * (1.0 + Ct3) - Psif[N - 1] * Ct3) * Rqd[I - N];
                    }
     
                N = 0;
     
                Psif[I] = Psif[I - 1];
     
                do
     
                    {
                        Psif[I] = (Psif[I] + 100.0 * PsiProv) /101.0; // Il faut calmer la convergence
                        PsiProv = Psif[I];
                        CsiPrime = Csi - Rau * Psif[I] - Beta * (1.0 - expl (-Tau/Beta));
                        Ct1 = expl (-Alpha * CsiPrime) / Lambda;
                        Ct2 = 4.0 / 3.0 * sqrtl (Dtau / Pi) * (1.0 + expl (-CsiPrime));
     
                        Psif[I] = (1.0 - Sum * Ct2 + (Ct1 * Ct3 + Ct2 * Ct3) * Psif[I - 1] ) / (Ct1 + Ct2 + Ct1 * Ct3 + Ct2 * Ct3);
     
                        N++;
     
                    }
     
                while ( fabs ( Psif[I] - PsiProv ) > 0.0001 );
     
                PsiTotal[I] = Psif[I] + Gamma * (1.0 - expl (-Tau / Beta));
     
                if ((I % 200) == 0)
     
                {
     
                    printf (" I= %d Csi= %Lf Psif= %Lf PsiTotal= %Lf N= %d \n", I, Csi, Psif[I], PsiTotal[I], N);
     
                }
     
                Csi = Csi + Dtau;
                Tau = Tau + Dtau;
            }
     
        Tauf = Tau - Dtau;
    }

    Cependant, je me heurte à un problème de boucle infinie dans mon programme python et qui n'existe pas en C++.
    Code Python : 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
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    import math
    import numpy as np
    import statistics
     
    Nmax = int(8000)
    I = int(0)
    J = int(0)
    N = int(0)
     
    Psif = [0,0]
    PsiTotal = [0,0]
    Rqd = [0,0]
    Pota = [0,0]
    for i in range(3*Nmax):
        Psif.append(np.longdouble(0e0))
        PsiTotal.append(np.longdouble(0e0))
        Rqd.append(np.longdouble(0e0))
        Pota.append(np.longdouble(0e0))
     
    Tau = float(0)
    Tauf = 0
    PsiProv = np.longdouble(0e0)
    F = float(96485.3329)
    R = float(8.3144621)
     
    Ei=-1
    Ef=1
    E0=0.39
    S=0.0706858
    S=S/10000
    C0=0.001
    C0=C0*1000
    T=294.65
    D=0.0001
    D=D/10000
    Ru=165
    Cd=0.000017
    ks=20
    ks=ks/100
    Alpha = 0.5
    v=0.05
     
    FRT=F/R/T
    FSC0=F*S*C0
    DFRT=D*FRT
    u=-(Ei-E0)*FRT
    CsiM=(Ef-E0)*FRT
    Lambda=ks/math.sqrt(v*DFRT)
    Rau=FRT*Ru*FSC0*math.sqrt(v*DFRT)
    Gamma=Cd*math.sqrt(v)/FSC0/math.sqrt(DFRT)
     
    "calculation of step time"
    Dtau = (CsiM+u)/Nmax
    Tau=Dtau
    Csi=-u
    Beta=Rau*Gamma
    Ct3=Beta/Dtau
     
    print(f"Lambda = {Lambda} ; Rau = {Rau} ; Gamma = {Gamma} ; Dtau = {Dtau}")
     
    for i in range(2*Nmax+1):
        Rqd[i]=math.pow((i+2),1.5)-(2*math.pow(i+1,1.5))+(math.pow((i),1.5))
     
    "Forward Scan"
     
    for i in range(Nmax+1):
        if i==0:
            pass
        else:
            Ct1=math.exp(-Alpha*Csi/Lambda)
            Ct2=np.longdouble(4/3*math.sqrt(Dtau/math.pi)*(1+math.exp(-Csi)))
            Pota[i]=Csi
            Sum=0
            for N in range(i):
                if N==0:
                    pass
                else :
                    Sum=Sum+(Psif[N]*(1+Ct3)-Psif[N-1]*Ct3)*Rqd[i-N]
            N=int(0)
            Psif[i]=Psif[i-1]
            while True:
                Psif
                Psif[i] = np.longdouble((Psif[i] + 100*PsiProv)/101)
                PsifProv = np.longdouble(Psif[i])
                CsiPrime = Csi - Rau*Psif[i] -Beta*(1+math.exp(-Tau/Beta))
                Ct1 = math.exp(-Alpha*CsiPrime)/Lambda
                Ct2 = np.longdouble(4/3*math.sqrt(Dtau/math.pi)*(1+math.exp(-CsiPrime)))
                Psif[i] = np.longdouble((1-Sum*Ct2+(Ct1*Ct3+Ct2*Ct3)*Psif[i-1])/(Ct1+Ct2+Ct1*Ct3+Ct2*Ct3))
                N=N+1
                if abs(Psif[i]-PsiProv)<=0.0001:
                    break
            PsiTotal[i]=Psif[i]+Gamma*(1-math.exp(-Tau/Beta))
            if i%200==0 and i!=0:
                print(f"I= {i}, Csi= {Csi}, Psif= {Psif[i]}, PsiTotal= {PsiTotal[i]}, N= {N}")
            Csi=Csi+Dtau
            Tau=Tau+Dtau
    Tauf=Tau-Dtau

    Quelqu'un aurait la solution svp ?
    le souci c'est que math.pi de python n'a pas la précision voulue ((~14 digits). Pour avoir une meilleure précision on peut utiliser le module mpmath (présent dans pypi).
    1 - On définit la précision pour les calculs ( ex:50 ):
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    from mpmath import mp
    mp.prec = 50
    2 - dans le code on remplace math.pi par mp.pi

    Ami calmant, J.P

  10. #10
    Membre Expert
    Homme Profil pro
    Enseignant
    Inscrit en
    Juin 2013
    Messages
    1 617
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Enseignant
    Secteur : Enseignement

    Informations forums :
    Inscription : Juin 2013
    Messages : 1 617
    Par défaut
    Citation Envoyé par jurassic pork Voir le message
    hello,


    le souci c'est que math.pi de python n'a pas la précision voulue ((~14 digits). Pour avoir une meilleure précision on peut utiliser le module mpmath (présent dans pypi).
    1 - On définit la précision pour les calculs ( ex:50 ):
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    from mpmath import mp
    mp.prec = 50
    2 - dans le code on remplace math.pi par mp.pi

    Ami calmant, J.P
    A priori, cela ne fonctionne pas mieux avec tes modifs.

  11. #11
    Expert confirmé
    Avatar de jurassic pork
    Homme Profil pro
    Bidouilleur
    Inscrit en
    Décembre 2008
    Messages
    4 210
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Bidouilleur
    Secteur : Industrie

    Informations forums :
    Inscription : Décembre 2008
    Messages : 4 210
    Par défaut
    heu c'est bizarre chez moi j'arrive à aller plus loin dans la boucle :
    D:\Dev\PycharmProjects\pythonProject\venv\Scripts\python.exe D:/Dev/PycharmProjects/pythonProject/binaryread/testC++.py
    Lambda = 1425.2289828701896 ; Rau = 0.6219328089590949 ; Gamma = 0.008881372761047751 ; Dtau = 0.009846019801588628
    I= 200, Csi= -52.78451215631721, Psif= 1.2543724463752129e-23, PsiTotal= 0.008881372761047751, N= 1
    I= 400, Csi= -50.81530819600007, Psif= 8.915267114499358e-23, PsiTotal= 0.008881372761047751, N= 1
    I= 600, Csi= -48.84610423568292, Psif= 6.3847702554760075e-22, PsiTotal= 0.008881372761047751, N= 1
    I= 800, Csi= -46.87690027536578, Psif= 4.574495692355078e-21, PsiTotal= 0.008881372761047751, N= 1
    I= 1000, Csi= -44.90769631504863, Psif= 3.277601068629375e-20, PsiTotal= 0.008881372761047751, N= 1
    I= 1200, Csi= -42.93849235473149, Psif= 2.348390683187784e-19, PsiTotal= 0.008881372761047751, N= 1
    I= 1400, Csi= -40.96928839441434, Psif= 1.6826149583944407e-18, PsiTotal= 0.008881372761047753, N= 1
    I= 1600, Csi= -39.0000844340972, Psif= 1.205588668328873e-17, PsiTotal= 0.008881372761047763, N= 1
    I= 1800, Csi= -37.03088047378005, Psif= 8.638007379958643e-17, PsiTotal= 0.008881372761047838, N= 1
    I= 2000, Csi= -35.061676513462906, Psif= 6.189106905936886e-16, PsiTotal= 0.00888137276104837, N= 1
    I= 2200, Csi= -33.09247255314576, Psif= 4.434476912244294e-15, PsiTotal= 0.008881372761052185, N= 1
    I= 2400, Csi= -31.123268592828296, Psif= 3.177289677331258e-14, PsiTotal= 0.008881372761079524, N= 1
    I= 2600, Csi= -29.15406463251044, Psif= 2.276518717263641e-13, PsiTotal= 0.008881372761275402, N= 1
    I= 2800, Csi= -27.184860672192585, Psif= 1.6311189705812705e-12, PsiTotal= 0.00888137276267887, N= 1
    I= 3000, Csi= -25.21565671187473, Psif= 1.168691947401768e-11, PsiTotal= 0.00888137277273467, N= 1
    I= 3200, Csi= -23.246452751556873, Psif= 8.373643419428793e-11, PsiTotal= 0.008881372844784186, N= 1
    I= 3400, Csi= -21.277248791239018, Psif= 5.999690819107269e-10, PsiTotal= 0.008881373361016832, N= 1
    I= 3600, Csi= -19.308044830921162, Psif= 4.298760663500684e-09, PsiTotal= 0.008881377059808414, N= 1
    I= 3800, Csi= -17.338840870603306, Psif= 3.080048970131414e-08, PsiTotal= 0.008881403561537453, N= 1
    I= 4000, Csi= -15.36963691028545, Psif= 2.20684488234085e-07, PsiTotal= 0.008881593445535985, N= 1
    I= 4200, Csi= -13.400432949967595, Psif= 1.5811938092893269e-06, PsiTotal= 0.00888295395485704, N= 1

  12. #12
    Membre Expert
    Homme Profil pro
    Enseignant
    Inscrit en
    Juin 2013
    Messages
    1 617
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Enseignant
    Secteur : Enseignement

    Informations forums :
    Inscription : Juin 2013
    Messages : 1 617
    Par défaut
    Citation Envoyé par jurassic pork Voir le message
    heu c'est bizarre chez moi j'arrive à aller plus loin dans la boucle :
    Le programme se termine chez toi ?

    Dans tous les cas, il faudrait connaître le résultat attendu...

  13. #13
    Expert confirmé
    Avatar de jurassic pork
    Homme Profil pro
    Bidouilleur
    Inscrit en
    Décembre 2008
    Messages
    4 210
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Bidouilleur
    Secteur : Industrie

    Informations forums :
    Inscription : Décembre 2008
    Messages : 4 210
    Par défaut
    Citation Envoyé par marco056 Voir le message
    Le programme se termine chez toi ?

    Dans tous les cas, il faudrait connaître le résultat attendu...
    ça coince à I=4600

  14. #14
    Membre Expert
    Homme Profil pro
    Enseignant
    Inscrit en
    Juin 2013
    Messages
    1 617
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations professionnelles :
    Activité : Enseignant
    Secteur : Enseignement

    Informations forums :
    Inscription : Juin 2013
    Messages : 1 617
    Par défaut
    Je n'y connais rien en C++ mais a priori :
    - erreur de signe à la ligne 85
    - modifier les boucles par exemple :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    for i in range(1,Nmax+1):
    virer la condition si i = 0 ensuite
    - es-tu sûr de cela ? Qu'est-ce que cela veut dire pour toi :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    for i in range(3*Nmax):
        Psif.append(np.longdouble(0e0))
        PsiTotal.append(np.longdouble(0e0))
        Rqd.append(np.longdouble(0e0))
        Pota.append(np.longdouble(0e0))

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. Boucle while infinie
    Par marcilles dans le forum Général JavaScript
    Réponses: 9
    Dernier message: 27/12/2018, 09h49
  2. Boucle while infini avec SESSION
    Par ZeeKinio dans le forum Langage
    Réponses: 8
    Dernier message: 29/01/2018, 22h28
  3. Boucle while infinie
    Par BlackoOSX dans le forum Arduino
    Réponses: 5
    Dernier message: 20/06/2016, 21h24
  4. boucle while infini!
    Par guizaniseifislam dans le forum Langage
    Réponses: 2
    Dernier message: 02/03/2012, 01h30
  5. Boucle while infinie
    Par god_enel dans le forum C
    Réponses: 1
    Dernier message: 16/04/2007, 11h35

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