Search

In this paper, we prove some results on uniqueness of meromorphic functions with three weighted sharing values. The results in this paper improve those given by H. X. Yi, I. Lahiri, T. C. Alzahary and H. X. Yi and other authors. [ABSTRACT FROM AUTHOR]

We consider the scalar conservation law with flux function discontinuous in the space variable, i.e., \begin{eqnarray} \label{eq1} u_t+(H(x)f(u)+(1-H(x))g(u))_{x} &=& 0 \quad \mbox{in } \R \times \R_{+}, \nonumber \\ u(0, x) &=& u_{0}(x) \quad \mbox{in } \R, \label{0.1} \end{eqnarray} where $H$ is the Heaviside function and $f$ and $g$ are smooth with the assumptions that either $f$ is convex and $g$ is concave or $f$ is concave and $g$ is convex. The existence of a weak solution of (\ref{eq1}) is proved by showing that upwind finite difference schemes of Godunov and Enquist--Osher type converge to a weak solution. Uniqueness follows from a Kruzkhov-type entropy condition. We also provide explicit solutions to the Riemann problem for (\ref{eq1}). At the level of numerics, we give easy-to-implement numerical schemes of Godunov and Enquist--Osher type. The central feature of this paper is the modification of the singular mapping technique (the main analytical tool for these types of equations) which allows us to show that the numerical schemes converge. Equations of type (\ref{eq1}) with the above hypothesis on the flux may occur when considering the following scalar conservation law with discontinuous flux: \begin{equation} \label{eq2} \begin{array}{r@{\;}l} u_t + (k (x) f (u))_x &= 0, \\ u (0, x) &= u_0 (x), \end{array} \end{equation} with $f$ convex and $k$ of indefinite sign. [ABSTRACT FROM AUTHOR]

A novel filled function is given in this paper to find a global minima for a nonsmooth constrained optimization problem. First, a modified concept of the filled function for nonsmooth constrained global optimization is introduced, and a filled function, which makes use of the idea of the filled function for unconstrained optimization and penalty function for constrained optimization, is proposed. Then, a solution algorithm based on the proposed filled function is developed. At last, some preliminary numerical results are reported. The results show that the proposed approach is promising. [ABSTRACT FROM AUTHOR]

In this paper, we study the relation between the error estimate of the bilinear interpolation on a general quadrilateral and the geometric characters of the quadrilateral. Some explicit bounds of the interpolation error are obtained based on some sharp estimates of the integral over 1/∣J∣p-1 for 1 ≤ p≤∞ on the reference element, where J is the Jacobian of the nonaffine mapping. This allows us to introduce weak geometric conditions (depending on p) leading to interpolation error estimates in the W1,p norm, for any p ϵ [1,∞), which can be regarded as a generalization of the regular decomposition property (RDP) condition introduced in [G. Acosta and R. G. Durán, SIAM J. Numer. Anal., 38 (2000), pp. 1073-1088] for p = 2 and new RDP conditions (NRDP) for p ≠ 2. We avoid the use of the reference family elements, which allows us to extend the results to a larger class of elements and to introduce the NRDP condition in a more unified way. As far as we know, the mesh condition presented in this paper is weaker than any other mesh conditions proposed in the literature for any p with 1 ≤ p≤∞. [ABSTRACT FROM AUTHOR]

In this paper, the frequency domain-based numerical methods for simulation of fractional order systems are studied in the sense of stability preservation. First, the stability boundary curve is exactly determined for these methods. Then, this boundary is analyzed and compared with an accurate (ideal) boundary in different frequency ranges. Also, the critical regions in which the stability does not preserve are determined. Finally, the analytical achievements are confirmed via some numerical illustrations. [ABSTRACT FROM AUTHOR]

