Unsymmetric positive definite linear systems

Is it necessary to pivot when solving an unsymmetric positive definite linear system Ax=b? Define $\text{T=}\text{(A+A}{}^{\mathrm{T}}\text{)}\text{2}$ and $\text{S=}\text{(A±A}{}^{\mathrm{T}}\text{)}\text{2}$. It is shown that pivoting is unnecessary if the quantity $\text{6ST}{}^{\mathrm{?1}}\text{S6}{}_2\text{6A6}{}_2$ is suitably small with respect to the working machine precision.     

Subproper splitting for rectangular matrices

In this paper iterative schemes for approximating a solution to a rectangular but consistent linear system Ax = b are studied. Let A?C^{m × n}_r. The splitting A = M ? N is called subproper if R(A) ? R(M) and $\text{R(A}{}^1\text{) ?R(M}{}^1\text{)}$. Consider the iteration $\text{x}{}_{\mathrm{i}}\text{= M}{}^2\text{Nx}{}_{\mathrm{i}}{}_{\mathrm{?1}}\text{+ M}{}^2\text{b}$. We characterize the convergence of this scheme to a solution of the linear system. When A?R^m×ⁿ_r, monotonicity and the concept of subproper regular splitting are used to determine a necessary and a sufficient condition for the scheme to converge to a solution.     

Lower bounds for the norm of the symmetric product

Suppose each of m, n, and k is a positive integer, k ? n, A is a (real-valued) symmetric n-linear function on E_m, and B is a k-linear symmetric function on E_m. The tensor and symmetric products of A and B are denoted, respectively, by A ?B and A?B. The identity

$\text{6A · B6}{}^2\text{=}\text{∑}\text{q=0}\text{n}\text{(}\text{n}\text{k}\text{)}\text{(}\text{n+k}\text{k}\text{)}\text{6A?}{}_{\mathrm{q}}\text{B6}{}^2$

is proven by Neuberger in [1]. An immediate consequence of this identity is the inequality

$\text{6A · B 6}{}^2\text{?}\text{n+k}\text{n}{}^{\mathrm{?1}}\text{6A · B 6}{}^2$

In this paper a necessary and sufficient condition for

$\text{6A · B 6}{}^2\text{=}\text{n+k}\text{n}{}^{\mathrm{?}}\text{6A · B 6}{}^2$

is given. It is also shown that under certain conditions the inequality can be considerably improved. This improvement results from an analysis of the terms 6A?_qB6, 1?q?n, appearing in the identity.     

Transcendental continued fractions

In a previous paper it was proven that given the continued fractions

$\text{A = a}{}_1\text{+}\text{1}\text{a}{}_2\text{+}\text{1}\text{a}{}_3\text{+… and B = b}{}_1\text{+}\text{1}\text{b}{}_2\text{+}\text{1}\text{b}{}_3\text{+…}$

where the a's and b's are positive integers, then A, B, A ± B, $\text{A}\text{B}$ and AB are irrational numbers if $\text{a}{}_{\mathrm{n}}\text{2}\text{> b}{}_{\mathrm{n}}\text{> a}{}_{\mathrm{n?1}}{}^{\mathrm{5n}}$ for all n sufficiently large, and transcendental numbers if for all n sufficiently large. Using a more direct approach it is proven in this paper that A, B, A ± B, $\text{A}\text{B}$ and AB are transcendental numbers if a_n > b_n > a_n?1^(n?1)2 for all n sufficiently large.     

Estimées Ck pour l'équation ∂̄b sur un convexe de type fini de Cn

