We use Newton divided differences for calculation of Greene sums -- the
rational functions determined by linear extensions of partially ordered sets.
Identities for Greene sums generate relations for Newton divided differences
and Arnold differential forms. Also generalizations of the Newton interpolation
series which are indexed by sequences of partially ordered sets are received.
lunes, 30 de noviembre de 2009
Calculus of linear extensions and Newton interpolation. (arXiv:0911.5620v1 [math.CO])
Coarse space over the ages. (arXiv:0911.5725v1 [math.NA])
The objective of this paper is to explain the principles of the design of a
coarse space in a simplified way and by pictures. The focus is on ideas rather
than on a more historically complete presentation. Also, space limitation does
not allow even to mention many important methods and papers that should be
rightfully included.
On Adaptive-Multilevel BDDC. (arXiv:0911.5730v1 [math.NA])
We combine the advantages of the adaptive and multilevel approaches, proposed
previously by the authors, to propose a new method that preserves both,
parallel scalability with increasing number of subdomains and excellent
convergence properties. Performance of the method is illustrated on a
two-dimensional problem of linear elasticity.
Numerical Study of Liquid Crystal Elastomer Using Mixed Finite Element Method. (arXiv:0911.5415v1 [math.NA])
In this paper, we tried to model the elastic behavior of liquid crystal
elastomer using mixed finite element method. We start from an energy functional
which includes Blandon's stored energy of LCE, penalization of change of
directors, and two Lagrangian terms enforcing incompressibility and the unity
of the directors. The resulting Euler-Lagrange equation is a nonlinear equation
of the displacement ${mathbf u}$, the director ${mathbf n}$, the pressure $p$
and the Lagrange multiplier $lambda$. Inf-sup conditions for the
well-posedness of the linearized system were proposed and some are proved. For
those inf-sup conditions that are not easy to prove, we suggest ways to do the
numerical verification (the inf-sup tests). Finally, some numerical results are
presented.
Quasi-Lp norm orthogonal Galerkin expansions in sums of Jacobi polynomials
Matroid Polytopes and their Volumes
Fast Multiresolution Algorithms and Their Related Variational Problems for Image Denoising
An Iterative Method for Solving Non-Linear Hydromagnetic Equations. (arXiv:0911.5214v1 [math.NA])
We propose an iterative finite element method for solving non-linear
hydromagnetic and steady Euler's equations. Some three-dimensional
computational tests are given to confirm the convergence and the high
efficiency of the method.
Non-convexly constrained linear inverse problems. (arXiv:0911.5098v1 [math.NA])
This paper considers the inversion of ill-posed linear operators. To
regularise the problem the solution is enforced to lie in a non-convex subset.
Theoretical properties for the stable inversion are derived and an iterative
algorithm akin to the projected Landweber algorithm is studied. This work
extends recent progress made on the efficient inversion of finite dimensional
linear systems under a sparsity constraint to the Hilbert space setting and to
more general non-convex constraints.
domingo, 29 de noviembre de 2009
A new efficient technique for finding the solution of initial-value problems using He's variational iteration method
The well-known He's variational iteration method has been implemented by researchers for solving various kinds of problems. Although it is a very efficient method for solving problems, the approximation is not accurate in large regions. In this work a new technique based on the successive use of the He's variational iteration method in smaller regions has been presented. The new method has been named as local He's variational iteration method. By applying the new technique, an accurate approximation provided in large domains in a shorter time and less computations. Copyright © 2009 John Wiley & Sons, Ltd.
domingo, 22 de noviembre de 2009
Coloring plane graphs with independent crossings
We show that every plane graph with maximum face size four in which all faces of size four are vertex-disjoint is cyclically 5-colorable. This answers a question of Albertson whether graphs drawn in the plane with all crossings independent are 5-colorable. © 2009 Wiley Periodicals, Inc. J Graph Theory
Balloons, cut-edges, matchings, and total domination in regular graphs of odd degree
A balloon in a graph G is a maximal 2-edge-connected subgraph incident to exactly one cut-edge of G. Let b(G) be the number of balloons, let c(G) be the number of cut-edges, and let [alpha][prime](G) be the maximum size of a matching. Let be the family of connected (2r+1)-regular graphs with n vertices, and let . For , we prove the sharp inequalities c(G)[les][r(n-2)-2]/(2r2+2r-1)-1 and [alpha][prime](G)[ges]n/2-rb/(2r+1). Using b[les][(2r-1)n+2]/(4r2+4r-2), we obtain a simple proof of the bound proved by Henning and Yeo. For each of these bounds and each r, the approach using balloons allows us to determine the infinite family where equality holds. For the total domination number [gamma]t(G) of a cubic graph, we prove [gamma]t(G)[les]n/2-b(G)/2 (except that [gamma]t(G) may be n/2-1 when b(G)=3 and the balloons cover all but one vertex). With [alpha][prime](G)[ges]n/2-b(G)/3 for cubic graphs, this improves the known inequality [gamma]t(G)[les][alpha][prime](G). © 2009 Wiley Periodicals, Inc. J Graph Theory
Choosability of toroidal graphs without short cycles
Let G be a toroidal graph without cycles of a fixed length k, and [chi]l(G) the list chromatic number of G. We establish tight upper bounds of [chi]l(G) for the following values of k: © 2009 Wiley Periodicals, Inc. J Graph Theory
The geometry of t-spreads in k-walk-regular graphs
A graph is walk-regular if the number of closed walks of length [ell] rooted at a given vertex is a constant through all the vertices for all [ell]. For a walk-regular graph G with d+1 different eigenvalues and spectrally maximum diameter D=d, we study the geometry of its d-spreads, that is, the sets of vertices which are mutually at distance d. When these vertices are projected onto an eigenspace of its adjacency matrix, we show that they form a simplex (or tetrahedron in a three-dimensional case) and we compute its parameters. Moreover, the results are generalized to the case of k-walk-regular graphs, a family which includes both walk-regular and distance-regular graphs, and their t-spreads or vertices at distance t from each other. © 2009 Wiley Periodicals, Inc. J Graph Theory
On Ramsey-type positional games
Beck introduced the concept of Ramsey games by studying the game versions of Ramsey and van der Waerden theorems. We contribute to this topic by investigating games corresponding to structural extensions of Ramsey and van der Waerden theorems - the theorem of Brauer, structural and restricted Ramsey theorems. © 2009 Wiley Periodicals, Inc. J Graph Theory
On recognizing graphs by numbers of homomorphisms
Let hom (G, H) be the number of homomorphisms from a graph G to a graph H. A well-known result of Lovász states that the function hom (·, H) from all graphs uniquely determines the graph H up to isomorphism. We study this function restricted to smaller classes of graphs. We show that several natural classes (2-degenerate graphs and graphs homomorphic to an arbitrary non-bipartite graph) are sufficient to recognize all graphs, and provide a description of graph properties that are recognizable by graphs with bounded tree-width. © 2009 Wiley Periodicals, Inc. J Graph Theory
Brown, IBM Unveil Multimillion-Dollar Supercomputer
Brown University and IBM today announced the opening of a multimillion-dollar supercomputer at Brown's Center for Computation and Visualization.
Preservation and destruction of periodic orbits by symplectic integrators
sábado, 21 de noviembre de 2009
Erratum to: MIR closures of polyhedral sets
Erratum to: MIR closures of polyhedral sets
- Content Type Journal Article
- Category Erratum
- DOI 10.1007/s10107-009-0328-z
- Authors
- Sanjeeb Dash, IBM, T.J. Watson Research Center P.O. Box 218 Yorktown Heights NY 10598 USA
- Oktay Günlük, IBM, T.J. Watson Research Center P.O. Box 218 Yorktown Heights NY 10598 USA
- Andrea Lodi, DEIS, University of Bologna viale Risorgimento 2 40136 Bologna Italy
- Journal Mathematical Programming
- Online ISSN 1436-4646
- Print ISSN 0025-5610
Numerical Studies of Adaptive Finite Element Methods for Two Dimensional Convection-Dominated Problems
Gaussian quadrature rules using function derivatives
For finite positive Borel measures supported on the real line we consider a new type of quadrature rule with maximal algebraic degree of exactness that involves function derivatives. We prove the existence of such quadrature rules and describe their basic properties. We also give an application of these quadrature rules to the solution of a Cauchy problem without solving it directly. Numerical examples are included as well.
viernes, 20 de noviembre de 2009
Numerical integration of an elastic-viscoplastic constitutive model for dry metamorphosed snow
A constitutive model for dry metamorphosed snow is proposed, within the framework of elasto-viscoplasticity, which is able to reproduce the most relevant features of the macroscopic behaviour of snow, particularly its time and rate dependency. The basic ideas for modelling stem from the conceptual forms proposed for bonded geomaterials, such as cemented soils or soft rocks. The high viscosity of snow is accounted for by adopting an overstress approach, suitably modified. An evolution law for the curvature-driven process of sintering, by which intergranular ice necks form and grow, is established. The system of constitutive equations is then numerically integrated via a fully implicit time stepping scheme. Selected results from finite element simulations of laboratory tests, available in the literature, are presented. Copyright © 2009 John Wiley & Sons, Ltd.
On explicit two-derivative Runge-Kutta methods
A stabilized Lagrange multiplier method for the finite element approximation of contact problems in elastostatics
Geometry of Configuration Spaces of Tensegrities
lunes, 16 de noviembre de 2009
Optimal Quadrature Formulas with Positive Coefficients in $L_2^{(m)}(0,1)$ Space. (arXiv:0911.2896v1 [math.NA])
In the Sobolev space $L_2^{(m)}(0,1)$ optimal quadrature formulas with the
nodes (1.5) are investigated. For optimal coefficients explicit form are
obtained and norm of the error functional is calculated. In particular, by
choosing parameter $eta_0$ in (1.5) the optimal quadrature formulas with
positive coefficients are obtained and compared with well known optimal
formulas.
miércoles, 11 de noviembre de 2009
The Penalized Lebesgue Constant for Surface Spline Interpolation. (arXiv:0911.1815v1 [math.CA])
Problems involving approximation from scattered data where data is arranged
quasi-uniformly have been treated by RBF methods for decades. Treating data
with spatially varying density has not been investigated with the same
intensity, and is far less well understood. In this article we consider the
stability of surface spline interpolation (a popular type of RBF interpolation)
for data with nonuniform arrangements. Using techniques similar to those
recently employed by Hangelbroek, Narcowich and Ward to demonstrate the
stability of interpolation from quasi-uniform data on manifolds, we show that
surface spline interpolation on R^d is stable, but in a stronger, local sense.
We also obtain pointwise estimates showing that the Lagrange function decays
very rapidly, and at a rate determined by the local spacing of datasites. These
results, in conjunction with a Lebesgue lemma, show that surface spline
interpolation enjoys the same rates of convergence as those of the local
approximation schemes recently developed by DeVore and Ron.
A convergent mixed method for the Stokes approximation of viscous compressible flow. (arXiv:0911.1870v1 [math.NA])
We propose a mixed finite element method for the motion of a strongly
viscous, ideal, and isentropic gas. At the boundary we impose a Navier-slip
condition such that the velocity equation can be posed in mixed form with the
vorticity as an auxiliary variable. In this formulation we design a finite
element method, where the velocity and vorticity is approximated with the div-
and curl- conforming Nedelec elements, respectively, of the first order and
first kind. The mixed scheme is coupled to a standard piecewise constant upwind
discontinuous Galerkin discretization of the continuity equation. For the time
discretization, implicit Euler time stepping is used. Our main result is that
the numerical solution converges to a weak solution as the discretization
parameters go to zero. The convergence analysis is inspired by the continuous
analysis of Feireisl and Lions for the compressible Navier-Stokes equations.
Tools used in the analysis include an equation for the effective viscous flux
and various renormalizations of the density scheme.
martes, 10 de noviembre de 2009
Symmetry-restrained molecular dynamics simulations improve homology models of potassium channels
Most crystallized homo-oligomeric ion channels are highly symmetric, which dramatically decreases conformational space and facilitates building homology models (HMs). However, in molecular dynamics (MD) simulations channels deviate from ideal symmetry and accumulate thermal defects, which complicate the refinement of HMs using MD. In this work we evaluate the ability of symmetry constrained MD simulations to improve HMs accuracy, using an approach conceptually similar to Critical Assessment of techniques for protein Structure Prediction (CASP) competition: build HMs of channels with known structure and evaluate the efficiency of proposed methods in improving HMs accuracy (measured as deviation from experimental structure). Results indicate that unrestrained MD does not improve the accuracy of HMs, instantaneous symmetrization improves accuracy but not stability of HMs during subsequent unrestrained MD, while gradually imposing symmetry constraints improves both accuracy (by 5-50%) and stability of HMs. Moreover, accuracy and stability are strongly correlated, making stability a reliable criterion in predicting the accuracy of new HMs. Proteins 2010. © 2009 Wiley-Liss, Inc.
lunes, 9 de noviembre de 2009
Waveform Transmission Method, a New Waveform-relaxation Based Algorithm to Solve Ordinary Differential Equations in Parallel. (arXiv:0911.1166v1 [math.NA])
Waveform Relaxation method (WR) is a beautiful algorithm to solve Ordinary
Differential Equations (ODEs). However, because of its poor convergence
capability, it was rarely used. Virtual Transmission Method (VTM) is a new
distributed algorithm to solve large sparse linear systems, with an appreciable
convergence speed. In this paper, by marrying WR with VTM, a new distributed
algorithm, named Waveform Transmission Method (WTM), was born. WTM has better
convergence capability than the traditional WR algorithms.
Selection of Corners for the BDDC Method. (arXiv:0911.1245v1 [math.NA])
The Balancing Domain Decomposition by Constraints (BDDC) method has evolved
quite fast since its introduction in 2003, as the primal counterpart to the
earlier FETI-DP method. Recent results have shown close connection of these
methods and theoretically supported equivalent rate of convergence. In both
methods, a fundamental role is played by the coarse space. Optimal choice of
constraints on continuity of the coarse space is still not a satisfactorily
solved problem. The usual basic choice is a 'minimal' set of coarse nodes
(sometimes called corners), that assures invertibility of local subdomain
problems and also of the global coarse problem. However, this set alone does
not suffice for optimal preconditioning in 3D. For this reason, continuity of
some generalized degrees of freedom, such as average values on edges or faces
of subdomains, have to be added. While theoretically correct, this approach
does not easily offer a flexible size of desired coarse problem. In our
contribution, we compare this approach with adding more coarse nodes into the
coarse problem, which is technically simpler and allows flexible setting of
desired approximation. Most attention is given to the basic set of corners, and
a new algorithm for its selection is proposed.