No description available.

PT-symmetric Hamiltonians proposed by Bender and Boettcher can have real energy spectra. As an extension of the Hermitian Hamiltonian, PT-symmetric systems have attracted a great interest in recent years. Understanding the underlying mathematical structure of these theories sheds insight on outstanding problems of physics. These problems include the nature of Higgs particles, the properties of dark matter, the matter-antimatter asymmetry in the universe, and neutrino oscillations. Furthermore, PT-phase transition has been observed in lasers, optical waveguides, microwave cavities, superconducting wires and circuits. The objective of this thesis is to extend the iterative method of solving Schrodinger equation used for an harmonic oscillator systems to Hamiltonians with PT-symmetric potentials. An important aspect of this approach is the high accuracy of eigenvalues and the fast convergence. Our method is a combination of Hill determinant method [8] and the power series expansion. eigenvalues and the fast convergence. One can transform the Schrodinger equation into a secular equation by using a trial wave function. A recursion structure can be obtained using the secular equation, which leads to accurate eigenvalues. Energy values approach to exact ones when the number of iterations is increased. We obtained eigenvalues for a set of PT-symmetric Hamiltonians.

Parallel magnetic resonance imaging (MRI) with non-Cartesian sampling pattern is a promising technique that increases the scan speed using multiple receiver coils with reduced samples. However, reconstruction is challenging due to the increased complexity. Three reconstruction methods were evaluated: gridding, blocked uniform resampling (BURS) and non-uniform FFT (NUFFT). Computer simulations of parallel reconstruction were performed. Root mean square error (RMSE) of the reconstructed images to the simulated phantom were used as image quality criterion. Gridding method showed best RMSE performance. Two type of a priori constraints to reduce noise and artifacts were evaluated: edge preserving penalty, which suppresses noise and aliasing artifact in image while preventing over-smoothness, and object support penalty, which reduces background noise amplification. A trust region based step-ratio method that iteratively calculates the penalty coefficient was proposed for the penalty functions. Two methods to alleviate computation burden were evaluated: smaller over sampling ratio, and interpolation coefficient matrix compression. The performance were individually tested using computer simulations. Edge preserving penalty and object support penalty were shown to have consistent improvement on RMSE. The performance of calculated penalty coefficients on the two penalties were close to the best RMSE. Oversampling ratio as low as 1.125 was shown to have impact of less than one percent on RMSE for the radial sampling pattern reconstruction. The value reduced the three dimensional data requirement to less than 1/5 of what the conventional 2x grid needed. Interpolation matrix compression with compression ratio up to 50 percent showed small impact on RMSE. The proposed method was validated on 25 MR data set from a GE MR scanner. Six image quality metrics were used to evaluate the performance. RMSE, normalized mutual information (NMI) and joint entropy (JE) relative to a reference image from a separate body coil scan were used to verify the fidelity of reconstruction to the reference. Region of interest (ROI) signal to noise ratio (SNR), two-data SNR and background noise were used to validate the quality of the reconstruction. The proposed method showed higher ROI SNR, two-data SNR, and lower background noise over conventional method with comparable RMSE, NMI and JE to the reference image at reduced computer resource requirement.

Fast and accurate solvers for capacitance extraction are needed by the VLSI industry in order to achieve good design quality in feasible time. With the development of technology, this demand is increasing dramatically. Three-dimensional capacitance extraction algorithms are desired due to their high accuracy. However, the present 3D algorithms are slow and thus their application is limited. In this dissertation, we present several novel techniques to significantly speed up capacitance extraction algorithms based on boundary element methods (BEM) and to compute the capacitance extraction in the presence of floating dummy conductors. We propose the PHiCap algorithm, which is based on a hierarchical refinement algorithm and the wavelet transform. Unlike traditional algorithms which result in dense linear systems, PHiCap converts the coefficient matrix in capacitance extraction problems to a sparse linear system. PHiCap solves the sparse linear system iteratively, with much faster convergence, using an efficient preconditioning technique. We also propose a variant of PHiCap in which the capacitances are solved for directly from a very small linear system. This small system is derived from the original large linear system by reordering the wavelet basis functions and computing an approximate LU factorization. We named the algorithm RedCap. To our knowledge, RedCap is the first capacitance extraction algorithm based on BEM that uses a direct method to solve a reduced linear system. In the presence of floating dummy conductors, the equivalent capacitances among regular conductors are required. For floating dummy conductors, the potential is unknown and the total charge is zero. We embed these requirements into the extraction linear system. Thus, the equivalent capacitance matrix is solved directly. The number of system solves needed is equal to the number of regular conductors. Based on a sensitivity analysis, we propose the selective coefficient enhancement method for increasing the accuracy of selected coupling or self-capacitances with only a small increase in the overall computation time. This method is desirable for applications, such as crosstalk and signal integrity analysis, where the coupling capacitances between some conductors needs high accuracy. We also propose the variable order multipole method which enhances the overall accuracy without raising the overall multipole expansion order. Finally, we apply the multigrid method to capacitance extraction to solve the linear system faster. We present experimental results to show that the techniques are significantly more efficient in comparison to existing techniques.