Let q=1,…,n?1 and D be a bounded convex domain in $\text{C}{}^{\mathrm{n}}$ of finite type m. We construct two integral operators T_q and $\text{T}{}_{\mathrm{q}}$ such that for all $\text{p∈}\text{N}\text{,}\text{T}{}_{\mathrm{q}}\text{,}\text{T}{}_{\mathrm{q}}\text{:C}{}^{\mathrm{p}}{}_{\mathrm{0,q}}\text{(bD)→C}{}^{\mathrm{p+1/m}}{}_{\mathrm{0,q?1}}\text{(bD)}$ are continuous, and for all (0,q)-forms h continuous on bD with $\text{?}{}_{\mathrm{b}}\text{h}$ continuous on bD too, with the additional hypothesis when q=n?1 that ∫_bDh∧φ=0 for all φ∈C^∞_n,0(bD) $\text{?}\text{?}{}_{\mathrm{b}}$-fermée, we show $\text{h=}\text{?}\text{?}{}_{\mathrm{b}}\text{(T}{}_{\mathrm{q}}\text{?}\text{T}{}_{\mathrm{q}}\text{)h+(T}{}_{\mathrm{q+1}}\text{?}\text{T}{}_{\mathrm{q+1}}\text{)}\text{?}\text{?}{}_{\mathrm{b}}\text{h}$. For this construction, we use the Diederich–Fornæss support function of Alexandre (Publ. IRMA Lille 54 (III) (2001)). To prove the continuity of T_q, we integrate by parts and take care of the tangential derivatives. The normal component in z of the kernel of $\text{T}{}_{\mathrm{q}}$ will have a bad behaviour, so, in order to find a good representative of its equivalence class, we isolate the tangential component of the kernel and then integrate by parts again. To cite this article: W. Alexandre, C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

Anti-Hadamard matrices

An anti-Hadamard matrix may be loosely defined as a real (0, 1) matrix which is invertible, but only just. Let A be an invertible (0, 1) matrix with eigenvalues λ_i, singular values σ_i, and inverse B = (b_ij). We are interested in the four closely related problems of finding λ(n) = min_{A, i}|λ_i|, σ(n) = min_{A, i}σ_i, χ(n) = max_{A, i, j} |b_ij|, and μ(n) = max_AΣ_ijb²_ij. Then A is an anti-Hadamard matrix if it attains μ(n). We show that λ(n), σ(n) are between $\text{(2n)}{}^{\mathrm{?1}}\text{(}\text{n}\text{4}\text{)}{}^{\mathrm{<mtext>?n</mtext><mtext>2</mtext>}}$ and c√n (2.274)^?n, where c is a constant, $\text{c}\text{(2.274)}{}^{\mathrm{n}}\text{?χ(n)?2(}\text{n}\text{4}\text{)}{}^{\mathrm{<mtext>n</mtext><mtext>2</mtext>}}$, and $\text{c}\text{(5.172)}{}^{\mathrm{n}}\text{?μ(n)?4n}{}^2\text{(}\text{n}\text{4}\text{)}{}^{\mathrm{n}}$. We also consider these problems when A is restricted to be a Toeplitz, triangular, circulant, or (+1, ?1) matrix. Besides the obvious application—to finding the most ill-conditioned (0, 1) matrices—there are connections with weighing designs, number theory, and geometry.     

Radial bounds for perturbations of elliptic operators

Elliptic operators $\text{A = ∑}{}_{\mathrm{¦\alpha ¦\; ?\; m}}\text{b}{}_{\mathrm{\alpha }}\text{(x) D}{}^{\mathrm{\alpha }}$, α a multi-index, with leading term positive and constant coefficient, and with lower order coefficients defined on $\text{R}{}^{\mathrm{<mi\; mathvariant="italic"\; is="true">n</mi>}}$ or a quotient space $\text{R}{}^{\mathrm{n}}\text{R}{}^{\mathrm{n}}\text{U}{}_{\mathrm{\alpha }}\text{, U}{}_{\mathrm{\alpha }}\text{? R}{}^{\mathrm{n}}$ are considered. It is shown that the L^p-spectrum of A is contained in a “parabolic region” Ω of the complex plane enclosing the positive real axis, uniformly in p. Outside Ω, the kernel of the resolvent of A is shown to be uniformly bounded by an L¹ radial convolution kernel. Some consequences are: A can be closed in all L^p (1 ? p ? ∞), and is essentially self-adjoint in L² if it is symmetric; A generates an analytic semigroup e^?tA in the right half plane, strongly L^p and pointwise continuous at t = 0. A priori estimates relating the leading term and remainder are obtained, and summability $\text{φ(εA)?→}{}_{\mathrm{\epsilon \; \to \; 0}}\text{φ(0) ?}$, with φ analytic, is proved for $\text{? ? L}{}^{\mathrm{p}}$, with convergence in L^p and on the Lebesgue set of ?. More comprehensive summability results are obtained when A has constant coefficients.     

Interpolation and inner maps that preserve measure

