No description available.

Probability-one homotopy methods are a class of methods for solving non-linear systems of equations that are globally convergent with probability one from an arbitrary starting point. The essence of these algorithms is the construction of an appropriate homotopy map and subsequent tracking of some smooth curve in the zero set of the homotopy map. Tracking a homotopy zero curve requires calculating the unit tangent vector at different points along the zero curve. Because of the way a homotopy map is constructed, the unit tangent vector at each point in the zero curve of a homotopy map ρ_α(λ,x) is in the one-dimensional kernel of the full rank n x (n + 1) Jacobian matrix Dρ_α(λ,x). Hence, tracking a zero curve of a homotopy map involves evaluating the Jacobian matrix and finding the one-dimensional kernel of the n x (n + 1) Jacobian matrix with rank n. Since accuracy is important, an orthogonal factorization of the Jacobian matrix is computed. The QR and LQ factorizations are considered here. Computational results are presented showing the performance of several different parallel orthogonal factorization/triangular system solving algorithms on a hypercube, in the context of parallel homotopy algorithms for problems with small, dense Jacobian matrices. This study also examines the effect of different component complexity distributions and the size of the Jacobian matrix on the different assignments of components to the processors, and determines in what context one assignment would perform better than others. / Ph. D.

The problem of finding a reduced order model, optimal in the H² sense, to a given system model is a fundamental one in control system analysis and design. The addition of a H^∞ constraint to the H² optimal model reduction problem results in a more practical yet computationally more difficult problem. Without the global convergence of homotopy methods, both the H² optimal and the combined H²/H^∞ model reduction problems are very difficult. For both problems homotopy algorithms based on several formulations—input normal form; Ly, Bryson, and Cannon's 2 X 2 block parametrization; a new nonminimal parametrization—are developed and compared here. For the H² optimal model order reduction problem, these numerical algorithms are also compared with that based on Hyland and Bernstein's optimal projection equations. Both the input normal form and Ly form are very efficient compared to the over parametrization formulation and the optimal projection equations approach, since they utilize the minimal number of possible degrees of freedom. However, they can fail to exist or be very ill conditioned. The conditions under which the input normal form and the Ly form become ill conditioned are examined. The over-parametrization formulation solves the ill conditioning issue, and usually is more efficient than the approach based on solving the optimal projection equations for the H² optimal model reduction problem. However, the over-parametrization formulation introduces a very high order singularity at the solution, and it is doubtful whether this singularity can be overcome by using interpolation or other existing methods. Although there are numerous algorithms for solving Riccati equations, there still remains a need for algorithms which can operate efficiently on large problems and on parallel machines and which can be generalized easily to solve variants of Riccati equations. This thesis gives a new homotopy-based algorithm for solving Riccati equations on a shared memory parallel computer. The central part of the algorithm is the computation of the kernel of the Jacobian matrix, which is essential for the corrector iterations along the homotopy zero curve. Using a Schur decomposition the tensor product structure of various matrices can be efficiently exploited. The algorithm allows for efficient parallelization on shared memory machines. The linear-quadratic-Gaussian (LQG) theory has engendered a systematic approach to synthesize high performance controllers for nominal models of complex, multi-input multioutput systems and hence it is a breakthrough in modern control theory. Homotopy algorithms for both full and reduced-order LQG controller design problems with an H^∞ constraint on disturbance attenuation are developed. The H^∞ constraint is enforced by replacing the covariance Lyapunov equation by a Riccati equation whose solution gives an upper bound on H² performance. The numerical algorithm, based on homotopy theory, solves the necessary conditions for a minimum of the upper bound on H² performance. The algorithms are based on two minimal parameter formulations: Ly, Bryson, and Cannon's 2 X 2 block parametrization and the input normal Riccati form parametrization. An over-parametrization formulation is also proposed. Numerical experiments suggest that the combination of a globally convergent homotopy method with a minimal parameter formulation applied to the upper bound minimization gives excellent results for mixed-norm synthesis. / Ph. D.