No description available.

The Randomized Gauss-Seidel Method (RGS) is an iterative algorithm that solves overdetermined systems of linear equations Ax = b. This paper studies an update on the RGS method, the Randomized Block Gauss-Seidel Method. At each step, the algorithm greedily minimizes the objective function L(x) = kAx bk2 with respect to a subset of coordinates. This paper describes a Randomized Block Gauss-Seidel Method (RBGS) which uses a randomized control method to choose a subset at each step. This algorithm is the first block RGS method with an expected linear convergence rate which can be described by the properties of the matrix A and its column submatrices. The analysis demonstrates that RBGS improves RGS more when given appropriate column-paving of the matrix, a partition of the columns into well-conditioned blocks. The main result yields a RBGS method that is more e cient than the simple RGS method.

Este trabalho propõe a aplicação de dois diferentes conceitos para a resolução do Problema de Estivagem de Unidades de Celulose - PEUC, que de acordo com Ribeiro e Lorena (2008) pode ser definido como um problema que busca alocar a máxima quantidade de unidades de celulose ao porão de cargas de um dado navio, respeitando as restrições físicas de dimensões, de posicionamento, de não-sobreposição das unidades e de capacidade máxima do porão do navio. Esse tipo de problema se encaixa, no contexto da Pesquisa Operacional, na classe de Corte e Empacotamento (Cutting and Packing - C&P) e pode ser classificado, de acordo com a tipologia de Wäscher, Haußner e Schumann (2007), como sendo um Single Large Object Placement Problem (SLOPP). Em última instância, o objetivo do PEUC é definir o melhor plano de estivagem para o carregamento de unidades de celulose em um dado porão de um navio, maximizando a área ocupada pelas unidades de celulose. Trata-se de um problema NP-Completo (DOWSLAND; DOWSLAND, 1992; BISCHOFF; WÄSCHER, 1995; MALAGUTI; DURáN; TOTH, 2013) e por isso foram propostas duas abordagens para buscar a melhoria das soluções encontradas e/ou redução do tempo computacional necessário. As abordagens propostas, o Modelo Matemático Modificado e o Método Iterativo de Solução, apresentaram bons resultados para instâncias experimentais, confirmando a efetividade de suas aplicações. Os resultados foram melhores tanto na qualidade das soluções (ocupação total do objeto), como no tempo computacional necessário. Também foram avaliadas quatro instâncias reais, com a comparação dos planos de estivagem resultantes da aplicação dos modelos matemáticos com os planos reais, elaborados manualmente por especialistas. Em três dos quatro casos os resultados das abordagens aqui propostas se mostraram melhores que os planos reais. / This work proposes the application of two different concepts to tackle the Woodpulp Stowage Problem - WSP, that according to Ribeiro e Lorena (2008) can be defined as a problem that seeks the allocation of the maximum quantity of woodpulp units inside the hold of a cargo vessel, always respecting the physical constraints, positioning constraints, non-overlapping of units and also the hold capacity. This kind of problem fits, in the context of Operational Research, into the class of Cutting & Packing and can be classified, according to Wäscher, Haußner e Schumann (2007) typology, as a Single Larga Object Placement Problem (SLOPP). Ultimately the objective of the WSP is to define the best stowage plan for the loading of woodpulp units inside a given hold of a given cargo vessel, maximizing the total area occupied by the woodpulp units. As it\'s a NP-Complete problem (DOWSLAND; DOWSLAND, 1992; BISCHOFF; WÄSCHER, 1995; MALAGUTI; DURáN; TOTH, 2013) two approaches were proposed to improve the quality of the resulting solutions and/or the reduction of the computational time needed. The proposed approaches, the Modified Mathematical Model and the Iterative Solution Method, showed good results for experimental instances, confirming the effectiveness of these approaches. The results were better regarding the quality of the solutions (total occupied area of the object) and also regarding the computational time needed. Also, four real instances were evaluated, comparing the results of the mathematical models with the real stowage plans, manually created by specialists. In three of the four instances, the proposed approaches showed better results than the real stowage plans.