It is shown, for n ? m ? 1, that there exist inner maps Φ: Bⁿ → B^m with boundary values $\text{Φ}{}_1\text{: B}{}^{\mathrm{n}}\text{→ B}{}^{\mathrm{m}}$ such that $\text{σ}{}_{\mathrm{m}}\text{(A) = σ}{}_{\mathrm{n}}\text{(Φ}{}_1{}^{\mathrm{?1}}\text{(A))}$. where σ_n and σ_m are the Haar measures on ?Bⁿ and ?B^m, respectively, and A ? Bⁿ is an arbitrary Borel set.     

An inequality involving permanents of certain direct products

Let A denote a decomposable symmetric complex valued n-linear function on C^m. We prove

$\text{6A·A6}{}^2\text{?2}{}^{\mathrm{n}}\text{2n}\text{n}{}^{\mathrm{?1}}\text{6A?A6}{}^2$

, where · denotes the symmetric product and ? the tensor product. As a consequence we have per

$\text{M}\text{M}\text{M}\text{M}\text{?2}{}^{\mathrm{n}}\text{[per(M)]}{}^2$

, where M is a positive semidefinite Hermitian matrix and per denotes the permanent function. A sufficient condition for equality in the matrix inequality is that M is a nonnegative diagonal matrix.     

On the fair coin tossing process

Let Ω = {1, 0} and for each integer n ≥ 1 let $\text{Ω}{}_{\mathrm{n}}\text{= Ω × Ω × … × Ω}$ (n-tuple) and $\text{Ω}{}_{\mathrm{n}}{}^{\mathrm{k}}\text{= {(a}{}_1\text{, a}{}_2\text{, …, a}{}_{\mathrm{n}}\text{)|(a}{}_1\text{, a}{}_2\text{, … , a}{}_{\mathrm{n}}\text{) ? Ω}{}_{\mathrm{n}}\text{and}\text{Σ}{}_{\mathrm{i=1}}{}^{\mathrm{n}}\text{a}{}_{\mathrm{i}}\text{= k}}$ for all k = 0,1,…,n. Let {Y_m}_m≥1 be a sequence of i.i.d. random variables such that $\text{P(Y}{}_1\text{= 0) = P(Y}{}_1\text{= 1) =}\text{1}\text{2}$. For each A in $\text{Ω}{}_{\mathrm{n}}$, let T_A be the first occurrence time of A with respect to the stochastic process {Y_m}_m≥1. R. Chen and A.Zame (1979, J. Multivariate Anal. 9, 150–157) prove that if n ≥ 3, then for each element A in $\text{Ω}{}_{\mathrm{n}}$, there is an element B in $\text{Ω}{}_{\mathrm{n}}$ such that the probability that T_B is less than T_A is greater than $\text{1}\text{2}$. This result is sharpened as follows: (I) for n ≥ 4 and 1 ≤ k ≤ n ? 1, each element A in $\text{Ω}{}_{\mathrm{n}}{}^{\mathrm{k}}$, there is an element B also in $\text{Ω}{}_{\mathrm{n}}{}^{\mathrm{k}}$ such that the probability that T_B is less than T_A is greater than $\text{1}\text{2}$; (II) for n ≥ 4 and 1 ≤ k ≤ n ? 1, each element A = (a₁, a₂,…,a_n) in $\text{Ω}{}_{\mathrm{n}}{}^{\mathrm{k}}$, there is an element C also in $\text{Ω}{}_{\mathrm{n}}{}^{\mathrm{k}}$ such that the probability that T_A is less than T_C is greater than $\text{1}\text{2}$ if n ≠ 2m or n = 2m but a_i = a_{i + 1} for some 1 ≤ i ≤ n?1. These new results provide us with a better and deeper understanding of the fair coin tossing process.     

A generalization of the fast LUP matrix decomposition algorithm and applications

We show that any m × n matrix A, over any field, can be written as a product, LSP, of three matrices, where L is a lower triangular matrix with l's on the main diagonal, S is an m × n matrix which reduces to an upper triangular matrix with nonzero diagonal elements when the zero rows are deleted, and P is an n × n permutation matrix. Moreover, L, S, and P can be found in O(m^α?1n) time, where the complexity of matrix multiplication is O(m^α). We use the LSP decomposition to construct fast algorithms for some important matrix problems. In particular, we develop O(m^α?1n) algorithms for the following problems, where A is any m × n matrix: (1) Determine if the system of equations $\text{A}\text{x}\text{=}\text{b}$ (where $\text{b}$ is a column vector) has a solution, and if so, find one such solution. (2) Find a generalized inverse, $\text{A}{}^1$, of A (i.e., $\text{AA}{}^1\text{A = A}$). (3) Find simultaneously a maximal independent set of rows and a maximal independent set of columns of A.     

