No abstract uploaded!

No abstract uploaded!

In this paper, a structure-preserving algorithm is developed for the computation of a semi-stabilizing solution of a Generalized Algebraic Riccati Equation (GARE). The semi-stabilizing solution of GAREs has been used to characterize the solvability of the (J, J ′ )-spectral factorization problem in control theory for general rational matrices which may have poles and zeros on the extended imaginary axis. The main difficulty in solving such a GARE lies in the fact that its associated Hamiltonian/skew-Hamiltonian pencil has eigenvalues on the extended imaginary axis. Consequently, it is not clear which eigenspace of the associated Hamiltonian/skew-Hamiltonian pencil can characterize the desired semi-stabilizing solution. That is, it is not clear which eigenvectors and principal vectors corresponding to the eigenvalues on the extended imaginary axis should be contained in the eigenspace that we wish to compute. Hence, the well-known generalized eigenspace approach for the classical algebraic Riccati equations cannot be employed directly. The proposed algorithm consists of a structure-preserving doubling algorithm (SDA) and a postprocessing procedure to determine the desired eigenvectors and principal vectors corresponding to the purely imaginary and infinite eigenvalues. Under mild assumptions, linear convergence of rate 1/2 for the SDA is proved. Numerical experiments illustrate that the proposed algorithm performs efficiently and reliably.

A natural metric in 2-manifold surfaces is to use geodesic distance. If a 2-manifold surface is represented by a triangle mesh T , the geodesic metric on T can be computed exactly using computational geometry methods. Previous work for establishing the geodesic metric on T only supports using half-edge data structures; i.e., each edge e in T is split into two halves (he1, he2) and each half-edge corresponds to one of two faces incident to e. In this paper, we prove that the exact-geodesic structures on two halfedges of e can be merged into one structure associated with e. Four merits are achieved based on the properties which are studied in this paper: (1) Existing CAD systems that use edge-based data structures can directly add the geodesic distance function without changing the kernel to a half-edge data structure; (2) To find the geodesic path from inquiry points to the source, the MMP algorithm can be run in an onthe- fly fashion such that the inquiry points are covered by correct wedges; (3) The MMP algorithm is sped up by pruning unnecessary wedges during the wedge propagation process; (4) The storage of the MMP algorithm is reduced since fewer wedges need to be stored in an edge-based data structure. Experimental results show that when compared to the classic half-edge data structure, the edge-based implementation of the MMP algorithm reduces running time by 44% and storage by 29% on average.

In this paper we study the Dirichlet problem for fully nonlinear second-order equations on a riemannian manifold. As in our previous paper [HL4], we define equations via closed subsets of the 2-jet bundle where each equation has a natural dual equation. Basic existence and uniqueness theorems are established in a wide variety of settings. However, the emphasis is on starting with a constant coefficient equation as a model, which then universally determines an equation on every riemannian manifold which is equipped with a topological reduction of the structure group to the invariance group of the model. For example, this covers all branches of the homogeneous complex Monge-Amp`ere equation on an almost complex hermitian manifold X. In general, for an equation F on a manifold X and a smooth domain  X, we give geometric conditions which imply that the Dirichlet problem on is uniquely solvable for all continuous boundary functions. We begin by introducing a weakened form of comparison which has the advantage that local implies global. We then introduce two fundamental concepts. The first is the notion of a monotonicity cone M for F. If X carries a global Msubharmonic function, then weak comparison implies full comparison. The second notion is that of boundary F-convexity, which is defined in terms of the asymptotics of F and is used to define barriers. In combining these notions the Dirichlet problem becomes uniquely solvable as claimed. This article also introduces the notion of local affine jet-equivalence for subequations. It is used in treating the cases above, but gives results for a much broader spectrum of equations on manifolds, including inhomogeneous equations and the Calabi-Yau equation on almost complex hermitian manifolds. A considerable portion of the paper is concerned with specific examples. They include a wide variety of equations which make sense on any riemannian manifold, and many which hold universally on almost complex or quaternionic hermitian manifolds, or topologically calibrated manifolds.