Search

In this paper, we consider a reaction-diffusion boundary value problem in a three-dimensional thin domain. The very different length scales in the geometry result in an anisotropy effect. Our study is motivated by a parabolic heat conduction problem in a thin foil leading to such anisotropic reaction-diffusion problems in each time step of an implicit time integration method [7]. The reaction-diffusion problem contains two important parameters, namely ε >0 which parameterizes the thickness of the domain and μ >0 denoting the measure for the size of the reaction term relative to that of the diffusion term. In this paper we analyze the convergence of a multigrid method with a robust (line) smoother. Both, for the W- and the V-cycle method we derive contraction number bounds smaller than one uniform with respect to the mesh size and the parameters ε and μ. [ABSTRACT FROM AUTHOR]

This paper deals with the numerical solution of optimal control problems for ODEs. The methods considered here rely on some standard optimization code to solve a discretized version of the control problem under consideration. We aim to make available to the optimization software not only the discrete objective functional, but also its gradient. The objective gradient can be computed either from forward (sensitivity) information or backward (adjoint) information. The purpose of this paper is to discuss various ways of adjoint computation. It will be shown both theoretically and numerically that methods based on the continuous adjoint equation require a careful choice of both the integrator and gradient assembly formulas in order to obtain a gradient consistent with the discretized control problem. Particular attention is given to automatic differentiation techniques which generate automatically a suitable integrator. [ABSTRACT FROM AUTHOR]

The main purpose of this paper is to introduce and study generalized hyperconnected spaces. Various characterizations of generalized hyperconnected spaces and preservation theorems are discussed. [ABSTRACT FROM AUTHOR]

The Fermat-Torricelli problem is an optimization problem associated with a finite subset $\{a_{j}\}_{j=1}^{q}$ of ℝ and a family $\{c_{j}\}_{j=1}^{q}$ of positive weights. The function F to be minimized is defined by $F(x)=\sum _{j=1}^{q}c_{j}\Vert x-a_{j}\Vert$. In this paper, we extend this problem to the case of volumes. [ABSTRACT FROM AUTHOR]

The notion of extended-well-posedness has been introduced by Zolezzi for scalar minimization problems and has been further generalized to vector minimization problems by Huang. In this paper, we study the extended well-posedness properties of vector minimization problems in which the objective function is C-quasiconvex. To achieve this task, we first study some stability properties of such problems. [ABSTRACT FROM AUTHOR]

Unique parametrizations of models are very important for parameter interpretation and consistency of estimators. In this paper we analyze the identifiability of a general class of finite mixtures of multinomial logits with varying and fixed effects, which includes the popular multinomial logit and conditional logit models. The application of the general identifiability conditions is demonstrated on several important special cases and relations to previously established results are discussed. The main results are illustrated with a simulation study using artificial data and a marketing dataset of brand choices. [ABSTRACT FROM AUTHOR]

Disjoint $\mathsf{NP}$ -pairs are a well studied complexity-theoretic concept with important applications in cryptography and propositional proof complexity. In this paper we introduce a natural generalization of the notion of disjoint $\mathsf{NP}$ -pairs to disjoint k-tuples of $\mathsf{NP}$ -sets for k≥2. We define subclasses of the class of all disjoint k-tuples of $\mathsf{NP}$ -sets. These subclasses are associated with a propositional proof system and possess complete tuples which are defined from the proof system. In our main result we show that complete disjoint $\mathsf{NP}$ -pairs exist if and only if complete disjoint k-tuples of $\mathsf{NP}$ -sets exist for all k≥2. Further, this is equivalent to the existence of a propositional proof system in which the disjointness of all k-tuples is shortly provable. We also show that a strengthening of this conditions characterizes the existence of optimal proof systems. [ABSTRACT FROM AUTHOR]

Abstract  The purpose of the paper is to apply monotone multigrid methods to static and dynamic biomechanical contact problems. In space, a finite element method involving a mortar discretization of the contact conditions is used. In time, a new contact-stabilized Newmark scheme is presented. Numerical experiments for a two body Hertzian contact problem and a biomechanical application are reported. [ABSTRACT FROM AUTHOR]

Abstract  In this paper, we propose monotonicity preserving and total variation diminishing (TVD) multigrid methods for solving scalar conservation laws. We generalize the upwind-biased residual restriction and interpolation operators for solving linear wave equations to nonlinear conservation laws. The idea is to define nonlinear restriction and interpolation based on local Riemann solutions. Theoretical analyses have been provided to analyze the monotonicity preserving and TVD properties of the resulting multigrid time stepping schemes. Numerical results are given to verify the theoretical results and demonstrate the effectiveness of the proposed schemes. Two dimensional extension is also discussed. [ABSTRACT FROM AUTHOR]