A generalization of doubly stochastic matrices over a field

Let X and Y be m×n matrices over a field F such that Y^TX is nonsingular, and let Λ and Λ′ be sets of n-square matrices over F. Solutions A to the simultaneous equations AX = XK and $\text{Y}{}^{\mathrm{T}}\text{A =}\text{K}\text{?}\text{Y}{}^{\mathrm{T}}$ where K?Λ and $\text{K}\text{?}\text{? Λ′}$ are considered. It is shown that many properties of doubly stochastic matrices over a field have a natural generalization in terms of the set Δ(Λ,Λ′) of all such solutions.     

Continuity of positive functionals on topological 1-algebras

It is shown that the set $\text{C}$_{m × n} of complex m × n matrices forms a lower semilattice under the partial ordering A ? B defined by $\text{A}{}^1\text{A = A}{}^1\text{B,}{}^1\text{AA}{}^1\text{= BA}{}^1\text{,}\text{where}\text{A}{}^1$ denotes the conjugate transpose of A. As a special case of a result for division rings, it is further shown that, over any field F, form = n = 2 and any proper involution 1 of F_{2 × 2}, the corresponding intersections A ∩ B all exist.     

On the convergence of the Neumann series in interval analysis

Let $\text{A}$ be a real or complex n × n interval matrix. Then it is shown that the Neumann series $\text{Σ}{}^{\mathrm{\infty }}{}_{\mathrm{k=0}}\text{A}{}^{\mathrm{k}}$ is convergent iff the sequence {$\text{A}$^k} converges to the null matrix $\text{O}$, i.e., iff the spectral radius of the real comparison matrix $\text{B}$ constructed in [2] is less than one.     

A family of quadratic forms associated to quadratic mappings of spheres

The general form of a real quadratic mapping of spheres can be determined by studying the diagonalization of each form in an associated family of quadratic forms. In particular, the eigenvalues provide a means for detecting maps which are of the Hopf type. When the eigenvalues are nonzero for every form in the family, the forms associated to $\text{?:}\text{S}{}^{\mathrm{n}}\text{→}\text{S}{}^{\mathrm{m}}$ give rise to a quadratic form on the tangent bundle of the unit sphere Sⁿ. If ? is of the Hopf type, nondegeneracy of each form occurs only when n=1,3,7,15.     

Extension de quelques theoremes sur les densites de series d'elements de n a des series de sous-ensembles finis de n