Globally convergent, probability-one homotopy methods have proven to be very effective for finding all the isolated solutions to polynomial systems of equations. After many years of development, homotopy path trackers based on probability-one homotopy methods are reliable and fast. Now, theoretical advances reducing the number of homotopy paths that must be tracked, and in the handling of singular solutions, have made probability-one homotopy methods even more practical. This thesis describes the theory behind and performance of the new code POLSYS_PLP, which consists of Fortran 90 modules for finding all isolated solutions of a complex coefficient polynomial system of equations by a probability-one homotopy method. The package is intended to be used in conjunction with HOMPACK90, and makes extensive use of Fortran 90 derived data types to support a partitioned linear product (PLP) polynomial system structure. PLP structure is a generalization of m-homogeneous structure, whereby each component of the system can have a different m-homogeneous structure. POLSYS_PLP employs a sophisticated power series end game for handling singular solutions, and provides support for problem definition both at a high level and via hand-crafted code. Different PLP structures and their corresponding Bezout numbers can be systematically explored before committing to root finding. / Master of Science

Probability-one homotopy algorithms have strong convergence characteristics under mild assumptions. Such algorithms for mixed complementarity problems (MCPs) have potentially wide impact because MCPs are pervasive in science and engineering. A probability-one homotopy algorithm for MCPs was developed earlier by Billups and Watson based on the default homotopy mapping. This algorithm had guaranteed global convergence under some mild conditions, and was able to solve most of the MCPs from the MCPLIB test library. This thesis extends that work by presenting some other homotopy mappings, enabling the solution of all the remaining problems from MCPLIB. The homotopy maps employed are the Newton homotopy and homotopy parameter embeddings. / Master of Science

There are algorithms for finding zeros or fixed points of nonlinear systems of (algebraic) equations that are globally convergent for almost all starting points, i.e., with probability one. The essence of all such algorithms is the construction of an appropriate homotopy map and then tracking some smooth curve in the zero set of this homotopy map. The augmented Jacobian matrix algorithm is part of the software package HOMPACK, and is based on an algorithm developed by W.C. Rheinboldt. The algorithm exists in two forms-one for dense Jacobian matrices, and the other for sparse Jacobian matrices. / M.S.

Yes / This article presents the homotopy perturbation method (HPM) employed to investigate the effects of inclination on the thermal behavior of a porous fin heat sink. The study aims to review the thermal characterization of heat sink with the inclined porous fin of rectangular geometry. The study establishes that heat sink of an inclined porous fin shows a higher thermal performance compared to a heat sink of equal dimension with a vertical porous fin. In addition, the study also shows that the performance of inclined or tilted fin increases with decrease in length–thickness aspect ratio. The study further reveals that increase in the internal heat generation variable decreases the fin temperature gradient, which invariably decreases the heat transfer of the fin. The obtained results using HPM highlights the accuracy of the present method for the analysis of nonlinear heat transfer problems, as it agrees well with the established results of Runge–Kutta. / Supported in part by the Tertiary Education Trust Fund of Federal Government of Nigeria, and the European Union’s Horizon 2020 research and innovation programme under grant agreement H2020-MSCA-ITN-2016SECRET-722424.

Fizemos uma apresentação detalhada, com um enfoque geométrico, da Teoria de Homotopia Simples e como aplicação, uma análise detalhada da classificação por homotopia e homotopia simples dos Espaços Lenticulares / We made a detailed presentation, with a geometric approach, of Simple Homotopy Theory and as a major application we present a detailed analysis of homotopy and simple homotopy classification of Lens Spaces

Fizemos uma apresentação detalhada, com um enfoque geométrico, da Teoria de Homotopia Simples e como aplicação, uma análise detalhada da classificação por homotopia e homotopia simples dos Espaços Lenticulares / We made a detailed presentation, with a geometric approach, of Simple Homotopy Theory and as a major application we present a detailed analysis of homotopy and simple homotopy classification of Lens Spaces

Masters of Science / The main purpose of this thesis is to study the interplay between relational structures and topology, and to portray pairings in terms of some finite poset models and order preserving maps. We show the interrelations between the categories of topological spaces, closure spaces and relational structures. We study the 4-point non-Hausdorff model S4 weakly homotopy equivalent to the circle s'. We study pairings of some objects in the category of relational structures, similar to the multiplication of Hopf spaces in topology. The multiplication S4 x S4 ---7 S4 fails to be order preserving for posets. Nevertheless, applying a single barycentric subdivision on S4 to get Ss, an 8-point model of the circle enables us to define an order preserving poset map Ss x Ss ---7 S4' Restricted to the axes, this map yields weak homotopy equivalences Ss ---7 S4' Hence it is a pairing. Further, using the non-Hausdorff join Ss ® Ss, we obtain a version of the Hopf map Ss ® Ss ---7 §S4. This model of the Hopf map is in fact a map of non-Hausdorff double mapping cylinders.