In this paper we develop new preconditioned iterative methods for solving monotone nonlinear eigenvalue problems. We investigate the convergence and derive grid-independent error estimates for these methods. Numerical experiments demonstrate the practical effectiveness of the proposed methods for a model problem.

Two solution strategies for large, sparse, and structured algebraic systems of equations are considered. The first strategy is to construct efficient preconditioners for iterative solvers. The second is to reduce the sparse algebraic system to a smaller, dense system of equations, which are called the boundary summation equations. The proposed preconditioners perform well when applied to equations that are discretizations of linear first order partial differential equations. Analysis shows that also very simple iterative methods converge in a number of iterations that is independent of the number of unknowns, if our preconditioners are applied to certain scalar model problems. Numerical experiments indicate that this property holds also for more complicated cases, and a flow problem modeled by the nonlinear Euler equations is treated successfully. The reduction process is applicable to a large class of difference equations. There is no approximation involved in the reduction, so the solution of the original algebraic equations is determined exactly if the reduced system is solved exactly. The reduced system is well suited for iterative solution, especially if the original system of equations is a discretization of a first order differential equation. The technique is used for several problems, ranging from scalar model problems to a semi-implicit discretization of the compressible Navier-Stokes equations. Both strategies use the concept of fundamental solutions, either of differential or difference operators. An algorithm for computing fundamental solutions of difference operators is also presented.

Este trabalho propõe a aplicação de dois diferentes conceitos para a resolução do Problema de Estivagem de Unidades de Celulose - PEUC, que de acordo com Ribeiro e Lorena (2008) pode ser definido como um problema que busca alocar a máxima quantidade de unidades de celulose ao porão de cargas de um dado navio, respeitando as restrições físicas de dimensões, de posicionamento, de não-sobreposição das unidades e de capacidade máxima do porão do navio. Esse tipo de problema se encaixa, no contexto da Pesquisa Operacional, na classe de Corte e Empacotamento (Cutting and Packing - C&P) e pode ser classificado, de acordo com a tipologia de Wäscher, Haußner e Schumann (2007), como sendo um Single Large Object Placement Problem (SLOPP). Em última instância, o objetivo do PEUC é definir o melhor plano de estivagem para o carregamento de unidades de celulose em um dado porão de um navio, maximizando a área ocupada pelas unidades de celulose. Trata-se de um problema NP-Completo (DOWSLAND; DOWSLAND, 1992; BISCHOFF; WÄSCHER, 1995; MALAGUTI; DURáN; TOTH, 2013) e por isso foram propostas duas abordagens para buscar a melhoria das soluções encontradas e/ou redução do tempo computacional necessário. As abordagens propostas, o Modelo Matemático Modificado e o Método Iterativo de Solução, apresentaram bons resultados para instâncias experimentais, confirmando a efetividade de suas aplicações. Os resultados foram melhores tanto na qualidade das soluções (ocupação total do objeto), como no tempo computacional necessário. Também foram avaliadas quatro instâncias reais, com a comparação dos planos de estivagem resultantes da aplicação dos modelos matemáticos com os planos reais, elaborados manualmente por especialistas. Em três dos quatro casos os resultados das abordagens aqui propostas se mostraram melhores que os planos reais. / This work proposes the application of two different concepts to tackle the Woodpulp Stowage Problem - WSP, that according to Ribeiro e Lorena (2008) can be defined as a problem that seeks the allocation of the maximum quantity of woodpulp units inside the hold of a cargo vessel, always respecting the physical constraints, positioning constraints, non-overlapping of units and also the hold capacity. This kind of problem fits, in the context of Operational Research, into the class of Cutting & Packing and can be classified, according to Wäscher, Haußner e Schumann (2007) typology, as a Single Larga Object Placement Problem (SLOPP). Ultimately the objective of the WSP is to define the best stowage plan for the loading of woodpulp units inside a given hold of a given cargo vessel, maximizing the total area occupied by the woodpulp units. As it\'s a NP-Complete problem (DOWSLAND; DOWSLAND, 1992; BISCHOFF; WÄSCHER, 1995; MALAGUTI; DURáN; TOTH, 2013) two approaches were proposed to improve the quality of the resulting solutions and/or the reduction of the computational time needed. The proposed approaches, the Modified Mathematical Model and the Iterative Solution Method, showed good results for experimental instances, confirming the effectiveness of these approaches. The results were better regarding the quality of the solutions (total occupied area of the object) and also regarding the computational time needed. Also, four real instances were evaluated, comparing the results of the mathematical models with the real stowage plans, manually created by specialists. In three of the four instances, the proposed approaches showed better results than the real stowage plans.