For a sequence A = {A_k} of finite subsets of N we introduce: $\text{δ(A) =}\text{inf}{}_{\mathrm{m?n}}\text{A(m)}\text{2}{}^{\mathrm{n}}$, $\text{d(A) =}\text{lim inf}{}_{\mathrm{n\to \infty }}\text{A(n)}\text{2}{}^{\mathrm{n}}$, where A(m) is the number of subsets A_k ? {1, 2, …, m}.The collection of all subsets of {1, …, n} together with the operation $\text{a ∪ b, (a ∩ b), (a 1 b = a ∪ b ? a ∩ b)}$ constitutes a finite semi-group N^∪ (semi-group N^∩) (group $\text{N}{}^1$). For N^∪, N^∩ we prove analogues of the Erdös-Landau theorem: δ(A+B) ? δ(A)(1+(2λ)^?1(1?δ(A>))), where B is a base of N of the average order λ. We prove for $\text{N}{}^{\mathrm{\cup }}\text{, N}{}^{\mathrm{\cap }}\text{, N}{}^1$ analogues of Schnirelmann's theorem (that δ(A) + δ(B) > 1 implies δ(A + B) = 1) and the inequalities λ ? 2h, where h is the order of the base.We introduce the concept of divisibility of subsets: a|b if b is a continuation of a. We prove an analog of the Davenport-Erdös theorem: if d(A) > 0, then there exists an infinite sequence {A_kr}, where A_kr | A_kr+1 for r = 1, 2, …. In Section 6 we consider for $\text{N∪, N∩, N1}$ analogues of Rohrbach inequality: $\text{2n}\text{? g(n) ? 2}\text{n}$, where g(n) = min k over the subsets {a₁ < … < a_k} ? {0, 1, 2, …, n}, such that every m? {0, 1, 2, …, n} can be expressed as m = a_i + a_j.Pour une série A = {A_k} de sous-ensembles finis de N on introduit les densités: $\text{δ(A) =}\text{inf}{}_{\mathrm{m?n}}\text{A(m)}\text{2}{}^{\mathrm{m}}$, $\text{d(A) =}\text{lim inf}{}_{\mathrm{n\to \infty }}\text{A(n)}\text{2}{}^{\mathrm{n}}$ où A(m) est le nombre d'ensembles A_k ? {1, 2, …, m}. L'ensemble de toutes les parties de {1, 2, …, n} devient, pour les opérations $\text{a ∪ b, a ∩ b, a 1 b = a ∪ b ? a ∩ b}$, un semi-groupe fini N^∪, N^∩ ou un groupe N¹ respectivement. Pour N^∪, N^∩ on démontre l'analogue du théorème de Erdös-Landau: δ(A + B) ? δ(A)(1 + (2λ)^?1(1?δ(A))), où B est une base de N d'ordre moyen λ. On démontre pour $\text{N}{}^{\mathrm{\cup }}\text{, N}{}^{\mathrm{\cap }}\text{, N}{}^1$ l'analogue du théorème de Schnirelmann (si δ(A) + δ(B) > 1, alors δ(A + B) = 1) et les inégalités λ ? 2h, où h est l'ordre de base. On introduit le rapport de divisibilité des enembles: a|b, si b est une continuation de a. On démontre l'analogue du théorème de Davenport-Erdös: si d(A) > 0, alors il existe une sous-série infinie {A_kr}, où A_kr|A_kr+1, pour r = 1, 2, … . Dans le Paragraphe 6 on envisage pour N^∪, $\text{N}{}^{\mathrm{\cap }}\text{, N}{}^1$ les analogues de l'inégalité de Rohrbach: $\text{2n}\text{? g(n) ? 2}\text{n}$, où g(n) = min k pour les ensembles {a₁ < … < a_k} ? {0, 1, 2, …, n} tels que pour tout m? {0, 1, 2, …, n} on a m = a_i + a_j.     

On the Schur multiplier norm of matrices