Abstract  This paper is devoted to multiresolution schemes that use a stencil selection procedure in order to obtain adaptation to the presence of edges in the images. Since non adapted schemes, based on a centered stencil, are less affected by the presence of texture, we propose the introduction of some weight that leads to a more frequent use of the centered stencil in regions without edges. In these regions the different stencils have similar weights and therefore the selection becomes an ill-posed problem with high risk of instabilities. In particular, numerical artifacts appear in the decompressed images. Our attention is centered in ENO schemes, but similar ideas can be developed for other multiresolution schemes. A nonlinear multiresolution scheme corresponding to a nonlinear interpolatory technique is analyzed. It is based on a modification of classical ENO schemes. As the original ENO stencil selection, our algorithm chooses the stencil within a region of smoothness of the interpolated function if the jump discontinuity is sufficiently big. The scheme is tested, allowing to compare its performances with other linear and nonlinear schemes. The algorithm gives results that are at least competitive in all the analyzed cases. The problems of the original ENO interpolation with the texture of real images seem solved in our numerical experiments. Our modified ENO multiresolution will lead to a reconstructed image free of numerical artifacts or blurred regions, obtaining similar results than WENO schemes. Similar ideas can be used in multiresolution schemes based in other stencil selection algorithms. [ABSTRACT FROM AUTHOR]

The idea of containment sets (csets) is due to Walster and Hansen, and the theory is mainly due to the first author. Now that floating point computation with infinities is widely accepted, it is necessary to achieve the same for interval computation. The cset of a function over a set in its domain space is the set of all limits of normal function values over that set. Csets form a sound basis for defining a number of practical models for interval arithmetic that handle division by zero and related operations in an intuitive and consistent manner. Cset-based systems offer new opportunities for compiler optimization by rearranging interval expressions, achieving tighter bounds by reducing dependencies within the expression. This paper presents basic theory. It discusses division by zero, the case for a global flag to support ``loose'' evaluation, performance, and semantics. It presents numerical examples using a trial Matlab implementation. [ABSTRACT FROM AUTHOR]

In this paper, we construct fast discrete transforms of wavelet type with an arbitrary number of zero moments for nonuniform grids. We obtain explicit formulas for the parameters defining the wavelet transform. [ABSTRACT FROM AUTHOR]

In this paper, we discuss the smoothing of an implicit function defined by a nonsmooth underdetermined system of equations F(y, z) =0. We apply a class of parametrized smoothing methods to smooth F and investigate the limiting behavior of the implicit function solving the smoothed equations. In particular, we discuss the approximation of the Clarke generalized Jacobian of the implicit function when F is piecewise smooth. As an application, we present an analysis of the generalized Karush-Kuhn-Tucker conditions of different forms for a piecewise-smooth equality-constrained minimization problem. [ABSTRACT FROM AUTHOR]

In this paper a Galerkin method is used to construct series solutions of a homogeneous Sturm-Liouville problem defined on [0, &pie;]. The series constructed are shown to converge to a specified du Bois-Reymond function f in L² [0,&Pie;]. It is then shown that the series solutions can be made to converge uniformly to the specified du Bois-Reymond function when averaged by the Cesaro-one summability method. Therefore, in the Cesaro-one sense, every continuous function / o n [0, &Pie;] is the uniform limit of solutions of non-homogeneous Sturm-Liouville problems. [ABSTRACT FROM AUTHOR]

In this paper, we consider a general primal-dual nonlinear rescaling (PDNR) method for convex optimization with inequality constraints. We prove the global convergence of the PDNR method and estimate the error bounds for the primal and dual sequences. In particular, we prove that, under the standard second-order optimality conditions, the error bounds for the primal and dual sequences converge to zero with linear rate. Moreover, for any given ratio 0 < γ < 1, there is a fixed scaling parameter kγ > 0 such that each PDNR step shrinks the primal-dual error bound by at least a factor 0 < γ < 1, for any k ≥ kγ. The PDNR solver was tested on a variety of NLP problems including the constrained optimization problems (COPS) set. The results obtained show that the PDNR solver is numerically stable and produces results with high accuracy. Moreover, for most of the problems solved, the number of Newton steps is practically independent of the problem size. [ABSTRACT FROM AUTHOR]

In this paper we develop techniques for computing elementwise conservative approximations of the flux on element boundaries for the continuous Galerkin method. The technique is based on computing a correction of the average normal flux on an edge or face. The correction is a jump in a piecewise constant or linear function. We derive a basic algorithm which is based on solving a global system of equations and a parallel algorithm based on solving local problems on stars. The methods work on meshes with different element types and hanging nodes. We prove existence, uniqueness, and optimal order error estimates. Lastly, we illustrate our results by a few numerical examples. [ABSTRACT FROM AUTHOR]

For the problem of minimizing the sum of Euclidean norms (MSN), most existing quadratically convergent algorithms require a strict complementarity assumption. However, this assumption is not satisfied for a number of MSN problems. In this paper, we present a globally and quadratically convergent algorithm for the MSN problem. In particular, the quadratic convergence result is obtained without assuming strict complementarity. Examples without strictly complementary solutions are given to show that our algorithm can indeed achieve quadratic convergence. Preliminary numerical results are reported. [ABSTRACT FROM AUTHOR]

In the present paper, we classify the finite subgroups $G\backslash subset\backslash operatorname\{GL\}\_4(\backslash Bbb\; C)$ such that the quotient $\backslash Bbb\; C^4$ by the action $G$ has only isolated terminal Gorenstein singularities. [ABSTRACT FROM AUTHOR]