WorldWideScience

Sample records for approximation

  1. Diophantine approximations

    CERN Document Server

    Niven, Ivan

    2008-01-01

    This self-contained treatment originated as a series of lectures delivered to the Mathematical Association of America. It covers basic results on homogeneous approximation of real numbers; the analogue for complex numbers; basic results for nonhomogeneous approximation in the real case; the analogue for complex numbers; and fundamental properties of the multiples of an irrational number, for both the fractional and integral parts.The author refrains from the use of continuous fractions and includes basic results in the complex case, a feature often neglected in favor of the real number discuss

  2. Approximation Techniques

    OpenAIRE

    2006-01-01

    This interactive tutorial presents the following concepts of Approximation Techniques: Methods of Weighted Residual (MWR), Weak Formulatioin, Piecewise Continuous Function, Galerkin Finite Element FormulationExplanations especially for mathematical statements are provided using mouseover the highlight equations. ME4613 Finite Element Methods

  3. Approximate Representations and Approximate Homomorphisms

    CERN Document Server

    Moore, Cristopher

    2010-01-01

    Approximate algebraic structures play a defining role in arithmetic combinatorics and have found remarkable applications to basic questions in number theory and pseudorandomness. Here we study approximate representations of finite groups: functions f:G -> U_d such that Pr[f(xy) = f(x) f(y)] is large, or more generally Exp_{x,y} ||f(xy) - f(x)f(y)||^2$ is small, where x and y are uniformly random elements of the group G and U_d denotes the unitary group of degree d. We bound these quantities in terms of the ratio d / d_min where d_min is the dimension of the smallest nontrivial representation of G. As an application, we bound the extent to which a function f : G -> H can be an approximate homomorphism where H is another finite group. We show that if H's representations are significantly smaller than G's, no such f can be much more homomorphic than a random function. We interpret these results as showing that if G is quasirandom, that is, if d_min is large, then G cannot be embedded in a small number of dimensi...

  4. Approximate Likelihood

    CERN Document Server

    CERN. Geneva

    2015-01-01

    Most physics results at the LHC end in a likelihood ratio test. This includes discovery and exclusion for searches as well as mass, cross-section, and coupling measurements. The use of Machine Learning (multivariate) algorithms in HEP is mainly restricted to searches, which can be reduced to classification between two fixed distributions: signal vs. background. I will show how we can extend the use of ML classifiers to distributions parameterized by physical quantities like masses and couplings as well as nuisance parameters associated to systematic uncertainties. This allows for one to approximate the likelihood ratio while still using a high dimensional feature vector for the data. Both the MEM and ABC approaches mentioned above aim to provide inference on model parameters (like cross-sections, masses, couplings, etc.). ABC is fundamentally tied Bayesian inference and focuses on the “likelihood free” setting where only a simulator is available and one cannot directly compute the likelihood for the dat...

  5. Diophantine approximation

    CERN Document Server

    Schmidt, Wolfgang M

    1980-01-01

    "In 1970, at the U. of Colorado, the author delivered a course of lectures on his famous generalization, then just established, relating to Roth's theorem on rational approxi- mations to algebraic numbers. The present volume is an ex- panded and up-dated version of the original mimeographed notes on the course. As an introduction to the author's own remarkable achievements relating to the Thue-Siegel-Roth theory, the text can hardly be bettered and the tract can already be regarded as a classic in its field."(Bull.LMS) "Schmidt's work on approximations by algebraic numbers belongs to the deepest and most satisfactory parts of number theory. These notes give the best accessible way to learn the subject. ... this book is highly recommended." (Mededelingen van het Wiskundig Genootschap)

  6. Diophantine approximation and badly approximable sets

    DEFF Research Database (Denmark)

    Kristensen, S.; Thorn, R.; Velani, S.

    2006-01-01

    . The classical set Bad of `badly approximable' numbers in the theory of Diophantine approximation falls within our framework as do the sets Bad(i,j) of simultaneously badly approximable numbers. Under various natural conditions we prove that the badly approximable subsets of Omega have full Hausdorff dimension....... Applications of our general framework include those from number theory (classical, complex, p-adic and formal power series) and dynamical systems (iterated function schemes, rational maps and Kleinian groups)....

  7. Optimal Belief Approximation

    CERN Document Server

    Leike, Reimar H

    2016-01-01

    In Bayesian statistics probability distributions express beliefs. However, for many problems the beliefs cannot be computed analytically and approximations of beliefs are needed. We seek a ranking function that quantifies how "embarrassing" it is to communicate a given approximation. We show that there is only one ranking under the requirements that (1) the best ranked approximation is the non-approximated belief and (2) that the ranking judges approximations only by their predictions for actual outcomes. We find that this ranking is equivalent to the Kullback-Leibler divergence that is frequently used in the literature. However, there seems to be confusion about the correct order in which its functional arguments, the approximated and non-approximated beliefs, should be used. We hope that our elementary derivation settles the apparent confusion. We show for example that when approximating beliefs with Gaussian distributions the optimal approximation is given by moment matching. This is in contrast to many su...

  8. Fuzzy Approximating Spaces

    OpenAIRE

    Bin Qin

    2014-01-01

    Relationships between fuzzy relations and fuzzy topologies are deeply researched. The concept of fuzzy approximating spaces is introduced and decision conditions that a fuzzy topological space is a fuzzy approximating space are obtained.

  9. Approximate flavor symmetries

    CERN Document Server

    Rasin, A

    1994-01-01

    We discuss the idea of approximate flavor symmetries. Relations between approximate flavor symmetries and natural flavor conservation and democracy models is explored. Implications for neutrino physics are also discussed.

  10. Approximate iterative algorithms

    CERN Document Server

    Almudevar, Anthony Louis

    2014-01-01

    Iterative algorithms often rely on approximate evaluation techniques, which may include statistical estimation, computer simulation or functional approximation. This volume presents methods for the study of approximate iterative algorithms, providing tools for the derivation of error bounds and convergence rates, and for the optimal design of such algorithms. Techniques of functional analysis are used to derive analytical relationships between approximation methods and convergence properties for general classes of algorithms. This work provides the necessary background in functional analysis a

  11. Approximation of distributed delays

    CERN Document Server

    Lu, Hao; Eberard, Damien; Simon, Jean-Pierre

    2010-01-01

    We address in this paper the approximation problem of distributed delays. Such elements are convolution operators with kernel having bounded support, and appear in the control of time-delay systems. From the rich literature on this topic, we propose a general methodology to achieve such an approximation. For this, we enclose the approximation problem in the graph topology, and work with the norm defined over the convolution Banach algebra. The class of rational approximates is described, and a constructive approximation is proposed. Analysis in time and frequency domains is provided. This methodology is illustrated on the stabilization control problem, for which simulations results show the effectiveness of the proposed methodology.

  12. Sparse approximation with bases

    CERN Document Server

    2015-01-01

    This book systematically presents recent fundamental results on greedy approximation with respect to bases. Motivated by numerous applications, the last decade has seen great successes in studying nonlinear sparse approximation. Recent findings have established that greedy-type algorithms are suitable methods of nonlinear approximation in both sparse approximation with respect to bases and sparse approximation with respect to redundant systems. These insights, combined with some previous fundamental results, form the basis for constructing the theory of greedy approximation. Taking into account the theoretical and practical demand for this kind of theory, the book systematically elaborates a theoretical framework for greedy approximation and its applications.  The book addresses the needs of researchers working in numerical mathematics, harmonic analysis, and functional analysis. It quickly takes the reader from classical results to the latest frontier, but is written at the level of a graduate course and do...

  13. AN OPTIMAL FUZZY APPROXIMATOR

    Institute of Scientific and Technical Information of China (English)

    YueShihong; ZhangKecun

    2002-01-01

    In a dot product space with the reproducing kernel (r. k. S. ) ,a fuzzy system with the estimation approximation errors is proposed ,which overcomes the defect that the existing fuzzy control system is difficult to estimate the errors of approximation for a desired function,and keeps the characteristics of fuzzy system as an inference approach. The structure of the new fuzzy approximator benefits a course got by other means.

  14. Approximation techniques for engineers

    CERN Document Server

    Komzsik, Louis

    2006-01-01

    Presenting numerous examples, algorithms, and industrial applications, Approximation Techniques for Engineers is your complete guide to the major techniques used in modern engineering practice. Whether you need approximations for discrete data of continuous functions, or you''re looking for approximate solutions to engineering problems, everything you need is nestled between the covers of this book. Now you can benefit from Louis Komzsik''s years of industrial experience to gain a working knowledge of a vast array of approximation techniques through this complete and self-contained resource.

  15. Expectation Consistent Approximate Inference

    DEFF Research Database (Denmark)

    Opper, Manfred; Winther, Ole

    2005-01-01

    We propose a novel framework for approximations to intractable probabilistic models which is based on a free energy formulation. The approximation can be understood from replacing an average over the original intractable distribution with a tractable one. It requires two tractable probability dis...

  16. Ordered cones and approximation

    CERN Document Server

    Keimel, Klaus

    1992-01-01

    This book presents a unified approach to Korovkin-type approximation theorems. It includes classical material on the approximation of real-valuedfunctions as well as recent and new results on set-valued functions and stochastic processes, and on weighted approximation. The results are notonly of qualitative nature, but include quantitative bounds on the order of approximation. The book is addressed to researchers in functional analysis and approximation theory as well as to those that want to applythese methods in other fields. It is largely self- contained, but the readershould have a solid background in abstract functional analysis. The unified approach is based on a new notion of locally convex ordered cones that are not embeddable in vector spaces but allow Hahn-Banach type separation and extension theorems. This concept seems to be of independent interest.

  17. Approximate Modified Policy Iteration

    CERN Document Server

    Scherrer, Bruno; Ghavamzadeh, Mohammad; Geist, Matthieu

    2012-01-01

    Modified policy iteration (MPI) is a dynamic programming (DP) algorithm that contains the two celebrated policy and value iteration methods. Despite its generality, MPI has not been thoroughly studied, especially its approximation form which is used when the state and/or action spaces are large or infinite. In this paper, we propose three approximate MPI (AMPI) algorithms that are extensions of the well-known approximate DP algorithms: fitted-value iteration, fitted-Q iteration, and classification-based policy iteration. We provide an error propagation analysis for AMPI that unifies those for approximate policy and value iteration. We also provide a finite-sample analysis for the classification-based implementation of AMPI (CBMPI), which is more general (and somehow contains) than the analysis of the other presented AMPI algorithms. An interesting observation is that the MPI's parameter allows us to control the balance of errors (in value function approximation and in estimating the greedy policy) in the fina...

  18. The Karlqvist approximation revisited

    OpenAIRE

    Tannous, C.

    2015-01-01

    The Karlqvist approximation signaling the historical beginning of magnetic recording head theory is reviewed and compared to various approaches progressing from Green, Fourier, Conformal mapping that obeys the Sommerfeld edge condition at angular points and leads to exact results.

  19. Approximations in Inspection Planning

    DEFF Research Database (Denmark)

    Engelund, S.; Sørensen, John Dalsgaard; Faber, M. H.;

    2000-01-01

    Planning of inspections of civil engineering structures may be performed within the framework of Bayesian decision analysis. The effort involved in a full Bayesian decision analysis is relatively large. Therefore, the actual inspection planning is usually performed using a number of approximations....... One of the more important of these approximations is the assumption that all inspections will reveal no defects. Using this approximation the optimal inspection plan may be determined on the basis of conditional probabilities, i.e. the probability of failure given no defects have been found...... by the inspection. In this paper the quality of this approximation is investigated. The inspection planning is formulated both as a full Bayesian decision problem and on the basis of the assumption that the inspection will reveal no defects....

  20. Approximation of irrationals

    Directory of Open Access Journals (Sweden)

    Malvina Baica

    1985-01-01

    Full Text Available The author uses a new modification of Jacobi-Perron Algorithm which holds for complex fields of any degree (abbr. ACF, and defines it as Generalized Euclidean Algorithm (abbr. GEA to approximate irrationals.

  1. Approximation Behooves Calibration

    DEFF Research Database (Denmark)

    da Silva Ribeiro, André Manuel; Poulsen, Rolf

    2013-01-01

    Calibration based on an expansion approximation for option prices in the Heston stochastic volatility model gives stable, accurate, and fast results for S&P500-index option data over the period 2005–2009.......Calibration based on an expansion approximation for option prices in the Heston stochastic volatility model gives stable, accurate, and fast results for S&P500-index option data over the period 2005–2009....

  2. Approximation and Computation

    CERN Document Server

    Gautschi, Walter; Rassias, Themistocles M

    2011-01-01

    Approximation theory and numerical analysis are central to the creation of accurate computer simulations and mathematical models. Research in these areas can influence the computational techniques used in a variety of mathematical and computational sciences. This collection of contributed chapters, dedicated to renowned mathematician Gradimir V. Milovanovia, represent the recent work of experts in the fields of approximation theory and numerical analysis. These invited contributions describe new trends in these important areas of research including theoretic developments, new computational alg

  3. Diophantine approximations on fractals

    CERN Document Server

    Einsiedler, Manfred; Shapira, Uri

    2009-01-01

    We exploit dynamical properties of diagonal actions to derive results in Diophantine approximations. In particular, we prove that the continued fraction expansion of almost any point on the middle third Cantor set (with respect to the natural measure) contains all finite patterns (hence is well approximable). Similarly, we show that for a variety of fractals in [0,1]^2, possessing some symmetry, almost any point is not Dirichlet improvable (hence is well approximable) and has property C (after Cassels). We then settle by similar methods a conjecture of M. Boshernitzan saying that there are no irrational numbers x in the unit interval such that the continued fraction expansions of {nx mod1 : n is a natural number} are uniformly eventually bounded.

  4. Approximate and Incomplete Factorizations

    NARCIS (Netherlands)

    Chan, T.F.; Vorst, H.A. van der

    2001-01-01

    In this chapter, we give a brief overview of a particular class of preconditioners known as incomplete factorizations. They can be thought of as approximating the exact LU factorization of a given matrix A (e.g. computed via Gaussian elimination) by disallowing certain ll-ins. As opposed to other PD

  5. Prestack wavefield approximations

    KAUST Repository

    Alkhalifah, Tariq

    2013-09-01

    The double-square-root (DSR) relation offers a platform to perform prestack imaging using an extended single wavefield that honors the geometrical configuration between sources, receivers, and the image point, or in other words, prestack wavefields. Extrapolating such wavefields, nevertheless, suffers from limitations. Chief among them is the singularity associated with horizontally propagating waves. I have devised highly accurate approximations free of such singularities which are highly accurate. Specifically, I use Padé expansions with denominators given by a power series that is an order lower than that of the numerator, and thus, introduce a free variable to balance the series order and normalize the singularity. For the higher-order Padé approximation, the errors are negligible. Additional simplifications, like recasting the DSR formula as a function of scattering angle, allow for a singularity free form that is useful for constant-angle-gather imaging. A dynamic form of this DSR formula can be supported by kinematic evaluations of the scattering angle to provide efficient prestack wavefield construction. Applying a similar approximation to the dip angle yields an efficient 1D wave equation with the scattering and dip angles extracted from, for example, DSR ray tracing. Application to the complex Marmousi data set demonstrates that these approximations, although they may provide less than optimal results, allow for efficient and flexible implementations. © 2013 Society of Exploration Geophysicists.

  6. The Zeldovich approximation

    CERN Document Server

    White, Martin

    2014-01-01

    This year marks the 100th anniversary of the birth of Yakov Zel'dovich. Amongst his many legacies is the Zel'dovich approximation for the growth of large-scale structure, which remains one of the most successful and insightful analytic models of structure formation. We use the Zel'dovich approximation to compute the two-point function of the matter and biased tracers, and compare to the results of N-body simulations and other Lagrangian perturbation theories. We show that Lagrangian perturbation theories converge well and that the Zel'dovich approximation provides a good fit to the N-body results except for the quadrupole moment of the halo correlation function. We extend the calculation of halo bias to 3rd order and also consider non-local biasing schemes, none of which remove the discrepancy. We argue that a part of the discrepancy owes to an incorrect prediction of inter-halo velocity correlations. We use the Zel'dovich approximation to compute the ingredients of the Gaussian streaming model and show that ...

  7. Approximating The DCM

    DEFF Research Database (Denmark)

    Madsen, Rasmus Elsborg

    2005-01-01

    The Dirichlet compound multinomial (DCM), which has recently been shown to be well suited for modeling for word burstiness in documents, is here investigated. A number of conceptual explanations that account for these recent results, are provided. An exponential family approximation of the DCM that...

  8. Prestack traveltime approximations

    KAUST Repository

    Alkhalifah, Tariq Ali

    2011-01-01

    Most prestack traveltime relations we tend work with are based on homogeneous (or semi-homogenous, possibly effective) media approximations. This includes the multi-focusing or double square-root (DSR) and the common reflection stack (CRS) equations. Using the DSR equation, I analyze the associated eikonal form in the general source-receiver domain. Like its wave-equation counterpart, it suffers from a critical singularity for horizontally traveling waves. As a result, I derive expansion based solutions of this eikonal based on polynomial expansions in terms of the reflection and dip angles in a generally inhomogenous background medium. These approximate solutions are free of singularities and can be used to estimate travetimes for small to moderate offsets (or reflection angles) in a generally inhomogeneous medium. A Marmousi example demonstrates the usefulness of the approach. © 2011 Society of Exploration Geophysicists.

  9. Approximate option pricing

    Energy Technology Data Exchange (ETDEWEB)

    Chalasani, P.; Saias, I. [Los Alamos National Lab., NM (United States); Jha, S. [Carnegie Mellon Univ., Pittsburgh, PA (United States)

    1996-04-08

    As increasingly large volumes of sophisticated options (called derivative securities) are traded in world financial markets, determining a fair price for these options has become an important and difficult computational problem. Many valuation codes use the binomial pricing model, in which the stock price is driven by a random walk. In this model, the value of an n-period option on a stock is the expected time-discounted value of the future cash flow on an n-period stock price path. Path-dependent options are particularly difficult to value since the future cash flow depends on the entire stock price path rather than on just the final stock price. Currently such options are approximately priced by Monte carlo methods with error bounds that hold only with high probability and which are reduced by increasing the number of simulation runs. In this paper the authors show that pricing an arbitrary path-dependent option is {number_sign}-P hard. They show that certain types f path-dependent options can be valued exactly in polynomial time. Asian options are path-dependent options that are particularly hard to price, and for these they design deterministic polynomial-time approximate algorithms. They show that the value of a perpetual American put option (which can be computed in constant time) is in many cases a good approximation to the value of an otherwise identical n-period American put option. In contrast to Monte Carlo methods, the algorithms have guaranteed error bounds that are polynormally small (and in some cases exponentially small) in the maturity n. For the error analysis they derive large-deviation results for random walks that may be of independent interest.

  10. Approximations to Euler's constant

    International Nuclear Information System (INIS)

    We study a problem of finding good approximations to Euler's constant γ=lim→∞ Sn, where Sn = Σk=Ln (1)/k-log(n+1), by linear forms in logarithms and harmonic numbers. In 1995, C. Elsner showed that slow convergence of the sequence Sn can be significantly improved if Sn is replaced by linear combinations of Sn with integer coefficients. In this paper, considering more general linear transformations of the sequence Sn we establish new accelerating convergence formulae for γ. Our estimates sharpen and generalize recent Elsner's, Rivoal's and author's results. (author)

  11. Finite elements and approximation

    CERN Document Server

    Zienkiewicz, O C

    2006-01-01

    A powerful tool for the approximate solution of differential equations, the finite element is extensively used in industry and research. This book offers students of engineering and physics a comprehensive view of the principles involved, with numerous illustrative examples and exercises.Starting with continuum boundary value problems and the need for numerical discretization, the text examines finite difference methods, weighted residual methods in the context of continuous trial functions, and piecewise defined trial functions and the finite element method. Additional topics include higher o

  12. The Compact Approximation Property does not imply the Approximation Property

    OpenAIRE

    Willis, George A.

    1992-01-01

    It is shown how to construct, given a Banach space which does not have the approximation property, another Banach space which does not have the approximation property but which does have the compact approximation property.

  13. Interacting boson approximation

    International Nuclear Information System (INIS)

    Lectures notes on the Interacting Boson Approximation are given. Topics include: angular momentum tensors; properties of T/sub i//sup (n)/ matrices; T/sub i//sup (n)/ matrices as Clebsch-Gordan coefficients; construction of higher rank tensors; normalization: trace of products of two s-rank tensors; completeness relation; algebra of U(N); eigenvalue of the quadratic Casimir operator for U(3); general result for U(N); angular momentum content of U(3) representation; p-Boson model; Hamiltonian; quadrupole transitions; S,P Boson model; expectation value of dipole operator; S-D model: U(6); quadratic Casimir operator; an O(5) subgroup; an O(6) subgroup; properties of O(5) representations; quadratic Casimir operator; quadratic Casimir operator for U(6); decomposition via SU(5) chain; a special O(3) decomposition of SU(3); useful identities; a useful property of D/sub αβγ/(α,β,γ = 4-8) as coupling coefficients; explicit construction of T/sub x//sup (2)/ and d/sub αβγ/; D-coefficients; eigenstates of T3; and summary of T = 2 states

  14. Prestack traveltime approximations

    KAUST Repository

    Alkhalifah, Tariq Ali

    2012-05-01

    Many of the explicit prestack traveltime relations used in practice are based on homogeneous (or semi-homogenous, possibly effective) media approximations. This includes the multifocusing, based on the double square-root (DSR) equation, and the common reflection stack (CRS) approaches. Using the DSR equation, I constructed the associated eikonal form in the general source-receiver domain. Like its wave-equation counterpart, it suffers from a critical singularity for horizontally traveling waves. As a result, I recasted the eikonal in terms of the reflection angle, and thus, derived expansion based solutions of this eikonal in terms of the difference between the source and receiver velocities in a generally inhomogenous background medium. The zero-order term solution, corresponding to ignoring the lateral velocity variation in estimating the prestack part, is free of singularities and can be used to estimate traveltimes for small to moderate offsets (or reflection angles) in a generally inhomogeneous medium. The higher-order terms include limitations for horizontally traveling waves, however, we can readily enforce stability constraints to avoid such singularities. In fact, another expansion over reflection angle can help us avoid these singularities by requiring the source and receiver velocities to be different. On the other hand, expansions in terms of reflection angles result in singularity free equations. For a homogenous background medium, as a test, the solutions are reasonably accurate to large reflection and dip angles. A Marmousi example demonstrated the usefulness and versatility of the formulation. © 2012 Society of Exploration Geophysicists.

  15. Operators of Approximations and Approximate Power Set Spaces

    Institute of Scientific and Technical Information of China (English)

    ZHANG Xian-yong; MO Zhi-wen; SHU Lan

    2004-01-01

    Boundary inner and outer operators are introduced; and union, intersection, complement operators of approximations are redefined. The approximation operators have a good property of maintaining union, intersection, complement operators, so the rough set theory has been enriched from the operator-oriented and set-oriented views. Approximate power set spaces are defined, and it is proved that the approximation operators are epimorphisms from power set space to approximate power set spaces. Some basic properties of approximate power set space are got by epimorphisms in contrast to power set space.

  16. Approximation algorithms and hardness of approximation for knapsack problems

    NARCIS (Netherlands)

    Buhrman, H.; Loff, B.; Torenvliet, L.

    2012-01-01

    We show various hardness of approximation algorithms for knapsack and related problems; in particular we will show that unless the Exponential-Time Hypothesis is false, then subset-sum cannot be approximated any better than with an FPTAS. We also give a simple new algorithm for approximating knapsac

  17. Approximation by planar elastic curves

    DEFF Research Database (Denmark)

    Brander, David; Gravesen, Jens; Nørbjerg, Toke Bjerge

    2015-01-01

    We give an algorithm for approximating a given plane curve segment by a planar elastic curve. The method depends on an analytic representation of the space of elastic curve segments, together with a geometric method for obtaining a good initial guess for the approximating curve. A gradient......-driven optimization is then used to find the approximating elastic curve....

  18. Approximate sine-Gordon solitons

    Energy Technology Data Exchange (ETDEWEB)

    Stratopoulos, G.N. (Dept. of Mathematical Sciences, Durham Univ. (United Kingdom)); Zakrzewski, W.J. (Dept. of Mathematical Sciences, Durham Univ. (United Kingdom))

    1993-08-01

    We look at the recently proposed scheme of approximating a sine-Gordon soliton by an expression derived from two dimensional instantons. We point out that the scheme of Sutcliffe in which he uses two dimensional instantons can be generalised to higher dimensions and that these generalisations produce even better approximations than the original approximation. We also comment on generalisations to other models. (orig.)

  19. Exact constants in approximation theory

    CERN Document Server

    Korneichuk, N

    1991-01-01

    This book is intended as a self-contained introduction for non-specialists, or as a reference work for experts, to the particular area of approximation theory that is concerned with exact constants. The results apply mainly to extremal problems in approximation theory, which in turn are closely related to numerical analysis and optimization. The book encompasses a wide range of questions and problems: best approximation by polynomials and splines; linear approximation methods, such as spline-approximation; optimal reconstruction of functions and linear functionals. Many of the results are base

  20. Approximate solutions for the skyrmion

    CERN Document Server

    Ponciano, J A; Fanchiotti, H; Canal-Garcia, C A

    2001-01-01

    We reconsider the Euler-Lagrange equation for the Skyrme model in the hedgehog ansatz and study the analytical properties of the solitonic solution. In view of the lack of a closed form solution to the problem, we work on approximate analytical solutions. We show that Pade approximants are well suited to continue analytically the asymptotic representation obtained in terms of a power series expansion near the origin, obtaining explicit approximate solutions for the Skyrme equations. We improve the approximations by applying the 2-point Pade approximant procedure whereby the exact behaviour at spatial infinity is incorporated. An even better convergence to the exact solution is obtained by introducing a modified form for the approximants. The new representations share the same analytical properties with the exact solution at both small and large values of the radial variable r.

  1. BDD Minimization for Approximate Computing

    OpenAIRE

    Soeken, Mathias; Grosse, Daniel; Chandrasekharan, Arun; Drechsler, Rolf

    2016-01-01

    We present Approximate BDD Minimization (ABM) as a problem that has application in approximate computing. Given a BDD representation of a multi-output Boolean function, ABM asks whether there exists another function that has a smaller BDD representation but meets a threshold w.r.t. an error metric. We present operators to derive approximated functions and present algorithms to exactly compute the error metrics directly on the BDD representation. An experimental evaluation demonstrates the app...

  2. Nonlinear approximation with redundant dictionaries

    DEFF Research Database (Denmark)

    Borup, Lasse; Nielsen, M.; Gribonval, R.

    2005-01-01

    In this paper we study nonlinear approximation and data representation with redundant function dictionaries. In particular, approximation with redundant wavelet bi-frame systems is studied in detail. Several results for orthonormal wavelets are generalized to the redundant case. In general...

  3. Approximate Inference in Probabilistic Models

    DEFF Research Database (Denmark)

    Opper, Manfred; Winther, Ole

    2004-01-01

    We present a framework for approximate inference in probabilistic data models which is based on free energies. The free energy is constructed from two approximating distributions which encode different aspects of the intractable model. Consistency between distributions is required on a chosen set...

  4. The Logic of Approximate Dependence

    OpenAIRE

    Väänänen, Jouko

    2014-01-01

    We extend the treatment of functional dependence, the basic concept of dependence logic, to include the possibility of dependence with a limited number of exceptions. We call this approximate dependence. The main result of the paper is a Completeness Theorem for approximate dependence atoms. We point out some problematic features of this which suggests that we should consider multi-teams, not just teams.

  5. Beyond the random phase approximation

    DEFF Research Database (Denmark)

    Olsen, Thomas; Thygesen, Kristian S.

    2013-01-01

    We assess the performance of a recently proposed renormalized adiabatic local density approximation (rALDA) for ab initio calculations of electronic correlation energies in solids and molecules. The method is an extension of the random phase approximation (RPA) derived from time-dependent density...

  6. Approximate circuits for increased reliability

    Energy Technology Data Exchange (ETDEWEB)

    Hamlet, Jason R.; Mayo, Jackson R.

    2015-12-22

    Embodiments of the invention describe a Boolean circuit having a voter circuit and a plurality of approximate circuits each based, at least in part, on a reference circuit. The approximate circuits are each to generate one or more output signals based on values of received input signals. The voter circuit is to receive the one or more output signals generated by each of the approximate circuits, and is to output one or more signals corresponding to a majority value of the received signals. At least some of the approximate circuits are to generate an output value different than the reference circuit for one or more input signal values; however, for each possible input signal value, the majority values of the one or more output signals generated by the approximate circuits and received by the voter circuit correspond to output signal result values of the reference circuit.

  7. Approximate circuits for increased reliability

    Energy Technology Data Exchange (ETDEWEB)

    Hamlet, Jason R.; Mayo, Jackson R.

    2015-08-18

    Embodiments of the invention describe a Boolean circuit having a voter circuit and a plurality of approximate circuits each based, at least in part, on a reference circuit. The approximate circuits are each to generate one or more output signals based on values of received input signals. The voter circuit is to receive the one or more output signals generated by each of the approximate circuits, and is to output one or more signals corresponding to a majority value of the received signals. At least some of the approximate circuits are to generate an output value different than the reference circuit for one or more input signal values; however, for each possible input signal value, the majority values of the one or more output signals generated by the approximate circuits and received by the voter circuit correspond to output signal result values of the reference circuit.

  8. N-variable rational approximants

    International Nuclear Information System (INIS)

    ''Desirable properties'' of a two-variable generalization of Pade approximants are laid down. The ''Chisholm approximants'' are defined and are shown to obey nearly all of these properties; the alternative ways of completing a unique definition are discussed, and the ''prong structure'' of the defining equations is elucidated. Several generalizations and variants of Chisholm approximants are described: N-variable diagonal, 2-variable simple off-diagonal, N-variable simple and general off-diagonal, and rotationally covariant 2-variable approximants. All of the 2-variable approximants are capable of representing singularities of functions of two variables, and of analytically continuing beyond the polycylinder of convergence of the double series. 8 figures

  9. The efficiency of Flory approximation

    International Nuclear Information System (INIS)

    The Flory approximation for the self-avoiding chain problem is compared with a conventional perturbation theory expansion. While in perturbation theory each term is averaged over the unperturbed set of configurations, the Flory approximation is equivalent to the perturbation theory with the averaging over the stretched set of configurations. This imposes restrictions on the integration domain in higher order terms and they can be treated self-consistently. The accuracy δν/ν of Flory approximation for self-avoiding chain problems is estimated to be 2-5% for 1 < d < 4. (orig.)

  10. Binary nucleation beyond capillarity approximation

    NARCIS (Netherlands)

    Kalikmanov, V.I.

    2010-01-01

    Large discrepancies between binary classical nucleation theory (BCNT) and experiments result from adsorption effects and inability of BCNT, based on the phenomenological capillarity approximation, to treat small clusters. We propose a model aimed at eliminating both of these deficiencies. Adsorption

  11. Approximate Implicitization Using Linear Algebra

    Directory of Open Access Journals (Sweden)

    Oliver J. D. Barrowclough

    2012-01-01

    Full Text Available We consider a family of algorithms for approximate implicitization of rational parametric curves and surfaces. The main approximation tool in all of the approaches is the singular value decomposition, and they are therefore well suited to floating-point implementation in computer-aided geometric design (CAGD systems. We unify the approaches under the names of commonly known polynomial basis functions and consider various theoretical and practical aspects of the algorithms. We offer new methods for a least squares approach to approximate implicitization using orthogonal polynomials, which tend to be faster and more numerically stable than some existing algorithms. We propose several simple propositions relating the properties of the polynomial bases to their implicit approximation properties.

  12. Weighted approximation with varying weight

    CERN Document Server

    Totik, Vilmos

    1994-01-01

    A new construction is given for approximating a logarithmic potential by a discrete one. This yields a new approach to approximation with weighted polynomials of the form w"n"(" "= uppercase)P"n"(" "= uppercase). The new technique settles several open problems, and it leads to a simple proof for the strong asymptotics on some L p(uppercase) extremal problems on the real line with exponential weights, which, for the case p=2, are equivalent to power- type asymptotics for the leading coefficients of the corresponding orthogonal polynomials. The method is also modified toyield (in a sense) uniformly good approximation on the whole support. This allows one to deduce strong asymptotics in some L p(uppercase) extremal problems with varying weights. Applications are given, relating to fast decreasing polynomials, asymptotic behavior of orthogonal polynomials and multipoint Pade approximation. The approach is potential-theoretic, but the text is self-contained.

  13. Approximate maximizers of intricacy functionals

    CERN Document Server

    Buzzi, Jerome

    2009-01-01

    G. Edelman, O. Sporns, and G. Tononi introduced in theoretical biology the neural complexity of a family of random variables. This functional is a special case of intricacy, i.e., an average of the mutual information of subsystems whose weights have good mathematical properties. Moreover, its maximum value grows at a definite speed with the size of the system. In this work, we compute exactly this speed of growth by building "approximate maximizers" subject to an entropy condition. These approximate maximizers work simultaneously for all intricacies. We also establish some properties of arbitrary approximate maximizers, in particular the existence of a threshold in the size of subsystems of approximate maximizers: most smaller subsystems are almost equidistributed, most larger subsystems determine the full system. The main ideas are a random construction of almost maximizers with a high statistical symmetry and the consideration of entropy profiles, i.e., the average entropies of sub-systems of a given size. ...

  14. Reinforcement Learning via AIXI Approximation

    OpenAIRE

    Veness, Joel; Ng, Kee Siong; Hutter, Marcus; Silver, David

    2010-01-01

    This paper introduces a principled approach for the design of a scalable general reinforcement learning agent. This approach is based on a direct approximation of AIXI, a Bayesian optimality notion for general reinforcement learning agents. Previously, it has been unclear whether the theory of AIXI could motivate the design of practical algorithms. We answer this hitherto open question in the affirmative, by providing the first computationally feasible approximation to the AIXI agent. To deve...

  15. Binary nucleation beyond capillarity approximation

    OpenAIRE

    Kalikmanov, V.I.

    2010-01-01

    Large discrepancies between binary classical nucleation theory (BCNT) and experiments result from adsorption effects and inability of BCNT, based on the phenomenological capillarity approximation, to treat small clusters. We propose a model aimed at eliminating both of these deficiencies. Adsorption is taken into account within Gibbsian approximation. Binary clusters are treated by means of statistical-mechanical considerations: tracing out the molecular degrees of freedom of the more volatil...

  16. Twisted inhomogeneous Diophantine approximation and badly approximable sets

    CERN Document Server

    Harrap, Stephen

    2010-01-01

    For any real pair i, j geq 0 with i+j=1 let Bad(i, j) denote the set of (i, j)-badly approximable pairs. That is, Bad(i, j) consists of irrational vectors x:=(x_1, x_2) in R^2 for which there exists a positive constant c(x) such that max {||qx_1||^(-i), ||qx_2||^(-j)} > c(x)/q for all q in N. Building on a result of Kurzweil, a new characterization of the set Bad(i, j) in terms of `well-approximable' vectors in the area of `twisted' inhomogeneous Diophantine approximation is established. In addition, it is shown that Bad^x(i, j), the `twisted' inhomogeneous analogue of Bad(i, j), has full Hausdorff dimension 2 when x is chosen from the set Bad(i, j).

  17. Wavelet Sparse Approximate Inverse Preconditioners

    Science.gov (United States)

    Chan, Tony F.; Tang, W.-P.; Wan, W. L.

    1996-01-01

    There is an increasing interest in using sparse approximate inverses as preconditioners for Krylov subspace iterative methods. Recent studies of Grote and Huckle and Chow and Saad also show that sparse approximate inverse preconditioner can be effective for a variety of matrices, e.g. Harwell-Boeing collections. Nonetheless a drawback is that it requires rapid decay of the inverse entries so that sparse approximate inverse is possible. However, for the class of matrices that, come from elliptic PDE problems, this assumption may not necessarily hold. Our main idea is to look for a basis, other than the standard one, such that a sparse representation of the inverse is feasible. A crucial observation is that the kind of matrices we are interested in typically have a piecewise smooth inverse. We exploit this fact, by applying wavelet techniques to construct a better sparse approximate inverse in the wavelet basis. We shall justify theoretically and numerically that our approach is effective for matrices with smooth inverse. We emphasize that in this paper we have only presented the idea of wavelet approximate inverses and demonstrated its potential but have not yet developed a highly refined and efficient algorithm.

  18. The structure of approximate groups

    CERN Document Server

    Breuillard, Emmanuel; Tao, Terence

    2011-01-01

    Let K >= 1 be a parameter. A K-approximate group is a finite set A in a (local) group which contains the identity, is symmetric, and such that A^2 is covered by K left translates of A. The main result of this paper is a qualitative description of approximate groups as being essentially finite-by-nilpotent, answering a conjecture of H. Helfgott and E. Lindenstrauss. This may be viewed as a generalisation of the Freiman-Ruzsa theorem on sets of small doubling in the integers to arbitrary groups. We begin by establishing a correspondence principle between approximate groups and locally compact (local) groups that allows us to recover many results recently established in a fundamental paper of Hrushovski. In particular we establish that approximate groups can be approximately modeled by Lie groups. To prove our main theorem we apply some additional arguments essentially due to Gleason. These arose in the solution of Hilbert's fifth problem in the 1950s. Applications of our main theorem include a finitary refineme...

  19. Relativistic regular approximations revisited: An infinite-order relativistic approximation

    International Nuclear Information System (INIS)

    The concept of the regular approximation is presented as the neglect of the energy dependence of the exact Foldy - Wouthuysen transformation of the Dirac Hamiltonian. Expansion of the normalization terms leads immediately to the zeroth-order regular approximation (ZORA) and first-order regular approximation (FORA) Hamiltonians as the zeroth- and first-order terms of the expansion. The expansion may be taken to infinite order by using an un-normalized Foldy - Wouthuysen transformation, which results in the ZORA Hamiltonian and a non-unit metric. This infinite-order regular approximation, IORA, has eigenvalues which differ from the Dirac eigenvalues by order E3/c4 for a hydrogen-like system, which is a considerable improvement over the ZORA eigenvalues, and similar to the non-variational FORA energies. A further perturbation analysis yields a third-order correction to the IORA energies, TIORA. Results are presented for several systems including the neutral U atom. The IORA eigenvalues for all but the 1s spinor of the neutral system are superior even to the scaled ZORA energies, which are exact for the hydrogenic system. The third-order correction reduces the IORA error for the inner orbitals to a very small fraction of the Dirac eigenvalue. copyright 1999 American Institute of Physics

  20. Approximating Graphic TSP by Matchings

    CERN Document Server

    Mömke, Tobias

    2011-01-01

    We present a framework for approximating the metric TSP based on a novel use of matchings. Traditionally, matchings have been used to add edges in order to make a given graph Eulerian, whereas our approach also allows for the removal of certain edges leading to a decreased cost. For the TSP on graphic metrics (graph-TSP), the approach yields a 1.461-approximation algorithm with respect to the Held-Karp lower bound. For graph-TSP restricted to a class of graphs that contains degree three bounded and claw-free graphs, we show that the integrality gap of the Held-Karp relaxation matches the conjectured ratio 4/3. The framework allows for generalizations in a natural way and also leads to a 1.586-approximation algorithm for the traveling salesman path problem on graphic metrics where the start and end vertices are prespecified.

  1. Approximation methods in probability theory

    CERN Document Server

    Čekanavičius, Vydas

    2016-01-01

    This book presents a wide range of well-known and less common methods used for estimating the accuracy of probabilistic approximations, including the Esseen type inversion formulas, the Stein method as well as the methods of convolutions and triangle function. Emphasising the correct usage of the methods presented, each step required for the proofs is examined in detail. As a result, this textbook provides valuable tools for proving approximation theorems. While Approximation Methods in Probability Theory will appeal to everyone interested in limit theorems of probability theory, the book is particularly aimed at graduate students who have completed a standard intermediate course in probability theory. Furthermore, experienced researchers wanting to enlarge their toolkit will also find this book useful.

  2. Reinforcement Learning via AIXI Approximation

    CERN Document Server

    Veness, Joel; Hutter, Marcus; Silver, David

    2010-01-01

    This paper introduces a principled approach for the design of a scalable general reinforcement learning agent. This approach is based on a direct approximation of AIXI, a Bayesian optimality notion for general reinforcement learning agents. Previously, it has been unclear whether the theory of AIXI could motivate the design of practical algorithms. We answer this hitherto open question in the affirmative, by providing the first computationally feasible approximation to the AIXI agent. To develop our approximation, we introduce a Monte Carlo Tree Search algorithm along with an agent-specific extension of the Context Tree Weighting algorithm. Empirically, we present a set of encouraging results on a number of stochastic, unknown, and partially observable domains.

  3. Concept Approximation between Fuzzy Ontologies

    Institute of Scientific and Technical Information of China (English)

    2006-01-01

    Fuzzy ontologies are efficient tools to handle fuzzy and uncertain knowledge on the semantic web; but there are heterogeneity problems when gaining interoperability among different fuzzy ontologies. This paper uses concept approximation between fuzzy ontologies based on instances to solve the heterogeneity problems. It firstly proposes an instance selection technology based on instance clustering and weighting to unify the fuzzy interpretation of different ontologies and reduce the number of instances to increase the efficiency. Then the paper resolves the problem of computing the approximations of concepts into the problem of computing the least upper approximations of atom concepts. It optimizes the search strategies by extending atom concept sets and defining the least upper bounds of concepts to reduce the searching space of the problem. An efficient algorithm for searching the least upper bounds of concept is given.

  4. An Approximation Ratio for Biclustering

    OpenAIRE

    Puolamäki, Kai; Hanhijärvi, Sami; Garriga, Gemma C

    2007-01-01

    The problem of biclustering consists of the simultaneous clustering of rows and columns of a matrix such that each of the submatrices induced by a pair of row and column clusters is as uniform as possible. In this paper we approximate the optimal biclustering by applying one-way clustering algorithms independently on the rows and on the columns of the input matrix. We show that such a solution yields a worst-case approximation ratio of 1+sqrt(2) under L1-norm for 0-1 valued matrices, and of 2...

  5. An Approximation Ratio for Biclustering

    CERN Document Server

    Puolamäki, Kai; Garriga, Gemma C

    2007-01-01

    The problem of biclustering consists of the simultaneous clustering of rows and columns of a matrix such that each of the submatrices induced by a pair of row and column clusters is as uniform as possible. In this paper we approximate the optimal biclustering by applying one-way clustering algorithms independently on the rows and on the columns of the input matrix. We show that such a solution yields a worst-case approximation ratio of 1+sqrt(2) under L1-norm for 0-1 valued matrices, and of 2 under L2-norm for real valued matrices.

  6. Tree wavelet approximations with applications

    Institute of Scientific and Technical Information of China (English)

    2005-01-01

    [1]Baraniuk, R. G., DeVore, R. A., Kyriazis, G., Yu, X. M., Near best tree approximation, Adv. Comput. Math.,2002, 16: 357-373.[2]Cohen, A., Dahmen, W., Daubechies, I., DeVore, R., Tree approximation and optimal encoding, Appl. Comput.Harmonic Anal., 2001, 11: 192-226.[3]Dahmen, W., Schneider, R., Xu, Y., Nonlinear functionals of wavelet expansions-adaptive reconstruction and fast evaluation, Numer. Math., 2000, 86: 49-101.[4]DeVore, R. A., Nonlinear approximation, Acta Numer., 1998, 7: 51-150.[5]Davis, G., Mallat, S., Avellaneda, M., Adaptive greedy approximations, Const. Approx., 1997, 13: 57-98.[6]DeVore, R. A., Temlyakov, V. N., Some remarks on greedy algorithms, Adv. Comput. Math., 1996, 5: 173-187.[7]Kashin, B. S., Temlyakov, V. N., Best m-term approximations and the entropy of sets in the space L1, Mat.Zametki (in Russian), 1994, 56: 57-86.[8]Temlyakov, V. N., The best m-term approximation and greedy algorithms, Adv. Comput. Math., 1998, 8:249-265.[9]Temlyakov, V. N., Greedy algorithm and m-term trigonometric approximation, Constr. Approx., 1998, 14:569-587.[10]Hutchinson, J. E., Fractals and self similarity, Indiana. Univ. Math. J., 1981, 30: 713-747.[11]Binev, P., Dahmen, W., DeVore, R. A., Petruchev, P., Approximation classes for adaptive methods, Serdica Math.J., 2002, 28: 1001-1026.[12]Gilbarg, D., Trudinger, N. S., Elliptic Partial Differential Equations of Second Order, Berlin: Springer-Verlag,1983.[13]Ciarlet, P. G., The Finite Element Method for Elliptic Problems, New York: North Holland, 1978.[14]Birman, M. S., Solomiak, M. Z., Piecewise polynomial approximation of functions of the class Wαp, Math. Sb.,1967, 73: 295-317.[15]DeVore, R. A., Lorentz, G. G., Constructive Approximation, New York: Springer-Verlag, 1993.[16]DeVore, R. A., Popov, V., Interpolation of Besov spaces, Trans. Amer. Math. Soc., 1988, 305: 397-414.[17]Devore, R., Jawerth, B., Popov, V., Compression of wavelet decompositions, Amer. J. Math., 1992, 114: 737-785.[18]Storozhenko, E

  7. Truthful approximations to range voting

    DEFF Research Database (Denmark)

    Filos-Ratsika, Aris; Miltersen, Peter Bro

    We consider the fundamental mechanism design problem of approximate social welfare maximization under general cardinal preferences on a finite number of alternatives and without money. The well-known range voting scheme can be thought of as a non-truthful mechanism for exact social welfare...

  8. Shearlets and Optimally Sparse Approximations

    DEFF Research Database (Denmark)

    Kutyniok, Gitta; Lemvig, Jakob; Lim, Wang-Q

    2012-01-01

    Multivariate functions are typically governed by anisotropic features such as edges in images or shock fronts in solutions of transport-dominated equations. One major goal both for the purpose of compression as well as for an efficient analysis is the provision of optimally sparse approximations...

  9. Approximate Reasoning with Fuzzy Booleans

    NARCIS (Netherlands)

    Broek, van den P.M.; Noppen, J.A.R.

    2004-01-01

    This paper introduces, in analogy to the concept of fuzzy numbers, the concept of fuzzy booleans, and examines approximate reasoning with the compositional rule of inference using fuzzy booleans. It is shown that each set of fuzzy rules is equivalent to a set of fuzzy rules with singleton crisp ante

  10. Analytical Approximations to Galaxy Clustering

    OpenAIRE

    Mo, H. J.

    1997-01-01

    We discuss some recent progress in constructing analytic approximations to the galaxy clustering. We show that successful models can be constructed for the clustering of both dark matter and dark matter haloes. Our understanding of galaxy clustering and galaxy biasing can be greatly enhanced by these models.

  11. Approximation properties of haplotype tagging

    Directory of Open Access Journals (Sweden)

    Dreiseitl Stephan

    2006-01-01

    Full Text Available Abstract Background Single nucleotide polymorphisms (SNPs are locations at which the genomic sequences of population members differ. Since these differences are known to follow patterns, disease association studies are facilitated by identifying SNPs that allow the unique identification of such patterns. This process, known as haplotype tagging, is formulated as a combinatorial optimization problem and analyzed in terms of complexity and approximation properties. Results It is shown that the tagging problem is NP-hard but approximable within 1 + ln((n2 - n/2 for n haplotypes but not approximable within (1 - ε ln(n/2 for any ε > 0 unless NP ⊂ DTIME(nlog log n. A simple, very easily implementable algorithm that exhibits the above upper bound on solution quality is presented. This algorithm has running time O((2m - p + 1 ≤ O(m(n2 - n/2 where p ≤ min(n, m for n haplotypes of size m. As we show that the approximation bound is asymptotically tight, the algorithm presented is optimal with respect to this asymptotic bound. Conclusion The haplotype tagging problem is hard, but approachable with a fast, practical, and surprisingly simple algorithm that cannot be significantly improved upon on a single processor machine. Hence, significant improvement in computatational efforts expended can only be expected if the computational effort is distributed and done in parallel.

  12. Ultrafast Approximation for Phylogenetic Bootstrap

    NARCIS (Netherlands)

    Bui Quang Minh, [No Value; Nguyen, Thi; von Haeseler, Arndt

    2013-01-01

    Nonparametric bootstrap has been a widely used tool in phylogenetic analysis to assess the clade support of phylogenetic trees. However, with the rapidly growing amount of data, this task remains a computational bottleneck. Recently, approximation methods such as the RAxML rapid bootstrap (RBS) and

  13. Approximation by Penultimate Stable Laws

    NARCIS (Netherlands)

    L.F.M. de Haan (Laurens); L. Peng (Liang); H. Iglesias Pereira

    1997-01-01

    textabstractIn certain cases partial sums of i.i.d. random variables with finite variance are better approximated by a sequence of stable distributions with indices \\\\alpha_n \\\\to 2 than by a normal distribution. We discuss when this happens and how much the convergence rate can be improved by using

  14. Approximate Reanalysis in Topology Optimization

    DEFF Research Database (Denmark)

    Amir, Oded; Bendsøe, Martin P.; Sigmund, Ole

    2009-01-01

    In the nested approach to structural optimization, most of the computational effort is invested in the solution of the finite element analysis equations. In this study, the integration of an approximate reanalysis procedure into the framework of topology optimization of continuum structures...

  15. Low Rank Approximation in $G_0W_0$ Approximation

    CERN Document Server

    Shao, Meiyue; Yang, Chao; Liu, Fang; da Jornada, Felipe H; Deslippe, Jack; Louie, Steven G

    2016-01-01

    The single particle energies obtained in a Kohn--Sham density functional theory (DFT) calculation are generally known to be poor approximations to electron excitation energies that are measured in transport, tunneling and spectroscopic experiments such as photo-emission spectroscopy. The correction to these energies can be obtained from the poles of a single particle Green's function derived from a many-body perturbation theory. From a computational perspective, the accuracy and efficiency of such an approach depends on how a self energy term that properly accounts for dynamic screening of electrons is approximated. The $G_0W_0$ approximation is a widely used technique in which the self energy is expressed as the convolution of a non-interacting Green's function ($G_0$) and a screened Coulomb interaction ($W_0$) in the frequency domain. The computational cost associated with such a convolution is high due to the high complexity of evaluating $W_0$ at multiple frequencies. In this paper, we discuss how the cos...

  16. Hydrogen Beyond the Classic Approximation

    CERN Document Server

    Scivetti, I

    2003-01-01

    The classical nucleus approximation is the most frequently used approach for the resolution of problems in condensed matter physics.However, there are systems in nature where it is necessary to introduce the nuclear degrees of freedom to obtain a correct description of the properties.Examples of this, are the systems with containing hydrogen.In this work, we have studied the resolution of the quantum nuclear problem for the particular case of the water molecule.The Hartree approximation has been used, i.e. we have considered that the nuclei are distinguishable particles.In addition, we have proposed a model to solve the tunneling process, which involves the resolution of the nuclear problem for configurations of the system away from its equilibrium position

  17. Approximate Matching of Hierarchial Data

    DEFF Research Database (Denmark)

    Augsten, Nikolaus

    formally proof that the pq-gram index can be incrementally updated based on the log of edit operations without reconstructing intermediate tree versions. The incremental update is independent of the data size and scales to a large number of changes in the data. We introduce windowed pq-grams for the......-gram based distance between streets, introduces a global greedy matching that guarantees stable pairs, and links addresses that are stored with different granularity. The connector has been successfully tested with public administration databases. Our extensive experiments on both synthetic and real world......The goal of this thesis is to design, develop, and evaluate new methods for the approximate matching of hierarchical data represented as labeled trees. In approximate matching scenarios two items should be matched if they are similar. Computing the similarity between labeled trees is hard as in...

  18. Validity of the eikonal approximation

    CERN Document Server

    Kabat, D

    1992-01-01

    We summarize results on the reliability of the eikonal approximation in obtaining the high energy behavior of a two particle forward scattering amplitude. Reliability depends on the spin of the exchanged field. For scalar fields the eikonal fails at eighth order in perturbation theory, when it misses the leading behavior of the exchange-type diagrams. In a vector theory the eikonal gets the exchange diagrams correctly, but fails by ignoring certain non-exchange graphs which dominate the asymptotic behavior of the full amplitude. For spin--2 tensor fields the eikonal captures the leading behavior of each order in perturbation theory, but the sum of eikonal terms is subdominant to graphs neglected by the approximation. We also comment on the eikonal for Yang-Mills vector exchange, where the additional complexities of the non-abelian theory may be absorbed into Regge-type modifications of the gauge boson propagators.

  19. Approximate Privacy: Foundations and Quantification

    CERN Document Server

    Feigenbaum, Joan; Schapira, Michael

    2009-01-01

    Increasing use of computers and networks in business, government, recreation, and almost all aspects of daily life has led to a proliferation of online sensitive data about individuals and organizations. Consequently, concern about the privacy of these data has become a top priority, particularly those data that are created and used in electronic commerce. There have been many formulations of privacy and, unfortunately, many negative results about the feasibility of maintaining privacy of sensitive data in realistic networked environments. We formulate communication-complexity-based definitions, both worst-case and average-case, of a problem's privacy-approximation ratio. We use our definitions to investigate the extent to which approximate privacy is achievable in two standard problems: the second-price Vickrey auction and the millionaires problem of Yao. For both the second-price Vickrey auction and the millionaires problem, we show that not only is perfect privacy impossible or infeasibly costly to achieve...

  20. Hydrogen: Beyond the Classic Approximation

    International Nuclear Information System (INIS)

    The classical nucleus approximation is the most frequently used approach for the resolution of problems in condensed matter physics.However, there are systems in nature where it is necessary to introduce the nuclear degrees of freedom to obtain a correct description of the properties.Examples of this, are the systems with containing hydrogen.In this work, we have studied the resolution of the quantum nuclear problem for the particular case of the water molecule.The Hartree approximation has been used, i.e. we have considered that the nuclei are distinguishable particles.In addition, we have proposed a model to solve the tunneling process, which involves the resolution of the nuclear problem for configurations of the system away from its equilibrium position

  1. Concentration Bounds for Stochastic Approximations

    CERN Document Server

    Frikha, Noufel

    2012-01-01

    We obtain non asymptotic concentration bounds for two kinds of stochastic approximations. We first consider the deviations between the expectation of a given function of the Euler scheme of some diffusion process at a fixed deterministic time and its empirical mean obtained by the Monte-Carlo procedure. We then give some estimates concerning the deviation between the value at a given time-step of a stochastic approximation algorithm and its target. Under suitable assumptions both concentration bounds turn out to be Gaussian. The key tool consists in exploiting accurately the concentration properties of the increments of the schemes. For the first case, as opposed to the previous work of Lemaire and Menozzi (EJP, 2010), we do not have any systematic bias in our estimates. Also, no specific non-degeneracy conditions are assumed.

  2. Variance approximation under balanced sampling

    OpenAIRE

    Deville, Jean-Claude; Tillé, Yves

    2016-01-01

    A balanced sampling design has the interesting property that Horvitz–Thompson estimators of totals for a set of balancing variables are equal to the totals we want to estimate, therefore the variance of Horvitz–Thompson estimators of variables of interest are reduced in function of their correlations with the balancing variables. Since it is hard to derive an analytic expression for the joint inclusion probabilities, we derive a general approximation of variance based on a residual technique....

  3. Approximate maximizers of intricacy functionals

    OpenAIRE

    Buzzi, Jerome; Zambotti, Lorenzo

    2009-01-01

    G. Edelman, O. Sporns, and G. Tononi introduced in theoretical biology the neural complexity of a family of random variables. This functional is a special case of intricacy, i.e., an average of the mutual information of subsystems whose weights have good mathematical properties. Moreover, its maximum value grows at a definite speed with the size of the system. In this work, we compute exactly this speed of growth by building "approximate maximizers" subject to an entropy condition. These appr...

  4. Approximating Metal-Insulator Transitions

    OpenAIRE

    Danieli, C.; Rayanov, K.; Pavlov, B.; Martin, G.; Flach, S

    2014-01-01

    We consider quantum wave propagation in one-dimensional quasiperiodic lattices. We propose an iterative construction of quasiperiodic potentials from sequences of potentials with increasing spatial period. At each finite iteration step the eigenstates reflect the properties of the limiting quasiperiodic potential properties up to a controlled maximum system size. We then observe approximate metal-insulator transitions (MIT) at the finite iteration steps. We also report evidence on mobility ed...

  5. Product Approximation of Grade and Precision

    Institute of Scientific and Technical Information of China (English)

    ZHANG Xian-yong; MO Zhi-wen

    2005-01-01

    The normal graded approximation and variable precision approximation are defined in approximate space. The relationship between graded approximation and variable precision approximation is studied, and an important formula of conversion between them is achieved. The product approximation of gradeand precision is defined and its basic properties are studied.

  6. Rollout Sampling Approximate Policy Iteration

    CERN Document Server

    Dimitrakakis, Christos

    2008-01-01

    Several researchers have recently investigated the connection between reinforcement learning and classification. We are motivated by proposals of approximate policy iteration schemes without value functions which focus on policy representation using classifiers and address policy learning as a supervised learning problem. This paper proposes variants of an improved policy iteration scheme which addresses the core sampling problem in evaluating a policy through simulation as a multi-armed bandit machine. The resulting algorithm offers comparable performance to the previous algorithm achieved, however, with significantly less computational effort. An order of magnitude improvement is demonstrated experimentally in two standard reinforcement learning domains: inverted pendulum and mountain-car.

  7. Quantum Tunneling Beyond Semiclassical Approximation

    CERN Document Server

    Banerjee, Rabin

    2008-01-01

    Hawking radiation as tunneling by Hamilton-Jacobi method beyond semiclassical approximation is analysed. We compute all quantum corrections in the single particle action revealing that these are proportional to the usual semiclassical contribution. We show that a simple choice of the proportionality constants reproduces the one loop back reaction effect in the spacetime, found by conformal field theory methods, which modifies the Hawking temperature of the black hole. Using the law of black hole mechanics we give the corrections to the Bekenstein-Hawking area law following from the modified Hawking temperature. Some examples are explicitly worked out.

  8. Quantum tunneling beyond semiclassical approximation

    Science.gov (United States)

    Banerjee, Rabin; Ranjan Majhi, Bibhas

    2008-06-01

    Hawking radiation as tunneling by Hamilton-Jacobi method beyond semiclassical approximation is analysed. We compute all quantum corrections in the single particle action revealing that these are proportional to the usual semiclassical contribution. We show that a simple choice of the proportionality constants reproduces the one loop back reaction effect in the spacetime, found by conformal field theory methods, which modifies the Hawking temperature of the black hole. Using the law of black hole mechanics we give the corrections to the Bekenstein-Hawking area law following from the modified Hawking temperature. Some examples are explicitly worked out.

  9. Fermion Tunneling Beyond Semiclassical Approximation

    CERN Document Server

    Majhi, Bibhas Ranjan

    2008-01-01

    Applying the Hamilton-Jacobi method beyond the semiclassical approximation prescribed in \\cite{Majhi3} for the scalar particle, Hawking radiation as tunneling of Dirac particle through an event horizon is analysed. We show that, as before, all quantum corrections in the single particle action are proportional to the usual semiclassical contribution. We also compute the modifications to the Hawking temperature and Bekenstein-Hawking entropy for the Schwarzschild black hole. Finally, the coefficient of the logarithmic correction to entropy is shown to be related with the trace anomaly.

  10. Fermion tunneling beyond semiclassical approximation

    Science.gov (United States)

    Majhi, Bibhas Ranjan

    2009-02-01

    Applying the Hamilton-Jacobi method beyond the semiclassical approximation prescribed in R. Banerjee and B. R. Majhi, J. High Energy Phys.JHEPFG1029-8479 06 (2008) 09510.1088/1126-6708/2008/06/095 for the scalar particle, Hawking radiation as tunneling of the Dirac particle through an event horizon is analyzed. We show that, as before, all quantum corrections in the single particle action are proportional to the usual semiclassical contribution. We also compute the modifications to the Hawking temperature and Bekenstein-Hawking entropy for the Schwarzschild black hole. Finally, the coefficient of the logarithmic correction to entropy is shown to be related with the trace anomaly.

  11. The distorted wave Glauber approximation

    International Nuclear Information System (INIS)

    A solution of the Pauli equation with non-zero potentials defines quantum scalar and vector potentials and magnetic fields and quantum trajectories. If a line integral of perturbing potentials and fields along these quantum trajectories is added to the phase of this solution, an approximate solution of the perturbed equation is found. Glauber theory is a special case and the conditions of applicability are similar. Applications given start from the harmonic oscillator and from a homogeneous magnetic field and add a perturbation. (author)

  12. Approximation of Surfaces by Cylinders

    DEFF Research Database (Denmark)

    Randrup, Thomas

    1998-01-01

    We present a new method for approximation of a given surface by a cylinder surface. It is a constructive geometric method, leading to a monorail representation of the cylinder surface. By use of a weighted Gaussian image of the given surface, we determine a projection plane. In the orthogonal...... projection of the surface onto this plane, a reference curve is determined by use of methods for thinning of binary images. Finally, the cylinder surface is constructed as follows: the directrix of the cylinder surface is determined by a least squares method minimizing the distance to the points...

  13. Wavelet Approximation in Data Assimilation

    Science.gov (United States)

    Tangborn, Andrew; Atlas, Robert (Technical Monitor)

    2002-01-01

    Estimation of the state of the atmosphere with the Kalman filter remains a distant goal because of high computational cost of evolving the error covariance for both linear and nonlinear systems. Wavelet approximation is presented here as a possible solution that efficiently compresses both global and local covariance information. We demonstrate the compression characteristics on the the error correlation field from a global two-dimensional chemical constituent assimilation, and implement an adaptive wavelet approximation scheme on the assimilation of the one-dimensional Burger's equation. In the former problem, we show that 99%, of the error correlation can be represented by just 3% of the wavelet coefficients, with good representation of localized features. In the Burger's equation assimilation, the discrete linearized equations (tangent linear model) and analysis covariance are projected onto a wavelet basis and truncated to just 6%, of the coefficients. A nearly optimal forecast is achieved and we show that errors due to truncation of the dynamics are no greater than the errors due to covariance truncation.

  14. Plasma Physics Approximations in Ares

    Energy Technology Data Exchange (ETDEWEB)

    Managan, R. A.

    2015-01-08

    Lee & More derived analytic forms for the transport properties of a plasma. Many hydro-codes use their formulae for electrical and thermal conductivity. The coefficients are complex functions of Fermi-Dirac integrals, Fn( μ/θ ), the chemical potential, μ or ζ = ln(1+e μ/θ ), and the temperature, θ = kT. Since these formulae are expensive to compute, rational function approximations were fit to them. Approximations are also used to find the chemical potential, either μ or ζ . The fits use ζ as the independent variable instead of μ/θ . New fits are provided for Aα (ζ ),Aβ (ζ ), ζ, f(ζ ) = (1 + e-μ/θ)F1/2(μ/θ), F1/2'/F1/2, Fcα, and Fcβ. In each case the relative error of the fit is minimized since the functions can vary by many orders of magnitude. The new fits are designed to exactly preserve the limiting values in the non-degenerate and highly degenerate limits or as ζ→ 0 or ∞. The original fits due to Lee & More and George Zimmerman are presented for comparison.

  15. Simple approximations for condensational growth

    Energy Technology Data Exchange (ETDEWEB)

    Kostinski, A B [Michigan Technological University, 1400 Townsend Drive, Houghton, MI 49931-1200 (United States)], E-mail: alex.kostinski@mtu.edu

    2009-01-15

    A simple geometric argument relating to the liquid water content of clouds is given. The phase relaxation time and the nature of the quasi-steady approximation for the diffusional growth of cloud drops are elucidated directly in terms of water vapor concentration. Spatial gradients of vapor concentration, inherent in the notion of quasi-steady growth, are discussed and we argue for an occasional reversal of the traditional point of view: rather than a drop growing in response to a given supersaturation, the observed values of the supersaturation in clouds are the result of a vapor field adjusting to droplet growth. Our perspective is illustrated by comparing the exponential decay of condensation trails with a quasi-steady regime of cirrus clouds. The role of aerosol loading in decreasing relaxation times and increasing the rate of growth of the liquid water content is also discussed.

  16. Strong shock implosion, approximate solution

    Science.gov (United States)

    Fujimoto, Y.; Mishkin, E. A.; Alejaldre, C.

    1983-01-01

    The self-similar, center-bound motion of a strong spherical, or cylindrical, shock wave moving through an ideal gas with a constant, γ= cp/ cv, is considered and a linearized, approximate solution is derived. An X, Y phase plane of the self-similar solution is defined and the representative curved of the system behind the shock front is replaced by a straight line connecting the mappings of the shock front with that of its tail. The reduced pressure P(ξ), density R(ξ) and velocity U1(ξ) are found in closed, quite accurate, form. Comparison with numerically obtained results, for γ= {5}/{3} and γ= {7}/{5}, is shown.

  17. Diophantine approximations and Diophantine equations

    CERN Document Server

    Schmidt, Wolfgang M

    1991-01-01

    "This book by a leading researcher and masterly expositor of the subject studies diophantine approximations to algebraic numbers and their applications to diophantine equations. The methods are classical, and the results stressed can be obtained without much background in algebraic geometry. In particular, Thue equations, norm form equations and S-unit equations, with emphasis on recent explicit bounds on the number of solutions, are included. The book will be useful for graduate students and researchers." (L'Enseignement Mathematique) "The rich Bibliography includes more than hundred references. The book is easy to read, it may be a useful piece of reading not only for experts but for students as well." Acta Scientiarum Mathematicarum

  18. Stochastic Approximation with Averaging Innovation

    CERN Document Server

    Laruelle, Sophie

    2010-01-01

    The aim of the paper is to establish a convergence theorem for multi-dimensional stochastic approximation in a setting with innovations satisfying some averaging properties and to study some applications. The averaging assumptions allow us to unify the framework where the innovations are generated (to solve problems from Numerical Probability) and the one with exogenous innovations (market data, output of "device" $e.g.$ an Euler scheme) with stationary or ergodic properties. We propose several fields of applications with random innovations or quasi-random numbers. In particular we provide in both setting a rule to tune the step of the algorithm. At last we illustrate our results on five examples notably in Finance.

  19. Benchmarking Declarative Approximate Selection Predicates

    CERN Document Server

    Hassanzadeh, Oktie

    2009-01-01

    Declarative data quality has been an active research topic. The fundamental principle behind a declarative approach to data quality is the use of declarative statements to realize data quality primitives on top of any relational data source. A primary advantage of such an approach is the ease of use and integration with existing applications. Several similarity predicates have been proposed in the past for common quality primitives (approximate selections, joins, etc.) and have been fully expressed using declarative SQL statements. In this thesis, new similarity predicates are proposed along with their declarative realization, based on notions of probabilistic information retrieval. Then, full declarative specifications of previously proposed similarity predicates in the literature are presented, grouped into classes according to their primary characteristics. Finally, a thorough performance and accuracy study comparing a large number of similarity predicates for data cleaning operations is performed.

  20. Approximating metal-insulator transitions

    Science.gov (United States)

    Danieli, Carlo; Rayanov, Kristian; Pavlov, Boris; Martin, Gaven; Flach, Sergej

    2015-12-01

    We consider quantum wave propagation in one-dimensional quasiperiodic lattices. We propose an iterative construction of quasiperiodic potentials from sequences of potentials with increasing spatial period. At each finite iteration step, the eigenstates reflect the properties of the limiting quasiperiodic potential properties up to a controlled maximum system size. We then observe approximate Metal-Insulator Transitions (MIT) at the finite iteration steps. We also report evidence on mobility edges, which are at variance to the celebrated Aubry-André model. The dynamics near the MIT shows a critical slowing down of the ballistic group velocity in the metallic phase, similar to the divergence of the localization length in the insulating phase.

  1. Multidimensional stochastic approximation Monte Carlo.

    Science.gov (United States)

    Zablotskiy, Sergey V; Ivanov, Victor A; Paul, Wolfgang

    2016-06-01

    Stochastic Approximation Monte Carlo (SAMC) has been established as a mathematically founded powerful flat-histogram Monte Carlo method, used to determine the density of states, g(E), of a model system. We show here how it can be generalized for the determination of multidimensional probability distributions (or equivalently densities of states) of macroscopic or mesoscopic variables defined on the space of microstates of a statistical mechanical system. This establishes this method as a systematic way for coarse graining a model system, or, in other words, for performing a renormalization group step on a model. We discuss the formulation of the Kadanoff block spin transformation and the coarse-graining procedure for polymer models in this language. We also apply it to a standard case in the literature of two-dimensional densities of states, where two competing energetic effects are present g(E_{1},E_{2}). We show when and why care has to be exercised when obtaining the microcanonical density of states g(E_{1}+E_{2}) from g(E_{1},E_{2}). PMID:27415383

  2. Randomized approximate nearest neighbors algorithm.

    Science.gov (United States)

    Jones, Peter Wilcox; Osipov, Andrei; Rokhlin, Vladimir

    2011-09-20

    We present a randomized algorithm for the approximate nearest neighbor problem in d-dimensional Euclidean space. Given N points {x(j)} in R(d), the algorithm attempts to find k nearest neighbors for each of x(j), where k is a user-specified integer parameter. The algorithm is iterative, and its running time requirements are proportional to T·N·(d·(log d) + k·(d + log k)·(log N)) + N·k(2)·(d + log k), with T the number of iterations performed. The memory requirements of the procedure are of the order N·(d + k). A by-product of the scheme is a data structure, permitting a rapid search for the k nearest neighbors among {x(j)} for an arbitrary point x ∈ R(d). The cost of each such query is proportional to T·(d·(log d) + log(N/k)·k·(d + log k)), and the memory requirements for the requisite data structure are of the order N·(d + k) + T·(d + N). The algorithm utilizes random rotations and a basic divide-and-conquer scheme, followed by a local graph search. We analyze the scheme's behavior for certain types of distributions of {x(j)} and illustrate its performance via several numerical examples.

  3. Fuzzy Set Approximations in Fuzzy Formal Contexts

    Institute of Scientific and Technical Information of China (English)

    Mingwen Shao; Shiqing Fan

    2006-01-01

    In this paper, a kind of multi-level formal concept is introduced. Based on the proposed multi-level formal concept, we present a pair of rough fuzzy set approximations within fuzzy formal contexts. By the proposed rough fuzzy set approximations, we can approximate a fuzzy set according to different precision level. We discuss the properties of the proposed approximation operators in detail.

  4. Obtaining exact value by approximate computations

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    Numerical approximate computations can solve large and complex problems fast.They have the advantage of high efficiency.However they only give approximate results,whereas we need exact results in some fields.There is a gap between approximate computations and exact results. In this paper,we build a bridge by which exact results can be obtained by numerical approximate computations.

  5. APPROXIMATE SAMPLING THEOREM FOR BIVARIATE CONTINUOUS FUNCTION

    Institute of Scientific and Technical Information of China (English)

    杨守志; 程正兴; 唐远炎

    2003-01-01

    An approximate solution of the refinement equation was given by its mask, and the approximate sampling theorem for bivariate continuous function was proved by applying the approximate solution. The approximate sampling function defined uniquely by the mask of the refinement equation is the approximate solution of the equation, a piece-wise linear function, and posseses an explicit computation formula. Therefore the mask of the refinement equation is selected according to one' s requirement, so that one may controll the decay speed of the approximate sampling function.

  6. Reduction of Linear Programming to Linear Approximation

    OpenAIRE

    Vaserstein, Leonid N.

    2006-01-01

    It is well known that every Chebyshev linear approximation problem can be reduced to a linear program. In this paper we show that conversely every linear program can be reduced to a Chebyshev linear approximation problem.

  7. Legendre rational approximation on the whole line

    Institute of Scientific and Technical Information of China (English)

    GUO; Benyu; WANG; Zhongqing

    2004-01-01

    The Legendre rational approximation is investigated. Some approximation results are established, which form the mathematical foundation of a new spectral method on the whole line. A model problem is considered. Numerical results show the efficiency of this new approach.

  8. An approximation technique for jet impingement flow

    Energy Technology Data Exchange (ETDEWEB)

    Najafi, Mahmoud; Fincher, Donald [Kent State University Ashtabula Department of Mathematical Sciences (United States); Rahni, Taeibi; Javadi, KH. [Department of Aerospace Engineering, Sharif University of Technology (Iran, Islamic Republic of); Massah, H. [Acoustic Research Center, Institute of Applied Physics (Iran, Islamic Republic of)

    2015-03-10

    The analytical approximate solution of a non-linear jet impingement flow model will be demonstrated. We will show that this is an improvement over the series approximation obtained via the Adomian decomposition method, which is itself, a powerful method for analysing non-linear differential equations. The results of these approximations will be compared to the Runge-Kutta approximation in order to demonstrate their validity.

  9. Diophantine approximation and special Liouville numbers

    OpenAIRE

    Schleischitz, Johannes

    2013-01-01

    This paper introduces some methods to determine the simultaneous approximation constants of a class of well approximable numbers $\\zeta_{1},\\zeta_{2},...,\\zeta_{k}$. The approach relies on results on the connection between the set of all $s$-adic expansions ($s\\geq 2$) of $\\zeta_{1},\\zeta_{2},...,\\zeta_{k}$ and their associated approximation constants. As an application, explicit construction of real numbers $\\zeta_{1},\\zeta_{2},...,\\zeta_{k}$ with prescribed approximation properties are dedu...

  10. Nonlinear approximation with bi-framelets

    DEFF Research Database (Denmark)

    Borup, Lasse; Nielsen, Morten; Gribonval, Rémi

    2005-01-01

    We study the approximation in Lebesgue spaces of wavelet bi-frame systems given by translations and dilations of a finite set of generators. A complete characterization of the approximation spaces associated with best m-term approximation of wavelet bi-framelet systems is given. The characterizat...

  11. Some relations between entropy and approximation numbers

    Institute of Scientific and Technical Information of China (English)

    郑志明

    1999-01-01

    A general result is obtained which relates the entropy numbers of compact maps on Hilbert space to its approximation numbers. Compared with previous works in this area, it is particularly convenient for dealing with the cases where the approximation numbers decay rapidly. A nice estimation between entropy and approximation numbers for noncompact maps is given.

  12. NONLINEAR APPROXIMATION WITH GENERAL WAVE PACKETS

    Institute of Scientific and Technical Information of China (English)

    L. Borup; M. Nielsen

    2005-01-01

    We study nonlinear approximation in the Triebel-Lizorkin spaces with dictionaries formed by dilating and translating one single function g. A general Jackson inequality is derived for best m-term approximation with such dictionaries. In some special cases where g has a special structure, a complete characterization of the approximation spaces is derived.

  13. Axiomatic Characterizations of IVF Rough Approximation Operators

    Directory of Open Access Journals (Sweden)

    Guangji Yu

    2014-01-01

    Full Text Available This paper is devoted to the study of axiomatic characterizations of IVF rough approximation operators. IVF approximation spaces are investigated. The fact that different IVF operators satisfy some axioms to guarantee the existence of different types of IVF relations which produce the same operators is proved and then IVF rough approximation operators are characterized by axioms.

  14. Approximation properties of fine hyperbolic graphs

    Indian Academy of Sciences (India)

    Benyin Fu

    2016-05-01

    In this paper, we propose a definition of approximation property which is called the metric invariant translation approximation property for a countable discrete metric space. Moreover, we use the techniques of Ozawa’s to prove that a fine hyperbolic graph has the metric invariant translation approximation property.

  15. Nonlinear approximation with dictionaries, I: Direct estimates

    DEFF Research Database (Denmark)

    Gribonval, Rémi; Nielsen, Morten

    $-term approximation with algorithmic constraints: thresholding and Chebychev approximation classes are studied respectively. We consider embeddings of the Jackson type (direct estimates) of sparsity spaces into the mentioned approximation classes. General direct estimates are based on the geometry of the Banach space...

  16. Approximate Nearest Neighbor Queries among Parallel Segments

    DEFF Research Database (Denmark)

    Emiris, Ioannis Z.; Malamatos, Theocharis; Tsigaridas, Elias

    2010-01-01

    We develop a data structure for answering efficiently approximate nearest neighbor queries over a set of parallel segments in three dimensions. We connect this problem to approximate nearest neighbor searching under weight constraints and approximate nearest neighbor searching on historical data...

  17. Truth Approximation, Social Epistemology, and Opinion Dynamics

    NARCIS (Netherlands)

    Douven, Igor; Kelp, Christoph

    2011-01-01

    This paper highlights some connections between work on truth approximation and work in social epistemology, in particular work on peer disagreement. In some of the literature on truth approximation, questions have been addressed concerning the efficiency of research strategies for approximating the

  18. Upper Bounds on Numerical Approximation Errors

    DEFF Research Database (Denmark)

    Raahauge, Peter

    2004-01-01

    This paper suggests a method for determining rigorous upper bounds on approximationerrors of numerical solutions to infinite horizon dynamic programming models.Bounds are provided for approximations of the value function and the policyfunction as well as the derivatives of the value function...... to approximations of a standard (strictly concave)growth model.KEYWORDS: Numerical approximation errors, Bellman contractions, Error bounds...

  19. APPROXIMATE AMENABILITY OF CERTAIN INVERSE SEMIGROUP ALGEBRAS

    Institute of Scientific and Technical Information of China (English)

    Mehdi ROSTAMI; Abdolrasoul POURABBAS; Morteza ESSMAILI

    2013-01-01

    In this article,the approximate amenability of semigroup algebra e1(S) is investigated,where S is a uniformly locally finite inverse semigroup.Indeed,we show that for a uniformly locally finite inverse semigroup S,the notions of amenability,approximate amenability and bounded approximate amenability of e1 (S) are equivalent.We use this to give a direct proof of the approximate amenability of e1(S) for a Brandt semigroup S.Moreover,we characterize the approximate amenability of e1(S),where S is a uniformly locally finite band semigroup.

  20. Mapping moveout approximations in TI media

    KAUST Repository

    Stovas, Alexey

    2013-11-21

    Moveout approximations play a very important role in seismic modeling, inversion, and scanning for parameters in complex media. We developed a scheme to map one-way moveout approximations for transversely isotropic media with a vertical axis of symmetry (VTI), which is widely available, to the tilted case (TTI) by introducing the effective tilt angle. As a result, we obtained highly accurate TTI moveout equations analogous with their VTI counterparts. Our analysis showed that the most accurate approximation is obtained from the mapping of generalized approximation. The new moveout approximations allow for, as the examples demonstrate, accurate description of moveout in the TTI case even for vertical heterogeneity. The proposed moveout approximations can be easily used for inversion in a layered TTI medium because the parameters of these approximations explicitly depend on corresponding effective parameters in a layered VTI medium.

  1. On Gakerkin approximations for the quasigeostrophic equations

    CERN Document Server

    Rocha, Cesar B; Grooms, Ian

    2015-01-01

    We study the representation of approximate solutions of the three-dimensional quasigeostrophic (QG) equations using Galerkin series with standard vertical modes. In particular, we show that standard modes are compatible with nonzero buoyancy at the surfaces and can be used to solve the Eady problem. We extend two existing Galerkin approaches (A and B) and develop a new Galerkin approximation (C). Approximation A, due to Flierl (1978), represents the streamfunction as a truncated Galerkin series and defines the potential vorticity (PV) that satisfies the inversion problem exactly. Approximation B, due to Tulloch and Smith (2009b), represents the PV as a truncated Galerkin series and calculates the streamfunction that satisfies the inversion problem exactly. Approximation C, the true Galerkin approximation for the QG equations, represents both streamfunction and PV as truncated Galerkin series, but does not satisfy the inversion equation exactly. The three approximations are fundamentally different unless the b...

  2. Frankenstein's Glue: Transition functions for approximate solutions

    CERN Document Server

    Yunes, N

    2006-01-01

    Approximations are commonly employed to find approximate solutions to the Einstein equations. These solutions, however, are usually only valid in some specific spacetime region. A global solution can be constructed by gluing approximate solutions together, but this procedure is difficult because discontinuities can arise, leading to large violations of the Einstein equations. In this paper, we provide an attempt to formalize this gluing scheme by studying transition functions that join approximate solutions together. In particular, we propose certain sufficient conditions on these functions and proof that these conditions guarantee that the joined solution still satisfies the Einstein equations to the same order as the approximate ones. An example is also provided for a binary system of non-spinning black holes, where the approximate solutions are taken to be given by a post-Newtonian expansion and a perturbed Schwarzschild solution. For this specific case, we show that if the transition functions satisfy the...

  3. A Conditional Saddlepoint Approximation for Testing Problems

    OpenAIRE

    Gatto, R.; Jammalamadaka, SR

    1999-01-01

    A saddlepoint approximation is provided for the distribution function of one M statistic conditional on another M statistic. Many interesting statistics based on dependent quantities (e.g., spacings, multinomial frequencies, rank differences) can be expressed in terms of independent identically distributed random variables conditioned on their sum, so that this conditional saddlepoint approximation yields accurate approximations for the distribution of such statistics. This saddlepoint approx...

  4. Approximation Resistant Predicates From Pairwise Independence

    CERN Document Server

    Austrin, Per

    2008-01-01

    We study the approximability of predicates on $k$ variables from a domain $[q]$, and give a new sufficient condition for such predicates to be approximation resistant under the Unique Games Conjecture. Specifically, we show that a predicate $P$ is approximation resistant if there exists a balanced pairwise independent distribution over $[q]^k$ whose support is contained in the set of satisfying assignments to $P$.

  5. Weak approximation of second-order BSDEs

    OpenAIRE

    Possamaï, Dylan; Tan, Xiaolu

    2013-01-01

    We study the weak approximation of the second-order backward SDEs (2BSDEs), when the continuous driving martingales are approximated by discrete time martingales. We establish a convergence result for a class of 2BSDEs, using both robustness properties of BSDEs, as proved in Briand, Delyon and M\\'{e}min [Stochastic Process. Appl. 97 (2002) 229-253], and tightness of solutions to discrete time BSDEs. In particular, when the approximating martingales are given by some particular controlled Mark...

  6. A Note on Generalized Approximation Property

    Directory of Open Access Journals (Sweden)

    Antara Bhar

    2013-01-01

    Full Text Available We introduce a notion of generalized approximation property, which we refer to as --AP possessed by a Banach space , corresponding to an arbitrary Banach sequence space and a convex subset of , the class of bounded linear operators on . This property includes approximation property studied by Grothendieck, -approximation property considered by Sinha and Karn and Delgado et al., and also approximation property studied by Lissitsin et al. We characterize a Banach space having --AP with the help of -compact operators, -nuclear operators, and quasi--nuclear operators. A particular case for ( has also been characterized.

  7. INVARIANT RANDOM APPROXIMATION IN NONCONVEX DOMAIN

    Directory of Open Access Journals (Sweden)

    R. Shrivastava

    2012-05-01

    Full Text Available Random fixed point results in the setup of compact and weakly compact domain of Banach spaces which is not necessary starshaped have been obtained in the present work. Invariant random approximation results have also been determined asits application. In this way, random version of invariant approximation results due toMukherjee and Som [13] and Singh [17] have been given.

  8. On approximating multi-criteria TSP

    NARCIS (Netherlands)

    Manthey, Bodo; Albers, S.; Marion, J.-Y.

    2009-01-01

    We present approximation algorithms for almost all variants of the multi-criteria traveling salesman problem (TSP), whose performances are independent of the number $k$ of criteria and come close to the approximation ratios obtained for TSP with a single objective function. We present randomized app

  9. A case where BO Approximation breaks down

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    @@ The Bom-Oppenheimer (BO)Approximation is ubiquitous in molecular physics,quantum physics and quantum chemistry. However, CAS researchers recently observed a breakdown of the Approximation in the reaction of fluorine with deuterium atoms.The result has been published in the August 24 issue of Science.

  10. Computing Functions by Approximating the Input

    Science.gov (United States)

    Goldberg, Mayer

    2012-01-01

    In computing real-valued functions, it is ordinarily assumed that the input to the function is known, and it is the output that we need to approximate. In this work, we take the opposite approach: we show how to compute the values of some transcendental functions by approximating the input to these functions, and obtaining exact answers for their…

  11. Approximation of the Inverse -Frame Operator

    Indian Academy of Sciences (India)

    M R Abdollahpour; A Najati

    2011-05-01

    In this paper, we introduce the concept of (strong) projection method for -frames which works for all conditional -Riesz frames. We also derive a method for approximation of the inverse -frame operator which is efficient for all -frames. We show how the inverse of -frame operator can be approximated as close as we like using finite-dimensional linear algebra.

  12. Two Point Pade Approximants and Duality

    CERN Document Server

    Banks, Tom

    2013-01-01

    We propose the use of two point Pade approximants to find expressions valid uniformly in coupling constant for theories with both weak and strong coupling expansions. In particular, one can use these approximants in models with a strong/weak duality, when the symmetries do not determine exact expressions for some quantity.

  13. Nonlinear approximation with dictionaries I. Direct estimates

    DEFF Research Database (Denmark)

    Gribonval, Rémi; Nielsen, Morten

    2004-01-01

    with algorithmic constraints: thresholding and Chebychev approximation classes are studied, respectively. We consider embeddings of the Jackson type (direct estimates) of sparsity spaces into the mentioned approximation classes. General direct estimates are based on the geometry of the Banach space, and we prove...

  14. Nonlinear approximation with general wave packets

    DEFF Research Database (Denmark)

    Borup, Lasse; Nielsen, Morten

    2005-01-01

    We study nonlinear approximation in the Triebel-Lizorkin spaces with dictionaries formed by dilating and translating one single function g. A general Jackson inequality is derived for best m-term approximation with such dictionaries. In some special cases where g has a special structure, a comple...

  15. Improved Approximation for the Directed Spanner Problem

    OpenAIRE

    Bhattacharyya, Arnab; Makarychev, Konstantin

    2010-01-01

    We prove that the size of the sparsest directed k-spanner of a graph can be approximated in polynomial time to within a factor of $\\tilde{O}(\\sqrt{n})$, for all k >= 3. This improves the $\\tilde{O}(n^{2/3})$-approximation recently shown by Dinitz and Krauthgamer.

  16. On the closedness of approximation spectra

    OpenAIRE

    Parkkonen, Jouni; Paulin, Frédéric

    2008-01-01

    Generalizing Cusick's theorem on the closedness of the classical Lagrange spectrum for the approximation of real numbers by rational ones, we prove that various approximation spectra are closed, using penetration properties of the geodesic flow in cusp neighbourhoods in negatively curved manifolds and a result of Maucourant.

  17. Inverse scattering problem in relativistic quasiclassical approximation

    International Nuclear Information System (INIS)

    Inverse scattering problem is solved on the basis of quasipotential approach in quantum field theory within the framework of relativistic quasiclassical approximation. Formulas of quasipotential restoration by phase shifts are derived. Cases of non-relativistic and ultra-relativistic approximations are investigated

  18. Approximations for stop-loss reinsurance premiums

    NARCIS (Netherlands)

    Reijnen, Rajko; Albers, Willem; Kallenberg, Wilbert C.M.

    2005-01-01

    Various approximations of stop-loss reinsurance premiums are described in literature. For a wide variety of claim size distributions and retention levels, such approximations are compared in this paper to each other, as well as to a quantitative criterion. For the aggregate claims two models are use

  19. Approximations for stop-loss reinsurance premiums

    NARCIS (Netherlands)

    Reijnen, R.; Albers, W.; Kallenberg, W.C.M.

    2003-01-01

    Various approximations of stop-loss reinsurance premiums are described in literature. For a wide variety of claim size distributions and retention levels, such approximations are compared in this paper to each other, as well as to a quantitative criterion. For the aggregate claims two models are use

  20. On approximating multi-criteria TSP

    NARCIS (Netherlands)

    Manthey, Bodo

    2012-01-01

    We present approximation algorithms for almost all variants of the multicriteria traveling salesman problem (TSP). First, we devise randomized approximation algorithms for multicriteria maximum traveling salesman problems (Max-TSP). For multicriteria Max-STSP where the edge weights have to be symmet

  1. Metadata for Approximate Query Answering Systems

    Directory of Open Access Journals (Sweden)

    Francesco Di Tria

    2012-01-01

    Full Text Available In business intelligence systems, data warehouse metadata management and representation are getting more and more attention by vendors and designers. The standard language for the data warehouse metadata representation is the Common Warehouse Metamodel. However, business intelligence systems include also approximate query answering systems, since these software tools provide fast responses for decision making on the basis of approximate query processing. Currently, the standard meta-model does not allow to represent the metadata needed by approximate query answering systems. In this paper, we propose an extension of the standard metamodel, in order to define the metadata to be used in online approximate analytical processing. These metadata have been successfully adopted in ADAP, a web-based approximate query answering system that creates and uses statistical data profiles.

  2. Trajectory averaging for stochastic approximation MCMC algorithms

    CERN Document Server

    Liang, Faming

    2010-01-01

    The subject of stochastic approximation was founded by Robbins and Monro [Ann. Math. Statist. 22 (1951) 400--407]. After five decades of continual development, it has developed into an important area in systems control and optimization, and it has also served as a prototype for the development of adaptive algorithms for on-line estimation and control of stochastic systems. Recently, it has been used in statistics with Markov chain Monte Carlo for solving maximum likelihood estimation problems and for general simulation and optimizations. In this paper, we first show that the trajectory averaging estimator is asymptotically efficient for the stochastic approximation MCMC (SAMCMC) algorithm under mild conditions, and then apply this result to the stochastic approximation Monte Carlo algorithm [Liang, Liu and Carroll J. Amer. Statist. Assoc. 102 (2007) 305--320]. The application of the trajectory averaging estimator to other stochastic approximation MCMC algorithms, for example, a stochastic approximation MLE al...

  3. An improved proximity force approximation for electrostatics

    CERN Document Server

    Fosco, C D; Mazzitelli, F D

    2012-01-01

    A quite straightforward approximation for the electrostatic interaction between two perfectly conducting surfaces suggests itself when the distance between them is much smaller than the characteristic lengths associated to their shapes. Indeed, in the so called "proximity force approximation" the electrostatic force is evaluated by first dividing each surface into a set of small flat patches, and then adding up the forces due two opposite pairs, the contribution of which are approximated as due to pairs of parallel planes. This approximation has been widely and successfully applied to different contexts, ranging from nuclear physics to Casimir effect calculations. We present here an improvement on this approximation, based on a derivative expansion for the electrostatic energy contained between the surfaces. The results obtained could be useful to discuss the geometric dependence of the electrostatic force, and also as a convenient benchmark for numerical analyses of the tip-sample electrostatic interaction i...

  4. Frankenstein's glue: transition functions for approximate solutions

    Science.gov (United States)

    Yunes, Nicolás

    2007-09-01

    Approximations are commonly employed to find approximate solutions to the Einstein equations. These solutions, however, are usually only valid in some specific spacetime region. A global solution can be constructed by gluing approximate solutions together, but this procedure is difficult because discontinuities can arise, leading to large violations of the Einstein equations. In this paper, we provide an attempt to formalize this gluing scheme by studying transition functions that join approximate analytic solutions together. In particular, we propose certain sufficient conditions on these functions and prove that these conditions guarantee that the joined solution still satisfies the Einstein equations analytically to the same order as the approximate ones. An example is also provided for a binary system of non-spinning black holes, where the approximate solutions are taken to be given by a post-Newtonian expansion and a perturbed Schwarzschild solution. For this specific case, we show that if the transition functions satisfy the proposed conditions, then the joined solution does not contain any violations to the Einstein equations larger than those already inherent in the approximations. We further show that if these functions violate the proposed conditions, then the matter content of the spacetime is modified by the introduction of a matter shell, whose stress energy tensor depends on derivatives of these functions.

  5. Orthorhombic rational approximants for decagonal quasicrystals

    Indian Academy of Sciences (India)

    S Ranganathan; Anandh Subramaniam

    2003-10-01

    An important exercise in the study of rational approximants is to derive their metric, especially in relation to the corresponding quasicrystal or the underlying clusters. Kuo’s model has been the widely accepted model to calculate the metric of the decagonal approximants. Using an alternate model, the metric of the approximants and other complex structures with the icosahedral cluster are explained elsewhere. In this work a comparison is made between the two models bringing out their equivalence. Further, using the concept of average lattices, a modified model is proposed.

  6. DIFFERENCE SCHEMES BASING ON COEFFICIENT APPROXIMATION

    Institute of Scientific and Technical Information of China (English)

    MOU Zong-ze; LONG Yong-xing; QU Wen-xiao

    2005-01-01

    In respect of variable coefficient differential equations, the equations of coefficient function approximation were more accurate than the coefficient to be frozen as a constant in every discrete subinterval. Usually, the difference schemes constructed based on Taylor expansion approximation of the solution do not suit the solution with sharp function.Introducing into local bases to be combined with coefficient function approximation, the difference can well depict more complex physical phenomena, for example, boundary layer as well as high oscillatory, with sharp behavior. The numerical test shows the method is more effective than the traditional one.

  7. Approximately Liner Phase IIR Digital Filter Banks

    Directory of Open Access Journals (Sweden)

    J. D. Ćertić

    2013-11-01

    Full Text Available In this paper, uniform and nonuniform digital filter banks based on approximately linear phase IIR filters and frequency response masking technique (FRM are presented. Both filter banks are realized as a connection of an interpolated half-band approximately linear phase IIR filter as a first stage of the FRM design and an appropriate number of masking filters. The masking filters are half-band IIR filters with an approximately linear phase. The resulting IIR filter banks are compared with linear-phase FIR filter banks exhibiting similar magnitude responses. The effects of coefficient quantization are analyzed.

  8. The Boussinesq approximation in rapidly rotating flows

    CERN Document Server

    Lopez, Jose M; Avila, Marc

    2013-01-01

    In the classical formulation of the Boussinesq approximation centrifugal buoyancy effects related to differential rotation, as well as strong vortices in the flow, are neglected. However, these may play an important role in rapidly rotating flows, such as in astrophysical and geophysical applications, and also in turbulent convection. We here provide a straightforward approach resulting in a Boussinesq-type approximation that consistently accounts for centrifugal effects. We further compare our new approach to the classical one in fluid flows confined between two differentially heated and rotating cylinders. The results justify the need of using the proposed approximation in rapidly rotating flows.

  9. Mathematical analysis, approximation theory and their applications

    CERN Document Server

    Gupta, Vijay

    2016-01-01

    Designed for graduate students, researchers, and engineers in mathematics, optimization, and economics, this self-contained volume presents theory, methods, and applications in mathematical analysis and approximation theory. Specific topics include: approximation of functions by linear positive operators with applications to computer aided geometric design, numerical analysis, optimization theory, and solutions of differential equations. Recent and significant developments in approximation theory, special functions and q-calculus along with their applications to mathematics, engineering, and social sciences are discussed and analyzed. Each chapter enriches the understanding of current research problems and theories in pure and applied research.

  10. Approximate equivalence in von Neumann algebras

    Institute of Scientific and Technical Information of China (English)

    DING Huiru; Don Hadwin

    2005-01-01

    One formulation of D. Voiculescu's theorem on approximate unitary equivalence is that two unital representations π and ρ of a separable C*-algebra are approximately unitarily equivalent if and only if rank o π = rank o ρ. We study the analog when the ranges of π and ρ are contained in a von Neumann algebra R, the unitaries inducing the approximate equivalence must come from R, and "rank" is replaced with "R-rank" (defined as the Murray-von Neumann equivalence of the range projection).

  11. Regression with Sparse Approximations of Data

    DEFF Research Database (Denmark)

    Noorzad, Pardis; Sturm, Bob L.

    2012-01-01

    We propose sparse approximation weighted regression (SPARROW), a method for local estimation of the regression function that uses sparse approximation with a dictionary of measurements. SPARROW estimates the regression function at a point with a linear combination of a few regressands selected by a...... sparse approximation of the point in terms of the regressors. We show SPARROW can be considered a variant of \\(k\\)-nearest neighbors regression (\\(k\\)-NNR), and more generally, local polynomial kernel regression. Unlike \\(k\\)-NNR, however, SPARROW can adapt the number of regressors to use based on the...

  12. Relativistic stellar pulsations in the Cowling approximation

    International Nuclear Information System (INIS)

    Much that is known about the general pulsational properties of non-rotating Newtonian stars is traceable to the fact that in the Cowling approximation, the stellar pulsation equations can be cast in a nearly Sturm-Liouville form. In this paper, the relativistic Cowling approximation is investigated, and it is shown that in this approximation the equations for non-radial relativistic stellar pulsations are also of nearly Sturm-Liouville character. The consequences of this are discussed as a series of theorems regarding the eigenfrequencies and eigenfunctions of g-, f- and p-modes in relativistic stars. (author)

  13. Approximation of the semi-infinite interval

    Directory of Open Access Journals (Sweden)

    A. McD. Mercer

    1980-01-01

    Full Text Available The approximation of a function f∈C[a,b] by Bernstein polynomials is well-known. It is based on the binomial distribution. O. Szasz has shown that there are analogous approximations on the interval [0,∞ based on the Poisson distribution. Recently R. Mohapatra has generalized Szasz' result to the case in which the approximating function is αe−ux∑k=N∞(uxkα+β−1Γ(kα+βf(kαuThe present note shows that these results are special cases of a Tauberian theorem for certain infinite series having positive coefficients.

  14. Low Rank Approximation Algorithms, Implementation, Applications

    CERN Document Server

    Markovsky, Ivan

    2012-01-01

    Matrix low-rank approximation is intimately related to data modelling; a problem that arises frequently in many different fields. Low Rank Approximation: Algorithms, Implementation, Applications is a comprehensive exposition of the theory, algorithms, and applications of structured low-rank approximation. Local optimization methods and effective suboptimal convex relaxations for Toeplitz, Hankel, and Sylvester structured problems are presented. A major part of the text is devoted to application of the theory. Applications described include: system and control theory: approximate realization, model reduction, output error, and errors-in-variables identification; signal processing: harmonic retrieval, sum-of-damped exponentials, finite impulse response modeling, and array processing; machine learning: multidimensional scaling and recommender system; computer vision: algebraic curve fitting and fundamental matrix estimation; bioinformatics for microarray data analysis; chemometrics for multivariate calibration; ...

  15. Methods of Fourier analysis and approximation theory

    CERN Document Server

    Tikhonov, Sergey

    2016-01-01

    Different facets of interplay between harmonic analysis and approximation theory are covered in this volume. The topics included are Fourier analysis, function spaces, optimization theory, partial differential equations, and their links to modern developments in the approximation theory. The articles of this collection were originated from two events. The first event took place during the 9th ISAAC Congress in Krakow, Poland, 5th-9th August 2013, at the section “Approximation Theory and Fourier Analysis”. The second event was the conference on Fourier Analysis and Approximation Theory in the Centre de Recerca Matemàtica (CRM), Barcelona, during 4th-8th November 2013, organized by the editors of this volume. All articles selected to be part of this collection were carefully reviewed.

  16. Seismic wave extrapolation using lowrank symbol approximation

    KAUST Repository

    Fomel, Sergey

    2012-04-30

    We consider the problem of constructing a wave extrapolation operator in a variable and possibly anisotropic medium. Our construction involves Fourier transforms in space combined with the help of a lowrank approximation of the space-wavenumber wave-propagator matrix. A lowrank approximation implies selecting a small set of representative spatial locations and a small set of representative wavenumbers. We present a mathematical derivation of this method, a description of the lowrank approximation algorithm and numerical examples that confirm the validity of the proposed approach. Wave extrapolation using lowrank approximation can be applied to seismic imaging by reverse-time migration in 3D heterogeneous isotropic or anisotropic media. © 2012 European Association of Geoscientists & Engineers.

  17. Approximation concepts for efficient structural synthesis

    Science.gov (United States)

    Schmit, L. A., Jr.; Miura, H.

    1976-01-01

    It is shown that efficient structural synthesis capabilities can be created by using approximation concepts to mesh finite element structural analysis methods with nonlinear mathematical programming techniques. The history of the application of mathematical programming techniques to structural design optimization problems is reviewed. Several rather general approximation concepts are described along with the technical foundations of the ACCESS 1 computer program, which implements several approximation concepts. A substantial collection of structural design problems involving truss and idealized wing structures is presented. It is concluded that since the basic ideas employed in creating the ACCESS 1 program are rather general, its successful development supports the contention that the introduction of approximation concepts will lead to the emergence of a new generation of practical and efficient, large scale, structural synthesis capabilities in which finite element analysis methods and mathematical programming algorithms will play a central role.

  18. Broadband Approximations for Doubly Curved Reflector Antenna

    Directory of Open Access Journals (Sweden)

    V. Schejbal

    2010-12-01

    Full Text Available The broadband approximations for shaped-beam doubly curved reflector antennas with primary feed (rectangular horn producing uniform amplitude and phase aperture distribution are derived and analyzed. They are very valuable for electromagnetic compatibility analyses both from electromagnetic interference and susceptibility point of view, because specialized more accurate methods such as physical optics are only used by antenna designers. To allow quick EMC analyses, typical values, beamwidth changes, sidelobe levels and aperture efficiencies are given for frequency changes approximately up to four times operating frequency. A comparison of approximated and measured patterns of doubly curved reflector antennas shows that the given approximation could be reliably used for analyses of pattern changes due to very broad frequency changes.

  19. Lattice quantum chromodynamics with approximately chiral fermions

    Energy Technology Data Exchange (ETDEWEB)

    Hierl, Dieter

    2008-05-15

    In this work we present Lattice QCD results obtained by approximately chiral fermions. We use the CI fermions in the quenched approximation to investigate the excited baryon spectrum and to search for the {theta}{sup +} pentaquark on the lattice. Furthermore we developed an algorithm for dynamical simulations using the FP action. Using FP fermions we calculate some LECs of chiral perturbation theory applying the epsilon expansion. (orig.)

  20. Selfconsistent approximations, symmetries and choice of representation

    OpenAIRE

    Leupold, Stefan

    2006-01-01

    In thermal field theory selfconsistent (Phi-derivable) approximations are used to improve (resum) propagators at the level of two-particle irreducible diagrams. At the same time vertices are treated at the bare level. Therefore such approximations typically violate the Ward identities connected to internal symmetries. Examples are presented how such violations can be tamed by a proper choice of representation for the fields which describe the system under consideration. These examples cover t...

  1. Adiabatic Approximation, Semiclassical Scattering, and Unidirectional Invisibility

    OpenAIRE

    Mostafazadeh, Ali

    2014-01-01

    arXiv:1401.4315v3 [quant-ph] 27 Feb 2014 Adiabatic Approximation, Semiclassical Scattering, and Unidirectional Invisibility Ali Mostafazadeh∗ Department of Mathematics, Ko¸c University, 34450 Sarıyer, Istanbul, Turkey Abstract The transfer matrix of a possibly complex and energy-dependent scattering potential can be identified with the S-matrix of a two-level time-dependent non-Hermitian Hamiltonian H( ). We show that the application of the adiabatic approximation ...

  2. Approximate Bayesian computation in population genetics.

    OpenAIRE

    Beaumont, Mark A; Zhang, Wenyang; Balding, David J.

    2002-01-01

    We propose a new method for approximate Bayesian statistical inference on the basis of summary statistics. The method is suited to complex problems that arise in population genetics, extending ideas developed in this setting by earlier authors. Properties of the posterior distribution of a parameter, such as its mean or density curve, are approximated without explicit likelihood calculations. This is achieved by fitting a local-linear regression of simulated parameter values on simulated summ...

  3. Approximation by Semigroups of Spherical Operators

    OpenAIRE

    Wang, Yuguang; Cao, Feilong

    2011-01-01

    This paper discusses the approximation by %semigroups of operators of class ($\\mathscr{C}_0$) on the sphere and focuses on a class of so called exponential-type multiplier operators. It is proved that such operators form a strongly continuous semigroup of contraction operators of class ($\\mathscr{C}_0$), from which the equivalence between approximation for these operators and $K$-functionals introduced by the operators is given. As examples, the constructed $r$-th Boolean of generalized spher...

  4. Polynomial approximation of functions in Sobolev spaces

    Science.gov (United States)

    Dupont, T.; Scott, R.

    1980-01-01

    Constructive proofs and several generalizations of approximation results of J. H. Bramble and S. R. Hilbert are presented. Using an averaged Taylor series, we represent a function as a polynomial plus a remainder. The remainder can be manipulated in many ways to give different types of bounds. Approximation of functions in fractional order Sobolev spaces is treated as well as the usual integer order spaces and several nonstandard Sobolev-like spaces.

  5. Introduction to the Maxwell Garnett approximation: tutorial.

    Science.gov (United States)

    Markel, Vadim A

    2016-07-01

    This tutorial is devoted to the Maxwell Garnett approximation and related theories. Topics covered in this first, introductory part of the tutorial include the Lorentz local field correction, the Clausius-Mossotti relation and its role in the modern numerical technique known as the discrete dipole approximation, the Maxwell Garnett mixing formula for isotropic and anisotropic media, multicomponent mixtures and the Bruggeman equation, the concept of smooth field, and Wiener and Bergman-Milton bounds. PMID:27409680

  6. A Ballistic Monte Carlo Approximation of {\\pi}

    CERN Document Server

    Dumoulin, Vincent

    2014-01-01

    We compute a Monte Carlo approximation of {\\pi} using importance sampling with shots coming out of a Mossberg 500 pump-action shotgun as the proposal distribution. An approximated value of 3.136 is obtained, corresponding to a 0.17% error on the exact value of {\\pi}. To our knowledge, this represents the first attempt at estimating {\\pi} using such method, thus opening up new perspectives towards computing mathematical constants using everyday tools.

  7. Time Stamps for Fixed-Point Approximation

    DEFF Research Database (Denmark)

    Damian, Daniela

    2001-01-01

    Time stamps were introduced in Shivers's PhD thesis for approximating the result of a control-flow analysis. We show them to be suitable for computing program analyses where the space of results (e.g., control-flow graphs) is large. We formalize time-stamping as a top-down, fixed......-point approximation algorithm which maintains a single copy of intermediate results. We then prove the correctness of this algorithm....

  8. Parallel local approximation MCMC for expensive models

    OpenAIRE

    Conrad, Patrick; Davis, Andrew; Marzouk, Youssef; Pillai, Natesh; Smith, Aaron

    2016-01-01

    Performing Bayesian inference via Markov chain Monte Carlo (MCMC) can be exceedingly expensive when posterior evaluations invoke the evaluation of a computationally expensive model, such as a system of partial differential equations. In recent work [Conrad et al. JASA 2015, arXiv:1402.1694] we described a framework for constructing and refining local approximations of such models during an MCMC simulation. These posterior--adapted approximations harness regularity of the model to reduce the c...

  9. The Actinide Transition Revisited by Gutzwiller Approximation

    Science.gov (United States)

    Xu, Wenhu; Lanata, Nicola; Yao, Yongxin; Kotliar, Gabriel

    2015-03-01

    We revisit the problem of the actinide transition using the Gutzwiller approximation (GA) in combination with the local density approximation (LDA). In particular, we compute the equilibrium volumes of the actinide series and reproduce the abrupt change of density found experimentally near plutonium as a function of the atomic number. We discuss how this behavior relates with the electron correlations in the 5 f states, the lattice structure, and the spin-orbit interaction. Our results are in good agreement with the experiments.

  10. Lattice quantum chromodynamics with approximately chiral fermions

    International Nuclear Information System (INIS)

    In this work we present Lattice QCD results obtained by approximately chiral fermions. We use the CI fermions in the quenched approximation to investigate the excited baryon spectrum and to search for the Θ+ pentaquark on the lattice. Furthermore we developed an algorithm for dynamical simulations using the FP action. Using FP fermions we calculate some LECs of chiral perturbation theory applying the epsilon expansion. (orig.)

  11. Staying Thermal with Hartree Ensemble Approximations

    CERN Document Server

    Salle, M; Vink, Jeroen C

    2000-01-01

    Using Hartree ensemble approximations to compute the real time dynamics of scalar fields in 1+1 dimension, we find that with suitable initial conditions, approximate thermalization is achieved much faster than found in our previous work. At large times, depending on the interaction strength and temperature, the particle distribution slowly changes: the Bose-Einstein distribution of the particle densities develops classical features. We also discuss variations of our method which are numerically more efficient.

  12. Approximations for the Erlang Loss Function

    DEFF Research Database (Denmark)

    Mejlbro, Leif

    1998-01-01

    Theoretically, at least three formulae are needed for arbitrarily good approximates of the Erlang Loss Function. In the paper, for convenience five formulae are presented guaranteeing a relative error <1E-2, and methods are indicated for improving this bound.......Theoretically, at least three formulae are needed for arbitrarily good approximates of the Erlang Loss Function. In the paper, for convenience five formulae are presented guaranteeing a relative error

  13. Approximate Assertional Reasoning Over Expressive Ontologies

    OpenAIRE

    Tserendorj, Tuvshintur

    2010-01-01

    In this thesis, approximate reasoning methods for scalable assertional reasoning are provided whose computational properties can be established in a well-understood way, namely in terms of soundness and completeness, and whose quality can be analyzed in terms of statistical measurements, namely recall and precision. The basic idea of these approximate reasoning methods is to speed up reasoning by trading off the quality of reasoning results against increased speed.

  14. Phase Transitions for Greedy Sparse Approximation Algorithms

    OpenAIRE

    Blanchard, Jeffrey D.; Cartis, Coralia; Tanner, Jared; Thompson, Andrew

    2010-01-01

    A major enterprise in compressed sensing and sparse approximation is the design and analysis of computationally tractable algorithms for recovering sparse, exact or approximate, solutions of underdetermined linear systems of equations. Many such algorithms have now been proven to have optimal-order uniform recovery guarantees using the ubiquitous Restricted Isometry Property (RIP) (Candes and Tao (2005) [11]). However, without specifying a matrix, or class of matrices, it is unclear when the ...

  15. Constrained Optimization via Stochastic approximation with a simultaneous perturbation gradient approximation

    DEFF Research Database (Denmark)

    Sadegh, Payman

    1997-01-01

    This paper deals with a projection algorithm for stochastic approximation using simultaneous perturbation gradient approximation for optimization under inequality constraints where no direct gradient of the loss function is available and the inequality constraints are given as explicit functions...

  16. Trajectory averaging for stochastic approximation MCMC algorithms

    KAUST Repository

    Liang, Faming

    2010-10-01

    The subject of stochastic approximation was founded by Robbins and Monro [Ann. Math. Statist. 22 (1951) 400-407]. After five decades of continual development, it has developed into an important area in systems control and optimization, and it has also served as a prototype for the development of adaptive algorithms for on-line estimation and control of stochastic systems. Recently, it has been used in statistics with Markov chain Monte Carlo for solving maximum likelihood estimation problems and for general simulation and optimizations. In this paper, we first show that the trajectory averaging estimator is asymptotically efficient for the stochastic approximation MCMC (SAMCMC) algorithm under mild conditions, and then apply this result to the stochastic approximation Monte Carlo algorithm [Liang, Liu and Carroll J. Amer. Statist. Assoc. 102 (2007) 305-320]. The application of the trajectory averaging estimator to other stochastic approximationMCMC algorithms, for example, a stochastic approximation MLE algorithm for missing data problems, is also considered in the paper. © Institute of Mathematical Statistics, 2010.

  17. Tree-fold loop approximation of AMD

    Energy Technology Data Exchange (ETDEWEB)

    Ono, Akira [Tohoku Univ., Sendai (Japan). Faculty of Science

    1997-05-01

    AMD (antisymmetrized molecular dynamics) is a frame work for describing a wave function of nucleon multi-body system by Slater determinant of Gaussian wave flux, and a theory for integrally describing a wide range of nuclear reactions such as intermittent energy heavy ion reaction, nucleon incident reaction and so forth. The aim of this study is induction on approximation equation of expected value, {nu}, in correlation capable of calculation with time proportional A (exp 3) (or lower), and to make AMD applicable to the heavier system such as Au+Au. As it must be avoided to break characteristics of AMD, it needs not to be anxious only by approximating the {nu}-value. However, in order to give this approximation any meaning, error of this approximation will have to be sufficiently small in comparison with bond energy of atomic nucleus and smaller than 1 MeV/nucleon. As the absolute expected value in correlation may be larger than 50 MeV/nucleon, the approximation is required to have a high accuracy within 2 percent. (G.K.)

  18. An Approximation Algorithm for #k-SAT

    CERN Document Server

    Thurley, Marc

    2011-01-01

    We present a simple randomized algorithm that approximates the number of satisfying assignments of Boolean formulas in conjunctive normal form. To the best of our knowledge this is the first algorithm which approximates #k-SAT for any k >= 3 within a running time that is not only non-trivial, but also significantly better than that of the currently fastest exact algorithms for the problem. More precisely, our algorithm is a randomized approximation scheme whose running time depends polynomially on the error tolerance and is mildly exponential in the number n of variables of the input formula. For example, even stipulating sub-exponentially small error tolerance, the number of solutions to 3-CNF input formulas can be approximated in time O(1.5366^n). For 4-CNF input the bound increases to O(1.6155^n). We further show how to obtain upper and lower bounds on the number of solutions to a CNF formula in a controllable way. Relaxing the requirements on the quality of the approximation, on k-CNF input we obtain sign...

  19. Fast wavelet based sparse approximate inverse preconditioner

    Energy Technology Data Exchange (ETDEWEB)

    Wan, W.L. [Univ. of California, Los Angeles, CA (United States)

    1996-12-31

    Incomplete LU factorization is a robust preconditioner for both general and PDE problems but unfortunately not easy to parallelize. Recent study of Huckle and Grote and Chow and Saad showed that sparse approximate inverse could be a potential alternative while readily parallelizable. However, for special class of matrix A that comes from elliptic PDE problems, their preconditioners are not optimal in the sense that independent of mesh size. A reason may be that no good sparse approximate inverse exists for the dense inverse matrix. Our observation is that for this kind of matrices, its inverse entries typically have piecewise smooth changes. We can take advantage of this fact and use wavelet compression techniques to construct a better sparse approximate inverse preconditioner. We shall show numerically that our approach is effective for this kind of matrices.

  20. Approximating light rays in the Schwarzschild field

    CERN Document Server

    Semerak, Oldrich

    2014-01-01

    A short formula is suggested which approximates photon trajectories in the Schwarzschild field better than other simple prescriptions from the literature. We compare it with various "low-order competitors", namely with those following from exact formulas for small $M$, with one of the results based on pseudo-Newtonian potentials, with a suitably adjusted hyperbola, and with the effective and often employed approximation by Beloborodov. Our main concern is the shape of the photon trajectories at finite radii, yet asymptotic behaviour is also discussed, important for lensing. An example is attached indicating that the newly suggested approximation is usable--and very accurate--for practical solving of the ray-deflection exercise.

  1. On the approximate zero of Newton method

    Institute of Scientific and Technical Information of China (English)

    黄正达

    2003-01-01

    A judgment criterion to guarantee a point to be a Chen' s approximate zero of Newton method for solving nonlinear equation is sought by dominating sequence techniques. The criterion is based on the fact that the dominating function may have only one simple positive zero, assuming that the operator is weak Lipschitz continuous, which is much more relaxed and can be checked much more easily than Lipschitz continuous in practice. It is demonstrated that a Chen' s approximate zero may not be a Smale' s approximate zero. The error estimate obtained indicated the convergent order when we use |f(x) | < ε to stop computation in software.The result can also be applied for solving partial derivative and integration equations.

  2. Large scale systems approximation: Analysis and control

    International Nuclear Information System (INIS)

    This work concerns the study of the approximation of high dimensional systems by low order models. This approximation is defined by aggregation methods, the method based on singular perturbations and a relatively recent method. This later one is formulated in a particular representation of the system and is called balanced realisation method. The application of the approximation is then studied in the case of suboptimal control theory for the several defined models. The results of numerical simulation for the analysis and control are carried on two examples defined by a model of a nuclear reactor core of order nine and a steam generator of a fuel station of order twenty and permitted to develop a comparative study of the performances of the different methods analysed in the case of open loop and closed loop models

  3. Numerical approximation of partial differential equations

    CERN Document Server

    Bartels, Sören

    2016-01-01

    Finite element methods for approximating partial differential equations have reached a high degree of maturity, and are an indispensible tool in science and technology. This textbook aims at providing a thorough introduction to the construction, analysis, and implementation of finite element methods for model problems arising in continuum mechanics. The first part of the book discusses elementary properties of linear partial differential equations along with their basic numerical approximation, the functional-analytical framework for rigorously establishing existence of solutions, and the construction and analysis of basic finite element methods. The second part is devoted to the optimal adaptive approximation of singularities and the fast iterative solution of linear systems of equations arising from finite element discretizations. In the third part, the mathematical framework for analyzing and discretizing saddle-point problems is formulated, corresponding finte element methods are analyzed, and particular ...

  4. Approximate Bayesian Computation: a nonparametric perspective

    CERN Document Server

    Blum, Michael

    2010-01-01

    Approximate Bayesian Computation is a family of likelihood-free inference techniques that are well-suited to models defined in terms of a stochastic generating mechanism. In a nutshell, Approximate Bayesian Computation proceeds by computing summary statistics s_obs from the data and simulating summary statistics for different values of the parameter theta. The posterior distribution is then approximated by an estimator of the conditional density g(theta|s_obs). In this paper, we derive the asymptotic bias and variance of the standard estimators of the posterior distribution which are based on rejection sampling and linear adjustment. Additionally, we introduce an original estimator of the posterior distribution based on quadratic adjustment and we show that its bias contains a fewer number of terms than the estimator with linear adjustment. Although we find that the estimators with adjustment are not universally superior to the estimator based on rejection sampling, we find that they can achieve better perfor...

  5. On approximation of Markov binomial distributions

    CERN Document Server

    Xia, Aihua; 10.3150/09-BEJ194

    2010-01-01

    For a Markov chain $\\mathbf{X}=\\{X_i,i=1,2,...,n\\}$ with the state space $\\{0,1\\}$, the random variable $S:=\\sum_{i=1}^nX_i$ is said to follow a Markov binomial distribution. The exact distribution of $S$, denoted $\\mathcal{L}S$, is very computationally intensive for large $n$ (see Gabriel [Biometrika 46 (1959) 454--460] and Bhat and Lal [Adv. in Appl. Probab. 20 (1988) 677--680]) and this paper concerns suitable approximate distributions for $\\mathcal{L}S$ when $\\mathbf{X}$ is stationary. We conclude that the negative binomial and binomial distributions are appropriate approximations for $\\mathcal{L}S$ when $\\operatorname {Var}S$ is greater than and less than $\\mathbb{E}S$, respectively. Also, due to the unique structure of the distribution, we are able to derive explicit error estimates for these approximations.

  6. The unitary convolution approximation for heavy ions

    CERN Document Server

    Grande, P L

    2002-01-01

    The convolution approximation for the impact-parameter dependent energy loss is reviewed with emphasis on the determination of the stopping force for heavy projectiles. In this method, the energy loss in different impact-parameter regions is well determined and interpolated smoothly. The physical inputs of the model are the projectile-screening function (in the case of dressed ions), the electron density and oscillators strengths of the target atoms. Moreover, the convolution approximation, in the perturbative mode (called PCA), yields remarkable agreement with full semi-classical-approximation (SCA) results for bare as well as for screened ions at all impact parameters. In the unitary mode (called UCA), the method contains some higher-order effects (yielding in some cases rather good agreement with full coupled-channel calculations) and approaches the classical regime similar as the Bohr model for large perturbations (Z/v>>1). The results are then used to compare with experimental values of the non-equilibri...

  7. Approximately-Balanced Drilling in Daqing Oilfield

    Institute of Scientific and Technical Information of China (English)

    Xia Bairu; Zheng Xiuhua; Li Guoqing; Tian Tuo

    2004-01-01

    The Daqing oilfield is a multilayered heterogeneous oil field where the pressure are different in the same vertical profile causing many troubles to the adjustment well drillings. The approximately-balanced drilling technique has been developed and proved to be efficient and successful in Daqing oilfield. This paper discusses the application of approximately-balanced drilling technique under the condition of multilayered pressure in Daqing oilfield, including the prediction of formation pressure, the pressure discharge technique for the drilling well and the control of the density of drilling fluid.

  8. Local density approximations from finite systems

    CERN Document Server

    Entwistle, Mike; Wetherell, Jack; Longstaff, Bradley; Ramsden, James; Godby, Rex

    2016-01-01

    The local density approximation (LDA) constructed through quantum Monte Carlo calculations of the homogeneous electron gas (HEG) is the most common approximation to the exchange-correlation functional in density functional theory. We introduce an alternative set of LDAs constructed from slab-like systems of one, two and three electrons that resemble the HEG within a finite region, and illustrate the concept in one dimension. Comparing with the exact densities and Kohn-Sham potentials for various test systems, we find that the LDAs give a good account of the self-interaction correction, but are less reliable when correlation is stronger or currents flow.

  9. Weisskopf--Wigner approximation in atomic physics

    International Nuclear Information System (INIS)

    Several approximations involved in the usual Weisskopf-Wigner treatment of the emission of light by an atom are investigated. The system considered is a recoilless, nonrelativistic hydrogen atom interacting with a quantized electromagnetic field, in dipole approximation (with a nonrelativistic cutoff in momentum space). Since only electric dipole waves interact with the atom, the Hamiltonian can be expressed in a simple one-dimensional form. The time evolution of the system is determined by resolvent operator techniques. The method goes beyond the analysis by Van Hove and Hugenholtz, allowing one to treat also fields of finite intensity in the infinite-volume limit. A comparison between this and other techniques is made

  10. Extending the Eikonal Approximation to Low Energy

    CERN Document Server

    Capel, Pierre; Ogata, Kazuyuki

    2014-01-01

    E-CDCC and DEA, two eikonal-based reaction models are compared to CDCC at low energy (e.g. 20AMeV) to study their behaviour in the regime at which the eikonal approximation is supposed to fail. We confirm that these models lack the Coulomb deflection of the projectile by the target. We show that a hybrid model, built on the CDCC framework at low angular momenta and the eikonal approximation at larger angular momenta gives a perfect agreement with CDCC. An empirical shift in impact parameter can also be used reliably to simulate this missing Coulomb deflection.

  11. Computing Nash Equilibria: Approximation and Smoothed Complexity

    OpenAIRE

    Chen, Xi; Deng, Xiaotie; Teng, Shang-Hua

    2006-01-01

    We show that the BIMATRIX game does not have a fully polynomial-time approximation scheme, unless PPAD is in P. In other words, no algorithm with time polynomial in n and 1/\\epsilon can compute an \\epsilon-approximate Nash equilibrium of an n by nbimatrix game, unless PPAD is in P. Instrumental to our proof, we introduce a new discrete fixed-point problem on a high-dimensional cube with a constant side-length, such as on an n-dimensional cube with side-length 7, and show that they are PPAD-co...

  12. Generalized companion matrix for approximate GCD

    CERN Document Server

    Boito, Paola

    2011-01-01

    We study a variant of the univariate approximate GCD problem, where the coe?- cients of one polynomial f(x)are known exactly, whereas the coe?cients of the second polynomial g(x)may be perturbed. Our approach relies on the properties of the matrix which describes the operator of multiplication by gin the quotient ring C[x]=(f). In particular, the structure of the null space of the multiplication matrix contains all the essential information about GCD(f; g). Moreover, the multiplication matrix exhibits a displacement structure that allows us to design a fast algorithm for approximate GCD computation with quadratic complexity w.r.t. polynomial degrees.

  13. Static correlation beyond the random phase approximation

    DEFF Research Database (Denmark)

    Olsen, Thomas; Thygesen, Kristian Sommer

    2014-01-01

    derived from Hedin's equations (Random Phase Approximation (RPA), Time-dependent Hartree-Fock (TDHF), Bethe-Salpeter equation (BSE), and Time-Dependent GW) all reproduce the correct dissociation limit. We also show that the BSE improves the correlation energies obtained within RPA and TDHF significantly...... for intermediate binding distances. A Hubbard model for the dimer allows us to obtain exact analytical results for the various approximations, which is readily compared with the exact diagonalization of the model. Moreover, the model is shown to reproduce all the qualitative results from the ab initio...

  14. BEST APPROXIMATION BY DOWNWARD SETS WITH APPLICATIONS

    Institute of Scientific and Technical Information of China (English)

    H.Mohebi; A. M. Rubinov

    2006-01-01

    We develop a theory of downward sets for a class of normed ordered spaces. We study best approximation in a normed ordered space X by elements of downward sets, and give necessary and sufficient conditions for any element of best approximation by a closed downward subset of X. We also characterize strictly downward subsets of X, and prove that a downward subset of X is strictly downward if and only if each its boundary point is Chebyshev. The results obtained are used for examination of some Chebyshev pairs (W,x), where x ∈ X and W is a closed downward subset of X.

  15. Nonlinear approximation in alpha-modulation spaces

    DEFF Research Database (Denmark)

    Borup, Lasse; Nielsen, Morten

    2006-01-01

    The α-modulation spaces are a family of spaces that contain the Besov and modulation spaces as special cases. In this paper we prove that brushlet bases can be constructed to form unconditional and even greedy bases for the α-modulation spaces. We study m -term nonlinear approximation with brushlet...

  16. An approximate classical unimolecular reaction rate theory

    Science.gov (United States)

    Zhao, Meishan; Rice, Stuart A.

    1992-05-01

    We describe a classical theory of unimolecular reaction rate which is derived from the analysis of Davis and Gray by use of simplifying approximations. These approximations concern the calculation of the locations of, and the fluxes of phase points across, the bottlenecks to fragmentation and to intramolecular energy transfer. The bottleneck to fragment separation is represented as a vibration-rotation state dependent separatrix, which approximation is similar to but extends and improves the approximations for the separatrix introduced by Gray, Rice, and Davis and by Zhao and Rice. The novel feature in our analysis is the representation of the bottlenecks to intramolecular energy transfer as dividing surfaces in phase space; the locations of these dividing surfaces are determined by the same conditions as locate the remnants of robust tori with frequency ratios related to the golden mean (in a two degree of freedom system these are the cantori). The flux of phase points across each dividing surface is calculated with an analytic representation instead of a stroboscopic mapping. The rate of unimolecular reaction is identified with the net rate at which phase points escape from the region of quasiperiodic bounded motion to the region of free fragment motion by consecutively crossing the dividing surfaces for intramolecular energy exchange and the separatrix. This new theory generates predictions of the rates of predissociation of the van der Waals molecules HeI2, NeI2 and ArI2 which are in very good agreement with available experimental data.

  17. Approximate Equilibrium Problems and Fixed Points

    Directory of Open Access Journals (Sweden)

    H. Mazaheri

    2013-01-01

    Full Text Available We find a common element of the set of fixed points of a map and the set of solutions of an approximate equilibrium problem in a Hilbert space. Then, we show that one of the sequences weakly converges. Also we obtain some theorems about equilibrium problems and fixed points.

  18. Median Approximations for Genomes Modeled as Matrices.

    Science.gov (United States)

    Zanetti, Joao Paulo Pereira; Biller, Priscila; Meidanis, Joao

    2016-04-01

    The genome median problem is an important problem in phylogenetic reconstruction under rearrangement models. It can be stated as follows: Given three genomes, find a fourth that minimizes the sum of the pairwise rearrangement distances between it and the three input genomes. In this paper, we model genomes as matrices and study the matrix median problem using the rank distance. It is known that, for any metric distance, at least one of the corners is a [Formula: see text]-approximation of the median. Our results allow us to compute up to three additional matrix median candidates, all of them with approximation ratios at least as good as the best corner, when the input matrices come from genomes. We also show a class of instances where our candidates are optimal. From the application point of view, it is usually more interesting to locate medians farther from the corners, and therefore, these new candidates are potentially more useful. In addition to the approximation algorithm, we suggest a heuristic to get a genome from an arbitrary square matrix. This is useful to translate the results of our median approximation algorithm back to genomes, and it has good results in our tests. To assess the relevance of our approach in the biological context, we ran simulated evolution tests and compared our solutions to those of an exact DCJ median solver. The results show that our method is capable of producing very good candidates. PMID:27072561

  19. WEIGHTED APPROXIMATION ON SZASZ-TYPE OPERATORS

    Institute of Scientific and Technical Information of China (English)

    Feng Guo

    2003-01-01

    In this paper, we use weighted modules ω2φλ (f,t)w to study the pointwise approximation on Szász-type operators, and obtain the direct and converse theorem, as well as characterizations of the pointwise approxi- mation of Jacobi-weighted Szász-type operators.

  20. Semi-classical approximation and microcanonical ensemble

    International Nuclear Information System (INIS)

    For quantum mechanical systems with spherically symmetric potential the improved W.K.B. approximation of Elworthy and Truman corresponds to the classical microcanonical ensemble in the limit where (h/2π) goes to zero, at least for small time. (orig.)

  1. Eignets for function approximation on manifolds

    CERN Document Server

    Mhaskar, H N

    2009-01-01

    Let $\\XX$ be a compact, smooth, connected, Riemannian manifold without boundary, $G:\\XX\\times\\XX\\to \\RR$ be a kernel. Analogous to a radial basis function network, an eignet is an expression of the form $\\sum_{j=1}^M a_jG(\\circ,y_j)$, where $a_j\\in\\RR$, $y_j\\in\\XX$, $1\\le j\\le M$. We describe a deterministic, universal algorithm for constructing an eignet for approximating functions in $L^p(\\mu;\\XX)$ for a general class of measures $\\mu$ and kernels $G$. Our algorithm yields linear operators. Using the minimal separation amongst the centers $y_j$ as the cost of approximation, we give modulus of smoothness estimates for the degree of approximation by our eignets, and show by means of a converse theorem that these are the best possible for every \\emph{individual function}. We also give estimates on the coefficients $a_j$ in terms of the norm of the eignet. Finally, we demonstrate that if any sequence of eignets satisfies the optimal estimates for the degree of approximation of a smooth function, measured in ter...

  2. Approximating the Balanced Minimum Evolution Problem

    CERN Document Server

    Fiorini, Samuel

    2011-01-01

    We prove a strong inapproximability result for the Balanced Minimum Evolution Problem. Our proof also implies that the problem remains NP-hard even when restricted to metric instances. Furthermore, we give a MST-based 2-approximation algorithm for the problem for such instances.

  3. Optical bistability without the rotating wave approximation

    Energy Technology Data Exchange (ETDEWEB)

    Sharaby, Yasser A., E-mail: Yasser_Sharaby@hotmail.co [Physics Department, Faculty of Applied Sciences, Suez Canal University, Suez (Egypt); Joshi, Amitabh, E-mail: ajoshi@eiu.ed [Department of Physics, Eastern Illinois University, Charleston, IL 61920 (United States); Hassan, Shoukry S., E-mail: Shoukryhassan@hotmail.co [Mathematics Department, College of Science, University of Bahrain, P.O. Box 32038 (Bahrain)

    2010-04-26

    Optical bistability for two-level atomic system in a ring cavity is investigated outside the rotating wave approximation (RWA) using non-autonomous Maxwell-Bloch equations with Fourier decomposition up to first harmonic. The first harmonic output field component exhibits reversed or closed loop bistability simultaneously with the usual (anti-clockwise) bistability in the fundamental field component.

  4. Kravchuk functions for the finite oscillator approximation

    Science.gov (United States)

    Atakishiyev, Natig M.; Wolf, Kurt Bernardo

    1995-01-01

    Kravchuk orthogonal functions - Kravchuk polynomials multiplied by the square root of the weight function - simplify the inversion algorithm for the analysis of discrete, finite signals in harmonic oscillator components. They can be regarded as the best approximation set. As the number of sampling points increases, the Kravchuk expansion becomes the standard oscillator expansion.

  5. An approximate analytical approach to resampling averages

    DEFF Research Database (Denmark)

    Malzahn, Dorthe; Opper, M.

    2004-01-01

    Using a novel reformulation, we develop a framework to compute approximate resampling data averages analytically. The method avoids multiple retraining of statistical models on the samples. Our approach uses a combination of the replica "trick" of statistical physics and the TAP approach for appr...

  6. Approximations in the PE-method

    DEFF Research Database (Denmark)

    Arranz, Marta Galindo

    1996-01-01

    Two differenct sources of errors may occur in the implementation of the PE methods; a phase error introduced in the approximation of a pseudo-differential operator and an amplitude error generated from the starting field. First, the inherent phase errors introduced in the solution are analyzed...

  7. Quantum electrodynamics in a classical approximation, 1

    International Nuclear Information System (INIS)

    Quantum electrodynamics is formulated in a classical approximation. A quantum mechanical proper-time is employed as a useful parameter, which enables us to elucidate the relationship between quantum electrodynamics and classical electrodynamics. The classical motion of a charged particle is realized as an asymptotic limit of quantum electrodynamics. (author)

  8. Revisiting Twomey's approximation for peak supersaturation

    Directory of Open Access Journals (Sweden)

    B. J. Shipway

    2014-10-01

    Full Text Available Twomey's seminal 1959 paper provided lower and upper bound approximations to the estimation of peak supersaturation within an updraft and thus provides the first closed expression for the number of nucleated cloud droplets. The form of this approximation is simple, but provides a surprisingly good estimate and has subsequently been employed in more sophisticated treatments of nucleation parametrization. In the current paper, we revisit the lower bound approximation of Twomey and make a small adjustment which can be used to obtain a more accurate calculation of peak supersaturation under all potential aerosol loadings and thermodynamic conditions. In order to make full use of this improved approximation, the underlying integro-differential equation for supersaturation evolution and the condition for calculating peak supersaturation are examined. A simple rearrangement of the algebra allows for an expression to be written down which can then be solved with a single lookup table with only one independent variable for an underlying lognormal aerosol population. Multimode aerosol with only N different dispersion characteristics require only N of these one-dimensional lookup tables. No additional information is required in the lookup table to deal with additional chemical, physical or thermodynamic properties. The resulting implementation provides a relatively simple, yet computationally cheap and very accurate physically-based parametrization of droplet nucleation for use in climate and NWP models.

  9. On operators with bounded approximation property

    OpenAIRE

    Reinov, Oleg

    2013-01-01

    It is known that any separable Banach space with BAP is a complemented subspace of a Banach space with a basis. We show that every operator with bounded approximation property, acting from a separable Banach space, can be factored through a Banach space with a basis.

  10. Approximation by Penultimate Extreme Value Distributions

    NARCIS (Netherlands)

    L.F.M. de Haan (Laurens)

    1998-01-01

    textabstractn certain cases the distribution of the normalized maximum of a sample can be better approximated by a sequence of different extreme value distributions than by the final one. We show that these cases are rather restricted and that the possible improvement is not spectacular.

  11. Approximate fixed point of Reich operator

    Directory of Open Access Journals (Sweden)

    M. Saha

    2013-01-01

    Full Text Available In the present paper, we study the existence of approximate fixed pointfor Reich operator together with the property that the ε-fixed points are concentrated in a set with the diameter tends to zero if ε $to$ > 0.

  12. An Approximation of Ultra-Parabolic Equations

    Directory of Open Access Journals (Sweden)

    Allaberen Ashyralyev

    2012-01-01

    Full Text Available The first and second order of accuracy difference schemes for the approximate solution of the initial boundary value problem for ultra-parabolic equations are presented. Stability of these difference schemes is established. Theoretical results are supported by the result of numerical examples.

  13. Approximation Algorithms for Model-Based Diagnosis

    NARCIS (Netherlands)

    Feldman, A.B.

    2010-01-01

    Model-based diagnosis is an area of abductive inference that uses a system model, together with observations about system behavior, to isolate sets of faulty components (diagnoses) that explain the observed behavior, according to some minimality criterion. This thesis presents greedy approximation a

  14. Neutral kaons without Weisskopf-Wigner approximation

    International Nuclear Information System (INIS)

    The model-independent formalism is constructed to describe decays of mixed particles without using the Weisskopf-Wigner approximation (WWA). Limitations due to various symmetries are traced for neutral K mesons. As an application we show that effects of CPT violation and going beyond WWA may be separated and studied independently. 16 refs

  15. Neutral Kaons without Weisskopf-Wigner Approximation

    OpenAIRE

    Azimov, Ya. I.

    1995-01-01

    The model-independent formalism is constructed to describe decays of mixed particles without using the Weisskopf-Wigner approximation. Limitations due to various symmetries are traced for neutral $K-$mesons. As an application we show that effects of $CPT-$violation and going beyond WWA may be separated and studied independently.

  16. Approximating the DGP of China's Quarterly GDP

    NARCIS (Netherlands)

    Ph.H.B.F. Franses (Philip Hans); H. Mees (Heleen)

    2010-01-01

    textabstractWe demonstrate that the data generating process (DGP) of China’s cumulated quarterly Gross Domestic Product (GDP, current prices), as it is reported by the National Bureau of Statistics of China, can be (very closely) approximated by a simple rule. This rule is that annual growth in any

  17. Empirical progress and nomic truth approximation revisited

    NARCIS (Netherlands)

    Kuipers, Theodorus

    2014-01-01

    In my From Instrumentalism to Constructive Realism (2000) I have shown how an instrumentalist account of empirical progress can be related to nomic truth approximation. However, it was assumed that a strong notion of nomic theories was needed for that analysis. In this paper it is shown, in terms of

  18. Image Compression Via a Fast DCT Approximation

    NARCIS (Netherlands)

    Bayer, F. M.; Cintra, R. J.

    2010-01-01

    Discrete transforms play an important role in digital signal processing. In particular, due to its transform domain energy compaction properties, the discrete cosine transform (DCT) is pivotal in many image processing problems. This paper introduces a numerical approximation method for the DCT based

  19. Double unresolved approximations to multiparton scattering amplitudes

    International Nuclear Information System (INIS)

    We present approximations to tree-level multiparton scattering amplitudes which are appropriate when two partons are unresolved. These approximations are required for the analytic isolation of infrared singularities of n+2 parton scattering processes contributing to the next-to-next-to-leading order corrections to n jet cross sections. In each case the colour ordered matrix elements factorise and yield a function containing the singular factors multiplying the n-parton amplitudes. When the unresolved particles are colour unconnected, the approximations are simple products of the familiar eikonal and Altarelli-Parisi splitting functions used to describe single unresolved emission. However, when the unresolved particles are colour connected the factorisation is more complicated and we introduce new and general functions to describe the triple collinear and soft/collinear limits in addition to the known double soft gluon limits of Berends and Giele. As expected the triple collinear splitting functions obey an N=1 SUSY identity. To illustrate the use of these double unresolved approximations, we have examined the singular limits of the tree-level matrix elements for e+e- →5 partons when only three partons are resolved. When integrated over the unresolved regions of phase space, these expressions will be of use in evaluating the O(αs3) corrections to the three-jet rate in electron-positron annihilation. (orig.)

  20. $\\Phi$-derivable approximations in gauge theories

    CERN Document Server

    Arrizabalaga, A

    2003-01-01

    We discuss the method of $\\Phi$-derivable approximations in gauge theories. There, two complications arise, namely the violation of Bose symmetry in correlation functions and the gauge dependence. For the latter we argue that the error introduced by the gauge dependent terms is controlled, therefore not invalidating the method.

  1. Approximate Furthest Neighbor in High Dimensions

    DEFF Research Database (Denmark)

    Pagh, Rasmus; Silvestri, Francesco; Sivertsen, Johan von Tangen;

    2015-01-01

    -dimensional Euclidean space. We build on the technique of Indyk (SODA 2003), storing random projections to provide sublinear query time for AFN. However, we introduce a different query algorithm, improving on Indyk’s approximation factor and reducing the running time by a logarithmic factor. We also present a variation...

  2. Virial expansion coefficients in the harmonic approximation

    DEFF Research Database (Denmark)

    R. Armstrong, J.; Zinner, Nikolaj Thomas; V. Fedorov, D.;

    2012-01-01

    The virial expansion method is applied within a harmonic approximation to an interacting N-body system of identical fermions. We compute the canonical partition functions for two and three particles to get the two lowest orders in the expansion. The energy spectrum is carefully interpolated...

  3. ON BEST SIMULTANEOUS APPROXIMATION IN QUOTIENT SPACES

    Institute of Scientific and Technical Information of China (English)

    M. Iranmanesh; H. Mohebi

    2007-01-01

    We assume that X is a normed linear space, W and M are subspaces of X.We develop a theory of best simultaneous approximation in quotient spaces and introduce equivalent assertions between the subspaces W and W + M and the quotient space W/M.

  4. Markov operators, positive semigroups and approximation processes

    CERN Document Server

    Altomare, Francesco; Leonessa, Vita; Rasa, Ioan

    2015-01-01

    In recent years several investigations have been devoted to the study of large classes of (mainly degenerate) initial-boundary value evolution problems in connection with the possibility to obtain a constructive approximation of the associated positive C_0-semigroups. In this research monograph we present the main lines of a theory which finds its root in the above-mentioned research field.

  5. Improved Approximations for Multiprocessor Scheduling Under Uncertainty

    CERN Document Server

    Crutchfield, Christopher; Fineman, Jeremy T; Karger, David R; Scott, Jacob

    2008-01-01

    This paper presents improved approximation algorithms for the problem of multiprocessor scheduling under uncertainty, or SUU, in which the execution of each job may fail probabilistically. This problem is motivated by the increasing use of distributed computing to handle large, computationally intensive tasks. In the SUU problem we are given n unit-length jobs and m machines, a directed acyclic graph G of precedence constraints among jobs, and unrelated failure probabilities q_{ij} for each job j when executed on machine i for a single timestep. Our goal is to find a schedule that minimizes the expected makespan, which is the expected time at which all jobs complete. Lin and Rajaraman gave the first approximations for this NP-hard problem for the special cases of independent jobs, precedence constraints forming disjoint chains, and precedence constraints forming trees. In this paper, we present asymptotically better approximation algorithms. In particular, we give an O(loglog min(m,n))-approximation for indep...

  6. On surface approximation using developable surfaces

    DEFF Research Database (Denmark)

    Chen, H. Y.; Lee, I. K.; Leopoldseder, S.;

    1998-01-01

    We introduce a method for approximating a given surface by a developable surface. It will be either a G_1 surface consisting of pieces of cones or cylinders of revolution or a G_r NURBS developable surface. Our algorithm will also deal properly with the problems of reverse engineering and produce...

  7. On surface approximation using developable surfaces

    DEFF Research Database (Denmark)

    Chen, H. Y.; Lee, I. K.; Leopoldseder, s.;

    1999-01-01

    We introduce a method for approximating a given surface by a developable surface. It will be either a G(1) surface consisting of pieces of cones or cylinders of revolution or a G(r) NURBS developable surface. Our algorithm will also deal properly with the problems of reverse engineering and produce...

  8. Nonlinear approximation with dictionaries. II. Inverse Estimates

    DEFF Research Database (Denmark)

    Gribonval, Rémi; Nielsen, Morten

    2006-01-01

    In this paper, which is the sequel to [16], we study inverse estimates of the Bernstein type for nonlinear approximation with structured redundant dictionaries in a Banach space. The main results are for blockwise incoherent dictionaries in Hilbert spaces, which generalize the notion of joint block-diagonal...

  9. Nonlinear approximation with dictionaries,.. II: Inverse estimates

    DEFF Research Database (Denmark)

    Gribonval, Rémi; Nielsen, Morten

    In this paper we study inverse estimates of the Bernstein type for nonlinear approximation with structured redundant dictionaries in a Banach space. The main results are for separated decomposable dictionaries in Hilbert spaces, which generalize the notion of joint block-diagonal mutually...

  10. Counting independent sets using the Bethe approximation

    Energy Technology Data Exchange (ETDEWEB)

    Chertkov, Michael [Los Alamos National Laboratory; Chandrasekaran, V [MIT; Gamarmik, D [MIT; Shah, D [MIT; Sin, J [MIT

    2009-01-01

    The authors consider the problem of counting the number of independent sets or the partition function of a hard-core model in a graph. The problem in general is computationally hard (P hard). They study the quality of the approximation provided by the Bethe free energy. Belief propagation (BP) is a message-passing algorithm can be used to compute fixed points of the Bethe approximation; however, BP is not always guarantee to converge. As the first result, they propose a simple message-passing algorithm that converges to a BP fixed pont for any grapy. They find that their algorithm converges within a multiplicative error 1 + {var_epsilon} of a fixed point in {Omicron}(n{sup 2}E{sup -4} log{sup 3}(nE{sup -1})) iterations for any bounded degree graph of n nodes. In a nutshell, the algorithm can be thought of as a modification of BP with 'time-varying' message-passing. Next, they analyze the resulting error to the number of independent sets provided by such a fixed point of the Bethe approximation. Using the recently developed loop calculus approach by Vhertkov and Chernyak, they establish that for any bounded graph with large enough girth, the error is {Omicron}(n{sup -{gamma}}) for some {gamma} > 0. As an application, they find that for random 3-regular graph, Bethe approximation of log-partition function (log of the number of independent sets) is within o(1) of corret log-partition - this is quite surprising as previous physics-based predictions were expecting an error of o(n). In sum, their results provide a systematic way to find Bethe fixed points for any graph quickly and allow for estimating error in Bethe approximation using novel combinatorial techniques.

  11. Hybrid diffusion approximation in highly absorbing media and its effects of source approximation

    Institute of Scientific and Technical Information of China (English)

    Huijuan Tian; Ying Liu; Lijun Wang; Yuhui Zhang; Lifeng Xiao

    2009-01-01

    A modified diffusion approximation model called the hybrid diffusion approximation that can be used for highly absorbing media is investigated.The analytic solution of the hybrid diffusion approximation for reflectance in two-source approximation and steady-state case with extrapolated boundary is obtained.The effects of source approximation on the analytic solution are investigated,and it is validated that two-source approximation in highly absorbing media to describe the optical properties of biological tissue is necessary.Monte Carlo simulation of recovering optical parameters from reflectant data is done with the use of this model.The errors of recovering μa and μ's are smaller than 15% for the reduced albedo between 0.77 and 0.5 with the source-detector separation of 0.4-3 ram.

  12. Approximate dynamic programming using fluid and diffusion approximations with applications to power management

    OpenAIRE

    Chen, Wei; Huang, Dayu; Kulkarni, Ankur A.; Unnikrishnan, Jayakrishnan; Zhu, Quanyan; Mehta, Prashant; Meyn, Sean; Wierman, Adam

    2013-01-01

    Neuro-dynamic programming is a class of powerful techniques for approximating the solution to dynamic programming equations. In their most computationally attractive formulations, these techniques provide the approximate solution only within a prescribed finite-dimensional function class. Thus, the question that always arises is how should the function class be chosen? The goal of this paper is to propose an approach using the solutions to associated fluid and diffusion approximations. In ord...

  13. Parallel computations and complex analytic approximations: From diophantine approximations to quantum mechanics

    Energy Technology Data Exchange (ETDEWEB)

    Chudnovsky, D.V.; Chudnovsky, G.V. [Columbia Univ., New York, NY (United States)

    1995-12-01

    High precision solution of extremal and (complex analytic) approximations problems that can be represented in terms of multiple integrals or integral equations involving hypergeornetric functions are examined. Fast algorithms of computations of (approximate) solutions are presented that are well suited for parallelization. Among problems considered are: WKB and adelic asymptotics of multidimensional hypergeometric Pade approximations to classical functions, and high accuracy computations of high order eigenvalues and eigenstates for 2D and 3D domains of complex geometry.

  14. Polynomial approximation and cubature at approximate Fekete and Leja points of the cylinder

    CERN Document Server

    De Marchi, Stefano

    2011-01-01

    The paper deals with polynomial interpolation, least-square approximation and cubature of functions defined on the rectangular cylinder, $K=D\\times [-1,1]$, with $D$ the unit disk. The nodes used for these processes are the {\\it Approximate Fekete Points} (AFP) and the {\\it Discrete Leja Points} (DLP) extracted from suitable {\\it Weakly Admissible Meshes} (WAMs) of the cylinder. From the analysis of the growth of the Lebesgue constants, approximation and cubature errors, we show that the AFP and the DLP extracted from WAM are good points for polynomial approximation and numerical integration of functions defined on the cylinder.

  15. Approximate inverse preconditioners for general sparse matrices

    Energy Technology Data Exchange (ETDEWEB)

    Chow, E.; Saad, Y. [Univ. of Minnesota, Minneapolis, MN (United States)

    1994-12-31

    Preconditioned Krylov subspace methods are often very efficient in solving sparse linear matrices that arise from the discretization of elliptic partial differential equations. However, for general sparse indifinite matrices, the usual ILU preconditioners fail, often because of the fact that the resulting factors L and U give rise to unstable forward and backward sweeps. In such cases, alternative preconditioners based on approximate inverses may be attractive. We are currently developing a number of such preconditioners based on iterating on each column to get the approximate inverse. For this approach to be efficient, the iteration must be done in sparse mode, i.e., we must use sparse-matrix by sparse-vector type operatoins. We will discuss a few options and compare their performance on standard problems from the Harwell-Boeing collection.

  16. Flow past a porous approximate spherical shell

    Science.gov (United States)

    Srinivasacharya, D.

    2007-07-01

    In this paper, the creeping flow of an incompressible viscous liquid past a porous approximate spherical shell is considered. The flow in the free fluid region outside the shell and in the cavity region of the shell is governed by the Navier Stokes equation. The flow within the porous annulus region of the shell is governed by Darcy’s Law. The boundary conditions used at the interface are continuity of the normal velocity, continuity of the pressure and Beavers and Joseph slip condition. An exact solution for the problem is obtained. An expression for the drag on the porous approximate spherical shell is obtained. The drag experienced by the shell is evaluated numerically for several values of the parameters governing the flow.

  17. Traveltime approximations for inhomogeneous HTI media

    KAUST Repository

    Alkhalifah, Tariq Ali

    2011-01-01

    Traveltimes information is convenient for parameter estimation especially if the medium is described by an anisotropic set of parameters. This is especially true if we could relate traveltimes analytically to these medium parameters, which is generally hard to do in inhomogeneous media. As a result, I develop traveltimes approximations for horizontaly transversely isotropic (HTI) media as simplified and even linear functions of the anisotropic parameters. This is accomplished by perturbing the solution of the HTI eikonal equation with respect to η and the azimuthal symmetry direction (usually used to describe the fracture direction) from a generally inhomogeneous elliptically anisotropic background medium. The resulting approximations can provide accurate analytical description of the traveltime in a homogenous background compared to other published moveout equations out there. These equations will allow us to readily extend the inhomogenous background elliptical anisotropic model to an HTI with a variable, but smoothly varying, η and horizontal symmetry direction values. © 2011 Society of Exploration Geophysicists.

  18. Adaptive approximate Bayesian computation for complex models

    CERN Document Server

    Lenormand, Maxime; Deffuant, Guillaume

    2011-01-01

    Approximate Bayesian computation (ABC) is a family of computational techniques in Bayesian statistics. These techniques allow to fit a model to data without relying on the computation of the model likelihood. They instead require to simulate a large number of times the model to be fitted. A number of refinements to the original rejection-based ABC scheme have been proposed, including the sequential improvement of posterior distributions. This technique allows to decrease the number of model simulations required, but it still presents several shortcomings which are particularly problematic for costly to simulate complex models. We here provide a new algorithm to perform adaptive approximate Bayesian computation, which is shown to perform better on both a toy example and a complex social model.

  19. Seismic modeling using the frozen Gaussian approximation

    CERN Document Server

    Yang, Xu; Fomel, Sergey

    2013-01-01

    We adopt the frozen Gaussian approximation (FGA) for modeling seismic waves. The method belongs to the category of ray-based beam methods. It decomposes seismic wavefield into a set of Gaussian functions and propagates these Gaussian functions along appropriate ray paths. As opposed to the classic Gaussian-beam method, FGA keeps the Gaussians frozen (at a fixed width) during the propagation process and adjusts their amplitudes to produce an accurate approximation after summation. We perform the initial decomposition of seismic data using a fast version of the Fourier-Bros-Iagolnitzer (FBI) transform and propagate the frozen Gaussian beams numerically using ray tracing. A test using a smoothed Marmousi model confirms the validity of FGA for accurate modeling of seismic wavefields.

  20. An Origami Approximation to the Cosmic Web

    CERN Document Server

    Neyrinck, Mark C

    2014-01-01

    The powerful Lagrangian view of structure formation was essentially introduced to cosmology by Zel'dovich. In the current cosmological paradigm, a dark-matter-sheet 3D manifold, inhabiting 6D position-velocity phase space, was flat (with vanishing velocity) at the big bang. Afterward, gravity stretched and bunched the sheet together in different places, forming a cosmic web when projected to the position coordinates. Here, I explain some properties of an origami approximation, in which the sheet does not stretch or contract (an assumption that is false in general), but is allowed to fold. Even without stretching, the sheet can form an idealized cosmic web, with convex polyhedral voids separated by straight walls and filaments, joined by convex polyhedral nodes. The nodes form in 'polygonal' or 'polyhedral' collapse, somewhat like spherical/ellipsoidal collapse, except incorporating simultaneous filament and wall formation. The origami approximation allows phase-space geometries of nodes, filaments, and walls ...

  1. Improved Approximations for Some Polymer Extension Models

    CERN Document Server

    Petrosyan, Rafayel

    2016-01-01

    We propose approximations for force-extension dependencies for the freely jointed chain (FJC) and worm-like chain (WLC) models as well as for extension-force dependence for the WLC model. Proposed expressions show less than 1% relative error in the useful range of the corresponding variables. These results can be applied for fitting force-extension curves obtained in molecular force spectroscopy experiments. Particularly they can be useful for cases where one has geometries of springs in series and/or in parallel where particular combination of expressions should be used for fitting the data. All approximations have been obtained following the same procedure of determining the asymptotes and then reducing the relative error of that expression by adding an appropriate term obtained from fitting its absolute error.

  2. On spline approximation of sliced inverse regression

    Institute of Scientific and Technical Information of China (English)

    Li-ping ZHU; Zhou YU

    2007-01-01

    The dimension reduction is helpful and often necessary in exploring the nonparametric regression structure. In this area, Sliced inverse regression (SIR) is a promising tool to estimate the central dimension reduction (CDR) space. To estimate the kernel matrix of the SIR, we herein suggest the spline approximation using the least squares regression. The heteroscedasticity can be incorporated well by introducing an appropriate weight function. The root-n asymptotic normality can be achieved for a wide range choice of knots. This is essentially analogous to the kernel estimation. Moreover,we also propose a modified Bayes information criterion (BIC) based on the eigenvalues of the SIR matrix. This modified BIC can be applied to any form of the SIR and other related methods. The methodology and some of the practical issues are illustrated through the horse mussel data. Empirical studies evidence the performance of our proposed spline approximation by comparison of the existing estimators.

  3. Approximated solutions to Born-Infeld dynamics

    Science.gov (United States)

    Ferraro, Rafael; Nigro, Mauro

    2016-02-01

    The Born-Infeld equation in the plane is usefully captured in complex language. The general exact solution can be written as a combination of holomorphic and anti-holomorphic functions. However, this solution only expresses the potential in an implicit way. We rework the formulation to obtain the complex potential in an explicit way, by means of a perturbative procedure. We take care of the secular behavior common to this kind of approach, by resorting to a symmetry the equation has at the considered order of approximation. We apply the method to build approximated solutions to Born-Infeld electrodynamics. We solve for BI electromagnetic waves traveling in opposite directions. We study the propagation at interfaces, with the aim of searching for effects susceptible to experimental detection. In particular, we show that a reflected wave is produced when a wave is incident on a semi-space containing a magnetostatic field.

  4. Approximate Lesion Localization in Dermoscopy Images

    CERN Document Server

    Celebi, M Emre; Schaefer, Gerald; Stoecker, William V; 10.1111/j.1600-0846.2009.00357.x

    2010-01-01

    Background: Dermoscopy is one of the major imaging modalities used in the diagnosis of melanoma and other pigmented skin lesions. Due to the difficulty and subjectivity of human interpretation, automated analysis of dermoscopy images has become an important research area. Border detection is often the first step in this analysis. Methods: In this article, we present an approximate lesion localization method that serves as a preprocessing step for detecting borders in dermoscopy images. In this method, first the black frame around the image is removed using an iterative algorithm. The approximate location of the lesion is then determined using an ensemble of thresholding algorithms. Results: The method is tested on a set of 428 dermoscopy images. The localization error is quantified by a metric that uses dermatologist determined borders as the ground truth. Conclusion: The results demonstrate that the method presented here achieves both fast and accurate localization of lesions in dermoscopy images.

  5. Approximate locality for quantum systems on graphs.

    Science.gov (United States)

    Osborne, Tobias J

    2008-10-01

    In this Letter we make progress on a long-standing open problem of Aaronson and Ambainis [Theory Comput. 1, 47 (2005)]: we show that if U is a sparse unitary operator with a gap Delta in its spectrum, then there exists an approximate logarithm H of U which is also sparse. The sparsity pattern of H gets more dense as 1/Delta increases. This result can be interpreted as a way to convert between local continuous-time and local discrete-time quantum processes. As an example we show that the discrete-time coined quantum walk can be realized stroboscopically from an approximately local continuous-time quantum walk. PMID:18851512

  6. Quasi-chemical approximation for polyatomic mixtures

    CERN Document Server

    Dávila, M V; Matoz-Fernandez, D A; Ramirez-Pastor, A J

    2016-01-01

    The statistical thermodynamics of binary mixtures of polyatomic species was developed on a generalization in the spirit of the lattice-gas model and the quasi-chemical approximation (QCA). The new theoretical framework is obtained by combining: (i) the exact analytical expression for the partition function of non-interacting mixtures of linear $k$-mers and $l$-mers (species occupying $k$ sites and $l$ sites, respectively) adsorbed in one dimension, and its extension to higher dimensions; and (ii) a generalization of the classical QCA for multicomponent adsorbates and multisite-occupancy adsorption. The process is analyzed through the partial adsorption isotherms corresponding to both species of the mixture. Comparisons with analytical data from Bragg-Williams approximation (BWA) and Monte Carlo simulations are performed in order to test the validity of the theoretical model. Even though a good fitting is obtained from BWA, it is found that QCA provides a more accurate description of the phenomenon of adsorpti...

  7. Analytical Ballistic Trajectories with Approximately Linear Drag

    Directory of Open Access Journals (Sweden)

    Giliam J. P. de Carpentier

    2014-01-01

    Full Text Available This paper introduces a practical analytical approximation of projectile trajectories in 2D and 3D roughly based on a linear drag model and explores a variety of different planning algorithms for these trajectories. Although the trajectories are only approximate, they still capture many of the characteristics of a real projectile in free fall under the influence of an invariant wind, gravitational pull, and terminal velocity, while the required math for these trajectories and planners is still simple enough to efficiently run on almost all modern hardware devices. Together, these properties make the proposed approach particularly useful for real-time applications where accuracy and performance need to be carefully balanced, such as in computer games.

  8. Numerical and approximate solutions for plume rise

    Science.gov (United States)

    Krishnamurthy, Ramesh; Gordon Hall, J.

    Numerical and approximate analytical solutions are compared for turbulent plume rise in a crosswind. The numerical solutions were calculated using the plume rise model of Hoult, Fay and Forney (1969, J. Air Pollut. Control Ass.19, 585-590), over a wide range of pertinent parameters. Some wind shear and elevated inversion effects are included. The numerical solutions are seen to agree with the approximate solutions over a fairly wide range of the parameters. For the conditions considered in the study, wind shear effects are seen to be quite small. A limited study was made of the penetration of elevated inversions by plumes. The results indicate the adequacy of a simple criterion proposed by Briggs (1969, AEC Critical Review Series, USAEC Division of Technical Information extension, Oak Ridge, Tennesse).

  9. Some approximation concepts for structural synthesis

    Science.gov (United States)

    Schmit, L. A., Jr.; Farshi, B.

    1974-01-01

    An efficient automated minimum weight design procedure is presented which is applicable to sizing structural systems that can be idealized by truss, shear panel, and constant strain triangles. Static stress and displacement constraints under alternative loading conditions are considered. The optimization algorithm is an adaptation of the method of inscribed hyperspheres and high efficiency is achieved by using several approximation concepts including temporary deletion of noncritical constraints, design variable linking, and Taylor series expansions for response variables in terms of design variables. Optimum designs for several planar and space truss examples problems are presented. The results reported support the contention that the innovative use of approximation concepts in structural synthesis can produce significant improvements in efficiency.

  10. Approximate Solutions in Planted 3-SAT

    Science.gov (United States)

    Hsu, Benjamin; Laumann, Christopher; Moessner, Roderich; Sondhi, Shivaji

    2013-03-01

    In many computational settings, there exists many instances where finding a solution requires a computing time that grows exponentially in the number of variables. Concrete examples occur in combinatorial optimization problems and cryptography in computer science or glassy systems in physics. However, while exact solutions are often known to require exponential time, a related and important question is the running time required to find approximate solutions. Treating this problem as a problem in statistical physics at finite temperature, we examine the computational running time in finding approximate solutions in 3-satisfiability for randomly generated 3-SAT instances which are guaranteed to have a solution. Analytic predictions are corroborated by numerical evidence using stochastic local search algorithms. A first order transition is found in the running time of these algorithms.

  11. Rough Sets in Approximate Solution Space

    Institute of Scientific and Technical Information of China (English)

    Hui Sun; Wei Tian; Qing Liu

    2006-01-01

    As a new mathematical theory, Rough sets have been applied to processing imprecise, uncertain and in complete data. It has been fruitful in finite and non-empty set. Rough sets, however, are only served as the theoretic tool to discretize the real function. As far as the real function research is concerned, the research to define rough sets in the real function is infrequent. In this paper, we exploit a new method to extend the rough set in normed linear space, in which we establish a rough set,put forward an upper and lower approximation definition, and make a preliminary research on the property of the rough set. A new tool is provided to study the approximation solutions of differential equation and functional variation in normed linear space. This research is significant in that it extends the application of rough sets to a new field.

  12. Pade approximants of random Stieltjes series

    CERN Document Server

    Marklof, Jens; Wolowski, Lech

    2007-01-01

    We consider the random continued fraction S(t) := 1/(s_1 + t/(s_2 + t/(s_3 + >...))) where the s_n are independent random variables with the same gamma distribution. For every realisation of the sequence, S(t) defines a Stieltjes function. We study the convergence of the finite truncations of the continued fraction or, equivalently, of the diagonal Pade approximants of the function S(t). By using the Dyson--Schmidt method for an equivalent one-dimensional disordered system, and the results of Marklof et al. (2005), we obtain explicit formulae (in terms of modified Bessel functions) for the almost-sure rate of convergence of these approximants, and for the almost-sure distribution of their poles.

  13. SOME CONVERSE RESULTS ON ONESIDED APPROXIMATION: JUSTIFICATIONS

    Institute of Scientific and Technical Information of China (English)

    Wang Jianli; Zhou Songping

    2003-01-01

    The present paper deals with best onesided approximation rate in Lp spaces ~En (f)Lp of f ∈ C2π. Although it is clear that the estimate ~En(f)Lp≤C ‖f‖ Lp cannot be correct for all f ∈ Lp2π in case p<∞, the question whether ~En (f)Lp ≤Cω (f, n-1 )Lp or ~En(f)Lp ≤CEn(f)Lp holds for f ∈ C2π remains totally untouched.Therefore it forms a basic problem to justify onesided approximation. The present paper will provide an answer to settle down the basis.

  14. Approximation of Marginal Abatement Cost Curve

    OpenAIRE

    Olga Kiuila; Rutherford, Thomas F.

    2011-01-01

    Top-down models usually include piecewise-smooth functions to describe marginal cost curves, while bottom-up models describe those curves with a step function. When a bottom-up cost curve is available, we can explicitly represent this curve with a top-down model in order to replicate its shape instead of arbitrary assumptions. We propose methods to approximate a piecewise function from a step function using constant elasticity of substitution technologies. Specifically, we consider a pollutio...

  15. Approximate Bayesian inference for complex ecosystems

    OpenAIRE

    Michael P H Stumpf

    2014-01-01

    Mathematical models have been central to ecology for nearly a century. Simple models of population dynamics have allowed us to understand fundamental aspects underlying the dynamics and stability of ecological systems. What has remained a challenge, however, is to meaningfully interpret experimental or observational data in light of mathematical models. Here, we review recent developments, notably in the growing field of approximate Bayesian computation (ABC), that allow us to calibrate mathe...

  16. Relativistic point interactions: Approximation by smooth potentials

    Science.gov (United States)

    Hughes, Rhonda J.

    1997-06-01

    We show that the four-parameter family of one-dimensional relativistic point interactions studied by Benvegnu and Dąbrowski may be approximated in the strong resolvent sense by smooth, local, short-range perturbations of the Dirac Hamiltonian. In addition, we prove that the nonrelativistic limits correspond to the Schrödinger point interactions studied extensively by the author and Paul Chernoff.

  17. Dynamic Approximate Vertex Cover and Maximum Matching

    OpenAIRE

    Onak, Krzysztof; Rubinfeld, Ronitt

    2010-01-01

    We consider the problem of maintaining a large matching or a small vertex cover in a dynamically changing graph. Each update to the graph is either an edge deletion or an edge insertion. We give the first randomized data structure that simultaneously achieves a constant approximation factor and handles a sequence of k updates in k. polylog(n) time. Previous data structures require a polynomial amount of computation per update. The starting point of our construction is a distributed algorit...

  18. Approximating viability kernels with support vector machines

    OpenAIRE

    Deffuant, G.; Chapel, L.; Martin, S.

    2007-01-01

    We propose an algorithm which performs a progressive approximation of a viability kernel, iteratively using a classification method. We establish the mathematical conditions that the classification method should fulfill to guarantee the convergence to the actual viability kernel. We study more particularly the use of support vector machines (SVMs) as classification techniques. We show that they make possible to use gradient optimisation techniques to find a viable control at each time step, a...

  19. Solving Math Problems Approximately: A Developmental Perspective.

    Directory of Open Access Journals (Sweden)

    Dana Ganor-Stern

    Full Text Available Although solving arithmetic problems approximately is an important skill in everyday life, little is known about the development of this skill. Past research has shown that when children are asked to solve multi-digit multiplication problems approximately, they provide estimates that are often very far from the exact answer. This is unfortunate as computation estimation is needed in many circumstances in daily life. The present study examined 4th graders, 6th graders and adults' ability to estimate the results of arithmetic problems relative to a reference number. A developmental pattern was observed in accuracy, speed and strategy use. With age there was a general increase in speed, and an increase in accuracy mainly for trials in which the reference number was close to the exact answer. The children tended to use the sense of magnitude strategy, which does not involve any calculation but relies mainly on an intuitive coarse sense of magnitude, while the adults used the approximated calculation strategy which involves rounding and multiplication procedures, and relies to a greater extent on calculation skills and working memory resources. Importantly, the children were less accurate than the adults, but were well above chance level. In all age groups performance was enhanced when the reference number was smaller (vs. larger than the exact answer and when it was far (vs. close from it, suggesting the involvement of an approximate number system. The results suggest the existence of an intuitive sense of magnitude for the results of arithmetic problems that might help children and even adults with difficulties in math. The present findings are discussed in the context of past research reporting poor estimation skills among children, and the conditions that might allow using children estimation skills in an effective manner.

  20. Gaussian Approximation Potentials: a brief tutorial introduction

    OpenAIRE

    Bart?k, Albert P.; Cs?nyi, G?bor

    2015-01-01

    We present a swift walk-through of our recent work that uses machine learning to t interatomic potentials based on quantum mechanical data. We describe our Gaussian Approximation Potentials (GAP) framework, discuss a variety of descriptors, how to train the model on total energies and derivatives and the simultaneous use of multiple models of di erent complexity. We also show a small example using QUIP, the software sandbox implementation of GAP that is available for non-comme...

  1. Dual Control for Approximate Bayesian Reinforcement Learning

    OpenAIRE

    Klenske, Edgar D.; Hennig, Philipp

    2015-01-01

    Control of non-episodic, finite-horizon dynamical systems with uncertain dynamics poses a tough and elementary case of the exploration-exploitation trade-off. Bayesian reinforcement learning, reasoning about the effect of actions and future observations, offers a principled solution, but is intractable. We review, then extend an old approximate approach from control theory---where the problem is known as dual control---in the context of modern regression methods, specifically generalized line...

  2. Space-Time Approximation with Sparse Grids

    Energy Technology Data Exchange (ETDEWEB)

    Griebel, M; Oeltz, D; Vassilevski, P S

    2005-04-14

    In this article we introduce approximation spaces for parabolic problems which are based on the tensor product construction of a multiscale basis in space and a multiscale basis in time. Proper truncation then leads to so-called space-time sparse grid spaces. For a uniform discretization of the spatial space of dimension d with O(N{sup d}) degrees of freedom, these spaces involve for d > 1 also only O(N{sup d}) degrees of freedom for the discretization of the whole space-time problem. But they provide the same approximation rate as classical space-time Finite Element spaces which need O(N{sup d+1}) degrees of freedoms. This makes these approximation spaces well suited for conventional parabolic and for time-dependent optimization problems. We analyze the approximation properties and the dimension of these sparse grid space-time spaces for general stable multiscale bases. We then restrict ourselves to an interpolatory multiscale basis, i.e. a hierarchical basis. Here, to be able to handle also complicated spatial domains {Omega}, we construct the hierarchical basis from a given spatial Finite Element basis as follows: First we determine coarse grid points recursively over the levels by the coarsening step of the algebraic multigrid method. Then, we derive interpolatory prolongation operators between the respective coarse and fine grid points by a least squares approach. This way we obtain an algebraic hierarchical basis for the spatial domain which we then use in our space-time sparse grid approach. We give numerical results on the convergence rate of the interpolation error of these spaces for various space-time problems with two spatial dimensions. Also implementational issues, data structures and questions of adaptivity are addressed to some extent.

  3. APPROXIMATION MULTIDIMENSION FUCTION WITH FUNCTIONAL NETWORK

    Institute of Scientific and Technical Information of China (English)

    Li Weibin; Liu Fang; Jiao Licheng; Zhang Shuling; Li Zongling

    2006-01-01

    The functional network was introduced by E.Catillo, which extended the neural network. Not only can it solve the problems solved, but also it can formulate the ones that cannot be solved by traditional network.This paper applies functional network to approximate the multidimension function under the ridgelet theory.The method performs more stable and faster than the traditional neural network. The numerical examples demonstrate the performance.

  4. Microscopic justification of the Equal Filling approximation

    CERN Document Server

    Perez-Martin, Sara

    2008-01-01

    The Equal Filling Approximation, a procedure widely used in mean field calculations to treat the dynamics of odd nuclei in a time reversal invariant way, is justified as the consequence of a variational principle over an average energy functional. The ideas of Statistical Quantum Mechanics are employed in the justification. As an illustration of the method, the ground and lowest lying states of some octupole deformed Radium isotopes are computed.

  5. On the diagonal approximation of full matrices

    OpenAIRE

    Lioen, W.M.

    1996-01-01

    In this paper the construction of diagonal matrices, in some sense approximating the inverse of a given square matrix, is described. The matrices are constructed using the well-known computer algebra system Maple. The techniques we show are applicable to square matrices in general. Results are given for use in Parallel diagonal-implicit Runge-Kutta (PDIRK) methods. For an s-stage Radau IIA corrector we conjecture $s!$ possibilities for the diagonal matrices.

  6. Fast approximate convex decomposition using relative concavity

    KAUST Repository

    Ghosh, Mukulika

    2013-02-01

    Approximate convex decomposition (ACD) is a technique that partitions an input object into approximately convex components. Decomposition into approximately convex pieces is both more efficient to compute than exact convex decomposition and can also generate a more manageable number of components. It can be used as a basis of divide-and-conquer algorithms for applications such as collision detection, skeleton extraction and mesh generation. In this paper, we propose a new method called Fast Approximate Convex Decomposition (FACD) that improves the quality of the decomposition and reduces the cost of computing it for both 2D and 3D models. In particular, we propose a new strategy for evaluating potential cuts that aims to reduce the relative concavity, rather than absolute concavity. As shown in our results, this leads to more natural and smaller decompositions that include components for small but important features such as toes or fingers while not decomposing larger components, such as the torso, that may have concavities due to surface texture. Second, instead of decomposing a component into two pieces at each step, as in the original ACD, we propose a new strategy that uses a dynamic programming approach to select a set of n c non-crossing (independent) cuts that can be simultaneously applied to decompose the component into n c+1 components. This reduces the depth of recursion and, together with a more efficient method for computing the concavity measure, leads to significant gains in efficiency. We provide comparative results for 2D and 3D models illustrating the improvements obtained by FACD over ACD and we compare with the segmentation methods in the Princeton Shape Benchmark by Chen et al. (2009) [31]. © 2012 Elsevier Ltd. All rights reserved.

  7. Broadband Approximations for Doubly Curved Reflector Antenna

    OpenAIRE

    V. Schejbal; J. Pidanic

    2010-01-01

    The broadband approximations for shaped-beam doubly curved reflector antennas with primary feed (rectangular horn) producing uniform amplitude and phase aperture distribution are derived and analyzed. They are very valuable for electromagnetic compatibility analyses both from electromagnetic interference and susceptibility point of view, because specialized more accurate methods such as physical optics are only used by antenna designers. To allow quick EMC analyses, typical values, beamwidth ...

  8. Neutrino Mass Matrix with Approximate Flavor Symmetry

    CERN Document Server

    Riazuddin, M

    2003-01-01

    Phenomenological implications of neutrino oscillations implied by recent experimental data on pattern of neutrino mass matrix are disscussed. It is shown that it is possible to have a neutrino mass matrix which shows approximate flavor symmetry; the neutrino mass differences arise from flavor violation in off-diagonal Yukawa couplings. Two modest extensions of the standard model, which can embed the resulting neutrino mass matix have also been discussed.

  9. Discrete least squares approximation with polynomial vectors

    OpenAIRE

    Van Barel, Marc; Bultheel, Adhemar

    1993-01-01

    We give a solution of a discrete least squares approximation problem in terms of orthogonal polynomial vectors. The degrees of the polynomial elements of these vectors can be different. An algorithm is constructed computing the coefficients of recurrence relations for the orthogonal polynomial vectors. In case the function values are prescribed in points on the real line or on the unit circle variants of the original algorithm can be designed which are an order of magnitude more efficient. Al...

  10. Relativistic impulse approximation for nuclear inelastic scattering

    International Nuclear Information System (INIS)

    The Relativistic Impulse Approximation (RIA) for proton-nucleus elastic and inelastic scattering is contrasted with its non-relativistic counterpart (the NRIA). Differences between the two approaches are examined with special emphasis on the nuclear convection current and its generalizations which may show signatures of strong relativistic nuclear potentials. A simple extension of the RIA to meson-nucleus scattering based on the linear, spin-zero Duffin-Kemmer wave equation is considered

  11. The Complexity of Approximately Counting Stable Matchings

    CERN Document Server

    Chebolu, Prasad; Martin, Russell

    2010-01-01

    We investigate the complexity of approximately counting stable matchings in the $k$-attribute model, where the preference lists are determined by dot products of "preference vectors" with "attribute vectors", or by Euclidean distances between "preference points" and "attribute points". Irving and Leather proved that counting the number of stable matchings in the general case is $#P$-complete. Counting the number of stable matchings is reducible to counting the number of downsets in a (related) partial order and is interreducible, in an approximation-preserving sense, to a class of problems that includes counting the number of independent sets in a bipartite graph ($#BIS$). It is conjectured that no FPRAS exists for this class of problems. We show this approximation-preserving interreducibilty remains even in the restricted $k$-attribute setting when $k \\geq 3$ (dot products) or $k \\geq 2$ (Euclidean distances). Finally, we show it is easy to count the number of stable matchings in the 1-attribute dot-product ...

  12. Development of the relativistic impulse approximation

    International Nuclear Information System (INIS)

    This talk contains three parts. Part I reviews the developments which led to the relativistic impulse approximation for proton-nucleus scattering. In Part II, problems with the impulse approximation in its original form - principally the low energy problem - are discussed and traced to pionic contributions. Use of pseudovector covariants in place of pseudoscalar ones in the NN amplitude provides more satisfactory low energy results, however, the difference between pseudovector and pseudoscalar results is ambiguous in the sense that it is not controlled by NN data. Only with further theoretical input can the ambiguity be removed. Part III of the talk presents a new development of the relativistic impulse approximation which is the result of work done in the past year and a half in collaboration with J.A. Tjon. A complete NN amplitude representation is developed and a complete set of Lorentz invariant amplitudes are calculated based on a one-meson exchange model and appropriate integral equations. A meson theoretical basis for the important pair contributions to proton-nucleus scattering is established by the new developments. 28 references

  13. Simultaneous perturbation stochastic approximation for tidal models

    KAUST Repository

    Altaf, M.U.

    2011-05-12

    The Dutch continental shelf model (DCSM) is a shallow sea model of entire continental shelf which is used operationally in the Netherlands to forecast the storm surges in the North Sea. The forecasts are necessary to support the decision of the timely closure of the moveable storm surge barriers to protect the land. In this study, an automated model calibration method, simultaneous perturbation stochastic approximation (SPSA) is implemented for tidal calibration of the DCSM. The method uses objective function evaluations to obtain the gradient approximations. The gradient approximation for the central difference method uses only two objective function evaluation independent of the number of parameters being optimized. The calibration parameter in this study is the model bathymetry. A number of calibration experiments is performed. The effectiveness of the algorithm is evaluated in terms of the accuracy of the final results as well as the computational costs required to produce these results. In doing so, comparison is made with a traditional steepest descent method and also with a newly developed proper orthogonal decompositionbased calibration method. The main findings are: (1) The SPSA method gives comparable results to steepest descent method with little computational cost. (2) The SPSA method with little computational cost can be used to estimate large number of parameters.

  14. CMB-lensing beyond the Born approximation

    Science.gov (United States)

    Marozzi, Giovanni; Fanizza, Giuseppe; Di Dio, Enea; Durrer, Ruth

    2016-09-01

    We investigate the weak lensing corrections to the cosmic microwave background temperature anisotropies considering effects beyond the Born approximation. To this aim, we use the small deflection angle approximation, to connect the lensed and unlensed power spectra, via expressions for the deflection angles up to third order in the gravitational potential. While the small deflection angle approximation has the drawback to be reliable only for multipoles l lesssim 2500, it allows us to consistently take into account the non-Gaussian nature of cosmological perturbation theory beyond the linear level. The contribution to the lensed temperature power spectrum coming from the non-Gaussian nature of the deflection angle at higher order is a new effect which has not been taken into account in the literature so far. It turns out to be the leading contribution among the post-Born lensing corrections. On the other hand, the effect is smaller than corrections coming from non-linearities in the matter power spectrum, and its imprint on CMB lensing is too small to be seen in present experiments.

  15. Cylindrical Helix Spline Approximation of Spatial Curves

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    In this paper, we present a new method for approximating spatial curves with a G1 cylindrical helix spline within a prescribed tolerance. We deduce the general formulation of a cylindrical helix,which has 11 freedoms. This means that it needs 11 restrictions to determine a cylindrical helix. Given a spatial parametric curve segment, including the start point and the end point of this segment, the tangent and the principal normal of the start point, we can always find a cylindrical segment to interpolate the given direction and position vectors. In order to approximate the known parametric curve within the prescribed tolerance, we adopt the trial method step by step. First, we must ensure the helix segment to interpolate the given two end points and match the principal normal and tangent of the start point, and then, we can keep the deviation between the cylindrical helix segment and the known curve segment within the prescribed tolerance everywhere. After the first segment had been formed, we can construct the next segment. Circularly, we can construct the G1 cylindrical helix spline to approximate the whole spatial parametric curve within the prescribed tolerance. Several examples are also given to show the efficiency of this method.

  16. Green-Ampt approximations: A comprehensive analysis

    Science.gov (United States)

    Ali, Shakir; Islam, Adlul; Mishra, P. K.; Sikka, Alok K.

    2016-04-01

    Green-Ampt (GA) model and its modifications are widely used for simulating infiltration process. Several explicit approximate solutions to the implicit GA model have been developed with varying degree of accuracy. In this study, performance of nine explicit approximations to the GA model is compared with the implicit GA model using the published data for broad range of soil classes and infiltration time. The explicit GA models considered are Li et al. (1976) (LI), Stone et al. (1994) (ST), Salvucci and Entekhabi (1994) (SE), Parlange et al. (2002) (PA), Barry et al. (2005) (BA), Swamee et al. (2012) (SW), Ali et al. (2013) (AL), Almedeij and Esen (2014) (AE), and Vatankhah (2015) (VA). Six statistical indicators (e.g., percent relative error, maximum absolute percent relative error, average absolute percent relative errors, percent bias, index of agreement, and Nash-Sutcliffe efficiency) and relative computer computation time are used for assessing the model performance. Models are ranked based on the overall performance index (OPI). The BA model is found to be the most accurate followed by the PA and VA models for variety of soil classes and infiltration periods. The AE, SW, SE, and LI model also performed comparatively better. Based on the overall performance index, the explicit models are ranked as BA > PA > VA > LI > AE > SE > SW > ST > AL. Results of this study will be helpful in selection of accurate and simple explicit approximate GA models for solving variety of hydrological problems.

  17. Generalized Quasilinear Approximation: Application to Zonal Jets

    Science.gov (United States)

    Marston, J. B.; Chini, G. P.; Tobias, S. M.

    2016-05-01

    Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. We present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through nonlocal spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta plane and show that it is accurate even for a small number of large-scale modes. As GQL is formally linear in the small zonal scales, it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems.

  18. A coastal ocean model with subgrid approximation

    Science.gov (United States)

    Walters, Roy A.

    2016-06-01

    A wide variety of coastal ocean models exist, each having attributes that reflect specific application areas. The model presented here is based on finite element methods with unstructured grids containing triangular and quadrilateral elements. The model optimizes robustness, accuracy, and efficiency by using semi-implicit methods in time in order to remove the most restrictive stability constraints, by using a semi-Lagrangian advection approximation to remove Courant number constraints, and by solving a wave equation at the discrete level for enhanced efficiency. An added feature is the approximation of the effects of subgrid objects. Here, the Reynolds-averaged Navier-Stokes equations and the incompressibility constraint are volume averaged over one or more computational cells. This procedure gives rise to new terms which must be approximated as a closure problem. A study of tidal power generation is presented as an example of this method. A problem that arises is specifying appropriate thrust and power coefficients for the volume averaged velocity when they are usually referenced to free stream velocity. A new contribution here is the evaluation of three approaches to this problem: an iteration procedure and two mapping formulations. All three sets of results for thrust (form drag) and power are in reasonable agreement.

  19. Approximation and universality of fuzzy Turing machines

    Institute of Scientific and Technical Information of China (English)

    LI YongMing

    2008-01-01

    Fuzzy Turing machines are the formal models of fuzzy algorithms or fuzzy computations.In this paper we give several different formulations of fuzzy Turing machine,which correspond to nondeterministic fuzzy Turing machine using max-★ composition for some t-norm ★ (or NFTM★,for short),nondeterministic fuzzy Turing machine (or NFTM),deterministic fuzzy Turing machine (or DFTM),and multi-tape versions of fuzzy Turing machines.Some distinct results compared to those of ordinary Turing machines are obtained.First,it is shown that NFTM★,NFTM,and DFTM are not necessarily equivalent in the power of recognizing fuzzy languages if the t-norm ★ does not satisfy finite generated condition,but are equivalent in the approximation sense.That is to say,we can approximate an NFTM★ by some NFTM with any given accuracy;the related constructions are also presented.The level characterization of fuzzy recursively enumerable languages and fuzzy recursive languages are exploited by ordinary r.e.languages and recursive languages.Second,we show that universal fuzzy Turing machine exists in the approximated sense.There is a universal fuzzy Turing machine that can simulate any NFTM★ on it with a given accuracy.

  20. CMB-lensing beyond the Born approximation

    CERN Document Server

    Marozzi, Giovanni; Di Dio, Enea; Durrer, Ruth

    2016-01-01

    We investigate the weak lensing corrections to the cosmic microwave background temperature anisotropies considering effects beyond the Born approximation. To this aim, we use the small deflection angle approximation, to connect the lensed and unlensed power spectra, via expressions for the deflection angles up to third order in the gravitational potential. While the small deflection angle approximation has the drawback to be reliable only for multipoles $\\ell\\lesssim 2500$, it allows us to consistently take into account the non-Gaussian nature of cosmological perturbation theory beyond the linear level. The contribution to the lensed temperature power spectrum coming from the non-Gaussian nature of the deflection angle at higher order is a new effect which has not been taken into account in the literature so far. It turns out to be the leading contribution among the post-Born lensing corrections. On the other hand, the effect is smaller than corrections coming from non-linearities in the matter power spectrum...

  1. Approximate Equalities on Rough Intuitionistic Fuzzy Sets and an Analysis of Approximate Equalities

    Directory of Open Access Journals (Sweden)

    B. K. Tripathy

    2012-03-01

    Full Text Available In order to involve user knowledge in determining equality of sets, which may not be equal in the mathematical sense, three types of approximate (rough equalities were introduced by Novotny and Pawlak ([8, 9, 10]. These notions were generalized by Tripathy, Mitra and Ojha ([13], who introduced the concepts of approximate (rough equivalences of sets. Rough equivalences capture equality of sets at a higher level than rough equalities. More properties of these concepts were established in [14]. Combining the conditions for the two types of approximate equalities, two more approximate equalities were introduced by Tripathy [12] and a comparative analysis of their relative efficiency was provided. In [15], the four types of approximate equalities were extended by considering rough fuzzy sets instead of only rough sets. In fact the concepts of leveled approximate equalities were introduced and properties were studied. In this paper we proceed further by introducing and studying the approximate equalities based on rough intuitionistic fuzzy sets instead of rough fuzzy sets. That is we introduce the concepts of approximate (rough equalities of intuitionistic fuzzy sets and study their properties. We provide some real life examples to show the applications of rough equalities of fuzzy sets and rough equalities of intuitionistic fuzzy sets.

  2. On the Use of Approximations in Statistical Physics

    CERN Document Server

    Hoffmann, C

    2003-01-01

    Two approximations are frequently used in statistical physics: the first one, which we shall name the mean values approximation, is generally (and improperly) named as "maximum term approximation". The second is the "Stirling approximation". In this paper we demonstrate that the error introduced by the first approximation is exactly compensated by the second approximation in the calculation of mean values of multinomial distributions.

  3. Approximate number and approximate time discrimination each correlate with school math abilities in young children.

    Science.gov (United States)

    Odic, Darko; Lisboa, Juan Valle; Eisinger, Robert; Olivera, Magdalena Gonzalez; Maiche, Alejandro; Halberda, Justin

    2016-01-01

    What is the relationship between our intuitive sense of number (e.g., when estimating how many marbles are in a jar), and our intuitive sense of other quantities, including time (e.g., when estimating how long it has been since we last ate breakfast)? Recent work in cognitive, developmental, comparative psychology, and computational neuroscience has suggested that our representations of approximate number, time, and spatial extent are fundamentally linked and constitute a "generalized magnitude system". But, the shared behavioral and neural signatures between number, time, and space may alternatively be due to similar encoding and decision-making processes, rather than due to shared domain-general representations. In this study, we investigate the relationship between approximate number and time in a large sample of 6-8 year-old children in Uruguay by examining how individual differences in the precision of number and time estimation correlate with school mathematics performance. Over four testing days, each child completed an approximate number discrimination task, an approximate time discrimination task, a digit span task, and a large battery of symbolic math tests. We replicate previous reports showing that symbolic math abilities correlate with approximate number precision and extend those findings by showing that math abilities also correlate with approximate time precision. But, contrary to approximate number and time sharing common representations, we find that each of these dimensions uniquely correlates with formal math: approximate number correlates more strongly with formal math compared to time and continues to correlate with math even when precision in time and individual differences in working memory are controlled for. These results suggest that there are important differences in the mental representations of approximate number and approximate time and further clarify the relationship between quantity representations and mathematics. PMID:26587963

  4. Approximate number and approximate time discrimination each correlate with school math abilities in young children.

    Science.gov (United States)

    Odic, Darko; Lisboa, Juan Valle; Eisinger, Robert; Olivera, Magdalena Gonzalez; Maiche, Alejandro; Halberda, Justin

    2016-01-01

    What is the relationship between our intuitive sense of number (e.g., when estimating how many marbles are in a jar), and our intuitive sense of other quantities, including time (e.g., when estimating how long it has been since we last ate breakfast)? Recent work in cognitive, developmental, comparative psychology, and computational neuroscience has suggested that our representations of approximate number, time, and spatial extent are fundamentally linked and constitute a "generalized magnitude system". But, the shared behavioral and neural signatures between number, time, and space may alternatively be due to similar encoding and decision-making processes, rather than due to shared domain-general representations. In this study, we investigate the relationship between approximate number and time in a large sample of 6-8 year-old children in Uruguay by examining how individual differences in the precision of number and time estimation correlate with school mathematics performance. Over four testing days, each child completed an approximate number discrimination task, an approximate time discrimination task, a digit span task, and a large battery of symbolic math tests. We replicate previous reports showing that symbolic math abilities correlate with approximate number precision and extend those findings by showing that math abilities also correlate with approximate time precision. But, contrary to approximate number and time sharing common representations, we find that each of these dimensions uniquely correlates with formal math: approximate number correlates more strongly with formal math compared to time and continues to correlate with math even when precision in time and individual differences in working memory are controlled for. These results suggest that there are important differences in the mental representations of approximate number and approximate time and further clarify the relationship between quantity representations and mathematics.

  5. Approximate protein structural alignment in polynomial time.

    Science.gov (United States)

    Kolodny, Rachel; Linial, Nathan

    2004-08-17

    Alignment of protein structures is a fundamental task in computational molecular biology. Good structural alignments can help detect distant evolutionary relationships that are hard or impossible to discern from protein sequences alone. Here, we study the structural alignment problem as a family of optimization problems and develop an approximate polynomial-time algorithm to solve them. For a commonly used scoring function, the algorithm runs in O(n(10)/epsilon(6)) time, for globular protein of length n, and it detects alignments that score within an additive error of epsilon from all optima. Thus, we prove that this task is computationally feasible, although the method that we introduce is too slow to be a useful everyday tool. We argue that such approximate solutions are, in fact, of greater interest than exact ones because of the noisy nature of experimentally determined protein coordinates. The measurement of similarity between a pair of protein structures used by our algorithm involves the Euclidean distance between the structures (appropriately rigidly transformed). We show that an alternative approach, which relies on internal distance matrices, must incorporate sophisticated geometric ingredients if it is to guarantee optimality and run in polynomial time. We use these observations to visualize the scoring function for several real instances of the problem. Our investigations yield insights on the computational complexity of protein alignment under various scoring functions. These insights can be used in the design of scoring functions for which the optimum can be approximated efficiently and perhaps in the development of efficient algorithms for the multiple structural alignment problem. PMID:15304646

  6. Photoelectron spectroscopy and the dipole approximation

    Energy Technology Data Exchange (ETDEWEB)

    Hemmers, O.; Hansen, D.L.; Wang, H. [Univ. of Nevada, Las Vegas, NV (United States)] [and others

    1997-04-01

    Photoelectron spectroscopy is a powerful technique because it directly probes, via the measurement of photoelectron kinetic energies, orbital and band structure in valence and core levels in a wide variety of samples. The technique becomes even more powerful when it is performed in an angle-resolved mode, where photoelectrons are distinguished not only by their kinetic energy, but by their direction of emission as well. Determining the probability of electron ejection as a function of angle probes the different quantum-mechanical channels available to a photoemission process, because it is sensitive to phase differences among the channels. As a result, angle-resolved photoemission has been used successfully for many years to provide stringent tests of the understanding of basic physical processes underlying gas-phase and solid-state interactions with radiation. One mainstay in the application of angle-resolved photoelectron spectroscopy is the well-known electric-dipole approximation for photon interactions. In this simplification, all higher-order terms, such as those due to electric-quadrupole and magnetic-dipole interactions, are neglected. As the photon energy increases, however, effects beyond the dipole approximation become important. To best determine the range of validity of the dipole approximation, photoemission measurements on a simple atomic system, neon, where extra-atomic effects cannot play a role, were performed at BL 8.0. The measurements show that deviations from {open_quotes}dipole{close_quotes} expectations in angle-resolved valence photoemission are observable for photon energies down to at least 0.25 keV, and are quite significant at energies around 1 keV. From these results, it is clear that non-dipole angular-distribution effects may need to be considered in any application of angle-resolved photoelectron spectroscopy that uses x-ray photons of energies as low as a few hundred eV.

  7. A linear approximation to black hole evaporation

    International Nuclear Information System (INIS)

    An evaporating Schwarzschild black hole is analysed including back reaction in a linear approximation. The analysis assumes a massless scalar field propagating in a spacetime consisting of two Vaidya metrics corresponding respectively to outgoing radiation and an infalling negative energy flux. For times late relative to the collapse but early relative to the lifetime of the hole, the standard rate is reproduced and has the correct time dependence. The event horizon shrinks at the expected rate. These results are independent of the exact location of the boundary between the regions. The magnitude of the quantum fluxes at various radii suggests that most of the pair production occurs far from the horizon

  8. TMB: Automatic differentiation and laplace approximation

    DEFF Research Database (Denmark)

    Kristensen, Kasper; Nielsen, Anders; Berg, Casper Willestofte;

    2016-01-01

    TMB is an open source R package that enables quick implementation of complex nonlinear random effects (latent variable) models in a manner similar to the established AD Model Builder package (ADMB, http://admb-project.org/; Fournier et al. 2011). In addition, it offers easy access to parallel...... computations. The user defines the joint likelihood for the data and the random effects as a C++ template function, while all the other operations are done in R; e.g., reading in the data. The package evaluates and maximizes the Laplace approximation of the marginal likelihood where the random effects...

  9. Casimir forces beyond the proximity approximation

    CERN Document Server

    Bimonte, G; Jaffe, R L; Kardar, M

    2011-01-01

    The proximity force approximation (PFA) relates the interaction between closely spaced, smoothly curved objects to the force between parallel plates. Precision experiments on Casimir forces necessitate, and spur research on, corrections to the PFA. We use a derivative expansion for gently curved surfaces to derive the leading curvature modifications to the PFA. Our methods apply to any homogeneous and isotropic materials; here we present results for Dirichlet and Neumann boundary conditions and for perfect conductors. A Pad\\'e extrapolation constrained by a multipole expansion at large distance and our improved expansion at short distances, provides an accurate expression for the sphere-plate Casimir force at all separations.

  10. Inaccurate approximation in the modelling of hyperinflations

    OpenAIRE

    Moffatt, Peter G.; Evens SALIES

    2006-01-01

    In time series macroeconometric models, the first difference in the logarithm of a variable is routinely used to represent the rate of change of that variable. It is often overlooked that the assumed approximation is accurate only if the rates of change are small. Models of hyper-inflation are a case in point, since in these models, by definition, changes in price are large. In this letter, Cagan's model is applied to Hungarian hyper-inflation data. It is then demonstrated that use of the app...

  11. Error Minimization of Polynomial Approximation of Delta

    Indian Academy of Sciences (India)

    Islam Sana; Sadiq Muhammad; Qureshi Muhammad Shahid

    2008-09-01

    The difference between Universal time (UT) and Dynamical time (TD), known as Delta ( ) is tabulated for the first day of each year in the Astronomical Almanac. During the last four centuries it is found that there are large differences between its values for two consecutive years. Polynomial approximations have been developed to obtain the values of for any time of a year for the period AD 1620 to AD 2000 (Meeu 2000) as no dynamical theories describe the variations in . In this work, a new set of polynomials for is obtained for the period AD 1620 to AD 2007 that is found to produce better results compared to previous attempts.

  12. Exponential Polynomial Approximation with Unrestricted Upper Density

    Institute of Scientific and Technical Information of China (English)

    Xiang Dong YANG

    2011-01-01

    We take a new approach to obtaining necessary and sufficient conditions for the incompleteness of exponential polynomials in Lp/α, where Lp/α is the weighted Banach space of complex continuous functions f defined on the real axis (R)satisfying (∫+∞/-∞|f(t)|pe-α(t)dt)1/p, 1 < p < ∞, and α(t) is a nonnegative continuous function defined on the real axis (R). In this paper, the upper density of the sequence which forms the exponential polynomials is not required to be finite. In the study of weighted polynomial approximation, consideration of the case is new.

  13. Approximate formulas for moderately small eikonal amplitudes

    CERN Document Server

    Kisselev, A V

    2015-01-01

    The eikonal approximation for moderately small scattering amplitudes is considered. With the purpose of using for their numerical estimations, the formulas are derived which contain no Bessel functions, and, hence, no rapidly oscillating integrands. To obtain these formulas, the improper integrals of the first kind which contain products of the Bessel functions J_0(z) are studied. The expression with four functions J_0(z) is generalized. The expressions for the integrals with the product of five and six Bessel functions J_0(z) are also found. The known formula for the improper integral with two functions J_nu(z) is generalized for non-integer nu.

  14. Partially coherent contrast-transfer-function approximation.

    Science.gov (United States)

    Nesterets, Yakov I; Gureyev, Timur E

    2016-04-01

    The contrast-transfer-function (CTF) approximation, widely used in various phase-contrast imaging techniques, is revisited. CTF validity conditions are extended to a wide class of strongly absorbing and refracting objects, as well as to nonuniform partially coherent incident illumination. Partially coherent free-space propagators, describing amplitude and phase in-line contrast, are introduced and their properties are investigated. The present results are relevant to the design of imaging experiments with partially coherent sources, as well as to the analysis and interpretation of the corresponding images. PMID:27140752

  15. High energy approximations in quantum field theory

    International Nuclear Information System (INIS)

    New theoretical methods in hadron physics based on a high-energy perturbation theory are discussed. The approximated solutions to quantum field theory obtained by this method appear to be sufficiently simple and rich in structure to encourage hadron dynamics studies. Operator eikonal form for field - theoretic Green's functions is derived and discussion is held on how the eikonal perturbation theory is to be renormalized. This method is extended to massive quantum electrodynamics of scalar charged bosons. Possible developments and applications of this theory are given

  16. On approximation of Markov binomial distributions

    OpenAIRE

    Xia, Aihua; Zhang, Mei

    2009-01-01

    For a Markov chain $\\mathbf{X}=\\{X_i,i=1,2,...,n\\}$ with the state space $\\{0,1\\}$, the random variable $S:=\\sum_{i=1}^nX_i$ is said to follow a Markov binomial distribution. The exact distribution of $S$, denoted $\\mathcal{L}S$, is very computationally intensive for large $n$ (see Gabriel [Biometrika 46 (1959) 454--460] and Bhat and Lal [Adv. in Appl. Probab. 20 (1988) 677--680]) and this paper concerns suitable approximate distributions for $\\mathcal{L}S$ when $\\mathbf{X}$ is stationary. We...

  17. Topics in multivariate approximation and interpolation

    CERN Document Server

    Jetter, Kurt

    2005-01-01

    This book is a collection of eleven articles, written by leading experts and dealing with special topics in Multivariate Approximation and Interpolation. The material discussed here has far-reaching applications in many areas of Applied Mathematics, such as in Computer Aided Geometric Design, in Mathematical Modelling, in Signal and Image Processing and in Machine Learning, to mention a few. The book aims at giving a comprehensive information leading the reader from the fundamental notions and results of each field to the forefront of research. It is an ideal and up-to-date introduction for gr

  18. An approximation algorithm for counting contingency tables

    OpenAIRE

    Barvinok, Alexander; Luria, Zur; Samorodnitsky, Alex; Yong, Alexander

    2008-01-01

    We present a randomized approximation algorithm for counting contingency tables, mxn non-negative integer matrices with given row sums R=(r_1, ..., r_m) and column sums C=(c_1, ..., c_n). We define smooth margins (R,C) in terms of the typical table and prove that for such margins the algorithm has quasi-polynomial N^{O(ln N)} complexity, where N=r_1+...+r_m=c_1+...+c_n. Various classes of margins are smooth, e.g., when m=O(n), n=O(m) and the ratios between the largest and the smallest row sum...

  19. Shape theory categorical methods of approximation

    CERN Document Server

    Cordier, J M

    2008-01-01

    This in-depth treatment uses shape theory as a ""case study"" to illustrate situations common to many areas of mathematics, including the use of archetypal models as a basis for systems of approximations. It offers students a unified and consolidated presentation of extensive research from category theory, shape theory, and the study of topological algebras.A short introduction to geometric shape explains specifics of the construction of the shape category and relates it to an abstract definition of shape theory. Upon returning to the geometric base, the text considers simplical complexes and

  20. A systematic approximate method for the study of evolution problem beyond rotating wave approximation

    Institute of Scientific and Technical Information of China (English)

    Zhang Zheng-Jie; Wang Ke-Lin; Qin Gan

    2005-01-01

    By a model of atwo-level particle coupled with boson field, we made it clear that an evolution problem can be solved beyond the rotating wave approximation. We have applied the coherent approximation method, which had been proved to be effective in dealing with stationary state problems of polaron, to the evolution problem of the system mentioned above. The results obtained showed that the coherent approximation method is effective to treat the evolution problem,and, in general cases, the non-rotating wave terms in Hamiltonian should not be ignored. Our results may provide a deep physical insight for further experiments to test the effects of non-rotating wave terms.

  1. Perturbation of Operators and Approximation of Spectrum

    Indian Academy of Sciences (India)

    K Kumar; M N N Namboodiri; S Serra-Capizzano

    2014-05-01

    Let $A(x)$ be a norm continuous family of bounded self-adjoint operators on a separable Hilbert space $\\mathbb{H}$ and let $A(x)_n$ be the orthogonal compressions of $A(x)$ to the span of first elements of an orthonormal basis of $\\mathbb{H}$. The problem considered here is to approximate the spectrum of $A(x)$ using the sequence of eigenvalues of $A(x)_n$. We show that the bounds of the essential spectrum and the discrete spectral values outside the bounds of essential spectrum of $A(x)$ can be approximated uniformly on all compact subsets by the sequence of eigenvalue functions of $A(x)_n$. The known results, for a bounded self-adjoint operator, are translated into the case of a norm continuous family of operators. Also an attempt is made to predict the existence of spectral gaps that may occur between the bounds of essential spectrum of $A(0)=A$ and study the effect of norm continuous perturbation of operators in the prediction of spectral gaps. As an example, gap issues of some block Toeplitz–Laurent operators are discussed. The pure linear algebraic approach is the main advantage of the results here.

  2. Refining Approximating Betweenness Centrality Based on Samplings

    CERN Document Server

    Ji, Shiyu

    2016-01-01

    Betweenness Centrality (BC) is an important measure used widely in complex network analysis, such as social network, web page search, etc. Computing the exact BC values is highly time consuming. Currently the fastest exact BC determining algorithm is given by Brandes, taking $O(nm)$ time for unweighted graphs and $O(nm+n^2\\log n)$ time for weighted graphs, where $n$ is the number of vertices and $m$ is the number of edges in the graph. Due to the extreme difficulty of reducing the time complexity of exact BC determining problem, many researchers have considered the possibility of any satisfactory BC approximation algorithms, especially those based on samplings. Bader et al. give the currently best BC approximation algorithm, with a high probability to successfully estimate the BC of one vertex within a factor of $1/\\varepsilon$ using $\\varepsilon t$ samples, where $t$ is the ratio between $n^2$ and the BC value of the vertex. However, some of the algorithmic parameters in Bader's work are not yet tightly boun...

  3. Function approximation using adaptive and overlapping intervals

    Energy Technology Data Exchange (ETDEWEB)

    Patil, R.B.

    1995-05-01

    A problem common to many disciplines is to approximate a function given only the values of the function at various points in input variable space. A method is proposed for approximating a function of several to one variable. The model takes the form of weighted averaging of overlapping basis functions defined over intervals. The number of such basis functions and their parameters (widths and centers) are automatically determined using given training data and a learning algorithm. The proposed algorithm can be seen as placing a nonuniform multidimensional grid in the input domain with overlapping cells. The non-uniformity and overlap of the cells is achieved by a learning algorithm to optimize a given objective function. This approach is motivated by the fuzzy modeling approach and a learning algorithms used for clustering and classification in pattern recognition. The basics of why and how the approach works are given. Few examples of nonlinear regression and classification are modeled. The relationship between the proposed technique, radial basis neural networks, kernel regression, probabilistic neural networks, and fuzzy modeling is explained. Finally advantages and disadvantages are discussed.

  4. Classical approximation to quantum cosmological correlations

    CERN Document Server

    van der Meulen, Meindert

    2007-01-01

    We investigate up to which order quantum effects can be neglected in calculating cosmological correlation functions after horizon exit. As a toy model, we study $\\phi^3$ theory on a de Sitter background for a massless minimally coupled scalar field $\\phi$. We find that for tree level and one loop contributions in the quantum theory, a good classical approximation can be constructed, but for higher loop corrections this is in general not expected to be possible. The reason is that loop corrections get non-negligible contributions from loop momenta with magnitude up to the Hubble scale $H$, at which scale classical physics is not expected to be a good approximation to the quantum theory. An explicit calculation of the one loop correction to the two point function, supports the argument that contributions from loop momenta of scale $H$ are not negligible. Generalization of the arguments for the toy model to derivative interactions and the curvature perturbation leads to the conclusion that the leading orders of ...

  5. On Approximating String Selection Problems with Outliers

    CERN Document Server

    Boucher, Christina; Levy, Avivit; Pritchard, David; Weimann, Oren

    2012-01-01

    Many problems in bioinformatics are about finding strings that approximately represent a collection of given strings. We look at more general problems where some input strings can be classified as outliers. The Close to Most Strings problem is, given a set S of same-length strings, and a parameter d, find a string x that maximizes the number of "non-outliers" within Hamming distance d of x. We prove this problem has no PTAS unless ZPP=NP, correcting a decade-old mistake. The Most Strings with Few Bad Columns problem is to find a maximum-size subset of input strings so that the number of non-identical positions is at most k; we show it has no PTAS unless P=NP. We also observe Closest to k Strings has no EPTAS unless W[1]=FPT. In sum, outliers help model problems associated with using biological data, but we show the problem of finding an approximate solution is computationally difficult.

  6. On spline approximation of sliced inverse regression

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    The dimension reduction is helpful and often necessary in exploring the nonparametric regression structure.In this area,Sliced inverse regression (SIR) is a promising tool to estimate the central dimension reduction (CDR) space.To estimate the kernel matrix of the SIR,we herein suggest the spline approximation using the least squares regression.The heteroscedasticity can be incorporated well by introducing an appropriate weight function.The root-n asymptotic normality can be achieved for a wide range choice of knots.This is essentially analogous to the kernel estimation.Moreover, we also propose a modified Bayes information criterion (BIC) based on the eigenvalues of the SIR matrix.This modified BIC can be applied to any form of the SIR and other related methods.The methodology and some of the practical issues are illustrated through the horse mussel data.Empirical studies evidence the performance of our proposed spline approximation by comparison of the existing estimators.

  7. Adaptive and Approximate Orthogonal Range Counting

    DEFF Research Database (Denmark)

    Chan, Timothy M.; Wilkinson, Bryan Thomas

    2013-01-01

    ]. •We give an O(n loglog n)-space data structure for approximate 2-D orthogonal range counting that can compute a (1+δ)-factor approximation to the count in O(loglog n) time for any fixed constant δ>0. Again, our bounds match the state of the art for the 2-D orthogonal range emptiness problem. •Lastly......Close Abstract We present three new results on one of the most basic problems in geometric data structures, 2-D orthogonal range counting. All the results are in the w-bit word RAM model. •It is well known that there are linear-space data structures for 2-D orthogonal range counting with worst......-case optimal query time O(log_w n). We give an O(n loglog n)-space adaptive data structure that improves the query time to O(loglog n + log_w k), where k is the output count. When k=O(1), our bounds match the state of the art for the 2-D orthogonal range emptiness problem [Chan, Larsen, and Pătraşcu, SoCG 2011...

  8. On some applications of diophantine approximations.

    Science.gov (United States)

    Chudnovsky, G V

    1984-03-01

    Siegel's results [Siegel, C. L. (1929) Abh. Preuss. Akad. Wiss. Phys.-Math. Kl. 1] on the transcendence and algebraic independence of values of E-functions are refined to obtain the best possible bound for the measures of irrationality and linear independence of values of arbitrary E-functions at rational points. Our results show that values of E-functions at rational points have measures of diophantine approximations typical to "almost all" numbers. In particular, any such number has the "2 + epsilon" exponent of irrationality: Theta - p/q > q(-2-epsilon) for relatively prime rational integers p,q, with q >/= q(0) (Theta, epsilon). These results answer some problems posed by Lang. The methods used here are based on the introduction of graded Padé approximations to systems of functions satisfying linear differential equations with rational function coefficients. The constructions and proofs of this paper were used in the functional (nonarithmetic case) in a previous paper [Chudnovsky, D. V. & Chudnovsky, G. V. (1983) Proc. Natl. Acad. Sci. USA 80, 5158-5162]. PMID:16593441

  9. Chiral baryon in the coherent pair approximation

    CERN Document Server

    Aly, T S T

    1999-01-01

    We revisit the work of K. Goeke, M. Harvey, F. Grümmer, and J. N. Urbano (Phys. Rev. {\\bf D37}, 754 (1988)) who considered a chiral model for the nucleon based on the linear sigma model with scalar-isoscalar scalar-isovector mesons coupled to quarks and solved using the coherent-pair approximation. In this way the quantum pion field can be treated in a non-perturbative fashion. In this work we review this model and the coherent pair approximation correcting several errors in the earlier work. We minimize the expectation value of the chiral hamiltonian in the ansatz coherent-pair ground state configuration and solve the resulting equations for nucleon quantum numbers. We calculate the canonical set of nucleon observables and compare with the Hedgehog model and experiment. Using the corrected equations yield slightly different values for nucleon observables but do not correct the large virial deviation in the $\\pi$-nucleon coupling. Our results therefore do not significantly alter the conclusions of Goeke, et ...

  10. The validity of the Background Field Approximation

    CERN Document Server

    Parentani, R

    1997-01-01

    In the absence of a tractable theory of quantum gravity, quantum matter field effects have been so far computed by treating gravity at the Background Field Approximation. The principle aim of this paper is to investigate the validity of this approximation which is not specific to gravity. To this end, for reasons of simplicity and clarity, we shall compare the descriptions of thermal processes induced by constant acceleration (i.e. the Unruh effect) in four dynamical frameworks. In this problem, the position of the ``heavy'' accelerated system plays the role of gravity. In the first framework, the trajectory is treated at the BFA: it is given from the outset and unaffected by radiative processes. In the second one, recoil effects induced by these emission processes are taken into account by describing the system's position by WKB wave functions. In the third one, the accelerated system is described by second quantized fields and in the fourth one, gravity is turned on. It is most interesting to see when and w...

  11. Analytic approximate radiation effects due to Bremsstrahlung

    Energy Technology Data Exchange (ETDEWEB)

    Ben-Zvi I.

    2012-02-01

    The purpose of this note is to provide analytic approximate expressions that can provide quick estimates of the various effects of the Bremsstrahlung radiation produced relatively low energy electrons, such as the dumping of the beam into the beam stop at the ERL or field emission in superconducting cavities. The purpose of this work is not to replace a dependable calculation or, better yet, a measurement under real conditions, but to provide a quick but approximate estimate for guidance purposes only. These effects include dose to personnel, ozone generation in the air volume exposed to the radiation, hydrogen generation in the beam dump water cooling system and radiation damage to near-by magnets. These expressions can be used for other purposes, but one should note that the electron beam energy range is limited. In these calculations the good range is from about 0.5 MeV to 10 MeV. To help in the application of this note, calculations are presented as a worked out example for the beam dump of the R&D Energy Recovery Linac.

  12. Chiral Magnetic Effect in Hydrodynamic Approximation

    CERN Document Server

    Zakharov, Valentin I

    2012-01-01

    We review derivations of the chiral magnetic effect (ChME) in hydrodynamic approximation. The reader is assumed to be familiar with the basics of the effect. The main challenge now is to account for the strong interactions between the constituents of the fluid. The main result is that the ChME is not renormalized: in the hydrodynamic approximation it remains the same as for non-interacting chiral fermions moving in an external magnetic field. The key ingredients in the proof are general laws of thermodynamics and the Adler-Bardeen theorem for the chiral anomaly in external electromagnetic fields. The chiral magnetic effect in hydrodynamics represents a macroscopic manifestation of a quantum phenomenon (chiral anomaly). Moreover, one can argue that the current induced by the magnetic field is dissipation free and talk about a kind of "chiral superconductivity". More precise description is a ballistic transport along magnetic field taking place in equilibrium and in absence of a driving force. The basic limitat...

  13. APPROXIMATING INNOVATION POTENTIAL WITH NEUROFUZZY ROBUST MODEL

    Directory of Open Access Journals (Sweden)

    Kasa, Richard

    2015-01-01

    Full Text Available In a remarkably short time, economic globalisation has changed the world’s economic order, bringing new challenges and opportunities to SMEs. These processes pushed the need to measure innovation capability, which has become a crucial issue for today’s economic and political decision makers. Companies cannot compete in this new environment unless they become more innovative and respond more effectively to consumers’ needs and preferences – as mentioned in the EU’s innovation strategy. Decision makers cannot make accurate and efficient decisions without knowing the capability for innovation of companies in a sector or a region. This need is forcing economists to develop an integrated, unified and complete method of measuring, approximating and even forecasting the innovation performance not only on a macro but also a micro level. In this recent article a critical analysis of the literature on innovation potential approximation and prediction is given, showing their weaknesses and a possible alternative that eliminates the limitations and disadvantages of classical measuring and predictive methods.

  14. Regularity and approximability of electronic wave functions

    CERN Document Server

    Yserentant, Harry

    2010-01-01

    The electronic Schrödinger equation describes the motion of N-electrons under Coulomb interaction forces in a field of clamped nuclei. The solutions of this equation, the electronic wave functions, depend on 3N variables, with three spatial dimensions for each electron. Approximating these solutions is thus inordinately challenging, and it is generally believed that a reduction to simplified models, such as those of the Hartree-Fock method or density functional theory, is the only tenable approach. This book seeks to show readers that this conventional wisdom need not be ironclad: the regularity of the solutions, which increases with the number of electrons, the decay behavior of their mixed derivatives, and the antisymmetry enforced by the Pauli principle contribute properties that allow these functions to be approximated with an order of complexity which comes arbitrarily close to that for a system of one or two electrons. The text is accessible to a mathematical audience at the beginning graduate level as...

  15. Proportional damping approximation using the energy gain and simultaneous perturbation stochastic approximation

    Science.gov (United States)

    Sultan, Cornel

    2010-10-01

    The design of vector second-order linear systems for accurate proportional damping approximation is addressed. For this purpose an error system is defined using the difference between the generalized coordinates of the non-proportionally damped system and its proportionally damped approximation in modal space. The accuracy of the approximation is characterized using the energy gain of the error system and the design problem is formulated as selecting parameters of the non-proportionally damped system to ensure that this gain is sufficiently small. An efficient algorithm that combines linear matrix inequalities and simultaneous perturbation stochastic approximation is developed to solve the problem and examples of its application to tensegrity structures design are presented.

  16. Supercriticality for Annealed Approximations of Boolean Networks

    CERN Document Server

    Mountford, Thomas

    2010-01-01

    We consider a model proposed by Derrida and Pomeau (1986) and recently studied by Chatterjee and Durrett (2009); it is defined as an approximation to S. Kauffman's boolean networks (1969). The model starts with the choice of a random directed graph on $n$ vertices; each node has $r$ input nodes pointing at it. A discrete time threshold contact process is then considered on this graph: at each instant, each site has probability $q$ of choosing to receive input; if it does, and if at least one of its inputs were occupied by a $1$ at the previous instant, then it is labeled with a $1$; in all other cases, it is labeled with a $0$. $r$ and $q$ are kept fixed and $n$ is taken to infinity. Improving a result of Chatterjee and Durrett, we show that if $qr > 1$, then the time of persistence of the dynamics is exponential in $n$.

  17. Localization in Sensor Network With Nystrom Approximation

    Directory of Open Access Journals (Sweden)

    Shailaja Patil

    2011-11-01

    Full Text Available The recent innovations in wireless technology and digital electronics have opened many areas of research in Wireless Sensor Networks. In the last few years these networks have been successfully usedin many applications such as localization, tracking, surveillance, battlefield monitoring, structural health monitoring, routing etc. Most of these applications need localization i.e. estimating location information either relative or absolute. In this paper we propose a computationally efficient algorithm namely, Light weight Multidimensional Scaling. This approach takes advantage of Nystrӧm approximation for estimating location of unknown sensor nodes, using the information of available distances between neighbours and anchors. Various node densities, noise factors and radio ranges are considered for simulation. The performance of the algorithm is obtained with Monte Carlo Simulation.

  18. Squashed entanglement and approximate private states

    Science.gov (United States)

    Wilde, Mark M.

    2016-09-01

    The squashed entanglement is a fundamental entanglement measure in quantum information theory, finding application as an upper bound on the distillable secret key or distillable entanglement of a quantum state or a quantum channel. This paper simplifies proofs that the squashed entanglement is an upper bound on distillable key for finite-dimensional quantum systems and solidifies such proofs for infinite-dimensional quantum systems. More specifically, this paper establishes that the logarithm of the dimension of the key system (call it log 2K ) in an ɛ -approximate private state is bounded from above by the squashed entanglement of that state plus a term that depends only ɛ and log 2K . Importantly, the extra term does not depend on the dimension of the shield systems of the private state. The result holds for the bipartite squashed entanglement, and an extension of this result is established for two different flavors of the multipartite squashed entanglement.

  19. The semiclassical approximation to supersymmetric quantum gravity

    CERN Document Server

    Kiefer, C; Moniz, P; Kiefer, Claus; Lueck, Tobias; Moniz, Paulo

    2005-01-01

    We develop a semiclassical approximation scheme for the constraint equations of supersymmetric canonical quantum gravity. This is achieved by a Born-Oppenheimer type of expansion, in analogy to the case of the usual Wheeler-DeWitt equation. The formalism is only consistent if the states at each order depend on the gravitino field. We recover at consecutive orders the Hamilton-Jacobi equation, the functional Schrodinger equation, and quantum gravitational correction terms to this Schrodinger equation. In particular, the following consequences are found: (i) the Hamilton-Jacobi equation and therefore the background spacetime must involve the gravitino, (ii) a (many fingered) local time parameter has to be present on $SuperRiem \\Sigma$ (the space of all possible tetrad and gravitino fields), (iii) quantum supersymmetric gravitational corrections affect the evolution of the very early universe. The physical meaning of these equations and results, in particular the similarities to and differences from the pure bos...

  20. Subset Selection by Local Convex Approximation

    DEFF Research Database (Denmark)

    Øjelund, Henrik; Sadegh, Payman; Madsen, Henrik;

    1999-01-01

    This paper concerns selection of the optimal subset of variables in a lenear regression setting. The posed problem is combinatiorial and the globally best subset can only be found in exponential time. We define a cost function for the subset selection problem by adding the penalty term to the usual...... least squares criterion. We propose an optimization technique for the posed probelm based on a modified version of the Newton-Raphson iterations, combined with a backward elimination type algorithm. THe Newton-Raphson modification concerns iterative approximations to the non-convex cost function...... of the subset selection problem so as to guarantee positive definiteness of the Hessian term, hence avoiding numerical instability. The backward Elemination type algorithm attempts to improve the results upon termination of the modified Newton-Raphson search by sing the current solution as an initial guess...

  1. An Approximate Model of Microchannel Cooling

    Institute of Scientific and Technical Information of China (English)

    ShipingYu; MingdaoXin

    1994-01-01

    Forced convective heat transfer in micro-rectangular channels can be described by a group of two-dimensional differential equations.These equations take the conduction in microchannel wall along the direction of flow of coolants into account,which are more generalized than those which neglect the conduction.For the same reason,they are suitable particularly for gases-cooled microchannels.With only numerical solution to the equations till today,an approximate analytic solution is derived here,From this solution,a rather simple formula can be introduced further,by which the differences between considering the conduction and neglecting it are easily found.In addition,the reasonableness of the classical fin method is also discussed.An experimental example of air-cooled microchannels is illustrated.

  2. Approximating Conditional Density Functions Using Dimension Reduction

    Institute of Scientific and Technical Information of China (English)

    Jian-qing Fan; Liang Peng; Qi-wei Yao; Wen-yang Zhang

    2009-01-01

    We propose to approximate the conditional density function of a random variable Y given a dependent random d-vector X by that of Y given θτX, where the unit vector θ is selected such that the average Kullback-Leibler discrepancy distance between the two conditional density functions obtains the minimum. Our approach is nonparametric as far as the estimation of the conditional density functions is concerned. We have shown that this nonparametric estimator is asymptotically adaptive to the unknown index θ in the sense that the first order asymptotic mean squared error of the estimator is the same as that when θ was known. The proposed method is illustrated using both simulated and real-data examples.

  3. A realistic formulation of approximate CP

    CERN Document Server

    Dent, T; Dent, Thomas; Silva-Marcos, Joaquim

    2003-01-01

    CP violation in the SM is naturally implemented as a small imaginary perturbation to real Yukawa couplings. For example, a large CP asymmetry in B_d decays can arise if the imaginary parts of quark mass matrices are of order 10^(-3)m_t,b or smaller. Applying the same principle of ``additive CP violation'' to soft SUSY-breaking terms, the electric dipole moments of the neutron and mercury atom are predicted near current experimental limits; for nonuniversal A-terms, EDM bounds can be satisfied given certain flavour structures. The proposal is conveniently formulated in a democratic basis, with Yukawas and soft terms of the form const x (1+eps+i zeta) where eps<<1, zeta<~10^(-3), motivated by approximate permutation x CP symmetry.

  4. Approximate Methods for State-Space Models

    CERN Document Server

    Koyama, Shinsuke; Shalizi, Cosma Rohilla; Kass, Robert E; 10.1198/jasa.2009.tm08326

    2010-01-01

    State-space models provide an important body of techniques for analyzing time-series, but their use requires estimating unobserved states. The optimal estimate of the state is its conditional expectation given the observation histories, and computing this expectation is hard when there are nonlinearities. Existing filtering methods, including sequential Monte Carlo, tend to be either inaccurate or slow. In this paper, we study a nonlinear filter for nonlinear/non-Gaussian state-space models, which uses Laplace's method, an asymptotic series expansion, to approximate the state's conditional mean and variance, together with a Gaussian conditional distribution. This {\\em Laplace-Gaussian filter} (LGF) gives fast, recursive, deterministic state estimates, with an error which is set by the stochastic characteristics of the model and is, we show, stable over time. We illustrate the estimation ability of the LGF by applying it to the problem of neural decoding and compare it to sequential Monte Carlo both in simulat...

  5. Approximation diffuse Hermite et ses applications

    OpenAIRE

    Savignat, Jean-Michel

    2000-01-01

    De nombreuses techniques de résolution d'équations aux dérivées partielles sans maillage ont été développées dans la dernière décennie, proposant une alternative attrayante lorsque les éléments finis atteignent leurs limites. Notre travail se concentre sur l'étude de l'approximation diffuse, de ses applications au lissage et a la résolution des équations différentielles : les éléments diffus. Cependant, les solutions proposées s'appliquent aussi à d'autres méthodes et de nombreux résultats nu...

  6. Approximate truncation robust computed tomography—ATRACT

    Science.gov (United States)

    Dennerlein, Frank; Maier, Andreas

    2013-09-01

    We present an approximate truncation robust algorithm to compute tomographic images (ATRACT). This algorithm targets at reconstructing volumetric images from cone-beam projections in scenarios where these projections are highly truncated in each dimension. It thus facilitates reconstructions of small subvolumes of interest, without involving prior knowledge about the object. Our method is readily applicable to medical C-arm imaging, where it may contribute to new clinical workflows together with a considerable reduction of x-ray dose. We give a detailed derivation of ATRACT that starts from the conventional Feldkamp filtered-backprojection algorithm and that involves, as one component, a novel original formula for the inversion of the two-dimensional Radon transform. Discretization and numerical implementation are discussed and reconstruction results from both, simulated projections and first clinical data sets are presented.

  7. Approximation by max-product type operators

    CERN Document Server

    Bede, Barnabás; Gal, Sorin G

    2016-01-01

    This monograph presents a broad treatment of developments in an area of constructive approximation involving the so-called "max-product" type operators. The exposition highlights the max-product operators as those which allow one to obtain, in many cases, more valuable estimates than those obtained by classical approaches. The text considers a wide variety of operators which are studied for a number of interesting problems such as quantitative estimates, convergence, saturation results, localization, to name several. Additionally, the book discusses the perfect analogies between the probabilistic approaches of the classical Bernstein type operators and of the classical convolution operators (non-periodic and periodic cases), and the possibilistic approaches of the max-product variants of these operators. These approaches allow for two natural interpretations of the max-product Bernstein type operators and convolution type operators: firstly, as possibilistic expectations of some fuzzy variables, and secondly,...

  8. Wave system and its approximate similarity solutions

    Institute of Scientific and Technical Information of China (English)

    Liu Ping; Fu Pei-Kai

    2011-01-01

    Recently,a new (2+1)-dimensional shallow water wave system,the (2+1)-dimensional displacement shallow water wave system (2DDSWWS),was constructed by applying the variational principle of the analytic mechanics in the Lagrange coordinates. The disadvantage is that fluid viscidity is not considered in the 2DDSWWS,which is the same as the famous Kadomtsev-Petviashvili equation and Korteweg-de Vries equation. Applying dimensional analysis,we modify the 2DDSWWS and add the term related to the fluid viscidity to the 2DDSWWS. The approximate similarity solutions of the modified 2DDSWWS (M2DDSWWS) is studied and four similarity solutions are obtained. For the perfect fluids,the coefficient of kinematic viscosity is zero,then the M2DDSWWS will degenerate to the 2DDSWWS.

  9. Approximation algorithm for multiprocessor parallel job scheduling

    Institute of Scientific and Technical Information of China (English)

    陈松乔; 黄金贵; 陈建二

    2002-01-01

    Pk|fix|Cmax problem is a new scheduling problem based on the multiprocessor parallel job, and it is proved to be NP-hard problem when k≥3. This paper focuses on the case of k=3. Some new observations and new techniques for P3|fix|Cmax problem are offered. The concept of semi-normal schedulings is introduced, and a very simple linear time algorithm Semi-normal Algorithm for constructing semi-normal schedulings is developed. With the method of the classical Graham List Scheduling, a thorough analysis of the optimal scheduling on a special instance is provided, which shows that the algorithm is an approximation algorithm of ratio of 9/8 for any instance of P3|fix|Cmax problem, and improves the previous best ratio of 7/6 by M.X.Goemans.

  10. The Electroweak Sudakov approximation in SHERPA

    CERN Document Server

    Thompson, Jennifer M

    2016-01-01

    As experimental particle physics becomes more and more precise, it is becoming increasingly important for Monte Carlo simulations to improve the precision of their predictions. In terms of the hard matrix element, this means calculating to a higher order in perturbation theory. To be consistent this requires both NNLO QCD corrections and NLO EW corrections to be included. There are also interference effects between these processes that are not simple to handle consistently. For a broad description of the behaviour of NLO EW corrections at high energies, the Sudakov logarithmic approach provides a good approximation, and is much less computationally expensive than the full calculation. The implementation of EW Sudakov logarithms within the SHERPA program are outlined here along with some initial results. As well as this, an overview of the status of full NLO EW computations with SHERPA is presented.

  11. Comparing numerical and analytic approximate gravitational waveforms

    Science.gov (United States)

    Afshari, Nousha; Lovelace, Geoffrey; SXS Collaboration

    2016-03-01

    A direct observation of gravitational waves will test Einstein's theory of general relativity under the most extreme conditions. The Laser Interferometer Gravitational-Wave Observatory, or LIGO, began searching for gravitational waves in September 2015 with three times the sensitivity of initial LIGO. To help Advanced LIGO detect as many gravitational waves as possible, a major research effort is underway to accurately predict the expected waves. In this poster, I will explore how the gravitational waveform produced by a long binary-black-hole inspiral, merger, and ringdown is affected by how fast the larger black hole spins. In particular, I will present results from simulations of merging black holes, completed using the Spectral Einstein Code (black-holes.org/SpEC.html), including some new, long simulations designed to mimic black hole-neutron star mergers. I will present comparisons of the numerical waveforms with analytic approximations.

  12. Approximation Algorithms for Secondary Spectrum Auctions

    CERN Document Server

    Hoefer, Martin; Vöcking, Berthold

    2010-01-01

    We study combinatorial auctions for the secondary spectrum market. In this market, short-term licenses shall be given to wireless nodes for communication in their local neighborhood. In contrast to the primary market, channels can be assigned to multiple bidders, provided that the corresponding devices are well separated such that the interference is sufficiently low. Interference conflicts are described in terms of a conflict graph in which the nodes represent the bidders and the edges represent conflicts such that the feasible allocations for a channel correspond to the independent sets in the conflict graph. In this paper, we suggest a novel LP formulation for combinatorial auctions with conflict graph using a non-standard graph parameter, the so-called inductive independence number. Taking into account this parameter enables us to bypass the well-known lower bound of \\Omega(n^{1-\\epsilon}) on the approximability of independent set in general graphs with n nodes (bidders). We achieve significantly better a...

  13. Adaptive Control with Approximated Policy Search Approach

    Directory of Open Access Journals (Sweden)

    Agus Naba

    2010-05-01

    Full Text Available Most of existing adaptive control schemes are designed to minimize error between plant state and goal state despite the fact that executing actions that are predicted to result in smaller errors only can mislead to non-goal states. We develop an adaptive control scheme that involves manipulating a controller of a general type to improve its performance as measured by an evaluation function. The developed method is closely related to a theory of Reinforcement Learning (RL but imposes a practical assumption made for faster learning. We assume that a value function of RL can be approximated by a function of Euclidean distance from a goal state and an action executed at the state. And, we propose to use it for the gradient search as an evaluation function. Simulation results provided through application of the proposed scheme to a pole-balancing problem using a linear state feedback controller and fuzzy controller verify the scheme’s efficacy.

  14. Random Phase Approximation in Relativistic Approach

    Institute of Scientific and Technical Information of China (English)

    MA Zhong-yu; YANG Ding; TIAN Yuan; CAO Li-gang

    2009-01-01

    Some special issues of the random phase approximation (RPA) in the relativistic approach are reviewed.A full consistency and proper treatment of coupling to the continuum are responsible for the successful application of the RPA in the description of dynamical properties of finite nuclei.The fully consistent relativistic RPA(RRPA) requires that the relativistic mean filed (RMF) wave function of the nucleus and the RRPA correlations are calculated in a same effective Lagrangian and the consistent treatment of the Dirac sea of negative energy states.The proper treatment of the single particle continuum with scattering asymptotic conditions in the RMF and RRPA is discussed.The full continuum spectrum can be described by the single particle Green's function and the relativistic continuum RPA is established.A separable form of the paring force is introduced in the relativistic quasi-particle RPA.

  15. Some Undecidable Problems on Approximability of NP Optimization Problems

    Institute of Scientific and Technical Information of China (English)

    黄雄

    1996-01-01

    In this paper some undecidable problems on approximability of NP optimization problems are investigated.In particular,the following problems are all undecidable:(1) Given an NP optimization problem,is it approximable in polynomial time?(2)For any polynomial-time computable function r(n),given a polynomial time approximable NP optimization problem,has it a polynomial-time approximation algorithm with approximation performance ratio r(n) (r(n)-approximable)?(3)For any polynomial-time computable functions r(n),r'(n),where r'(n)approximable NP optimization problem,is it r'(n)-approximable?

  16. Stochastic approximation with state-dependent noise

    Institute of Scientific and Technical Information of China (English)

    2000-01-01

    The purpose of stochastic approximation (SA) is to find the roots of f(·) or the maximizer (minimizer) of L(·) when the unknown function f(·) or L(·) can be observed but with noise. SA is an important tool in dealing with many problems arising from systems and control, whose solutions often rely on convergence of the SA algorithm applied. Here the pathwise convergence of SA algorithms is considered, when the observation noise may depend on state by which we mean those x at which f(x) or L(x) are observed. The conditions imposed on the observation noise are the weakest in comparison with the existing ones. When the algorithm is to find the roots of f(·), the superiority of the condition given in the paper over those used in literature consists in the fact that the present condition is directly verifiable, needless to see the behaviour of the algorithm. When the algorithm is to find the maximizer (minimizer) of L(·), the present conditioin allows the observation noise to depend on the state. The conditions imposed on f(·) and L(·) are truly general: f(·) is required to be measurable and locally bounded if the roots of f(·) are sought, and the gradient of L(·) is required to be locally Lipschitz continuous if the maximizer (minimizer) of L(·) is searched.

  17. Configuring Airspace Sectors with Approximate Dynamic Programming

    Science.gov (United States)

    Bloem, Michael; Gupta, Pramod

    2010-01-01

    In response to changing traffic and staffing conditions, supervisors dynamically configure airspace sectors by assigning them to control positions. A finite horizon airspace sector configuration problem models this supervisor decision. The problem is to select an airspace configuration at each time step while considering a workload cost, a reconfiguration cost, and a constraint on the number of control positions at each time step. Three algorithms for this problem are proposed and evaluated: a myopic heuristic, an exact dynamic programming algorithm, and a rollouts approximate dynamic programming algorithm. On problem instances from current operations with only dozens of possible configurations, an exact dynamic programming solution gives the optimal cost value. The rollouts algorithm achieves costs within 2% of optimal for these instances, on average. For larger problem instances that are representative of future operations and have thousands of possible configurations, excessive computation time prohibits the use of exact dynamic programming. On such problem instances, the rollouts algorithm reduces the cost achieved by the heuristic by more than 15% on average with an acceptable computation time.

  18. Precision variational approximations in statistical data assimilation

    Directory of Open Access Journals (Sweden)

    J. Ye

    2014-10-01

    Full Text Available Data assimilation transfers information from observations of a complex system to physically-based system models with state variables x(t. Typically, the observations are noisy, the model has errors, and the initial state of the model is uncertain, so the data assimilation is statistical. One can thus ask questions about expected values of functions ⟨G(X⟩ on the path X = {x(t0, ..., x(tm} of the model as it moves through an observation window where measurements are made at times {t0, ..., tm}. The probability distribution on the path P(X = exp[−A0(X] determines these expected values. Variational methods seeking extrema of the "action" A0(X, widely known as 4DVar (Talagrand and Courtier, 1987; Evensen, 2009,, are widespread for estimating ⟨G(X ⟩ in many fields of science. In a path integral formulation of statistical data assimilation, we consider variational approximations in a standard realization of the action where measurement and model errors are Gaussian. We (a discuss an annealing method for locating the path X0 giving a consistent global minimum of the action A0(X0, (b consider the explicit role of the number of measurements at each measurement time in determining A0(X0, and (c identify a parameter regime for the scale of model errors which allows X0 to give a precise estimate of ⟨G(X0⟩ with computable, small higher order corrections.

  19. Semiclassical approximation to supersymmetric quantum gravity

    Science.gov (United States)

    Kiefer, Claus; Lück, Tobias; Moniz, Paulo

    2005-08-01

    We develop a semiclassical approximation scheme for the constraint equations of supersymmetric canonical quantum gravity. This is achieved by a Born-Oppenheimer type of expansion, in analogy to the case of the usual Wheeler-DeWitt equation. The formalism is only consistent if the states at each order depend on the gravitino field. We recover at consecutive orders the Hamilton-Jacobi equation, the functional Schrödinger equation, and quantum gravitational correction terms to this Schrödinger equation. In particular, the following consequences are found: (i) the Hamilton-Jacobi equation and therefore the background spacetime must involve the gravitino, (ii) a (many-fingered) local time parameter has to be present on super Riem Σ (the space of all possible tetrad and gravitino fields), (iii) quantum supersymmetric gravitational corrections affect the evolution of the very early Universe. The physical meaning of these equations and results, in particular, the similarities to and differences from the pure bosonic case, are discussed.

  20. Coulomb crystals in the harmonic lattice approximation

    CERN Document Server

    Baiko, D A; De Witt, H E; Slattery, W L

    2000-01-01

    The dynamic structure factor ${\\tilde S}({\\bf k},\\omega)$ and the two-particle distribution function $g({\\bf r},t)$ of ions in a Coulomb crystal are obtained in a closed analytic form using the harmonic lattice (HL) approximation which takes into account all processes of multi-phonon excitation and absorption. The static radial two-particle distribution function $g(r)$ is calculated for classical ($T \\gtrsim \\hbar \\omega_p$, where $\\omega_p$ is the ion plasma frequency) and quantum ($T \\ll \\hbar \\omega_p$) body-centered cubic (bcc) crystals. The results for the classical crystal are in a very good agreement with extensive Monte Carlo (MC) calculations at $1.5 \\lesssim r/a calculated for classical and quantum bcc and face-centered cubic crystals, and anharmonic corrections are discussed. The inelastic part of the HL static structure factor $S''(k)$, averaged over orientations of wave-vector {\\bf k}, is shown to contain pronounced singularities at Bragg diffraction positions. The type of the singularities is di...

  1. Distributed Verification and Hardness of Distributed Approximation

    CERN Document Server

    Sarma, Atish Das; Kor, Liah; Korman, Amos; Nanongkai, Danupon; Pandurangan, Gopal; Peleg, David; Wattenhofer, Roger

    2010-01-01

    We study the {\\em verification} problem in distributed networks, stated as follows. Let $H$ be a subgraph of a network $G$ where each vertex of $G$ knows which edges incident on it are in $H$. We would like to verify whether $H$ has some properties, e.g., if it is a tree or if it is connected. We would like to perform this verification in a decentralized fashion via a distributed algorithm. The time complexity of verification is measured as the number of rounds of distributed communication. In this paper we initiate a systematic study of distributed verification, and give almost tight lower bounds on the running time of distributed verification algorithms for many fundamental problems such as connectivity, spanning connected subgraph, and $s-t$ cut verification. We then show applications of these results in deriving strong unconditional time lower bounds on the {\\em hardness of distributed approximation} for many classical optimization problems including minimum spanning tree, shortest paths, and minimum cut....

  2. Approximate Model for Turbulent Stagnation Point Flow.

    Energy Technology Data Exchange (ETDEWEB)

    Dechant, Lawrence [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

    2016-01-01

    Here we derive an approximate turbulent self-similar model for a class of favorable pressure gradient wedge-like flows, focusing on the stagnation point limit. While the self-similar model provides a useful gross flow field estimate this approach must be combined with a near wall model is to determine skin friction and by Reynolds analogy the heat transfer coefficient. The combined approach is developed in detail for the stagnation point flow problem where turbulent skin friction and Nusselt number results are obtained. Comparison to the classical Van Driest (1958) result suggests overall reasonable agreement. Though the model is only valid near the stagnation region of cylinders and spheres it nonetheless provides a reasonable model for overall cylinder and sphere heat transfer. The enhancement effect of free stream turbulence upon the laminar flow is used to derive a similar expression which is valid for turbulent flow. Examination of free stream enhanced laminar flow suggests that the rather than enhancement of a laminar flow behavior free stream disturbance results in early transition to turbulent stagnation point behavior. Excellent agreement is shown between enhanced laminar flow and turbulent flow behavior for high levels, e.g. 5% of free stream turbulence. Finally the blunt body turbulent stagnation results are shown to provide realistic heat transfer results for turbulent jet impingement problems.

  3. Dynamical Vertex Approximation for the Hubbard Model

    Science.gov (United States)

    Toschi, Alessandro

    A full understanding of correlated electron systems in the physically relevant situations of three and two dimensions represents a challenge for the contemporary condensed matter theory. However, in the last years considerable progress has been achieved by means of increasingly more powerful quantum many-body algorithms, applied to the basic model for correlated electrons, the Hubbard Hamiltonian. Here, I will review the physics emerging from studies performed with the dynamical vertex approximation, which includes diagrammatic corrections to the local description of the dynamical mean field theory (DMFT). In particular, I will first discuss the phase diagram in three dimensions with a special focus on the commensurate and incommensurate magnetic phases, their (quantum) critical properties, and the impact of fluctuations on electronic lifetimes and spectral functions. In two dimensions, the effects of non-local fluctuations beyond DMFT grow enormously, determining the appearance of a low-temperature insulating behavior for all values of the interaction in the unfrustrated model: Here the prototypical features of the Mott-Hubbard metal-insulator transition, as well as the existence of magnetically ordered phases, are completely overwhelmed by antiferromagnetic fluctuations of exponentially large extension, in accordance with the Mermin-Wagner theorem. Eventually, by a fluctuation diagnostics analysis of cluster DMFT self-energies, the same magnetic fluctuations are identified as responsible for the pseudogap regime in the holed-doped frustrated case, with important implications for the theoretical modeling of the cuprate physics.

  4. The time-dependent Gutzwiller approximation

    Science.gov (United States)

    Fabrizio, Michele

    2015-03-01

    The time-dependent Gutzwiller Approximation (t-GA) is shown to be capable of tracking the off-equilibrium evolution both of coherent quasiparticles and of incoherent Hubbard bands. The method is used to demonstrate that the sharp dynamical crossover observed by time-dependent DMFT in the quench-dynamics of a half-filled Hubbard model can be identified within the t-GA as a genuine dynamical transition separating two distinct physical phases. This result, strictly variational for lattices of infinite coordination number, is intriguing as it actually questions the occurrence of thermalization. Next, we shall present how t-GA works in a multi-band model for V2O3 that displays a first-order Mott transition. We shall show that a physically accessible excitation pathway is able to collapse the Mott gap down and drive off-equilibrium the insulator into a metastable metal phase. Work supported by the European Union, Seventh Framework Programme, under the project GO FAST, Grant Agreement No. 280555.

  5. Adaptive approximation of higher order posterior statistics

    KAUST Repository

    Lee, Wonjung

    2014-02-01

    Filtering is an approach for incorporating observed data into time-evolving systems. Instead of a family of Dirac delta masses that is widely used in Monte Carlo methods, we here use the Wiener chaos expansion for the parametrization of the conditioned probability distribution to solve the nonlinear filtering problem. The Wiener chaos expansion is not the best method for uncertainty propagation without observations. Nevertheless, the projection of the system variables in a fixed polynomial basis spanning the probability space might be a competitive representation in the presence of relatively frequent observations because the Wiener chaos approach not only leads to an accurate and efficient prediction for short time uncertainty quantification, but it also allows to apply several data assimilation methods that can be used to yield a better approximate filtering solution. The aim of the present paper is to investigate this hypothesis. We answer in the affirmative for the (stochastic) Lorenz-63 system based on numerical simulations in which the uncertainty quantification method and the data assimilation method are adaptively selected by whether the dynamics is driven by Brownian motion and the near-Gaussianity of the measure to be updated, respectively. © 2013 Elsevier Inc.

  6. Network histograms and universality of blockmodel approximation

    Science.gov (United States)

    Olhede, Sofia C.; Wolfe, Patrick J.

    2014-01-01

    In this paper we introduce the network histogram, a statistical summary of network interactions to be used as a tool for exploratory data analysis. A network histogram is obtained by fitting a stochastic blockmodel to a single observation of a network dataset. Blocks of edges play the role of histogram bins and community sizes that of histogram bandwidths or bin sizes. Just as standard histograms allow for varying bandwidths, different blockmodel estimates can all be considered valid representations of an underlying probability model, subject to bandwidth constraints. Here we provide methods for automatic bandwidth selection, by which the network histogram approximates the generating mechanism that gives rise to exchangeable random graphs. This makes the blockmodel a universal network representation for unlabeled graphs. With this insight, we discuss the interpretation of network communities in light of the fact that many different community assignments can all give an equally valid representation of such a network. To demonstrate the fidelity-versus-interpretability tradeoff inherent in considering different numbers and sizes of communities, we analyze two publicly available networks—political weblogs and student friendships—and discuss how to interpret the network histogram when additional information related to node and edge labeling is present. PMID:25275010

  7. Extending the random-phase approximation for electronic correlation energies: the renormalized adiabatic local density approximation

    DEFF Research Database (Denmark)

    Olsen, Thomas; Thygesen, Kristian S.

    2012-01-01

    The adiabatic connection fluctuation-dissipation theorem with the random phase approximation (RPA) has recently been applied with success to obtain correlation energies of a variety of chemical and solid state systems. The main merit of this approach is the improved description of dispersive forces...

  8. Bilayer graphene spectral function in the random phase approximation and self-consistent GW approximation

    Science.gov (United States)

    Sabashvili, Andro; Östlund, Stellan; Granath, Mats

    2013-08-01

    We calculate the single-particle spectral function for doped bilayer graphene in the low energy limit, described by two parabolic bands with zero band gap and long range Coulomb interaction. Calculations are done using thermal Green's functions in both the random phase approximation (RPA) and the fully self-consistent GW approximation. Consistent with previous studies RPA yields a spectral function which, apart from the Landau quasiparticle peaks, shows additional coherent features interpreted as plasmarons, i.e., composite electron-plasmon excitations. In the GW approximation the plasmaron becomes incoherent and peaks are replaced by much broader features. The deviation of the quasiparticle weight and mass renormalization from their noninteracting values is small which indicates that bilayer graphene is a weakly interacting system. The electron energy loss function, Im[-ɛq-1(ω)] shows a sharp plasmon mode in RPA which in the GW approximation becomes less coherent and thus consistent with the weaker plasmaron features in the corresponding single-particle spectral function.

  9. Collisionless magnetic reconnection under anisotropic MHD approximation

    Science.gov (United States)

    Hirabayashi, Kota; Hoshino, Masahiro

    We study the formation of slow-mode shocks in collisionless magnetic reconnection by using one- and two-dimensional collisionless magneto-hydro-dynamic (MHD) simulations based on the double adiabatic approximation, which is an important step to bridge the gap between the Petschek-type MHD reconnection model accompanied by a pair of slow shocks and the observational evidence of the rare occasion of in-situ slow shock observation. According to our results, a pair of slow shocks does form in the reconnection layer. The resultant shock waves, however, are quite weak compared with those in an isotropic MHD from the point of view of the plasma compression and the amount of the magnetic energy released across the shock. Once the slow shock forms, the downstream plasma are heated in highly anisotropic manner and a firehose-sense (P_{||}>P_{⊥}) pressure anisotropy arises. The maximum anisotropy is limited by the marginal firehose criterion, 1-(P_{||}-P_{⊥})/B(2) =0. In spite of the weakness of the shocks, the resultant reconnection rate is kept at the same level compared with that in the corresponding ordinary MHD simulations. It is also revealed that the sequential order of propagation of the slow shock and the rotational discontinuity, which appears when the guide field component exists, changes depending on the magnitude of the guide field. Especially, when no guide field exists, the rotational discontinuity degenerates with the contact discontinuity remaining at the position of the initial current sheet, while with the slow shock in the isotropic MHD. Our result implies that the slow shock does not necessarily play an important role in the energy conversion in the reconnection system and is consistent with the satellite observation in the Earth's magnetosphere.

  10. Generalized stationary phase approximations for mountain waves

    Science.gov (United States)

    Knight, H.; Broutman, D.; Eckermann, S. D.

    2016-04-01

    Large altitude asymptotic approximations are derived for vertical displacements due to mountain waves generated by hydrostatic wind flow over arbitrary topography. This leads to new asymptotic analytic expressions for wave-induced vertical displacement for mountains with an elliptical Gaussian shape and with the major axis oriented at any angle relative to the background wind. The motivation is to understand local maxima in vertical displacement amplitude at a given height for elliptical mountains aligned at oblique angles to the wind direction, as identified in Eckermann et al. ["Effects of horizontal geometrical spreading on the parameterization of orographic gravity-wave drag. Part 1: Numerical transform solutions," J. Atmos. Sci. 72, 2330-2347 (2015)]. The standard stationary phase method reproduces one type of local amplitude maximum that migrates downwind with increasing altitude. Another type of local amplitude maximum stays close to the vertical axis over the center of the mountain, and a new generalized stationary phase method is developed to describe this other type of local amplitude maximum and the horizontal variation of wave-induced vertical displacement near the vertical axis of the mountain in the large altitude limit. The new generalized stationary phase method describes the asymptotic behavior of integrals where the asymptotic parameter is raised to two different powers (1/2 and 1) rather than just one power as in the standard stationary phase method. The vertical displacement formulas are initially derived assuming a uniform background wind but are extended to accommodate both vertical shear with a fixed wind direction and vertical variations in the buoyancy frequency.

  11. Hydration thermodynamics beyond the linear response approximation.

    Science.gov (United States)

    Raineri, Fernando O

    2016-10-19

    The solvation energetics associated with the transformation of a solute molecule at infinite dilution in water from an initial state A to a final state B is reconsidered. The two solute states have different potentials energies of interaction, [Formula: see text] and [Formula: see text], with the solvent environment. Throughout the A [Formula: see text] B transformation of the solute, the solvation system is described by a Hamiltonian [Formula: see text] that changes linearly with the coupling parameter ξ. By focusing on the characterization of the probability density [Formula: see text] that the dimensionless perturbational solute-solvent interaction energy [Formula: see text] has numerical value y when the coupling parameter is ξ, we derive a hierarchy of differential equation relations between the ξ-dependent cumulant functions of various orders in the expansion of the appropriate cumulant generating function. On the basis of this theoretical framework we then introduce an inherently nonlinear solvation model for which we are able to find analytical results for both [Formula: see text] and for the solvation thermodynamic functions. The solvation model is based on the premise that there is an upper or a lower bound (depending on the nature of the interactions considered) to the amplitude of the fluctuations of Y in the solution system at equilibrium. The results reveal essential differences in behavior for the model when compared with the linear response approximation to solvation, particularly with regards to the probability density [Formula: see text]. The analytical expressions for the solvation properties show, however, that the linear response behavior is recovered from the new model when the room for the thermal fluctuations in Y is not restricted by the existence of a nearby bound. We compare the predictions of the model with the results from molecular dynamics computer simulations for aqueous solvation, in which either (1) the solute

  12. Coronal Loops: Evolving Beyond the Isothermal Approximation

    Science.gov (United States)

    Schmelz, J. T.; Cirtain, J. W.; Allen, J. D.

    2002-05-01

    Are coronal loops isothermal? A controversy over this question has arisen recently because different investigators using different techniques have obtained very different answers. Analysis of SOHO-EIT and TRACE data using narrowband filter ratios to obtain temperature maps has produced several key publications that suggest that coronal loops may be isothermal. We have constructed a multi-thermal distribution for several pixels along a relatively isolated coronal loop on the southwest limb of the solar disk using spectral line data from SOHO-CDS taken on 1998 Apr 20. These distributions are clearly inconsistent with isothermal plasma along either the line of sight or the length of the loop, and suggested rather that the temperature increases from the footpoints to the loop top. We speculated originally that these differences could be attributed to pixel size -- CDS pixels are larger, and more `contaminating' material would be expected along the line of sight. To test this idea, we used CDS iron line ratios from our data set to mimic the isothermal results from the narrowband filter instruments. These ratios indicated that the temperature gradient along the loop was flat, despite the fact that a more complete analysis of the same data showed this result to be false! The CDS pixel size was not the cause of the discrepancy; rather, the problem lies with the isothermal approximation used in EIT and TRACE analysis. These results should serve as a strong warning to anyone using this simplistic method to obtain temperature. This warning is echoed on the EIT web page: ``Danger! Enter at your own risk!'' In other words, values for temperature may be found, but they may have nothing to do with physical reality. Solar physics research at the University of Memphis is supported by NASA grant NAG5-9783. This research was funded in part by the NASA/TRACE MODA grant for Montana State University.

  13. Bond selective chemistry beyond the adiabatic approximation

    Energy Technology Data Exchange (ETDEWEB)

    Butler, L.J. [Univ. of Chicago, IL (United States)

    1993-12-01

    One of the most important challenges in chemistry is to develop predictive ability for the branching between energetically allowed chemical reaction pathways. Such predictive capability, coupled with a fundamental understanding of the important molecular interactions, is essential to the development and utilization of new fuels and the design of efficient combustion processes. Existing transition state and exact quantum theories successfully predict the branching between available product channels for systems in which each reaction coordinate can be adequately described by different paths along a single adiabatic potential energy surface. In particular, unimolecular dissociation following thermal, infrared multiphoton, or overtone excitation in the ground state yields a branching between energetically allowed product channels which can be successfully predicted by the application of statistical theories, i.e. the weakest bond breaks. (The predictions are particularly good for competing reactions in which when there is no saddle point along the reaction coordinates, as in simple bond fission reactions.) The predicted lack of bond selectivity results from the assumption of rapid internal vibrational energy redistribution and the implicit use of a single adiabatic Born-Oppenheimer potential energy surface for the reaction. However, the adiabatic approximation is not valid for the reaction of a wide variety of energetic materials and organic fuels; coupling between the electronic states of the reacting species play a a key role in determining the selectivity of the chemical reactions induced. The work described below investigated the central role played by coupling between electronic states in polyatomic molecules in determining the selective branching between energetically allowed fragmentation pathways in two key systems.

  14. Hydration thermodynamics beyond the linear response approximation

    Science.gov (United States)

    Raineri, Fernando O.

    2016-10-01

    The solvation energetics associated with the transformation of a solute molecule at infinite dilution in water from an initial state A to a final state B is reconsidered. The two solute states have different potentials energies of interaction, {{\\Psi}\\text{A}} and {{\\Psi}\\text{B}} , with the solvent environment. Throughout the A \\to B transformation of the solute, the solvation system is described by a Hamiltonian H≤ft(ξ \\right) that changes linearly with the coupling parameter ξ. By focusing on the characterization of the probability density {{\\wp}ξ}≤ft( y\\right) that the dimensionless perturbational solute-solvent interaction energy Y=β ≤ft({{\\Psi}\\text{B}}-{{\\Psi}\\text{A}}\\right) has numerical value y when the coupling parameter is ξ, we derive a hierarchy of differential equation relations between the ξ-dependent cumulant functions of various orders in the expansion of the appropriate cumulant generating function. On the basis of this theoretical framework we then introduce an inherently nonlinear solvation model for which we are able to find analytical results for both {{\\wp}ξ} ≤ft( y\\right) and for the solvation thermodynamic functions. The solvation model is based on the premise that there is an upper or a lower bound (depending on the nature of the interactions considered) to the amplitude of the fluctuations of Y in the solution system at equilibrium. The results reveal essential differences in behavior for the model when compared with the linear response approximation to solvation, particularly with regards to the probability density {{\\wp}ξ} ≤ft( y\\right) . The analytical expressions for the solvation properties show, however, that the linear response behavior is recovered from the new model when the room for the thermal fluctuations in Y is not restricted by the existence of a nearby bound. We compare the predictions of the model with the results from molecular dynamics computer simulations for aqueous solvation, in

  15. Approximate nearest neighbors via dictionary learning

    Science.gov (United States)

    Cherian, Anoop; Morellas, Vassilios; Papanikolopoulos, Nikolaos

    2011-06-01

    Approximate Nearest Neighbors (ANN) in high dimensional vector spaces is a fundamental, yet challenging problem in many areas of computer science, including computer vision, data mining and robotics. In this work, we investigate this problem from the perspective of compressive sensing, especially the dictionary learning aspect. High dimensional feature vectors are seldom seen to be sparse in the feature domain; examples include, but not limited to Scale Invariant Feature Transform (SIFT) descriptors, Histogram Of Gradients, Shape Contexts, etc. Compressive sensing advocates that if a given vector has a dense support in a feature space, then there should exist an alternative high dimensional subspace where the features are sparse. This idea is leveraged by dictionary learning techniques through learning an overcomplete projection from the feature space so that the vectors are sparse in the new space. The learned dictionary aids in refining the search for the nearest neighbors to a query feature vector into the most likely subspace combination indexed by its non-zero active basis elements. Since the size of the dictionary is generally very large, distinct feature vectors are most likely to have distinct non-zero basis. Utilizing this observation, we propose a novel representation of the feature vectors as tuples of non-zero dictionary indices, which then reduces the ANN search problem into hashing the tuples to an index table; thereby dramatically improving the speed of the search. A drawback of this naive approach is that it is very sensitive to feature perturbations. This can be due to two possibilities: (i) the feature vectors are corrupted by noise, (ii) the true data vectors undergo perturbations themselves. Existing dictionary learning methods address the first possibility. In this work we investigate the second possibility and approach it from a robust optimization perspective. This boils down to the problem of learning a dictionary robust to feature

  16. A Lattice-Theoretic Approach to Multigranulation Approximation Space

    OpenAIRE

    Xiaoli He; Yanhong She

    2014-01-01

    In this paper, we mainly investigate the equivalence between multigranulation approximation space and single-granulation approximation space from the lattice-theoretic viewpoint. It is proved that multigranulation approximation space is equivalent to single-granulation approximation space if and only if the pair of multigranulation rough approximation operators ( Σ i = 1 n R i ¯ , Σ i = 1 n R i _ ) forms an order-preserving Galois connection, if and only if the collection of lower (resp., upp...

  17. Rational offset approximation of rational Bézier curves

    Institute of Scientific and Technical Information of China (English)

    CHENG Min; WANG Guo-jin

    2006-01-01

    The problem of parametric speed approximation of a rational curve is raised in this paper. Offset curves are widely used in various applications. As for the reason that in most cases the offset curves do not preserve the same polynomial or rational polynomial representations, it arouses difficulty in applications. Thus approximation methods have been introduced to solve this problem. In this paper, it has been pointed out that the crux of offset curve approximation lies in the approximation of parametric speed. Based on the Jacobi polynomial approximation theory with endpoints interpolation, an algebraic rational approximation algorithm of offset curve, which preserves the direction of normal, is presented.

  18. APPROXIMATE DUALITY OF g-FRAMES IN HILBERT SPACES

    Institute of Scientific and Technical Information of China (English)

    Amir KHOSRAVI; Morteza MIRZAEE AZANDARYANI

    2014-01-01

    In this article, we introduce and characterize approximate duality for g-frames. We get some important properties and applications of approximate duals. We also obtain some new results in approximate duality of frames, and generalize some of the known results in approximate duality of frames to g-frames. We also get some results for fusion frames, and perturbation of approximately dual g-frames. We show that approximate duals are stable under small perturbations and they are useful for erasures and reconstruction.

  19. Towards a Philosophy of Approximations in the 'Exact' Sciences

    Directory of Open Access Journals (Sweden)

    Valentin N. Ostrovsky

    2005-11-01

    Full Text Available The issue of approximations is mostly neglected in the philosophy of science, and sometimes misinterpreted. The paper demonstrates that approximations are in fact in the core of some recent discussions in the philosophy of chemistry: on the shape of molecules, the Born-Oppenheimer approximation, the role of orbitals, and the physical explanation of the Periodic Table of Elements. The ontological and epistemological significance of approximations in the exact sciences is analyzed. The crucial role of approximations in generating qualitative images and comprehensible models is emphasized. A complementarity relation between numerically 'exact' theories and explanatory approximate approaches is claimed.

  20. Approximation law for discrete-time variable structure control systems

    Institute of Scientific and Technical Information of China (English)

    Yan ZHENG; Yuanwei JING

    2006-01-01

    Two approximation laws of sliding mode for discrete-time variable structure control systems are proposed to overcome the limitations of the exponential approximation law and the variable rate approximation law. By applying the proposed approximation laws of sliding mode to discrete-time variable structure control systems, the stability of origin can be guaranteed, and the chattering along the switching surface caused by discrete-time variable structure control can be restrained effectively. In designing of approximation laws, the problem that the system control input is restricted is also considered, which is very important in practical systems. Finally a simulation example shows the effectiveness of the two approximation laws proposed.

  1. A working-set framework for sequential convex approximation methods

    DEFF Research Database (Denmark)

    Stolpe, Mathias

    2008-01-01

    We present an active-set algorithmic framework intended as an extension to existing implementations of sequential convex approximation methods for solving nonlinear inequality constrained programs. The framework is independent of the choice of approximations and the stabilization technique used...

  2. Approximate Generalized Conditional Symmetries for Perturbed Evolution Equations

    Institute of Scientific and Technical Information of China (English)

    ZHANG Shun-Li; WANG Yong; LOU Sen-Yue

    2007-01-01

    The concept of approximate generalized conditional symmetry (AGCS) for the perturbed evolution equations is introduced, and how to derive approximate conditional invariant solutions to the perturbed equations via their AGCSs is illustrated with examples.

  3. Stochastic approximation Monte Carlo importance sampling for approximating exact conditional probabilities

    KAUST Repository

    Cheon, Sooyoung

    2013-02-16

    Importance sampling and Markov chain Monte Carlo methods have been used in exact inference for contingency tables for a long time, however, their performances are not always very satisfactory. In this paper, we propose a stochastic approximation Monte Carlo importance sampling (SAMCIS) method for tackling this problem. SAMCIS is a combination of adaptive Markov chain Monte Carlo and importance sampling, which employs the stochastic approximation Monte Carlo algorithm (Liang et al., J. Am. Stat. Assoc., 102(477):305-320, 2007) to draw samples from an enlarged reference set with a known Markov basis. Compared to the existing importance sampling and Markov chain Monte Carlo methods, SAMCIS has a few advantages, such as fast convergence, ergodicity, and the ability to achieve a desired proportion of valid tables. The numerical results indicate that SAMCIS can outperform the existing importance sampling and Markov chain Monte Carlo methods: It can produce much more accurate estimates in much shorter CPU time than the existing methods, especially for the tables with high degrees of freedom. © 2013 Springer Science+Business Media New York.

  4. On the Consistency of the Ladder Approximation and the Rainbow Approximation of Dyson-Schwinger Equation of QCD

    CERN Document Server

    Chang, L; Chang, Lei; Liu, Yu-xin

    2006-01-01

    We study the consistency of the ladder approximation and the rainbow approximation of the Dyson-Schwinger equation of QCD. By considering the non-Abelian property of QCD, we show that the QED-type Ward-Takahashi identity is not required for the rainbow-ladder approximation of QCD. It indicates that there does not exists any internal inconsistency in the usual rainbow-ladder approximation of QCD. In addition, we propose an modified ladder approximation which guarantees the Slavnov-Taylor identity for the quark-gluon vertex omitting the ghost effect in the approximation.

  5. Pawlak algebra and approximate structure on fuzzy lattice.

    Science.gov (United States)

    Zhuang, Ying; Liu, Wenqi; Wu, Chin-Chia; Li, Jinhai

    2014-01-01

    The aim of this paper is to investigate the general approximation structure, weak approximation operators, and Pawlak algebra in the framework of fuzzy lattice, lattice topology, and auxiliary ordering. First, we prove that the weak approximation operator space forms a complete distributive lattice. Then we study the properties of transitive closure of approximation operators and apply them to rough set theory. We also investigate molecule Pawlak algebra and obtain some related properties.

  6. SOME NONLINEAR APPROXIMATIONS FOR MATRIX-VALUED FUNCTIONS

    Institute of Scientific and Technical Information of China (English)

    Guo-liang Xu

    2003-01-01

    Some nonlinear approximants, i.e., exponential-sum interpolation with equal distance or at origin, (0,1)-type, (0,2)-type and (1,2)-type fraction-sum approximations, for matrixvalued functions are introduced. All these approximation problems lead to a same form system of nonlinear equations. Solving methods for the nonlinear system are discussed.Conclusions on uniqueness and convergence of the approximants for certain class of functions are given.

  7. ON APPROXIMATION BY SPHERICAL REPRODUCING KERNEL HILBERT SPACES

    Institute of Scientific and Technical Information of China (English)

    2007-01-01

    The spherical approximation between two nested reproducing kernels Hilbert spaces generated from different smooth kernels is investigated. It is shown that the functions of a space can be approximated by that of the subspace with better smoothness. Furthermore, the upper bound of approximation error is given.

  8. A UNIFIED APPROACH TO CERTAIN PROBLEMS OF BEST LOCAL APPROXIMATION

    Institute of Scientific and Technical Information of China (English)

    H.H.Cuenya; M.D.Lorenzo; C.N.Rodriguez

    2007-01-01

    In this paper we study best local quasi-rational approximation and best local approximation from finite dimensional subspaces of vectorial functions of several variables.Our approach extends and unifies several problems concerning best local multi-point approximation in different norms.

  9. Approximate inference for spatial functional data on massively parallel processors

    DEFF Research Database (Denmark)

    Rakêt, Lars Lau; Markussen, Bo

    2014-01-01

    is considered, and so-called operator approximations for doing inference in the resulting models are presented. These approximations embed observations in function space, transferring likelihood calculations to the functional domain. The resulting approximated problems are naturally parallel and can be solved...

  10. Real space Dynamical Super Cell Approximation for interacting disordered systems

    OpenAIRE

    Moradian, Rostam

    2004-01-01

    Effective medium super-cell approximation method which is introduced for disordered systems is extended to a general case of interacting disordered systems. We found that the dynamical cluster approximation (DCA) and also the non local coherent potential approximation (NLCPA) are two simple case of this technique. Whole equations of this formalism derived by using the effective medium theory in real space.

  11. The Use of Approximations in a High School Chemistry Course

    Science.gov (United States)

    Matsumoto, Paul S.; Tong, Gary; Lee, Stephanie; Kam, Bonita

    2009-01-01

    While approximations are used frequently in science, high school students may be unaware of the use of approximations in science, the motivation for their use, and the limitations of their use. In the article, we consider the use of approximations in a high school chemistry class as opportunities to increase student understanding of the use of…

  12. Bioluminescence tomography based on the phase approximation model

    OpenAIRE

    Cong, W; Wang, G.

    2010-01-01

    A reconstruction method of bioluminescence sources is proposed based on a phase approximation model. Compared with the diffuse approximation, this phase approximation model more correctly predicts bioluminescence photon propagation in biological tissues, so that bioluminescence tomography can accurately locate and quantify the distribution of bioluminescence sources. The compressive sensing (CS) technique is applied to regularize the inverse source reconstruction to enhance numerical stabilit...

  13. Soft Rough Approximation Operators on a Complete Atomic Boolean Lattice

    Directory of Open Access Journals (Sweden)

    Heba I. Mustafa

    2013-01-01

    Full Text Available The concept of soft sets based on complete atomic Boolean lattice, which can be seen as a generalization of soft sets, is introduced. Some operations on these soft sets are discussed, and new types of soft sets such as full, keeping infimum, and keeping supremum are defined and supported by some illustrative examples. Two pairs of new soft rough approximation operators are proposed and the relationship among soft set is investigated, and their related properties are given. We show that Järvinen's approximations can be viewed as a special case of our approximation. If , then our soft approximations coincide with crisp soft rough approximations (Feng et al. 2011.

  14. Approximating the spin distributions in capture reactions between spherical nuclei

    Energy Technology Data Exchange (ETDEWEB)

    Chushnyakova, M.V., E-mail: maria.chushnyakova@gmail.com [Physics Department, Omsk State Technical University, 644050 Omsk (Russian Federation); Applied Physics Department, Tomsk Polytechnic University, 634028 Tomsk (Russian Federation); Gontchar, I.I. [Physics and Chemistry Department, Omsk State Transport University, 644046 Omsk (Russian Federation)

    2015-09-15

    Twenty years ago an approximation for the spin distribution of the dinuclear systems formed in capture reactions with heavy ions was proposed. This approximation is used nowadays. However, since that time the experimental errors of the measured capture cross sections were reduced drastically. We show that the accuracy of the old spin distribution approximation is certainly out of date and propose a new approximation built on the dynamical modeling of the capture process. Results suggest that this new approximation might be useful especially for performing modeling of decay of excited dinuclear systems (compound nuclei) formed during heavy-ion collisions.

  15. Approximating the spin distributions in capture reactions between spherical nuclei

    Science.gov (United States)

    Chushnyakova, M. V.; Gontchar, I. I.

    2015-09-01

    Twenty years ago an approximation for the spin distribution of the dinuclear systems formed in capture reactions with heavy ions was proposed. This approximation is used nowadays. However, since that time the experimental errors of the measured capture cross sections were reduced drastically. We show that the accuracy of the old spin distribution approximation is certainly out of date and propose a new approximation built on the dynamical modeling of the capture process. Results suggest that this new approximation might be useful especially for performing modeling of decay of excited dinuclear systems (compound nuclei) formed during heavy-ion collisions.

  16. Approximate Augmented Lagrangian Functions and Nonlinear Semidefinite Programs

    Institute of Scientific and Technical Information of China (English)

    X. X. HUANG; K. L. TEO; X. Q. YANG

    2006-01-01

    In this paper, an approximate augmented Lagrangian function for nonlinear semidefinite programs is introduced. Some basic properties of the approximate augmented Lagrange function such as monotonicity and convexity are discussed. Necessary and sufficient conditions for approximate strong duality results are derived. Conditions for an approximate exact penalty representation in the framework of augmented Lagrangian are given. Under certain conditions, it is shown that any limit point of a sequence of stationary points of approximate augmented Lagrangian problems is a KKT point of the original semidefinite program and that a sequence of optimal solutions to augmented Lagrangian problems converges to a solution of the original semidefinite program.

  17. Approximate Counting for Complex-Weighted Boolean Constraint Satisfaction Problems

    CERN Document Server

    Yamakami, Tomoyuki

    2010-01-01

    Constraint satisfaction problems (or CSPs) have been extensively studied in AI, database theory, graph theory, etc. From an approximation viewpoint, it has been important to approximate the total number of assignments that satisfy all given Boolean constraints. There is a trichotomy theorem for such approximate counting for (non-weighted) Boolean CSPs; namely, all such counting problems are neatly classified into three categories under polynomial-time approximation-preserving reductions [Dyer, Goldberg, and Jerrum, 2010]. We extend this result to approximate counting for complex-weighted Boolean CSPs, provided that all arity-1 constraints are freely available to use. This makes a significant progress in the quest for the approximation classification of all counting Boolean CSPs in the most general form. To deal with complex weights, we employ proof techniques along the line of solving Holant problems [Valiant, 2002, 2008]. Our result also gives an approximation version of the dichotomy theorem of the complexi...

  18. Methods of Approximation Theory in Complex Analysis and Mathematical Physics

    CERN Document Server

    Saff, Edward

    1993-01-01

    The book incorporates research papers and surveys written by participants ofan International Scientific Programme on Approximation Theory jointly supervised by Institute for Constructive Mathematics of University of South Florida at Tampa, USA and the Euler International Mathematical Instituteat St. Petersburg, Russia. The aim of the Programme was to present new developments in Constructive Approximation Theory. The topics of the papers are: asymptotic behaviour of orthogonal polynomials, rational approximation of classical functions, quadrature formulas, theory of n-widths, nonlinear approximation in Hardy algebras,numerical results on best polynomial approximations, wavelet analysis. FROM THE CONTENTS: E.A. Rakhmanov: Strong asymptotics for orthogonal polynomials associated with exponential weights on R.- A.L. Levin, E.B. Saff: Exact Convergence Rates for Best Lp Rational Approximation to the Signum Function and for Optimal Quadrature in Hp.- H. Stahl: Uniform Rational Approximation of x .- M. Rahman, S.K. ...

  19. Approximate trace and singleton failures equivalences for transition systems

    Institute of Scientific and Technical Information of China (English)

    Chao Wang; Jinzhao Wu; Hongyan Tan

    2015-01-01

    Established system equivalences for transition systems, such as trace equivalence and failures equivalence, require the ob-servations to be exactly identical. However, an accurate measure-ment is impossible when interacting with the physical world, hence exact equivalence is restrictive and not robust. Using Baire met-ric, a generalized framework of transition system approximation is proposed by developing the notions of approximate language equivalence and approximate singleton failures (SF) equivalence. The framework takes the traditional exact equivalence as a special case. The approximate language equivalence is coarser than the approximate SF equivalence, just like the hierarchy of the exact ones. The main conclusion is that the two approximate equiva-lences satisfy the transitive property, consequently, they can be successively used in transition system approximation.

  20. Density functional approximation for van der Waals fluids: based on hard sphere density functional approximation

    Institute of Scientific and Technical Information of China (English)

    Zhou Shi-Qi

    2007-01-01

    A universal theoretical approach is proposed which enables all hard sphere density functional approximations(DFAs) applicable to van der Waals fluids. The resultant DFA obtained by combining the universal theoretical approach with any hard sphere DFAs only needs as input a second-order direct correlation function (DCF) of a coexistence bulk fluid, and is applicable in both supercritical and subcritical temperature regions. The associated effective hard sphere density can be specified by a hard wall sum rule. It is indicated that the value of the effective hard sphere density so determined can be universal, i.e. can be applied to any external potentials different from the single hard wall. As an illustrating example, the universal theoretical approach is combined with a hard sphere bridge DFA to predict the density profile of a hard core attractive Yukawa model fluid influenced by diverse external fields; agreement between the present formalism's predictions and the corresponding simulation data is good or at least comparable to several previous DFT approaches. The primary advantage of the present theoretical approach combined with other hard sphere DFAs is discussed.

  1. High order source approximation for the EFEN method

    International Nuclear Information System (INIS)

    The flat source approximation in one dimensional Exponential Function Expansion Nodal (EFEN) method is extended to a high order polynomial approximation while maintaining the simplicity of the nodal response matrix. By applying the new method to a one dimensional PWR pin-by-pin problem, it has been observed that quadratic source approximation is good enough for PWR pin-by-pin calculation, while the flat source approximation causes about 5% of relative error to the thermal flux. By applying the new method to a one dimensional assembly homogenized problem, it has been found that the EFEN method with cubic source approximation can be employed to handle PWR core diffusion problems. Numerical results suggest the optimization of source approximation order for different energy groups and different spacial locations to achieve more accurate results with less computing effort. (author)

  2. Energetics of a fluid under the Boussinesq approximation

    CERN Document Server

    Maruyama, Kiyoshi

    2014-01-01

    This paper presents a theory describing the energy budget of a fluid under the Boussinesq approximation: the theory is developed in a manner consistent with the conservation law of mass. It shows that no potential energy is available under the Boussinesq approximation, and also reveals that the work done by the buoyancy force due to changes in temperature corresponds to the conversion between kinetic and internal energy. This energy conversion, however, makes only an ignorable contribution to the distribution of temperature under the approximation. The Boussinesq approximation is, in physical oceanography, extended so that the motion of seawater can be studied. This paper considers this extended approximation as well. Under the extended approximation, the work done by the buoyancy force due to changes in salinity corresponds to the conversion between kinetic and potential energy. It also turns out that the conservation law of mass does not allow the condition $\

  3. Upper bounds on minimum cardinality of exact and approximate reducts

    KAUST Repository

    Chikalov, Igor

    2010-01-01

    In the paper, we consider the notions of exact and approximate decision reducts for binary decision tables. We present upper bounds on minimum cardinality of exact and approximate reducts depending on the number of rows (objects) in the decision table. We show that the bound for exact reducts is unimprovable in the general case, and the bound for approximate reducts is almost unimprovable in the general case. © 2010 Springer-Verlag Berlin Heidelberg.

  4. Energy of Bardeen Model Using Approximate Symmetry Method

    OpenAIRE

    Sharif, M.; Waheed, Saira

    2010-01-01

    In this paper, we investigate the energy problem in general relativity using approximate Lie symmetry methods for differential equations. This procedure is applied to Bardeen model (the regular black hole solution). Here we are forced to evaluate the third-order approximate symmetries of the orbital and geodesic equations. It is shown that energy must be re-scaled by some factor in the third-order approximation. We discuss the insights of this re-scaling factor.

  5. Approximate analytical calculations of photon geodesics in the Schwarzschild metric

    CERN Document Server

    De Falco, Vittorio; Stella, Luigi

    2016-01-01

    We develop a method for deriving approximate analytical formulae to integrate photon geodesics in a Schwarzschild spacetime. Based on this, we derive the approximate equations for light bending and propagation delay that have been introduced empirically. We then derive for the first time an approximate analytical equation for the solid angle. We discuss the accuracy and range of applicability of the new equations and present a few simple applications of them to known astrophysical problems.

  6. Adaptive Algorithms of Nonlinear Approximation with Finite Terms

    Institute of Scientific and Technical Information of China (English)

    Wen Bin WEI; Yue Sheng XU; Pei Xin YE

    2007-01-01

    This paper deals with realizable adaptive algorithms of the nonlinear approximation with finite terms based on wavelets. We present a concrete algorithm by which we may find the required index set Am for the greedy algorithm GPm(·,ψ). This makes the greedy algorithm realize the near best approximation in practice. Moreover, we study the efficiency of the finite-term approximation of another algorithm introduced by Birge and Massart.

  7. Approximations of continuous Newton's method: An extension of Cayley's problem

    Directory of Open Access Journals (Sweden)

    Jon Jacobsen

    2007-02-01

    Full Text Available Continuous Newton's Method refers to a certain dynamical system whose associated flow generically tends to the roots of a given polynomial. An Euler approximation of this system, with step size $h=1$, yields the discrete Newton's method algorithm for finding roots. In this note we contrast Euler approximations with several different approximations of the continuous ODE system and, using computer experiments, consider their impact on the associated fractal basin boundaries of the roots.

  8. The cost of approximate controllability for semilinear heat equations

    Institute of Scientific and Technical Information of China (English)

    Yuqing YAN; Yi ZHAO; Yu HUANG

    2009-01-01

    We consider the semilinear heat equation with globally Lipschitz non-linearity involving gradient terms in a bounded domain of Rn.In this paper,we obtain explicit bounds of the cost of approximate controllability,i.e.,of the minimal norm of a control needed to control the system approximately.The methods we used combine global Carleman estimates,the variational approach to approximate controllability and Schauder's fixed point theorem.

  9. Closed-Form Approximations for Spread Option Prices and Greeks

    OpenAIRE

    Li, Minqiang

    2008-01-01

    We develop a new closed-form approximation method for pricing spread options. Numerical analysis shows that our method is more accurate than existing analytical approximations. Our method is also extremely fast, with computing time more than two orders of magnitude shorter than one-dimensional numerical integration. We also develop closed-form approximations for the greeks of spread options. In addition, we analyze the price sensitivities of spread options and provide lower and upper bounds f...

  10. APPROXIMATION LAWS OF DISCRETE-TIME VARIABLE STRUCTURE CONTROL SYSTEMS

    Institute of Scientific and Technical Information of China (English)

    2008-01-01

    Two new approximation laws of sliding mode for discrete-time variable structure control systems are proposed in this paper. By applying the proposed approximation laws of sliding mode to discrete-time variable structure control systems,the stability of origin can be guaranteed,and the chattering along the switching surface caused by discrete-time variable structure control can be restrained effectively. In designing of these approximation laws,the problem that the system control input is restricted is also ...

  11. An Approximation Method of NURBS Curves in NC Machining

    Institute of Scientific and Technical Information of China (English)

    YUE Ying; HAN Qingyao; WANG Zhangqi

    2006-01-01

    An algorithm for approximating arbitrary NURBS curve with straight line is presented. Firstly, NURBS curve is acquired according to data points on the curve. Secondly, Approximating arbitrary NURBS curve is based on dichotomy. The resulting straight line approaches to the original curve with relatively fewer segments within the required tolerance. The example shows that the algorithm is simple and its approximation precision is high. The method is most useful in numerical control to drive the cutter along straight line or circular paths.

  12. CHARACTERIZATION OF BEST APPROXIMATIONS IN METRIC LINEAR SPACES

    Institute of Scientific and Technical Information of China (English)

    Sizwe Mabizela

    2003-01-01

    Let (X,d) be a real metric linear space, with translation-invariant metric d and G a linear subspace of X. In this paper we use functionals in the Lipschitz dual of X to characterize those elements of G which are best approximations to elements of X. We also give simultaneous characterization of elements of best approximation and also consider elements of e-approximation.

  13. ON APPROXIMATION BY REPRODUCING KERNEL SPACES IN WEIGHTED Lp SPACES

    Institute of Scientific and Technical Information of China (English)

    Baohuai SHENG

    2007-01-01

    In this paper, we investigate the order of approximation by reproducing kernel spaces on (-1, 1) in weighted Lp spaces. We first restate the translation network from the view of reproducing kernel spaces and then construct a sequence of approximating operators with the help of Jacobi orthogonal polynomials, with which we establish a kind of Jackson inequality to describe the error estimate.Finally, The results are used to discuss an approximation problem arising from learning theory.

  14. Approximate Lie Group Analysis of Finite-difference Equations

    OpenAIRE

    Latypov, Azat M.

    1995-01-01

    Approximate group analysis technique, that is, the technique combining the methodology of group analysis and theory of small perturbations, is applied to finite-difference equations approximating ordinary differential equations. Finite-difference equations are viewed as a system of algebraic equations with a small parameter, introduced through the definitions of finite-difference derivatives. It is shown that application of the approximate invariance criterion to this algebraic system results...

  15. Generalized Lorentzian approximations for the Voigt line shape.

    Science.gov (United States)

    Martin, P; Puerta, J

    1981-01-15

    The object of the work reported in this paper was to find a simple and easy to calculate approximation to the Voigt function using the Padé method. To do this we calculated the multipole approximation to the complex function as the error function or as the plasma dispersion function. This generalized Lorentzian approximation can be used instead of the exact function in experiments that do not require great accuracy. PMID:20309100

  16. Proving acceptability properties of relaxed nondeterministic approximate programs

    OpenAIRE

    Carbin, Michael James; Kim, Deokhwan; Misailovic, Sasa; Rinard, Martin C

    2012-01-01

    Approximate program transformations such as skipping tasks [29, 30], loop perforation [21, 22, 35], reduction sampling [38], multiple selectable implementations [3, 4, 16, 38], dynamic knobs [16], synchronization elimination [20, 32], approximate function memoization [11],and approximate data types [34] produce programs that can execute at a variety of points in an underlying performance versus accuracy tradeoff space. These transformed programs have the ability to trade accuracy of their res...

  17. Generalized eikonal approximation for strong-field ionization

    OpenAIRE

    Vélez, F. Cajiao; Krajewska, K.; Kamiński, J. Z.

    2015-01-01

    We develop the eikonal perturbation theory to describe the strong-field ionization by finite laser pulses. This approach in the first order with respect to the binding potential (the so-called generalized eikonal approximation) avoids a singularity at the potential center. Thus, in contrast to the ordinary eikonal approximation, it allows to treat rescattering phenomena in terms of quantum trajectories. We demonstrate how the first Born approximation and its domain of validity follow from eik...

  18. A Linear-Programming Approximation of AC Power Flows

    OpenAIRE

    Coffrin, Carleton; Van Hentenryck, Pascal

    2012-01-01

    Linear active-power-only DC power flow approximations are pervasive in the planning and control of power systems. However, these approximations fail to capture reactive power and voltage magnitudes, both of which are necessary in many applications to ensure voltage stability and AC power flow feasibility. This paper proposes linear-programming models (the LPAC models) that incorporate reactive power and voltage magnitudes in a linear power flow approximation. The LPAC models are built on a co...

  19. Approximation for a large-angle simple pendulum period

    Energy Technology Data Exchange (ETDEWEB)

    Belendez, A; Rodes, J J; Belendez, T; Hernandez, A [Departamento de Fisica, IngenierIa de Sistemas y TeorIa de la Senal, Universidad de Alicante, Apartado 99, E-03080 Alicante (Spain)], E-mail: a.belendez@ua.es

    2009-03-15

    An approximation scheme to obtain the period for large amplitude oscillations of a simple pendulum is analysed and discussed. The analytical approximate formula for the period is the same as that suggested by Hite (2005 Phys. Teach. 43 290), but it is now obtained analytically by means of a term-by-term comparison of the power-series expansion for the approximate period with the corresponding series for the exact period. (letters and comments)

  20. Kernel Approximate Bayesian Computation for Population Genetic Inferences

    OpenAIRE

    Nakagome, Shigeki; Fukumizu, Kenji; Mano, Shuhei

    2012-01-01

    Approximate Bayesian computation (ABC) is a likelihood-free approach for Bayesian inferences based on a rejection algorithm method that applies a tolerance of dissimilarity between summary statistics from observed and simulated data. Although several improvements to the algorithm have been proposed, none of these improvements avoid the following two sources of approximation: 1) lack of sufficient statistics: sampling is not from the true posterior density given data but from an approximate po...

  1. Accuracy of the relativistic Cowling approximation in slowly rotating stars

    OpenAIRE

    Yoshida, Shijun; Kojima, Yasufumi

    1997-01-01

    We have calculated the non-radial oscillation in slowly rotating relativistic stars with the Cowling approximation. The frequencies are compared with those based on the complete linearized equations of general relativity. It is found that the results with the approximation differ by less than about $20 %$ for typical relativistic stellar models. The approximation is more accurate for higher-order modes as in the Newtonian case.

  2. APPROXIMATION OF RANDOM SUMS OF RANDOM VARIABLES IN INSURANCE

    Directory of Open Access Journals (Sweden)

    Paweł Wlaź

    2014-09-01

    Full Text Available The paper deals with approximations of random sums. By random sum we mean a sum of random number of independent and identically distributed random variables. Distribution of this sum is called a compound distribution. The model is especially important in non-life insurance. There are many methods for approximating compound distributions, one of the most popular one is approximation with shifted gamma distribution. In this work we show an alternative way – using kernel density, Fast Fourier Transform and numerical optimization methods – for finding shifted gamma approximations and show results suggesting its superiority over classical method.

  3. Legendre-tau approximations for functional differential equations

    Science.gov (United States)

    Ito, K.; Teglas, R.

    1986-01-01

    The numerical approximation of solutions to linear retarded functional differential equations are considered using the so-called Legendre-tau method. The functional differential equation is first reformulated as a partial differential equation with a nonlocal boundary condition involving time-differentiation. The approximate solution is then represented as a truncated Legendre series with time-varying coefficients which satisfy a certain system of ordinary differential equations. The method is very easy to code and yields very accurate approximations. Convergence is established, various numerical examples are presented, and comparison between the latter and cubic spline approximation is made.

  4. Divergence measures and a general framework for local variational approximation.

    Science.gov (United States)

    Watanabe, Kazuho; Okada, Masato; Ikeda, Kazushi

    2011-12-01

    The local variational method is a technique to approximate an intractable posterior distribution in Bayesian learning. This article formulates a general framework for local variational approximation and shows that its objective function is decomposable into the sum of the Kullback information and the expected Bregman divergence from the approximating posterior distribution to the Bayesian posterior distribution. Based on a geometrical argument in the space of approximating posteriors, we propose an efficient method to evaluate an upper bound of the marginal likelihood. Moreover, we demonstrate that the variational Bayesian approach for the latent variable models can be viewed as a special case of this general framework.

  5. Two-potential eikonal approximation for electron-atom collisions

    International Nuclear Information System (INIS)

    The Glauber approximation is known to be in appreciable error at all angles when applied to the elastic electron-atom scattering at medium and lower energies. It is shown that this is not due to the frozen-target approximation but mainly a result of the inadequate semiclassical treatment of close-encounter collisions in the Glauber approximation. A simple method is proposed to correct this inadequacy and is applied to e-H elastic scattering at energies from 20 to 100 eV. A remarkable improvement over the Glauber approximation is obtained, and the results agree with experiments very well at all angles where measurements are available

  6. A Quotient Space Approximation Model of Multiresolution Signal Analysis

    Institute of Scientific and Technical Information of China (English)

    Ling Zhang; Bo Zhang

    2005-01-01

    In this paper, we present a quotient space approximation model of multiresolution signal analysis and discuss the properties and characteristics of the model. Then the comparison between wavelet transform and the quotient space approximation is made. First, when wavelet transform is viewed from the new quotient space approximation perspective, it may help us to gain an insight into the essence of multiresolution signal analysis. Second, from the similarity between wavelet and quotient space approximations, it is possible to transfer the rich wavelet techniques into the latter so that a new way for multiresolution analysis may be found.

  7. On the Approximation and Smoothed Complexity of Leontief Market Equilibria

    OpenAIRE

    Huang, Li-Sha; Teng, Shang-Hua

    2006-01-01

    We show that the problem of finding an \\epsilon-approximate Nash equilibrium of an n by n two-person games can be reduced to the computation of an (\\epsilon/n)^2-approximate market equilibrium of a Leontief economy. Together with a recent result of Chen, Deng and Teng, this polynomial reduction implies that the Leontief market exchange problem does not have a fully polynomial-time approximation scheme, that is, there is no algorithm that can compute an \\epsilon-approximate market equilibrium ...

  8. Scattering from rough thin films: discrete-dipole-approximation simulations.

    Science.gov (United States)

    Parviainen, Hannu; Lumme, Kari

    2008-01-01

    We investigate the wave-optical light scattering properties of deformed thin circular films of constant thickness using the discrete-dipole approximation. Effects on the intensity distribution of the scattered light due to different statistical roughness models, model dependent roughness parameters, and uncorrelated, random, small-scale porosity of the inhomogeneous medium are studied. The suitability of the discrete-dipole approximation for rough-surface scattering problems is evaluated by considering thin films as computationally feasible rough-surface analogs. The effects due to small-scale inhomogeneity of the scattering medium are compared with the analytic approximation by Maxwell Garnett, and the results are found to agree with the approximation.

  9. Accuracy of the non-relativistic approximation for momentum diffusion

    Science.gov (United States)

    Liang, Shiuan-Ni; Lan, Boon Leong

    2016-06-01

    The accuracy of the non-relativistic approximation, which is calculated using the same parameter and the same initial ensemble of trajectories, to relativistic momentum diffusion at low speed is studied numerically for a prototypical nonlinear Hamiltonian system -the periodically delta-kicked particle. We find that if the initial ensemble is a non-localized semi-uniform ensemble, the non-relativistic approximation to the relativistic mean square momentum displacement is always accurate. However, if the initial ensemble is a localized Gaussian, the non-relativistic approximation may not always be accurate and the approximation can break down rapidly.

  10. Impact of inflow transport approximation on light water reactor analysis

    Science.gov (United States)

    Choi, Sooyoung; Smith, Kord; Lee, Hyun Chul; Lee, Deokjung

    2015-10-01

    The impact of the inflow transport approximation on light water reactor analysis is investigated, and it is verified that the inflow transport approximation significantly improves the accuracy of the transport and transport/diffusion solutions. A methodology for an inflow transport approximation is implemented in order to generate an accurate transport cross section. The inflow transport approximation is compared to the conventional methods, which are the consistent-PN and the outflow transport approximations. The three transport approximations are implemented in the lattice physics code STREAM, and verification is performed for various verification problems in order to investigate their effects and accuracy. From the verification, it is noted that the consistent-PN and the outflow transport approximations cause significant error in calculating the eigenvalue and the power distribution. The inflow transport approximation shows very accurate and precise results for the verification problems. The inflow transport approximation shows significant improvements not only for the high leakage problem but also for practical large core problem analyses.

  11. The Log-Linear Return Approximation, Bubbles, and Predictability

    DEFF Research Database (Denmark)

    Engsted, Tom; Pedersen, Thomas Quistgaard; Tanggaard, Carsten

    2012-01-01

    We study in detail the log-linear return approximation introduced by Campbell and Shiller (1988a). First, we derive an upper bound for the mean approximation error, given stationarity of the log dividend-price ratio. Next, we simulate various rational bubbles which have explosive conditional...

  12. The log-linear return approximation, bubbles, and predictability

    DEFF Research Database (Denmark)

    Engsted, Tom; Pedersen, Thomas Quistgaard; Tanggaard, Carsten

    We study in detail the log-linear return approximation introduced by Campbell and Shiller (1988a). First, we derive an upper bound for the mean approximation error, given stationarity of the log dividendprice ratio. Next, we simulate various rational bubbles which have explosive conditional...

  13. The best simultaneous approximation in linear 2-normed spaces

    CERN Document Server

    Acikgoz, Mehmet

    2012-01-01

    In this paper, we shall investigate and analyse a new study on the best simultaneous approximation in the context of linear 2-normed spaces inspired by Elumalai and his coworkers in Elumalai. The basis of this investigation is to extend and refinement the definition of the classical aproximation, best approximation and some related concepts to linear 2-normed spaces.

  14. Nonrelativistic atom-photon interaction beyond the multipole approximation

    NARCIS (Netherlands)

    Boers, D.J.; Wijers, C.M.J.

    2003-01-01

    We investigate the interaction between the hydrogen atomic orbitals and the quantized modes of the electromagnetic field within the domain of nonrelativistic quantum electrodynamics in the Coulomb gauge. Contrary to the conventional dipole approximation and higher-order multipole approximations, whi

  15. Sharp Bounds for Symmetric and Asymmetric Diophantine Approximation

    Institute of Scientific and Technical Information of China (English)

    Cornelis KRAAIKAMP; Ionica SMEETS

    2011-01-01

    In 2004,Tong found bounds for the approximation quality of a regular continued fraction convergent to a rational number,expressed in bounds for both the previous and next approximation.The authors sharpen his results with a geometric method and give both sharp upper and lower bounds.The asymptotic frequencies that these bounds occur are also calculated.

  16. Democracy functions and optimal embeddings for approximation spaces

    CERN Document Server

    Garrigós, Gustavo; de Natividade, Maria

    2009-01-01

    We prove optimal embeddings for nonlinear approximation spaces in terms of weighted Lorentz sequence spaces, with the weights depending on the democracy functions of the basis. As applications we recover known embeddings for $N$-term wavelet approximation in Lebesgue, Orlicz, and Lorentz norms. We also study the "greedy classes" introduced by Gribonval and Nielsen.

  17. Continuum Random Phase Approximation with finite-range interactions

    Energy Technology Data Exchange (ETDEWEB)

    Co' , Giampaolo [Universita del Salento, Dipartimento di Fisica ' ' E. De Giorgi' ' , Lecce (Italy); INFN, Sezione di Lecce, Lecce (Italy); De Donno, Viviana [Universita del Salento, Dipartimento di Fisica ' ' E. De Giorgi' ' , Lecce (Italy); Anguiano, Marta; Lallena, Antonio M. [Universidad de Granada, Departamento de Fisica Atomica, Molecular y Nuclear, Granada (Spain)

    2016-05-15

    We rewrite the Random Phase Approximation secular equations in a form which allows the treatment of the continuum part of the single-particle spectrum without approximations. Within this formalism finite-range interactions can be used without restrictions. We present some results, obtained with Gogny interactions, where the role of the continuum is relevant. (orig.)

  18. Saddlepoint approximations: a review and some new applications

    NARCIS (Netherlands)

    S.A. Broda; M.S. Paolella

    2012-01-01

    The saddlepoint method of approximation is attributed to Daniels (1954), and can be described in basic terms as yielding an accurate and usually fast and very numerically reliable approximation to the mass or density function (hereafter pdf), and the cumulative distribution function (cdf), of a rand

  19. An Approximation Algorithm for the Capacitated Arc Routing Problem

    DEFF Research Database (Denmark)

    Wøhlk, Sanne

    2008-01-01

    In this paper we consider approximation of the Capacitated Arc Routing Problem, which is the problem of servicing a set of edges in a graph using a fleet of capacity constrained vehicles. We present a 7/2 - 3/W-approximation algorithm for the problem and prove that this algorithm outperforms...

  20. Evaluation of the average ion approximation for a tokamak plasma

    International Nuclear Information System (INIS)

    The average ion approximation, sometimes used to calculated atomic processes in plasmas, is assessed by computing deviations in various rates over a set of conditions representative of tokamak edge plasmas. Conditions are identified under which the rates are primarily a function of the average ion charge and plasma parameters, as assumed in the average ion approximation. (Author) 19 refs., tab., 5 figs

  1. lambda-fuzzy approximate fixed point in fuzzy metric spaces

    OpenAIRE

    H. Mazaheri; S. A. M. Mohsenalhosseini

    2014-01-01

    In this paper $\\lambda$-fuzzy approximate fixed point for a map and $\\lambda$-fuzzy pair approximate fixed points for two maps in fuzzy metric spaces are defined. Instead of fuzzy numbers or real numbers are used to define fuzzy points. We also prove the existence theorems.

  2. Approximate Inertial Manifolds for Chemotaxis-Growth System

    Institute of Scientific and Technical Information of China (English)

    Hong LUO; Zhilin PU

    2012-01-01

    The long-time behaviour of solution to chemotaxis-growth system with Neumann condition is considered in this paper.The approximate inertial manifolds of such equations are constructed based on the contraction principle,and the orders of approximations of the manifolds to the global attractor are derived.

  3. Neutrino masses from an approximate mixing matrix with $\\theta_{13}\

    CERN Document Server

    Damanik, Asan

    2016-01-01

    An approximate neutrino mixing matrix is formutated by using the standard neutrino mixing matrix as a basis and experimental data of neutrino oscillations as inputs. By using the resulted approximate neutrino mixing matrix to proceed the neutrino mass matrix and constraining the resulted neutrino mass matrix with zero texture: $M_{\

  4. Fifth International Conference on "Approximation and Optimization in the Caribbean"

    CERN Document Server

    Approximation, Optimization and Mathematical Economic

    2001-01-01

    The articles in this proceedings volume reflect the current trends in the theory of approximation, optimization and mathematical economics, and include numerous applications. The book will be of interest to researchers and graduate students involved in functional analysis, approximation theory, mathematical programming and optimization, game theory, mathematical finance and economics.

  5. Accurate Period Approximation for Any Simple Pendulum Amplitude

    Institute of Scientific and Technical Information of China (English)

    XUE De-Sheng; ZHOU Zhao; GAO Mei-Zhen

    2012-01-01

    Accurate approximate analytical formulae of the pendulum period composed of a few elementary functions for any amplitude are constructed.Based on an approximation of the elliptic integral,two new logarithmic formulae for large amplitude close to 180° are obtained.Considering the trigonometric function modulation results from the dependence of relative error on the amplitude,we realize accurate approximation period expressions for any amplitude between 0 and 180°.A relative error less than 0.02% is achieved for any amplitude.This kind of modulation is also effective for other large-amplitude logarithmic approximation expressions.%Accurate approximate analytical formulae of the pendulum period composed of a few elementary functions for any amplitude are constructed. Based on an approximation of the elliptic integral, two new logarithmic formulae for large amplitude close to 180° are obtained. Considering the trigonometric function modulation results from the dependence of relative error on the amplitude, we realize accurate approximation period expressions for any amplitude between 0 and 180°. A relative error less than 0.02% is achieved for any amplitude. This kind of modulation is also effective for other large-amplitude logarithmic approximation expressions.

  6. Merging Belief Propagation and the Mean Field Approximation

    DEFF Research Database (Denmark)

    Riegler, Erwin; Kirkelund, Gunvor Elisabeth; Manchón, Carles Navarro;

    2010-01-01

    We present a joint message passing approach that combines belief propagation and the mean field approximation. Our analysis is based on the region-based free energy approximation method proposed by Yedidia et al., which allows to use the same objective function (Kullback-Leibler divergence) as a ...

  7. Quantum random walk approximation on locally compact quantum groups

    OpenAIRE

    Lindsay, J Martin; Skalski, Adam G.

    2011-01-01

    A natural scheme is established for the approximation of quantum Levy processes on locally compact quantum groups by quantum random walks. We work in the somewhat broader context of discrete approximations of completely positive quantum stochastic convolution cocycles on C*-bialgebras.

  8. Monte Carlo Greeks for financial products via approximative transition densities

    OpenAIRE

    Joerg Kampen; Anastasia Kolodko; John Schoenmakers

    2008-01-01

    In this paper we introduce efficient Monte Carlo estimators for the valuation of high-dimensional derivatives and their sensitivities (''Greeks''). These estimators are based on an analytical, usually approximative representation of the underlying density. We study approximative densities obtained by the WKB method. The results are applied in the context of a Libor market model.

  9. Approximation of Hardy space on the unit sphere

    Institute of Scientific and Technical Information of China (English)

    余纯武; 陈莘萌; 王昆扬; 戴峰

    2003-01-01

    The authors discuss the boundedness and approximation properties of translation and mean operator on H1(∑) by the estimates of high degree difference on ultraspherical polynomials, atom de-composition and construct properties on sphere. Also the boundedness and approximation of linear means at all kinds of indexes on Hp(0 < p < 1) and the almost everywhere convergence of Cesaro means are established.

  10. SOME GENERALIZATIONS OF KY FAN'S BEST APPROXIMATION THEOREM

    Institute of Scientific and Technical Information of China (English)

    A.R.Khan; N.Hussain; A.B.Thaheem

    2004-01-01

    We pressent new Ky Fan type best approximation theorems for a discontinuous multivalued map on metrizable topological vector spaces and hyperconvex spaces. In addition, fixed point results are derived for the map studied. Our work generalizes severl results in approximation theory.

  11. A MARKOVIAN APPROXIMATED SOLUTION TO A PORTFOLIO MANAGEMENT PROBLEM

    OpenAIRE

    Krawczyk, Jacek B.

    2000-01-01

    A portfolio management problem is approximated through a Markov decision chain. The weak Euler scheme is applied to discretise the time evolution of a portfolio and an inverse distance method is used to describe the transition probabilities. The approximating Markov decision problem is solved by value iteration. Numerical solutions of varying degrees of accuracy to a few specific portfolio problems are obtained.

  12. An Analysis of the Morris Loe Angle Trisection Approximation.

    Science.gov (United States)

    Aslan, Farhad,; And Others

    1992-01-01

    Presents the Morris Loe Angle Trisection Approximation Method to introduce students to areas of mathematics where approximations are used when exact answers are difficult or impossible to obtain. Examines the accuracy of the method using the laws of sines and cosines and a BASIC computer program that is provided. (MDH)

  13. The modified signed likelihood statistic and saddlepoint approximations

    DEFF Research Database (Denmark)

    Jensen, Jens Ledet

    1992-01-01

    SUMMARY: For a number of tests in exponential families we show that the use of a normal approximation to the modified signed likelihood ratio statistic r * is equivalent to the use of a saddlepoint approximation. This is also true in a large deviation region where the signed likelihood ratio...... statistic r is of order √ n. © 1992 Biometrika Trust....

  14. Neutral kaons beyond the Weisskopf-Wigner approximation

    International Nuclear Information System (INIS)

    A model-independent formalism for describing mixed states is constructed without recourse to the Weisskopf-Wigner approximation. Constraints imposed by symmetries are considered for K mesons. The formalism is used to show that effects associated with CPT violation and with going beyond the Weisskopf-Wigner approximation can be separated and studied independently

  15. Nucleon Properties from Approximating Chiral Quark Sigma Model

    CERN Document Server

    Abu-Shady, M

    2009-01-01

    We apply the approximating chiral quark model. This chiral quark model is based on an effective Lagrangian which the interactions between quarks via sigma and pions mesons. The field equations have been solved in the mean field approximation for the hedgehog baryon state. Good results are obtained for nucleon properties in comparison with original model.

  16. Examination of the adiabatic approximation in open systems

    International Nuclear Information System (INIS)

    We examine the notion of the adiabatic approximation in open systems by applying it to closed systems. Our results shows that the notion is equivalent to the standard adiabatic approximation if the systems are initially in eigenstates, and it leads to a more general expression if the systems are in mixed states

  17. Approximations in fusion and breakup reactions induced by radioactive beams

    Energy Technology Data Exchange (ETDEWEB)

    Cardenas, W.H.Z.; Carlin Filho, N.; Hussein, M.S. [Sao Paulo Univ., SP (Brazil). Inst. de Fisica; Canto, L.F.; Donangelo, R. [Universidade Federal, Rio de Janeiro, RJ (Brazil). Inst. de Fisica; Lubian, J. [Universidade Federal Fluminense, Niteroi, RJ (Brazil). Inst. de Fisica; Centro de Aplicaciones Tecnologicas y Desarrollo Nuclear (CEADEN), Havana (Cuba); Romanelli, A. [Facultad de Ingenieria, Montevideo (Uruguay). Inst. de Fisica

    2000-07-01

    Some commonly used approximations for complete fusion and breakup transmission coefficients in collisions of weakly bound projectiles at near barrier energies are assessed. We show that they strongly depend on the adopted classical trajectory and can be significantly improved with proper treatment of the incident and emergent currents in the WKB approximation. (author)

  18. p-optimal robust designs for multiresponse approximately linear regression

    Institute of Scientific and Technical Information of China (English)

    LIU Xin; YUE Rong-xian

    2008-01-01

    This paper deals with the problem of p-optimal robust designs for multiresponse approximately linear regression models. Each response is assumed to be only approximately linear in the regressors, and the bias function varies over a given ??2-neighbourhood. A kind of bivariate models with two responses is taken as an example to illustrate how to get the expression of the design measure.

  19. -Statistical Extension of the Korovkin Type Approximation Theorem

    Indian Academy of Sciences (India)

    Esra Erkus; Oktay Duman

    2005-11-01

    In this paper, using the concept of -statistical convergence which is a regular (non-matrix) summability method, we obtain a general Korovkin type approximation theorem which concerns the problem of approximating a function by means of a sequence $\\{L_nf\\}$ of positive linear operators.

  20. Transport coefficients in the P3 and P5 approximations

    International Nuclear Information System (INIS)

    Several procedures to account for space dependence of diffusion coefficients can be found in the literature. In this work, a relationship between neutron current and flux in the P3 and P5 approximations (generalized Fick's law) is obtained. Thermal neutron diffusion coefficient and other transport coefficients are found in P2 and P3 approximations. (author)

  1. Born-Oppenheimer approximation for a harmonic molecule

    OpenAIRE

    Fernandez, Francisco M.

    2008-01-01

    We apply the Born-Oppenheimer approximation to a harmonic diatomic molecule with one electron. We compare the exact and approximate results not only for the internal degrees of freedom but also for the motion of the center of mass. We address the problem of identical nuclei and discuss other applications of the model and its limitations.

  2. The use of neural networks for approximation of nuclear data

    Energy Technology Data Exchange (ETDEWEB)

    Korovin, Yu. A.; Maksimushkina, A. V., E-mail: AVMaksimushkina@mephi.ru [National Research Nuclear University MEPhI (Moscow Engineering Physics Institute) (Russian Federation)

    2015-12-15

    The article discusses the possibility of using neural networks for approximation or reconstruction of data such as the reaction cross sections. The quality of the approximation using fitting criteria is also evaluated. The activity of materials under irradiation is calculated from data obtained using neural networks.

  3. Approximations for strongly-coupled supersymmetric quantum mechanics

    CERN Document Server

    Kabat, D; Kabat, Daniel; Lifschytz, Gilad

    2000-01-01

    We advocate a set of approximations for studying the finite temperature behavior of strongly-coupled theories in 0+1 dimensions. The approximation consists of expanding about a Gaussian action, with the width of the Gaussian determined by a set of gap equations. The approximation can be applied to supersymmetric systems, provided that the gap equations are formulated in superspace. It can be applied to large-N theories, by keeping just the planar contribution to the gap equations. We analyze several models of scalar supersymmetric quantum mechanics, and show that the Gaussian approximation correctly distinguishes between a moduli space, mass gap, and supersymmetry breaking at strong coupling. Then we apply the approximation to a bosonic large-N gauge theory, and argue that a Gross-Witten transition separates the weak-coupling and strong-coupling regimes. A similar transition should occur in a generic large-N gauge theory, in particular in 0-brane quantum mechanics.

  4. Discretely sampled variance and volatility swaps versus their continuous approximations

    CERN Document Server

    Jarrow, Robert; Larsson, Martin; Protter, Philip

    2011-01-01

    Discretely sampled variance and volatility swaps trade actively in OTC markets. To price these swaps, the continuously sampled approximation is often used to simplify the computations. The purpose of this paper is to study the conditions under which this approximation is valid. Our first set of theorems characterize the conditions under which the discretely sampled swap values are finite, given the values of the continuous approximations exist. Surprisingly, for some otherwise reasonable price processes, the discretely sampled swap prices do not exist, thereby invalidating the approximation. Examples are provided. Assuming further that both swap values exist, we study sufficient conditions under which the discretely sampled values converge to their continuous counterparts. Because of its popularity in the literature, we apply our theorems to the 3/2 stochastic volatility model. Although we can show finiteness of all swap values, we can prove convergence of the approximation only for some parameter values.

  5. Approximate solutions of theWei Hua oscillator using the Pekeris approximation and Nikiforov–Uvarov method

    Indian Academy of Sciences (India)

    P K Bera

    2012-01-01

    The approximate analytical bound-state solutions of the Schrödinger equation for the Wei Hua oscillator are carried out in N-dimensional space by taking Pekeris approximation scheme to the orbital centrifugal term. Solutions of the corresponding hyper-radial equation are obtained using the conventional Nikiforov–Uvarov (NU) method.

  6. Frankenstein's glue: transition functions for approximate solutions

    Energy Technology Data Exchange (ETDEWEB)

    Yunes, Nicolas [Center for Gravitational Wave Physics, Institute for Gravitational Physics and Geometry, Department of Physics, Pennsylvania State University, University Park, PA 16802-6300 (United States)

    2007-09-07

    Approximations are commonly employed to find approximate solutions to the Einstein equations. These solutions, however, are usually only valid in some specific spacetime region. A global solution can be constructed by gluing approximate solutions together, but this procedure is difficult because discontinuities can arise, leading to large violations of the Einstein equations. In this paper, we provide an attempt to formalize this gluing scheme by studying transition functions that join approximate analytic solutions together. In particular, we propose certain sufficient conditions on these functions and prove that these conditions guarantee that the joined solution still satisfies the Einstein equations analytically to the same order as the approximate ones. An example is also provided for a binary system of non-spinning black holes, where the approximate solutions are taken to be given by a post-Newtonian expansion and a perturbed Schwarzschild solution. For this specific case, we show that if the transition functions satisfy the proposed conditions, then the joined solution does not contain any violations to the Einstein equations larger than those already inherent in the approximations. We further show that if these functions violate the proposed conditions, then the matter content of the spacetime is modified by the introduction of a matter shell, whose stress-energy tensor depends on derivatives of these functions.

  7. Approximate bayesian parameter inference for dynamical systems in systems biology

    International Nuclear Information System (INIS)

    This paper proposes to use approximate instead of exact stochastic simulation algorithms for approximate Bayesian parameter inference of dynamical systems in systems biology. It first presents the mathematical framework for the description of systems biology models, especially from the aspect of a stochastic formulation as opposed to deterministic model formulations based on the law of mass action. In contrast to maximum likelihood methods for parameter inference, approximate inference method- share presented which are based on sampling parameters from a known prior probability distribution, which gradually evolves toward a posterior distribution, through the comparison of simulated data from the model to a given data set of measurements. The paper then discusses the simulation process, where an over- view is given of the different exact and approximate methods for stochastic simulation and their improvements that we propose. The exact and approximate simulators are implemented and used within approximate Bayesian parameter inference methods. Our evaluation of these methods on two tasks of parameter estimation in two different models shows that equally good results are obtained much faster when using approximate simulation as compared to using exact simulation. (Author)

  8. Minimax rational approximation of the Fermi-Dirac distribution

    CERN Document Server

    Moussa, Jonathan E

    2016-01-01

    Accurate rational approximations of the Fermi-Dirac distribution are a useful component in many numerical algorithms for electronic structure calculations. The best known approximations use $O( \\log (\\beta \\Delta) \\log (\\epsilon^{-1}))$ poles to achieve an error tolerance $\\epsilon$ at temperature $\\beta^{-1}$ over an energy interval $\\Delta$. We apply minimax approximation to reduce the number of poles by a factor of four and replace $\\Delta$ with $\\Delta_{\\mathrm{occ}}$, the occupied energy interval. This is particularly beneficial when $\\Delta \\gg \\Delta_{\\mathrm{occ}}$, such as in electronic structure calculations that use a large basis set.

  9. Geometrical optics approximation of light scattering by large air bubbles

    Institute of Scientific and Technical Information of China (English)

    Haitao Yu; Jianqi Shen; Yuehuan Wei

    2008-01-01

    For large spherical bubbles in water,geometrical optics approximation is considered a better method for calculating light scattering patterns.In this paper,the basic theory of geometrical optics approximation is clarified.The change of phase for bubbles is calculated when total reflection occurs,which is different from particles with relative refractive indices larger than 1.Verification of the method was achieved by assuming a spherical particle and comparing present results to Mie scattering and Debye calculation.Agreement with the Mie theory was excellent in all directions when the dimensionless size parameter is larger than 50.Limitations of the geometrical optics approximation are also discussed.

  10. Approximation algorithms for Euler genus and related problems

    OpenAIRE

    Chekuri, Chandra; Sidiropoulos, Anastasios

    2013-01-01

    The Euler genus of a graph is a fundamental and well-studied parameter in graph theory and topology. Computing it has been shown to be NP-hard by [Thomassen '89 & '93], and it is known to be fixed-parameter tractable. However, the approximability of the Euler genus is wide open. While the existence of an O(1)-approximation is not ruled out, only an O(sqrt(n))-approximation [Chen, Kanchi, Kanevsky '97] is known even in bounded degree graphs. In this paper we give a polynomial-time algorithm wh...

  11. Weak Approximation of SDEs by Discrete-Time Processes

    Directory of Open Access Journals (Sweden)

    Henryk Zähle

    2008-01-01

    Full Text Available We consider the martingale problem related to the solution of an SDE on the line. It is shown that the solution of this martingale problem can be approximated by solutions of the corresponding time-discrete martingale problems under some conditions. This criterion is especially expedient for establishing the convergence of population processes to SDEs. We also show that the criterion yields a weak Euler scheme approximation of SDEs under fairly weak assumptions on the driving force of the approximating processes.

  12. Study on High-Speed Magnitude Approximation for Complex Vectors

    Institute of Scientific and Technical Information of China (English)

    陈建春; 杨万海; 许少英

    2003-01-01

    High-speed magnitude approximation algorithms for complex vectors are discussed intensively. The performance and the convergence speed of these approximation algorithms are analyzed. For the polygon fitting algorithms, the approximation formula under the least mean square error criterion is derived. For the iterative algorithms, a modified CORDIC (coordinate rotation digital computer) algorithm is developed. This modified CORDIC algorithm is proved to be with a maximum relative error about one half that of the original CORDIC algorithm. Finally, the effects of the finite register length on these algorithms are also concerned, which shows that 9 to 12-bit coefficients are sufficient for practical applications.

  13. An effective algorithm for approximating adaptive behavior in seasonal environments

    DEFF Research Database (Denmark)

    Sainmont, Julie; Andersen, Ken Haste; Thygesen, Uffe Høgsbro;

    2015-01-01

    behavior; the so-called "myopic approximation", "short sighted", or "static optimization". We explore the performance of the myopic approximation with diel vertical migration (DVM) as an example of a daily routine, a behavior with seasonal dependence that trades off predation risk with foraging......-annual variations, aspects that can only be accessed in dynamic programming approaches with escalating computational costs. Furthermore, the explanatory power of the myopic approximation is notably higher than when behavior is not implemented, highlighting the importance for adaptive DVM behavior in ecological...

  14. Revisiting the envelope approximation: gravitational waves from bubble collisions

    CERN Document Server

    Weir, David J

    2016-01-01

    We study the envelope approximation and its applicability to first-order phase transitions in the early universe. We demonstrate that the power laws seen in previous studies exist independent of the nucleation rate. We also compare the envelope approximation prediction to results from large-scale phase transition simulations. For phase transitions where the contribution to gravitational waves from scalar fields dominates over that from the coupled plasma of light particles, the envelope approximation is in agreement, giving a power spectrum of the same form and order of magnitude. In all other cases the form and amplitude of the gravitational wave power spectrum is markedly different and new techniques are required.

  15. Revisiting the envelope approximation: Gravitational waves from bubble collisions

    Science.gov (United States)

    Weir, David J.

    2016-06-01

    We study the envelope approximation and its applicability to first-order phase transitions in the early Universe. We demonstrate that the power laws seen in previous studies exist independently of the nucleation rate. We also compare the envelope approximation prediction to results from large-scale phase transition simulations. For phase transitions where the contribution to gravitational waves from scalar fields dominates over that from the coupled plasma of light particles, the envelope approximation is in agreement, giving a power spectrum of the same form and order of magnitude. In all other cases the form and amplitude of the gravitational wave power spectrum is markedly different and new techniques are required.

  16. On the Beebe-Linderberg two-electron integral approximation

    Science.gov (United States)

    Røeggen, I.; Wisløff-Nilssen, E.

    1986-12-01

    The Beebe-Linderberg two-electron integral approximation, which is generated by a Cholesky decomposition of the two-electron integral matrix ([μν|λσ]), is slightly modified. On the basis of test calculations, two key questions concerning this approximation are discussed: The numerical rank of the two-electron integral matrix and the relationship between the integral threshold and electronic properties. The numerical results presented in this work suggest that the modified Beebe-Linderberg approximation might be considered as an alternative to effective core potential methods.

  17. Some approximations in the linear dynamic equations of thin cylinders

    Science.gov (United States)

    El-Raheb, M.; Babcock, C. D., Jr.

    1981-01-01

    Theoretical analysis is performed on the linear dynamic equations of thin cylindrical shells to find the error committed by making the Donnell assumption and the neglect of in-plane inertia. At first, the effect of these approximations is studied on a shell with classical simply supported boundary condition. The same approximations are then investigated for other boundary conditions from a consistent approximate solution of the eigenvalue problem. The Donnell assumption is valid at frequencies high compared with the ring frequencies, for finite length thin shells. The error in the eigenfrequencies from omitting tangential inertia is appreciable for modes with large circumferential and axial wavelengths, independent of shell thickness and boundary conditions.

  18. Approximating Edit Distance in Near-Linear Time

    OpenAIRE

    Andoni, Alexandr; Onak, Krzysztof

    2011-01-01

    We show how to compute the edit distance between two strings of length n up to a factor of 2^{\\~O(sqrt(log n))} in n^(1+o(1)) time. This is the first sub-polynomial approximation algorithm for this problem that runs in near-linear time, improving on the state-of-the-art n^(1/3+o(1)) approximation. Previously, approximation of 2^{\\~O(sqrt(log n))} was known only for embedding edit distance into l_1, and it is not known if that embedding can be computed in less than quadratic time.

  19. Approximation to real numbers by cubic algebraic integers I

    OpenAIRE

    Roy, Damien

    2002-01-01

    In 1969, H. Davenport and W. M. Schmidt studied the problem of approximation to a real number \\xi by algebraic integers of degree at most three. They did so, using geometry of numbers, by resorting to the dual problem of finding simultaneous approximations to \\xi and \\xi^2 by rational numbers with the same denominator. In this paper, we show that their measure of approximation for the dual problem is optimal and that it is realized for a countable set of real numbers \\xi. We give several prop...

  20. WAVELET KERNEL SUPPORT VECTOR MACHINES FOR SPARSE APPROXIMATION

    Institute of Scientific and Technical Information of China (English)

    Tong Yubing; Yang Dongkai; Zhang Qishan

    2006-01-01

    Wavelet, a powerful tool for signal processing, can be used to approximate the target function. For enhancing the sparse property of wavelet approximation, a new algorithm was proposed by using wavelet kernel Support Vector Machines (SVM), which can converge to minimum error with better sparsity. Here, wavelet functions would be firstly used to construct the admitted kernel for SVM according to Mercy theory; then new SVM with this kernel can be used to approximate the target funciton with better sparsity than wavelet approxiamtion itself. The results obtained by our simulation experiment show the feasibility and validity of wavelet kernel support vector machines.

  1. SFU-Driven Transparent Approximation Acceleration on GPUs

    Energy Technology Data Exchange (ETDEWEB)

    Li, Ang; Song, Shuaiwen; Wijtvliet, Mark; Kumar, Akash; Corporaal, Henk

    2016-06-01

    Approximate computing, the technique that sacrifices certain amount of accuracy in exchange for substantial performance boost or power reduction, is one of the most promising solutions to enable power control and performance scaling towards exascale. Although most existing approximation designs target the emerging data-intensive applications that are comparatively more error-tolerable, there is still high demand for the acceleration of traditional scientific applications (e.g., weather and nuclear simulation), which often comprise intensive transcendental function calls and are very sensitive to accuracy loss. To address this challenge, we focus on a very important but often ignored approximation unit on GPUs.

  2. New approximation for the effective energy of nonlinear conducting composites

    Science.gov (United States)

    Gibiansky, Leonid; Torquato, Salvatore

    1998-07-01

    Approximations for the effective energy and, thus, effective conductivity of nonlinear, isotropic conducting dispersions are developed. This is accomplished by using the Ponte Castaneda variational principles [Philos. Trans. R. Soc. London Ser. A 340, 1321 (1992)] and the Torquato approximation [J. Appl. Phys. 58, 3790 (1985)] of the effective conductivity of corresponding linear composites. The results are obtained for dispersions with superconducting or insulating inclusions, and, more generally, for phases with a power-law energy. It is shown that the new approximations lie within the best available rigorous upper and lower bounds on the effective energy.

  3. Baby Skyrme model, near-BPS approximations, and supersymmetric extensions

    Science.gov (United States)

    Bolognesi, S.; Zakrzewski, W.

    2015-02-01

    We study the baby Skyrme model as a theory that interpolates between two distinct BPS systems. For this, a near-BPS approximation can be used when there is a small deviation from each of the two BPS limits. We provide analytical explanation and numerical support for the validity of this approximation. We then study the set of all possible supersymmetric extensions of the baby Skyrme model with N =1 and the particular ones with extended N =2 supersymmetries and relate this to the above mentioned almost-BPS approximation.

  4. Baby Skyrme Model, Near-BPS Approximations and Supersymmetric Extensions

    CERN Document Server

    Bolognesi, S

    2014-01-01

    We study the baby Skyrme model as a theory that interpolates between two distinct BPS systems. For this a near-BPS approximation can be used which, however, involves a small deviation from each of the two BPS limits. We provide analytical explanation and numerical support for the validity of this approximation. We then study the set of all possible supersymmetric extensions of the baby Skyrme model with ${\\cal N}=1$ and the particular ones with extended ${\\cal N}=2$ supersymmetries and relate this to the above mentioned almost-BPS approximation.

  5. Generalised quasi-linear approximation of the HMRI

    CERN Document Server

    Child, Adam; Marston, Brad; Tobias, Steven

    2016-01-01

    Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in descri...

  6. Preservation of Linear Constraints in Approximation of Tensors

    Institute of Scientific and Technical Information of China (English)

    Eugene Tyrtyshnikov

    2009-01-01

    For an arbitrary tensor (multi-index array) with linear constraints at each direction, it is proved that the factors of any minimal canonical tensor approximation to this tensor satisfy the same linear constraints for the corresponding directions.

  7. Kullback-Leibler divergence and the Pareto-Exponential approximation.

    Science.gov (United States)

    Weinberg, G V

    2016-01-01

    Recent radar research interests in the Pareto distribution as a model for X-band maritime surveillance radar clutter returns have resulted in analysis of the asymptotic behaviour of this clutter model. In particular, it is of interest to understand when the Pareto distribution is well approximated by an Exponential distribution. The justification for this is that under the latter clutter model assumption, simpler radar detection schemes can be applied. An information theory approach is introduced to investigate the Pareto-Exponential approximation. By analysing the Kullback-Leibler divergence between the two distributions it is possible to not only assess when the approximation is valid, but to determine, for a given Pareto model, the optimal Exponential approximation.

  8. Reply to Steele & Ferrer : Modeling Oscillation, Approximately or Exactly?

    NARCIS (Netherlands)

    Oud, Johan H. L.; Folmer, Henk

    2011-01-01

    This article addresses modeling oscillation in continuous time. It criticizes Steele and Ferrer's article "Latent Differential Equation Modeling of Self-Regulatory and Coregulatory Affective Processes" (2011), particularly the approximate estimation procedure applied. This procedure is the latent ve

  9. Approximate inverse preconditioning of iterative methods for nonsymmetric linear systems

    Energy Technology Data Exchange (ETDEWEB)

    Benzi, M. [Universita di Bologna (Italy); Tuma, M. [Inst. of Computer Sciences, Prague (Czech Republic)

    1996-12-31

    A method for computing an incomplete factorization of the inverse of a nonsymmetric matrix A is presented. The resulting factorized sparse approximate inverse is used as a preconditioner in the iterative solution of Ax = b by Krylov subspace methods.

  10. Linear approximation of the first eigenvalue on compact manifolds

    Institute of Scientific and Technical Information of China (English)

    CHEN; Mufa(陈木法); E.; Scacciatelli; YAO; Liang(姚亮)

    2002-01-01

    For compact, connected Riemannian manifolds with Ricci curvature bounded below by a constant, what is the linear approximation of the first eigenvalue of Laplacian? The answer is presented with computer assisted proof and the result is optimal in certain sense.

  11. Inner approximations for polynomial matrix inequalities and robust stability regions

    CERN Document Server

    Henrion, Didier

    2011-01-01

    Following a polynomial approach, many robust fixed-order controller design problems can be formulated as optimization problems whose set of feasible solutions is modelled by parametrized polynomial matrix inequalities (PMI). These feasibility sets are typically nonconvex. Given a parametrized PMI set, we provide a hierarchy of linear matrix inequality (LMI) problems whose optimal solutions generate inner approximations modelled by a single polynomial sublevel set. Those inner approximations converge in a strong analytic sense to the nonconvex original feasible set, with asymptotically vanishing conservatism. One may also impose the hierarchy of inner approximations to be nested or convex. In the latter case they do not converge any more to the feasible set, but they can be used in a convex optimization framework at the price of some conservatism. Finally, we show that the specific geometry of nonconvex polynomial stability regions can be exploited to improve convergence of the hierarchy of inner approximation...

  12. Numerical Approximation of Asymptotically Disappearing Solutions of Maxwell's Equations

    CERN Document Server

    Adler, J H; Zikatanov, L T

    2012-01-01

    This work is on the numerical approximation of incoming solutions to Maxwell's equations with dissipative boundary conditions whose energy decays exponentially with time. Such solutions are called asymptotically disappearing (ADS) and they play an importarnt role in inverse back-scatering problems. The existence of ADS is a difficult mathematical problem. For the exterior of a sphere, such solutions have been constructed analytically by Colombini, Petkov and Rauch [7] by specifying appropriate initial conditions. However, for general domains of practical interest (such as Lipschitz polyhedra), the existence of such solutions is not evident. This paper considers a finite-element approximation of Maxwell's equations in the exterior of a polyhedron, whose boundary approximates the sphere. Standard Nedelec-Raviart-Thomas elements are used with a Crank-Nicholson scheme to approximate the electric and magnetic fields. Discrete initial conditions interpolating the ones chosen in [7] are modified so that they are (we...

  13. APPROXIMATION BY GENERALIZED MKZ-OPERATORS IN POLYNOMIAL WEIGHTED SPACES

    Institute of Scientific and Technical Information of China (English)

    Lucyna Rempulska; Mariola Skorupka

    2007-01-01

    We prove some approximation properties of generalized Meyer-K(o)nig and Zeller operators for differentiable functions in polynomial weighted spaces. The results extend some results proved in [ 1-3,7-16].

  14. An approximate procedure for solving base-isolated structures

    International Nuclear Information System (INIS)

    Dynamic analysis of several shear type structures with base isolation indicates that the response of these structures follows their fundamental mode shape. Based on this observation, this paper uses an approximate procedure for computing the response of base-isolated structures. The procedure consists of modeling the structure and its base by a two-degree of freedom system, one representing the base and the other the structure. The response from the two-degree of freedom model and mode shapes of the structure are used to compute the response of the structure to earthquake excitation. The approximate procedure is simple, requires substantially less computational time than other methods, and gives results that are in excellent agreement with those from direct integration. Nonlinear properties and nonproportional damping are easily included in the model. Savings of approximately 54--77 percent in computational time result by using the approximate model

  15. Born approximation in linear-time invariant system

    CERN Document Server

    Gumjudpai, Burin

    2015-01-01

    Linear-time invariant (LTI) oscillation systems such as forced mechanical vibration, series RLC and parallel RLC circuits can be solved by using simplest initial conditions or employing of Green's function of which knowledge of initial condition of the force term is needed. Here we show a mathematical connection of the LTI system and the Helmholtz equation form of the time-independent Schr\\"{o}dinger equation in quantum mechanical scattering problem. We apply Born approximation in quantum mechanics to obtain LTI general solution in form of infinite Born series which can be expressed as a series of one-dimensional Feynman graphs. Conditions corresponding to the approximation are given for the case of harmonic driving force. The Born series of the harmonic forced oscillation case are derived by directly applying the approximation to the LTI system or by transforming the LTI system to Helmholtz equation prior to doing the approximation.

  16. The Sobolev approximation for line formation with partial frequency redistribution

    Science.gov (United States)

    Hummer, D. G.; Rybicki, G. B.

    1992-01-01

    Attention is given to the formation of a spectral line in a uniformly expanding infinite medium in the Sobolev approximation, with emphasis on the various mechanisms for frequency redistribution. Numerical and analytic solutions of the transfer equation are presented of a number of redistribution functions and their approximations, including type I and type II partial redistribution, coherent scattering and complete redistribution, and the Fokker-Planck and uncorrelated approximation to the R sub II function. The solutions for the mean intensity are shown to depend very much on the type of redistribution mechanism, while for the frequency-weighted mean intensity, which enters the rate equations, this dependence is weak. It is inferred that use of Sobolev escape probabilities based on complete redistribution can be an adequate approximation for many calculations for which only the radiative excitation rates are needed.

  17. Approximate supernova remnant dynamics with cosmic ray production

    International Nuclear Information System (INIS)

    Supernova explosions are the most violent and energetic events in the galaxy and have long been considered probable sources of cosmic rays. Recent shock acceleration models treating the cosmic rays (CR's) as test particles nb a prescribed supernova remnant (SNR) evolution, indeed indicate an approximate power law momentum distribution f sub source (p) approximation p(-a) for the particles ultimately injected into the interstellar medium (ISM). This spectrum extends almost to the momentum p = 1 million GeV/c, where the break in the observed spectrum occurs. The calculated power law index approximately less than 4.2 agrees with that inferred for the galactic CR sources. The absolute CR intensity can however not be well determined in such a test particle approximation

  18. Static NLO susceptibilities testing approximation schemes against exact results

    CERN Document Server

    Del Freo, L; Painelli, A; Freo, Luca Del; Terenziani, Francesca; Painelli, Anna

    2001-01-01

    The reliability of the approximations commonly adopted in the calculation of static optical (hyper)polarizabilities is tested against exact results obtained for an interesting toy-model. The model accounts for the principal features of typical nonlinear organic materials with mobile electrons strongly coupled to molecular vibrations. The approximations introduced in sum over states and finite field schemes are analyzed in detail. Both the Born-Oppenheimer and the clamped nucleus approximations turn out to be safe for molecules, whereas for donor-acceptor charge transfer complexes deviations from adiabaticity are expected. In the regime of low vibrational frequency, static susceptibilities are strongly dominated by the successive derivatives of the potential energy and large vibrational contributions to hyperpolarizabilities are found. In this regime anharmonic corrections to hyperpolarizabilities are very large, and the harmonic approximation, exact for the linear polarizability, turns out totally inadequate ...

  19. Eigenvalue approximation from below using non-conforming finite elements

    Institute of Scientific and Technical Information of China (English)

    2010-01-01

    This is a survey article about using non-conforming finite elements in solving eigenvalue problems of elliptic operators,with emphasis on obtaining lower bounds. In addition,this article also contains some new materials for eigenvalue approximations of the Laplace operator,which include:1) the proof of the fact that the non-conforming Crouzeix-Raviart element approximates eigenvalues associated with smooth eigenfunctions from below;2) the proof of the fact that the non-conforming EQ rot1 element approximates eigenvalues from below on polygonal domains that can be decomposed into rectangular elements;3) the explanation of the phenomena that numerical eigenvalues λ 1,h and λ 3,h of the non-conforming Q rot1 element approximate the true eigenvalues from below for the L-shaped domain. Finally,we list several unsolved problems.

  20. ON COPOSITIVE APPROXIMATION IN SOME CLASSICAL SPACES OF SEQUENCES

    Institute of Scientific and Technical Information of China (English)

    Aref Kamal

    2003-01-01

    In this paper the author writes a simple characterization for the best copositive approximation in c; the space of convergent sequences, by elements of finite dimensional Chebyshev subspaces, and shows that it is unique.

  1. Optimal approximations of power-laws with exponentials

    OpenAIRE

    Thierry Bochud; Damien Challet

    2006-01-01

    We propose an explicit recursive method to approximate a power-law with a finite sum of weighted exponentials. Applications to moving averages with long memory are discussed in relationship with stochastic volatility models.

  2. Approximate controllability and regularity for nonlinear differential equations

    OpenAIRE

    Kim Jin-Ran; Jeong Jin-Mun; Ju Eun-Young

    2011-01-01

    Abstract In this article, we deal with the existence, uniqueness, and a variation of solutions of the nonlinear control system with nonlinear monotone hemicontinuous and coercive operator. Moreover, the approximate controllability for the given nonlinear control system is studied.

  3. Best Simultaneous Approximation to Totally Bounded Sequences in Banach Spaces

    Institute of Scientific and Technical Information of China (English)

    Xian Fa LUO; Chong LI; Genaro LOPEZ

    2008-01-01

    This paper is concerned with the problem of best weighted simultaneous approximations to totally bounded sequences in Banach spaces. Characterization results from convex sets in Banach spaces are established under the assumption that the Banach space is uniformly smooth.

  4. RESTRICTED NONLINEAR APPROXIMATION AND SINGULAR SOLUTIONS OF BOUNDARY INTEGRAL EQUATIONS

    Institute of Scientific and Technical Information of China (English)

    Reinhard Hochmuth

    2002-01-01

    This paper studies several problems, which are potentially relevant for the construction of adaptive numerical schemes. First, biorthogonal spline wavelets on [0,1 ] are chosen as a starting point for characterizations of functions in Besov spaces B , (0,1) with 0<σ<∞ and (1+σ)-1<τ<∞. Such function spaces are known to be related to nonlinear approximation. Then so called restricted nonlinear approximation procedures with respect to Sobolev space norms are considered. Besides characterization results Jackson type estimates for various tree-type and tresholding algorithms are investigated. Finally known approximation results for geometry induced singularity functions of boundary integeral equations are combined with the characterization results for restricted nonlinear approximation to show Besov space regularity results.

  5. Hamilton's Principle and Approximate Solutions to Problems in Classical Mechanics

    Science.gov (United States)

    Schlitt, D. W.

    1977-01-01

    Shows how to use the Ritz method for obtaining approximate solutions to problems expressed in variational form directly from the variational equation. Application of this method to classical mechanics is given. (MLH)

  6. Approximating the Stieltjes integral by using the generalized trapezoid rule

    Directory of Open Access Journals (Sweden)

    Mustafa Haluk Chelik

    2011-12-01

    Full Text Available Accurate approximations for the Stieltjes integral by the generalized trapzoid rule. The generalized trapezoid rule is established on the notion of the derivative of function with respect to an increasing function, defined in [9].

  7. Vapor Intrusion Facility Boundaries - Approximate Extent of Contamination

    Data.gov (United States)

    U.S. Environmental Protection Agency — Theses polygons represent the approximate extent of contamination for Vapor Intrusion facilities. Polygons for non-Federal Facility sites were updated in January,...

  8. Approximate Solutions of Klein-Gordon Equation with Kratzer Potential

    Directory of Open Access Journals (Sweden)

    H. Hassanabadi

    2011-01-01

    Full Text Available Approximate solutions of the D-dimensional Klein-Gordon equation are obtained for the scalar and vector general Kratzer potential for any l by using the ansatz method. The energy behavior is numerically discussed.

  9. Badly approximable systems of linear forms in absolute value

    DEFF Research Database (Denmark)

    Hussain, M.; Kristensen, Simon

    In this paper we show that the set of mixed type badly approximable simultaneously small linear forms is of maximal dimension. As a consequence of this theorem we settle the conjecture stated in [9]....

  10. Numerical Stability and Convergence of Approximate Methods for Conservation Laws

    Science.gov (United States)

    Galkin, V. A.

    We present the new approach to background of approximate methods convergence based on functional solutions theory for conservation laws. The applications to physical kinetics, gas and fluid dynamics are considered.

  11. Envelope induced ionization dynamics beyond the dipole approximation

    CERN Document Server

    Simonsen, Aleksander Skjerlie; Førre, Morten; Lindroth, Eva; Selstø, Sølve

    2015-01-01

    When atoms and molecules are ionized by laser pulses of finite duration and increasingly high intensities, the validity of the much used dipole approximation, in which the spatial dependence and magnetic component of the external field are neglected, eventually breaks down. We report that when going beyond the dipole approximation for the description of atoms exposed to ultraviolet light, the spatial dependence of the pulse shape, the envelope, provides the dominant correction, while the spatial dependence of the carrier may safely be neglected in the general case. We present a first order beyond-dipole correction to the Hamiltonian which accounts exclusively for effects stemming from the carrier-envelope of the pulse. This much simpler form of the correction is further discussed in connection with various descriptions of the light-matter interaction. We demonstrate by ab initio calculations that this approximation, which we will refer to as the envelope approximation, reproduces the full interaction beyond t...

  12. Techniques for correcting approximate finite difference solutions. [considering transonic flow

    Science.gov (United States)

    Nixon, D.

    1978-01-01

    A method of correcting finite-difference solutions for the effect of truncation error or the use of an approximate basic equation is presented. Applications to transonic flow problems are described and examples are given.

  13. Approximate Riemann solvers for the Godunov SPH (GSPH)

    Science.gov (United States)

    Puri, Kunal; Ramachandran, Prabhu

    2014-08-01

    The Godunov Smoothed Particle Hydrodynamics (GSPH) method is coupled with non-iterative, approximate Riemann solvers for solutions to the compressible Euler equations. The use of approximate solvers avoids the expensive solution of the non-linear Riemann problem for every interacting particle pair, as required by GSPH. In addition, we establish an equivalence between the dissipative terms of GSPH and the signal based SPH artificial viscosity, under the restriction of a class of approximate Riemann solvers. This equivalence is used to explain the anomalous “wall heating” experienced by GSPH and we provide some suggestions to overcome it. Numerical tests in one and two dimensions are used to validate the proposed Riemann solvers. A general SPH pairing instability is observed for two-dimensional problems when using unequal mass particles. In general, Ducowicz Roe's and HLLC approximate Riemann solvers are found to be suitable replacements for the iterative Riemann solver in the original GSPH scheme.

  14. Rotational and Helical Surface Approximation for Reverse Engineering

    DEFF Research Database (Denmark)

    Randrup, Thomas; Pottmann, Helmut

    1997-01-01

    to basic shapes used in computer aided design. The algorithms apply methods of line geometry to the set of surface normals in combination with techniques of numerical approximation. The presented results possess applications in reverse engineering and computer aided manufacturing....

  15. Parallel Approximation Algorithms for Facility-Location Problems

    CERN Document Server

    Blelloch, Guy E

    2010-01-01

    This paper presents the design and analysis of parallel approximation algorithms for facility-location problems, including $\\NC$ and $\\RNC$ algorithms for (metric) facility location, $k$-center, $k$-median, and $k$-means. These problems have received considerable attention during the past decades from the approximation algorithms community, concentrating primarily on improving the approximation guarantees. In this paper, we ask, is it possible to parallelize some of the beautiful results from the sequential setting? Our starting point is a small, but diverse, subset of results in approximation algorithms for facility-location problems, with a primary goal of developing techniques for devising their efficient parallel counterparts. We focus on giving algorithms with low depth, near work efficiency (compared to the sequential versions), and low cache complexity. Common in algorithms we present is the idea that instead of picking only the most cost-effective element, we make room for parallelism by allowing a sm...

  16. Approximate penetration factors for nuclear reactions of astrophysical interest

    Science.gov (United States)

    Humblet, J.; Fowler, W. A.; Zimmerman, B. A.

    1987-01-01

    The ranges of validity of approximations of P(l), the penetration factor which appears in the parameterization of nuclear-reaction cross sections at low energies and is employed in the extrapolation of laboratory data to even lower energies of astrophysical interest, are investigated analytically. Consideration is given to the WKB approximation, P(l) at the energy of the total barrier, approximations derived from the asymptotic expansion of G(l) for large eta, approximations for small values of the parameter x, applications of P(l) to nuclear reactions, and the dependence of P(l) on channel radius. Numerical results are presented in tables and graphs, and parameter ranges where the danger of serious errors is high are identified.

  17. Approximate equations at breaking for nearshore wave transformation coefficients

    Digital Repository Service at National Institute of Oceanography (India)

    Chandramohan, P.; Nayak, B.U.; SanilKumar, V.

    Based on small amplitude wave theory approximate equations are evaluated for determining the coefficients of shoaling, refraction, bottom friction, bottom percolation and viscous dissipation at breaking. The results obtainEd. by these equations...

  18. A coefficient average approximation towards Gutzwiller wavefunction formalism

    Science.gov (United States)

    Liu, Jun; Yao, Yongxin; Wang, Cai-Zhuang; Ho, Kai-Ming

    2015-06-01

    Gutzwiller wavefunction is a physically well-motivated trial wavefunction for describing correlated electron systems. In this work, a new approximation is introduced to facilitate the evaluation of the expectation value of any operator within the Gutzwiller wavefunction formalism. The basic idea is to make use of a specially designed average over Gutzwiller wavefunction coefficients expanded in the many-body Fock space to approximate the ratio of expectation values between a Gutzwiller wavefunction and its underlying noninteracting wavefunction. To check with the standard Gutzwiller approximation (GA), we test its performance on single band systems and find quite interesting properties. On finite systems, we noticed that it gives superior performance over GA, while on infinite systems it asymptotically approaches GA. Analytic analysis together with numerical tests are provided to support this claimed asymptotical behavior. Finally, possible improvements on the approximation and its generalization towards multiband systems are illustrated and discussed.

  19. Isoscalar Giant Monopole Resonance in Relativistic Continuum Random Phase Approximation

    Institute of Scientific and Technical Information of China (English)

    2008-01-01

    <正>The isoscalar giant monopole resonance (ISGMR) in nuclei is studied in the framework of a fully consistent relativistic continuum random phase approximation (RCRPA). In this method the contribution

  20. An approximation for homogeneous freezing temperature of water droplets

    Science.gov (United States)

    O, K.-T.; Wood, R.

    2015-11-01

    In this work, based on the well-known formulae of classical nucleation theory (CNT), the temperature TNc = 1 at which the mean number of critical embryos inside a droplet is unity is derived and proposed as a new approximation for homogeneous freezing temperature of water droplets. Without consideration of time dependence and stochastic nature of the ice nucleation process, the approximation TNc = 1 is able to reproduce the dependence of homogeneous freezing temperature on drop size and water activity of aqueous drops observed in a wide range of experimental studies. We use the TNc = 1 approximation to argue that the distribution of homogeneous freezing temperatures observed in the experiments may largely be explained by the spread in the size distribution of droplets used in the particular experiment. It thus appears that this approximation is useful for predicting homogeneous freezing temperatures of water droplets in the atmosphere.