A parallel algorithm for discrete least squares rational approximation.
We discuss the formulation of a simulator in three spatial dimensions for a multicomponent, two phase (air, water) system of groundwater flow and transport with biodegradation kinetics and wells with multiple screens. The simulator has been developed for parallel, distributed memory, message passing machines. The numerical procedures employed are a fully implicit expanded mixed finite element method for flow and either a characteristics-mixed method or a Godunov method for transport and reactions...
In the present work we describe HPEC (High Performance Eigenvalues Computation), a parallel software package for the evaluation of some eigenvalues of a large sparse symmetric matrix. It implements an efficient and portable Block Lanczos algorithm for distributed memory multicomputers. HPEC is based on basic linear algebra operations for sparse and dense matrices, some of which have been derived by ScaLAPACK library modules. Numerical experiments have been carried out to evaluate HPEC performance...
In this paper, we present a parallel scheme to solve the population balance equations based on the method of characteristics and the finite element discretization. The application of the method of characteristics transform the higher dimensional population balance equation into a series of lower dimensional convection-diffusion-reaction equations which can be solved in a parallel way. Some numerical results are presented to show the accuracy and efficiency.
A direct projection method for solving systems of linear algebraic equations is described. The algorithm is equivalent to the algorithm for minimization of the corresponding quadratic function and can be generalized for the minimization of a strictly convex function.
We give a derivation of an a-posteriori strategy for choosing the regularization parameter in Tikhonov regularization for solving nonlinear ill-posed problems, which leads to optimal convergence rates. This strategy requires a special stability estimate for the regularized solutions. A new proof fot this stability estimate is given.
This paper presents a postprocessing technique for estimating the local regularity of numerical solutions in high-resolution finite element schemes. A derivative of degree p ≥ 0 is considered to be smooth if a discontinuous linear reconstruction does not create new maxima or minima. The intended use of this criterion is the identification of smooth cells in the context of p-adaptation or selective flux limiting. As a model problem, we consider a 2D convection equation discretized with bilinear finite...