Initialisation

[jbcaillau]

@alesiagr

• initialisation vectorielle (pour état et / ou contrôle), pas matricielle mais vecteur de vecteur, je pense
• oui, on peut donner (pour état et / ou contrôle) soit des fonctions du temps, soit des vecteurs (version discrétisée)
• fonction du temps donc $u(t)$, en effet, pas $u(x)$ (mais si tu connais un u(x) et un x_init, il suffit de passer u_init = u.(x_init))
• le tuto pour initialisations

[jbcaillau]

@alesiagr dans ton pb, tu ne déclares pas de variable (= variable d'optimisation, par ex. temps final libre...) alors que tu en passes une pour l'init : de quoi s'agit-il ?

https://github.com/mctao-inria/sail-therm/blob/c44976797e6ef084b2f631dc37695abfbf0a92bb/interstellar_sail.jl#L187

[alesiagr]

@alesiagr dans ton pb, tu ne déclares pas de variable (= variable d'optimisation, par ex. temps final libre...) alors que tu en passes une pour l'init : de quoi s'agit-il ?

https://github.com/mctao-inria/sail-therm/blob/c44976797e6ef084b2f631dc37695abfbf0a92bb/interstellar_sail.jl#L187

Il s'agit d'une erreur, je voulais le supprimer, mais pas encore fait car le code avant ne marche pas. Je supprime!

[jbcaillau]

@alesiagr

j'ai dit du 💩 :

* initialisation vectorielle (pour état et / ou contrôle), pas matricielle mais vecteur de vecteur, je pense

pour l'instant : soit init constante, soit init par fonction. pas de vecteur (de vecteurs). voir https://github.com/control-toolbox/CTDirect.jl/issues/112 par contre, assez facile de générer une fonction (du temps) à partir des valeurs (état / contrôle) samplées sur une grille (ce que tu fournis) ; par interpolation, tu vois comment faire ? sachant que le mieux est d'essayer avec des valeurs constantes (sur chaque composante) au cas où ça marcherait (pas clair, en effet)

[alesiagr]

@alesiagr

j'ai dit du 💩 :

~* initialisation vectorielle (pour état et / ou contrôle), pas matricielle mais vecteur de vecteur, je pense~

pour l'instant : soit init constante, soit init par fonction. pas de vecteur (de vecteurs). voir control-toolbox/CTDirect.jl#112 par contre, assez facile de générer une fonction (du temps) à partir des valeurs (état / contrôle) samplées sur une grille (ce que tu fournis) ; par interpolation, tu vois comment faire ? sachant que le mieux est d'essayer avec des valeurs constantes (sur chaque composante) au cas où ça marcherait (pas clair, en effet)

Oui, j'ai voulu faire ça hier encore, je parle de l'interpolation, mais je me suis dit que je vais attendre ta réponse. Bon bah du coup je vais le faire :D

[jbcaillau]

@alesiagr Voir Interpolations.jl. Une interpolation linéaire (composante à composante pour l'état, qui est vectoriel) fait le job.

Oui, j'ai voulu faire ça hier encore, je parle de l'interpolation, mais je me suis dit que je vais attendre ta réponse. Bon bah du coup je vais le faire :D

[jbcaillau]

@alesiagr pour synchro Github, tu es bien membre du repo avec les droits ad hoc, donc ça doit marcher. si pas déjà fait, installe le client gh : https://cli.github.com

et re-clone le repo avec le lien web (https plutôt que ssh). Depuis vscode, il suffit d'ailleurs que tu te connectes à ton compte GitHub, e basta.

[jbcaillau]

@alesiagr je regarde le code, je vois pas de push aujourd'hui : c'est normal ou tu as toujours le pb d'accès ?

[alesiagr]

@jbcaillau je viens de faire un push parce que j'arrive pas à le faire fonctionner. J'obtiens des résultats un peu étranges sans initial guess, j'essaie de voir comment les visualiser en regardant https://github.com/control-toolbox/CTBase.jl/blob/b73f56d6447e1967bf56c59f49dc8bcdb09dd71c/src/plot.jl

[jbcaillau]

OK : je pulle et regarde

[jbcaillau]

@alesiagr dans la branche develop :

• quelques corrections de style 🙂 (j'espère n'avoir rien cassé)
• j'ai vu après coup que tu avais extrait plusieurs fonctions auxiliaires (et fait des include), je vais mettre à jour
• surtout : l'opti ne converge pas (erreurs numériques...) À suivre

[alesiagr]

J'ai vu, merci! Je pense avoir trouvé le pb de git, je viens de faire le request pour que tu m'autorises s'il te plaît :)

[alesiagr]

J'ai réussi à faire marcher GitHub après 4 jours! Il fallait faire un personal access token apparemment.... 🥲 Enfin ça marche... Je regarde le code pour comprendre pourquoi il veut pas converger.

C'est possible que j'ai écrasé un de tes push @jbcaillau , désolée pour ça. Mais c'est parce que jusqu'à présent je ne pouvais faire commit qu'en uploadant le fichier dans le browser.

[alesiagr]

Salut @jbcaillau, il y a qqch que je n'arrive pas à comprendre en Julia, c'est la manière de créer des vecteurs/ matrices. Ma fonction F0 ne fonctionne pas quand j'essaie de la propager avec ODEProblem (voir fichier test.jl), donc j'ai dû la re-écrire sous la forme non-matricielle, ligne par ligne, comme on le faisait avant. Si je fais ça avec OptimalControl.jl, ça fonctionne pas non plus. J'ai aussi l'erreur de ForwardDiff.jl qui apparait lorsque je réécris la fonction ligne par ligne, mais sinon le code ne converge pas tout court...

[jbcaillau]

coucou @alesiagr dispo demain 17:00 pour un point ?

[alesiagr]

hello @jbcaillau , 18h fonctionne pour toi?

[jbcaillau]

hello @jbcaillau , 18h fonctionne pour toi? yes: webex