In this paper we propose and analyze a new, fast method for the numerical solution of time domain boundary integral formulations of the wave equation. We employ Lubich's convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The coefficient matrix of the arising system of linear equations is a triangular block Toeplitz matrix. Possible choices for solving the linear system arising from the above discretization include the use of fast Fourier transform (FFT) techniques and the use of data-sparse approximations. By using FFT techniques, the computational complexity can be reduced substantially while the storage cost remains unchanged and is, typically, high. Using data-sparse approximations, the gain is reversed; i.e., the computational cost is (approximately) unchanged while the storage cost is substantially reduced. The method proposed in this paper combines the advantages of these two approaches. First, the discrete convolution (related to the block Toeplitz system) is transformed into the (discrete) Fourier image, thereby arriving at a decoupled system of discretized Helmholtz equations with complex wave numbers. A fast data-sparse (e.g., fast multipole or panel-clustering) method can then be applied to the transformed system. Additionally, significant savings can be achieved if the boundary data are smooth and time-limited. In this case the right-hand sides of many of the Helmholtz problems are almost zero, and hence can be disregarded. Finally, the proposed method is inherently parallel. We analyze the stability and convergence of these methods, thereby deriving the choice of parameters that preserves the convergence rates of the unperturbed convolution quadrature. We also present numerical results which illustrate the predicted convergence behavior. [ABSTRACT FROM AUTHOR]

In this paper, we consider the singularly perturbed reaction-diffusion problem in one and two dimensions. The boundary value problem is decomposed into a first-order system to which a suitable weighted least-squares formulation is proposed. A robust, stable, and efficient approach is developed based on local discontinuous Galerkin (LDG) discretization for the weak form. Uniform error estimates are derived. Numerical examples are presented to illustrate the method and the theoretical results. Comparison studies are made between the proposed method and other methods. [ABSTRACT FROM AUTHOR]

This paper studies the boundedness, global asymptotic stability, and periodicity of positive solutions of the equation xn = f(xn-2)/g(xn-1), n ϵ ℕ0, where f, g ϵ C[(0,∞), (0,∞)]. It is shown that if f and g are nondecreasing, then for every solution of the equation the subsequences {x2n} and {x2n-1} are eventually monotone. For the case when f(x) = a + βx and g satisfies the conditions g(0) = 1, g is nondecreasing, and x/g(x) is increasing, we prove that every prime periodic solution of the equation has period equal to one or two. We also investigate the global periodicity of the equation, showing that if all solutions of the equation are periodic with period three, then f(x) = c1/x and g(x) = c2x, for some positive c1 and c2. [ABSTRACT FROM AUTHOR]

Kupershmidt and Tuenter have introduced reflection symmetries for the q-Bernoulli numbers and the Bernoulli polynomials in (2005), (2001), respectively. However, they have not dealt with congruence properties for these numbers entirely. Kupershmidt gave a quantization of the reflection symmetry for the classical Bernoulli polynomials. Tuenter derived a symmetry of power sum polynomials and the classical Bernoulli numbers. In this paper, we study the new symmetries of the q-Bernoulli numbers and polynomials, which are different from Kupershmidt's and Tuenter's results. By using our symmetries for the q-Bernoulli polynomials, we can obtain some interesting relationships between q-Bernoulli numbers and polynomials. [ABSTRACT FROM AUTHOR]

The main purpose of this paper is to study the distribution of Genocchi polynomials. Finally, we construct the Genocchi zeta function which interpolates Genocchi polynomials at negative integers. [ABSTRACT FROM AUTHOR]

A new evolving network based on the scale-free network of Barabási and Albert (BA) is studied, and the accelerated attachment of new edges is considered in its evolving process. The accelerated attachment is different from the previous accelerated growth of edges and has two particular meanings in this paper. One is that a new vertex with the edges is inserted into the network with acceleration at each time step; the other is that, with a given probability, some additional edges are linked with the vertices in proportion to the number of their obtained edges in the latest evolving periods. The new model describes the cases of those complex networks with a few exceptional vertices. The attachment mechanism of the new adding edges for these vertices does not follow the preferential attachment rule. Comparing with the linear edge growth model, the characteristics of the accelerated growth model are studied theoretically and numerically. We show that the degree distributions of these models have a power law decay and the exponents are larger than that of the BA model. We point out that the characteristics of the exceptional vertices and the aging vertices in an aging network are not identical. The reasons for neglecting this attachment in most of evolving networks are also summarized. [ABSTRACT FROM AUTHOR]