The Schur product of two n×n complex matrices A=(a_ij), B=(b_ij) is defined by A°B=(a_ijb_ij. By a result of Schur [2], the algebra of n×n matrices with Schur product and the usual addition is a commutative Banach algebra under the operator norm (the norm of the operator defined on $\text{C}$ⁿ by the matrix). For a fixed matrix A, the norm of the operator B?A°B on this Banach algebra is called the Schur multiplier norm of A, and is denoted by ∥A∥^m. It is proved here that $\text{∥A∥=∥U}{}^1\text{AU∥}{}_{\mathrm{m}}$ for all unitary U (where ∥·∥ denotes the operator norm) iff A is a scalar multiple of a unitary matrix; and that ∥A∥_m=∥A∥ iff there exist two permutations P, Q, a p×p (1?p?n) unitary U, an (n?p)×(n?p)1 contraction C, and a nonnegative number λ such that

$\text{A=λP}\text{U}\text{0}\text{0}\text{C}\text{Q;}$

and this is so iff $\text{∥A°}\text{A}\text{?}\text{∥=∥A∥}{}^2$, where ā is the matrix obtained by taking entrywise conjugates of A.     

Series inversion of some convolution transforms

We study degeneration for ? → + 0 of the two-point boundary value problems

$\text{τ?}{}^{\mathrm{\pm }}\text{u := ?((au′)′ + bu′ + cu) ± xu′ ? κu = h, u(±1) = A ± B}$

, and convergence of the operators T_?⁺ and T_?^? on $\text{L}$²(?1, 1) connected with them, T_?^±u := τ_?^±u for all

$\text{u?D(T}{}_{\mathrm{?}}{}^{\mathrm{\pm }}\text{, D(T}{}_{\mathrm{?}}{}^{\mathrm{\pm }}\text{) := {u ? L}{}^2\text{(?1, 1) ∣ u″ ? L}{}^2\text{(?1, 1) \&; u(?1) = u(1) = O}, T}{}_0{}^{\mathrm{+}}\text{u: = xu′}$

for all

$\text{u?D(T}{}_{\mathrm{O}}{}^{\mathrm{+}}\text{), D(T}{}_{\mathrm{O}}{}^{\mathrm{+}}\text{) := {u ? L}{}^2\text{(?1, 1) ∣ xu′ ? L}{}^2\text{(?1, 1) \&; u(?1) = u(1) = O}}$

. Here ? is a small positive parameter, λ a complex “spectral” parameter; a, b and c are real $\text{b}$^∞-functions, a(x) ? γ > 0 for all x? [?1, 1] and h is a sufficiently smooth complex function. We prove that the limits of the eigenvalues of T_?⁺ and of T_?^? are the negative and nonpositive integers respectively by comparison of the general case to the special case in which a  1 and b  c  0 and in which we can compute the limits exactly. We show that (T_?⁺ ? λ)^?1 converges for ? → +0 strongly to (T₀⁺ ? λ)^?1 if $\text{R e λ > ?}\text{1}\text{2}$. In an analogous way, we define the operator T_?⁺, n (n ? $\text{N}$ in the Sobolev space H₀^?n(? 1, 1) as a restriction of τ_?⁺ and prove strong convergence of (T⁺_?,n ? λ)^?1 for ? → +0 in this space of distributions if $\text{R e λ > ?n ?}\text{1}\text{2}$. With aid of the maximum principle we infer from this that, if h?$\text{C}$¹, the solution of τ_?⁺u ? λu = h, u(±1) = A ± B converges for ? → +0 uniformly on [?1, ? ?] ∪ [?, 1] to the solution of xu′ ? λu = h, u(±1) = A ± B for each p > 0 and for each λ ? $\text{C}$ if ? ?$\text{N}$.Finally we prove by duality that the solution of τ_?^?u ? λu = h converges to a definite solution of the reduced equation uniformly on each compact subset of (?1, 0) ∪ (0, 1) if h is sufficiently smooth and if 1 ? ?$\text{N}$.     

On derivative-dependent infinitesimal deformations of differentiable maps

We study, by means of flows in jet bundles, infinitesimal deformations of germs of $\text{b}$^∞ maps $\text{?:R}{}^{\mathrm{n}}\text{→R}{}^{\mathrm{m}}$ that depend on a finite number of derivatives of these maps. We show that for m = 1 such deformations cannot essentially depend on derivatives higher than the first and that the deformations generalize the notion of an infinitesimal contact transformation. Analogous results hold for m > 1. Some applications to differential equations are given.     

Estimées Ck pour les couronnes de convexes de type fini de Cn

C^k estimates for convex domains of finite type in $\text{C}{}^{\mathrm{n}}$ are known from Alexandre (C. R. Acad. Paris, Ser. I 335 (2002) 23–26). We now want to show the same result for annuli. Precisely, we show that for all convex domains D and D′ relatively compact of $\text{C}{}^{\mathrm{n}}$, of finite type m and m′ such that $\text{D}\text{?D′}$, for all q=1,…,n?2, there exists a linear operator $\text{T}{}^1{}_{\mathrm{q}}$ from $\text{C}{}_{\mathrm{0,q}}\text{(}\text{D′?D}\text{)}$ to $\text{C}{}_{\mathrm{0,q?1}}\text{(}\text{D′?D}\text{)}$ such that for all $\text{k∈}\text{N}$ and all (0,q)-form f, $\text{?}\text{?}$-closed of regularity C^k up to the boundary, $\text{T}{}^1{}_{\mathrm{q}}\text{f}$ is of regularity C^{k+1/max(m,m′)} up to the boundary and $\text{?}\text{?}\text{T}{}_{\mathrm{q}}{}^1\text{f=f}$. We fit the method of Diederich, Fisher and Fornaess to the annuli by switching z and ζ. However, the integration kernel will not have the same behavior on the frontier as in the Diederich–Fischer–Fornaess case and we have to alter the Diederich–Fornaess support function which will not be holomorphic anymore. Also, we take care of the so generated residual term in the homotopy formula and show that it is extremely regular so that solve the $\text{?}\text{?}$ problem for it will not be difficult. To cite this article: W. Alexandre, C. R. Acad. Sci. Paris, Ser. I 336 (2003).