Comptes Rendus
Équations aux dérivées partielles
Version continue de l'algorithme d'Uzawa
Comptes Rendus. Mathématique, Volume 337 (2003) no. 1, pp. 31-36.

Nous avons proposé dans Carlier et al. (ESAIM Proceedings, CEMRACS 1999) un algorithme permettant d'approximer la projection d'une fonction fH 0 1 (Ω) (où Ω est un domaine convexe) sur le cône des fonctions convexes. Cet algorithme est basé sur une expression duale de la contrainte de convexité, qui conduit à un problème de point-selle qui n'a pas de solution en général. Nous montrons ici que l'algorithme d'Uzawa appliqué à cette situation peut être vu comme une discrétisation semi-implicite d'une équation d'évolution du type

dλ dt+Ψ(λ)0,
Ψ est une fonction convexe, propre, et s.c.i. Dans le cas où le problème de point-selle n'admet pas de solution, on a
0R(Ψ) ¯ mais Ψ -1 (0)=.
Nous établissons que λ(t) diverge alors, mais qu'une sous-suite de la composante primale de la trajectoire converge faiblement vers la solution du problème de projection initial.

In Carlier et al. (ESAIM Proceedings, CEMRACS 1999), an algorithm was proposed to approximate the projection of a function fH 0 1 (Ω) (where Ω is a convex domain) onto the cone of convex functions. This algorithm is based on a dual expression of the constraint, which leads to a saddle-point problem which has no solution in general. We show here that the Uzawa algorithm for this saddle-point problem can be seen as the semi-discretization of an evolution equation

dλ dt+Ψ(λ)0,
where Ψ is a convex, l.s.c., proper function. In case the saddle-point problem has no solution, one has 0R(Ψ) ¯ but ∂Ψ−1(0)=∅. We establish that λ(t) is then divergent, and that a subsequence of the associated trajectory in the primal space converges weakly to the solution of the initial projection problem.

Accepté le :
Publié le :
DOI : 10.1016/S1631-073X(03)00267-X

Bertrand Maury 1

1 Laboratoire Jacques-Louis Lions, Université Pierre et Marie Curie, boîte courrier 187, 75252 Paris cedex 05, France
@article{CRMATH_2003__337_1_31_0,
     author = {Bertrand Maury},
     title = {Version continue de l'algorithme {d'Uzawa}},
     journal = {Comptes Rendus. Math\'ematique},
     pages = {31--36},
     publisher = {Elsevier},
     volume = {337},
     number = {1},
     year = {2003},
     doi = {10.1016/S1631-073X(03)00267-X},
     language = {fr},
}
TY  - JOUR
AU  - Bertrand Maury
TI  - Version continue de l'algorithme d'Uzawa
JO  - Comptes Rendus. Mathématique
PY  - 2003
SP  - 31
EP  - 36
VL  - 337
IS  - 1
PB  - Elsevier
DO  - 10.1016/S1631-073X(03)00267-X
LA  - fr
ID  - CRMATH_2003__337_1_31_0
ER  - 
%0 Journal Article
%A Bertrand Maury
%T Version continue de l'algorithme d'Uzawa
%J Comptes Rendus. Mathématique
%D 2003
%P 31-36
%V 337
%N 1
%I Elsevier
%R 10.1016/S1631-073X(03)00267-X
%G fr
%F CRMATH_2003__337_1_31_0
Bertrand Maury. Version continue de l'algorithme d'Uzawa. Comptes Rendus. Mathématique, Volume 337 (2003) no. 1, pp. 31-36. doi : 10.1016/S1631-073X(03)00267-X. https://comptes-rendus.academie-sciences.fr/mathematique/articles/10.1016/S1631-073X(03)00267-X/

[1] H. Brezis Opérateurs maximaux monotones et semi–groupes de contraction dans les espaces de Hilbert, North-Holland, 1973

[2] R.E. Bruck Asymptotic convergence of nonlinear contraction semigroups in Hilbert space, J. Funct. Anal., Volume 18 (1975), pp. 15-26

[3] G. Carlier, T. Lachand-Robert, B. Maury, H1-projection into the set of convex functions: a saddle-point formulation, in: ESAIM Proceedings, CEMRACS 1999

[4] P.G. Ciarlet, Introduction à l'Analyse Numérique Matricielle et à l'Optimisation, Masson, Paris

[5] A. Haraux Nonlinear Evolution Equations – Global Behaviour of Solutions, Lecture Notes in Math., Springer, 1981

[6] J.J. Moreau Décomposition orthogonale d'un espace Hilbertien selon deux cônes mutuellement polaires, C. R. Acad. Sci. Paris, Sér. I, Volume 255 (1962), pp. 238-240

[7] J.-C. Rochet; P. Choné Ironing, sweeping and multidimensional screening, Econometrica, Volume 66 (1998), pp. 783-826

Cité par Sources :

Commentaires - Politique