The filter generator is an important building block in many stream ciphers. The generator consists of a linear feedback shift register of length n that generates an m-sequence of period 2′ - 1 filtered through a Boolean function of degree d that combines bits from the shift register and creates an output bit zt at any time t. The previous best attacks aimed at reconstructing the initial state from an observed keystream, have essentially reduced the problem to solving a nonlinear system of D = (Multiple line equation(s) cannot be represented in ASCII text) (i) equations in n unknowns using techniques based on linear algebra. This attack needs about D bits of keystream and the system can be solved in complexity O (Dω), where ω can be taken to be Strassen's reduction exponent ω = log2 (7) ≈ 2.807. This paper describes a new algorithm that recovers the initial state of most filter generators after observing O(D) keystream bits with complexity O((D - n)/2) ≈ O(D), after a pre-computation with complexity O(D(log2 D)³). [ABSTRACT FROM AUTHOR]

This paper introduces a new combinatorial construction for q-ary constant-weight codes which yields several families of optimal codes and asymptotically optimal codes. The construction reveals intimate connection between q-ary constant-weight codes and sets of pairwise disjoint combinatorial designs of various types. [ABSTRACT FROM AUTHOR]

In this paper, we analyze convergence rates of wavelet schemes for time-dependent convection-reaction equations within the framework of the Eulerian--Lagrangian localized adjoint method (ELLAM). Under certain minimal assumptions that guarantee $ H^1 $-regularity of exact solutions, we show that a generic ELLAM scheme has a convergence rate $ \mathcal{O}(h/\sqrt{\Delta t} + \Delta t) $ in $ L^2 $-norm. Then, applying the theory of operator interpolation, we obtain error estimates for initial data with even lower regularity. Namely, it is shown that the error of such a scheme is $ \mathcal{O}((h/\sqrt{\Delta t})^\theta + (\Delta t)^\theta) $ for initial data in a Besov space $ \displaystyle B^\theta_{2,q} (0 < \theta < 1, 0 < q <= infinity) $. The error estimates are {a priori} and optimal in some cases. Numerical experiments using orthogonal wavelets are presented to illustrate the theoretical estimates. [ABSTRACT FROM AUTHOR]

The aim of this work is to solve hyperbolic conservation laws by means of a finite volume method for both spatial and time discretization. We extend the ideas developed in [X.-D. Liu and S. Osher, SIAM J. Numer. Anal., 33 (1996), pp. 760--779; X.-D. Liu and E. Tadmor, Numer. Math., 79 (1998), pp. 397--425] to fourth-order upwind and central schemes. In order to do this, once we know the cell-averages of the solution, $\overline {u}_j ^n$, in cells $I_{j}$ at time $T=t^n$, we define a new three-degree reconstruction polynomial that in each cell, $I_{j}$, presents the same shape as the cell-averages $\{ {\overline {u}_{j-1} ^n,\overline {u}_j ^n,\overline {u}_{j+1} ^n}\}$. By combining this reconstruction with the nonoscillatory property and the maximum principle requirement described in [X.-D. Liu and S. Osher, SIAM J. Numer. Anal., 33 (1996), pp. 760--779] we obtain a fourth-order scheme that satisfies the total variation bounded (TVB) property. Extension to systems is carried out by componentwise application of the scalar framework. Numerical experiments confirm the order of the schemes presented in this paper and their nonoscillatory behavior in different test problems. [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 study option pricing under a regime-switching exponential Lévy model. Assuming that the coefficients are time-dependent and modulated by a finite state Markov chain, we generalise the work in Momeya and Morales (Method Comput Appl Probab, 2014, doi:), and Siu and Yang (Acta Mathe Appl Sin 2:369-388, 2009), that is, we use a pricing method based on the Esscher transform conditional on the information available on the Markov chain. We also carry out numerical analysis, to show the impact of the risk induced by the underlying Markov chain on the price of the option. [ABSTRACT FROM AUTHOR]

We present in this paper numerical simulations of reactive flows interacting with shock waves. We argue that spectral methods are suitable for these problems and review the recent developments in spectral methods that have made them a powerful numerical tool appropriate for long-term integrations of complicated flows, even in the presence of shock waves. A spectral code is described in detail, and the theory that leads to each of its components is explained. Results of interactions of hydrogen jets with shock waves are presented and analyzed, and comparisons with ENO finite difference schemes are carried out. [ABSTRACT FROM AUTHOR]

This paper develops a least-squares approach to the solution of the incompressible Navier­Stokes equations in primitive variables. As with our earlier work on Stokes equations, we recast the Navier­Stokes equations as a first-order system by introducing a velocity-flux variable and associated curl and trace equations. We show that a least-squares principle based on L2 norms applied to this system yields optimal discretization error estimates in the H1 norm in each variable, including the velocity flux. An analogous principle based on the use of an H-1 norm for the reduced system (with no curl or trace constraints) is shown to yield similar estimates, but now in the L2 norm for velocity-flux and pressure. Although the H-1 least-squares principle does not allow practical implementation, these results are critical to the analysis of a practical least-squares method for the reduced system based on a discrete equivalent of the negative norm. A practical method of this type is the subject of a companion paper. Finally, we establish optimal multigrid convergence estimates for the algebraic system resulting from the L2 norm approach. [ABSTRACT FROM AUTHOR]

Mesh matching is an effective way to convert the non-conforming interfaces between two hexahedral meshes into conforming ones, which is very important for achieving high-quality finite element analysis. However, the existing mesh matching algorithm is neither efficient nor effective enough to handle complex interfaces and self-intersecting sheets. In this paper, the algorithm is improved in three aspects: (1) by introducing a more precise criteria for chord matching and the concept of partition chord set, complex interfaces with internal loops can be handled more effectively; (2) by proposing a new solution, self-intersecting sheet can be inflated and extracted locally; and (3) by putting forward a mesh quality evaluation method, the sheet extraction operation during mesh matching can be done more efficiently. Our improved mesh matching algorithm is fully automatic, and its effectiveness is demonstrated by several examples in different matching situations. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we apply a complete pivoting strategy to compute the left-looking version of AINV preconditioner for linear systems. As the preprocessing, the MultiLevel Nested Dissection reordering has also been applied. We have used this preconditioner as the right preconditioner for several linear systems where the coefficient matrices have been downloaded from the University of Florida Sparse Matrix Collection. Numerical experiments presented in this paper indicate the effectiveness of such a complete pivoting on the quality of left-looking version of AINV preconditioner. [Copyright &y& Elsevier]

In this paper we obtain sufficient conditions for the asymptotic properties of solutions of two dimensional neutral difference systems. Our result extends some existing results in the literature. An example is given to illustrate the result. [ABSTRACT FROM AUTHOR]

In this correspondence, we point out an error in equation (4) in the above paper. We also correct equation (4) by a detailed mathematical deduction. [ABSTRACT FROM AUTHOR]

Abstract Given a finite group G , we say that G has property P n if for every prime integer p , G has at most n − 1 irreducible characters whose degrees are multiples of p. In this paper, we classify all finite groups that have property P 4. We show that the groups satisfying property P 4 are exactly the finite groups with at most three nonlinear irreducible characters, one solvable group of order 168, SL 2 (3) , A 5 , S 5 , PSL 2 (7) and A 6. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we use the exp-function method to construct some new soliton solutions of the Benjamin–Bona–Mahony and modified Benjamin–Bona–Mahony equations. These equations have important and fundamental applications in mathematical physics and engineering sciences. The exp-function method is used to find the soliton solution of a wide class of nonlinear evolution equations with symbolic computation. This method provides the concise and straightforward solution in a very easier way. The results obtained in this paper can be viewed as a refinement and improvement of the previously known results. [Copyright &y& Elsevier]

This paper presents a general and comprehensive description of Optimization Methods, and Algorithms from a novel viewpoint. It is shown, in particular, that Direct Methods, Iterative Methods, and Computer Science Algorithms belong to a well-defined general class of both Finite and Infinite Procedures, characterized by suitable descent directions. [ABSTRACT FROM AUTHOR]

A strong indication about the existence of a (7 p, 4, 1) difference family with p ≡ 7 (mod 12) a prime has been given in [11]. Here, developing some ideas of that paper, we give, much more generally, a strong indication about the existence of a cyclic ( pq, 4, 1) difference family whenever p and q are primes congruent to 7 (mod 12) and of a cyclic ( pq, 5, 1) difference family whenever p and q are primes congruent to 11 (mod 20). Indeed we give an algorithm for their construction that seems to be always successful and we have checked it works whenever both primes p and q do not exceed 1,000. All our ( pq, 4, 1) and ( pq, 5, 1) difference families have the nice property of admitting a multiplier of order 3 or 5, respectively, that fixes almost all base blocks. As an intermediate result we also find an optimal ( p, 5, 1) optical orthogonal code for every prime p ≡ 11 (mod 20) not exceeding 10,000. [ABSTRACT FROM AUTHOR]

This paper presents a new method to isolate all configurations that a multiloop linkage can adopt. The problem is tackled by means of formulation and resolution techniques that fit particularly well together. The adopted formulation yields a system of simple equations (only containing linear, bilinear, and quadratic monomials, and trivial trigonometric terms for the helical pair only) whose structure is later exploited by a branch-and-prune method based on linear relaxations. The method is general, as it can be applied to linkages with single or multiple loops with arbitrary topology, involving lower pairs of any kind, and complete, as all possible solutions get accurately bounded, irrespective of whether the linkage is rigid or mobile. [ABSTRACT FROM AUTHOR]

This paper deals with superconvergence phenomena in general grids when projectionbased approximations are used for solving Fredholm integral equations of the second kind with weakly singular kernels. Four variants of the Galerkin method are considered. They are the classical Galerkin method, the iterated Galerkin method, the Kantorovich method, and the iterated Kantorovich method. It is proved that the iterated Kantorovich approximation exhibits the best superconvergence rate if the right-hand side of the integral equation is nonsmooth. All error estimates are derived for an arbitrary grid without any uniformity or quasi-uniformity condition on it, and are formulated in terms of the data without any additional assumption on the solution. Numerical examples concern the equation governing transfer of photons in stellar atmospheres. The numerical results illustrate the fact that the error estimates proposed in the different theorems are quite sharp, and confirm the superiority of the iterated Kantorovich scheme. [ABSTRACT FROM AUTHOR]

A recent study [D. N. Arnold, D. Boffi, and R. S. Falk, SIAM J. Numer. Anal., 42 (2005), pp. 2429-2451] reveals that convergence of finite element methods using H(div , O)-compatible finite element spaces deteriorates on nonaffine quadrilateral grids. This phenomena is particularly troublesome for the lowest-order Raviart-Thomas elements, because it implies loss of convergence in some norms for finite element solutions of mixed and least-squares methods. In this paper we propose reformulation of finite element methods, based on the natural mimetic divergence operator [M. Shashkov, Conservative Finite Difference Methods on General Grids, CRC Press, Boca Raton, FL, 1996], which restores the order of convergence. Reformulations of mixed Galerkin and leastsquares methods for the Darcy equation illustrate our approach. We prove that reformulated methods converge optimally with respect to a norm involving the mimetic divergence operator. Furthermore, we prove that standard and reformulated versions of the mixed Galerkin method lead to identical linear systems, but the two versions of the least-squares method are veritably different. The surprising conclusion is that the degradation of convergence in the mixed method on nonaffine quadrilateral grids is superficial, and that the lowest-order Raviart-Thomas elements are safe to use in this method. However, the breakdown in the least-squares method is real, and there one should use our proposed reformulation. [ABSTRACT FROM AUTHOR]

We develop a family of new interior penalty discontinuous Galerkin methods for the Keller-Segel chemotaxis model. This model is described by a system of two nonlinear PDEs: a convection-diffusion equation for the cell density coupled with a reaction-diffusion equation for the chemoattractant concentration. It has been recently shown that the convective part of this system is of a mixed hyperbolic-elliptic-type, which may cause severe instabilities when the studied system is solved by straightforward numerical methods. Therefore, the first step in the derivation of our new methods is made by introducing the new variable for the gradient of the chemoattractant concentration and by reformulating the original Keller-Segel model in the form of a convection-diffusionreaction system with a hyperbolic convective part. We then design interior penalty discontinuous Galerkin methods for the rewritten Keller-Segel system. Our methods employ the central-upwind numerical fluxes, originally developed in the context of finite-volume methods for hyperbolic systems of conservation laws. In this paper, we consider Cartesian grids and prove error estimates for the proposed high-order discontinuous Galerkin methods. Our proof is valid for pre-blow-up times since we assume boundedness of the exact solution. We also show that the blow-up time of the exact solution is bounded from above by the blow-up time of our numerical solution. In the numerical tests presented below, we demonstrate that the obtained numerical solutions have no negative values and are oscillation-free, even though no slope-limiting technique has been implemented. [ABSTRACT FROM AUTHOR]

In a paper by J. V. Pereira and the author it was proved that any pencil of plane curves of degree $d>1$ with irreducible generic fiber can have at most five completely reducible fibers although no examples with five such fibers had ever been found. Recently Janis Stipins has proved that if a pencil has a base of $d^2$ points, then it cannot have five completely reducible fibers. In this paper we generalize Stipins' result to arbitrary pencils. We also include into consideration more general special fibers that are the unions of lines and non-reduced curves. These fibers are important for characteristic varieties of hyperplane complements. [ABSTRACT FROM AUTHOR]

Judy and Motta developed a customizable electrodeposition process for fabrication of very small metal structures on a substrate. In this process, layers of metal of various shapes are placed on the substrate, then the substrate is inserted in an electroplating solution. Some of the metal layers have power applied to them, while the rest of the metal layers are not connected to the power initially. Metal ions in the plating solution start depositing on the powered layers and a surface grows from the powered layers. As the surface grows, it will touch metal layers that were initially unpowered, causing them to become powered and to start growing with the rest of the surface. The metal layers on the substrate are known as seed layer patterns, and different seed layer patterns can produce different shapes. This paper presents a mathematical model, a forward simulation method, and an inverse problem solution for the growth of a surface from a seed layer pattern. The model describes the surface evolution as uniform growth in the direction normal to the surface. This growth is simulated in two and three dimensions using the level set method. The inverse problem is to design a seed layer pattern that produces a desired surface shape. Some surface shapes are not attainable by any seed layer pattern. For smooth attainable shapes, we present a computational method that solves this inverse problem. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, high-resolution finite volume schemes are combined with an adaptive mesh technique inspired by multiresolution analysis to improve the computational efficiency for two-dimensional hyperbolic conservation laws. The method is conservative. Moreover, it is stable which is proven numerically in this paper. The computational grid is dynamically adapted so that higher spatial resolution is automatically allocated to regions where strong gradients are observed. Using this proposed scheme, we compute several two-dimensional model problems and a compressive rate ranging from about 5–10 is observed in all simulations. [Copyright &y& Elsevier]

The aim of this paper, firstly, is to construct generating functions of q-Euler numbers and polynomials of higher order by applying the fermionic p-adic q-Volkenborn integral, secondly, to define multivariate q-Euler zeta function (Barnes-type Hurwitz q-Euler zeta function) and l-function which interpolate these numbers and polynomials at negative integers, respectively. We give relation between Barnes-type Hurwitz q-Euler zeta function and multivariate q-Euler l-function. Moreover, complete sums of products of these numbers and polynomials are found. We give some applications related to these numbers and functions as well. [ABSTRACT FROM AUTHOR]

Abstract: The sub-elliptic obstacle problem arises in various branches of the applied sciences, e.g., in mechanical engineering and robotics, mathematical finance, image reconstruction and neurophysiology. In the recent paper [Donatella Danielli, Nicola Garofalo, Sandro Salsa, Variational inequalities with lack of ellipticity. I. Optimal interior regularity and non-degeneracy of the free boundary, Indiana Univ. Math. J. 52 (2) (2003) 361–398; MR1976081 (2004c:35424)] it was proved that weak solutions to the sub-elliptic obstacle problem in a Carnot group belong to the Folland–Stein (optimal) Lipschitz class (the analogue of the well-known interior local regularity for the classical obstacle problem). However, the regularity of the free boundary remained a challenging open problem. In this paper we prove that, in Carnot groups of step , the free boundary is (Euclidean) near points satisfying a certain thickness condition. This constitutes the sub-elliptic counterpart of a celebrated result due to Caffarelli [Luis A. Caffarelli, The regularity of free boundaries in higher dimensions, Acta Math. 139 (3–4) (1977) 155–184; MR0454350 (56 #12601)]. [Copyright &y& Elsevier]

In this paper we present sufficient conditions for global optimality of non-convex quadratic programs involving linear matrix inequality (LMI) constraints. Our approach makes use of the concept of a quadratic subgradient. We develop optimality conditions for quadratic programs with LMI constraints by using Lagrangian function and by examining conditions which minimizes a quadratic subgradient of the Lagrangian function over simple bounding constraints. As applications, we obtain sufficient optimality condition for quadratic programs with LMI and box constraints by minimizing a quadrtic subgradient over box constraints. We also give optimality conditions for quadratic minimization involving LMI and binary constraints. [ABSTRACT FROM AUTHOR]

Abstract: This paper is devoted to the study of the internal stabilization of the Bernoulli–Euler plate equation in a square. The continuous and the space semi-discretized problems are successively considered and analyzed using a frequency domain approach. For the infinite-dimensional problem, we provide a new proof of the exponential stability result, based on a two-dimensional Ingham''s type result. In the second and main part of the paper, we propose a finite difference space semi-discretization scheme and we prove that this scheme yields a uniform exponential decay rate (with respect to the mesh size). [Copyright &y& Elsevier]

We study two space-time mesh refinement methods as the one introduced in [F. Collino, T. Fouquet, and P. Joly, Numer. Math., 95 (2003), pp. 197--221]. The stability of such methods is guaranteed by construction through the conservation of a discrete energy. In this paper, we show the L2 convergence of these schemes and provide optimal error estimates. The proof is based on energy techniques and bootstrap arguments. Our results are validated with numerical simulations and compared with results from plane wave analysis [F. Collino, T. Fouquet, and P. Joly, Numer. Math., 95 (2003), pp. 223--251]. [ABSTRACT FROM AUTHOR]

A finite element discretization is proposed and analyzed for a linear parabolic problems with discontinuous coefficients. Due to low global regularity of the solution, it seems difficult to achieve optimal order of convergence with classical finite element methods [Numer. Math., 79 (1998), pp. 175--202]. In this paper, we have used a finite element discretization, where interface triangles are assumed to be curved triangles instead of straight triangles as in classical finite element methods. Optimal order error estimates in L2 and H1 norms are shown to hold even if the regularity of the solution is low on the whole domain. While the continuous time Galerkin method is discussed for the spatially discrete scheme, the discontinuous Galerkin method is analyzed for the fully discrete scheme. The interfaces and boundaries of the domains are assumed to be smooth for our purpose. [ABSTRACT FROM AUTHOR]

A complete error analysis is performed for the spectral-Galerkin approximation of a model Helmholtz equation with high wave numbers. The analysis presented in this paper does not rely on the explicit knowledge of continuous/discrete Green's functions and does not require any mesh condition to be satisfied. Furthermore, new error estimates are also established for multi-dimensional radial and spherical symmetric domains. Illustrative numerical results in agreement with the theoretical analysis are presented. [ABSTRACT FROM AUTHOR]

This paper deals with the upwind finite volume method applied to the linear advection equation on a bounded domain and with natural boundary conditions. We introduce what we call the geometric corrector, which is a sequence associated with every finite volume mesh in $\mathbf{R}^{nd}$ and every nonvanishing vector $\mathbf{a}$ of $\mathbf{R}^{nd}$. First we show that if the continuous solution is regular enough and if the norm of this corrector is bounded by the mesh size, then an order one error estimate for the finite volume scheme occurs. Afterwards we prove that this norm is indeed bounded by the mesh size in several cases, including the one where an arbitrary coarse conformal triangular mesh is uniformly refined in two dimensions. Computing numerically exactly this corrector allows us to state that this result might be extended under conditions to more general cases, such as the one with independent refined meshes. [ABSTRACT FROM AUTHOR]

In this paper, we consider a multilithology diffusion model used in the field of stratigraphic basin simulations to simulate large scale depositional transport processes of sediments described as a mixture of L lithologies. This model is a simplified one for which the surficial fluxes are proportional to the slope of the topography and to a lithology fraction with unitary diffusion coefficients. The main variables of the system are the sediment thickness h, the L surface concentrations cis in lithology i of the sediments at the top of the basin, and the L concentrations ci in lithology i in the sediments inside the basin. For this simplified model, the sediment thickness decouples from the other unknowns and satisfies a linear parabolic equation. The remaining equations account for the mass conservation of the lithologies, and couple, for each lithology, a first order linear equation for cis with a linear advection equation for ci for which cis appears as an input boundary condition. For this coupled system, a weak formulation is introduced. The system is discretized by an implicit time integration and a cell centered finite volume method. This numerical scheme is shown to satisfy stability estimates and to converge, up to a subsequence, to a weak solution of the problem. [ABSTRACT FROM AUTHOR]

This paper presents an approach to transform asymmetric systems into symmetric systems by equivalence transformation and discusses what forms of restrictions should be imposed on the system matrices so that they can be simultaneously transformed into symmetric matrices. Conditions of symmetrizability obtained here are more "liberal" and numerical calculations of this transformation are more straightforward. Several examples are provided to illustrate the new approach. [ABSTRACT FROM AUTHOR]

A Markov property associates a set of conditional independencies to a graph. Two alternative Markov properties are available for chain graphs (CGs), the Lauritzen–Wermuth–Frydenberg (LWF) and the Andersson–Madigan– Perlman (AMP) Markov properties, which are different in general but coincide for the subclass ofCGs with no flags. Markov equivalence induces a partition of the class of CGs into equivalence classes and every equivalence class contains a, possibly empty, subclass of CGs with no flags itself containing a, possibly empty, subclass of directed acyclic graphs (DAGs). LWF-Markov equivalence classes of CGs can be naturally characterized by means of the so-calledlargest CGs, whereas a graphical characterization of equivalence classes of DAGs is provided by theessential graphs. In this paper, we show the existence of largest CGs with no flags that provide a natural characterization of equivalence classes of CGs of this kind, with respect to both the LWF- and the AMP-Markov properties. We propose a procedure for the construction of the largest CGs, the largest CGs with no flags and the essential graphs, thereby providing a unified approach to the problem. As by-products we obtain a characterization of graphs that are largest CGs with no flags and an alternative characterization of graphs which are largest CGs. Furthermore, a known characterization of the essential graphs is shown to be a special case of our more general framework. The three graphical characterizations have a common structure: they use two versions of a locally verifiable graphical rule. Moreover, in case of DAGs, an immediate comparison of three characterizing graphs is possible. [ABSTRACT FROM AUTHOR]

Stochastic discrimination (SD) depends on a discriminant function for classification. In this paper, an improved SD is introduced to reduce the error rate of the standard SD In the context of a two-class classification problem. The learning procedure of the improved SD consists of two stages. Initially a standard SD, but with shorter learning period is carried out to identify an important space where all the misclassified samples are located. Then the standard SD is modified by 1) restricting sampling in the important space, and 2) introducing a new discriminant function for samples in the important space. It is shown by mathematical derivation that the new discriminant function has the same mean, but with a smaller variance than that of the standard SD for samples in the important space. It is also analyzed that the smaller the variance of the discriminant function, the lower the error rate of the classifier. Consequently, the proposed improved SD improves standard SD by its capability of achieving higher classification accuracy. Illustrative examples are provided to demonstrate the effectiveness of the proposed improved SD. [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]