Milgram, A
2011-02-21
This comment addresses critics on the claimed stability of solution to the accelerated-predator-satiety Lotka-Volterra predator-prey problem, proposed by Dubey al. (2010. A solution to the accelerated-predator-satiety Lotka-Volterra predator-prey problem using Boubaker polynomial expansion scheme. Journal of Theoretical Biology 264, 154-160). Critics are based on incompatibilities between the claimed asymptotic behavior and the presumed Malthusian growth of prey population in absence of predator.
Directory of Open Access Journals (Sweden)
B. Kafash
2014-04-01
Full Text Available In this paper, we present a computational method for solving optimal control problems and the controlled Duffing oscillator. This method is based on state parametrization. In fact, the state variable is approximated by Boubaker polynomials with unknown coefficients. The equation of motion, performance index and boundary conditions are converted into some algebraic equations. Thus, an optimal control problem converts to a optimization problem, which can then be solved easily. By this method, the numerical value of the performance index is obtained. Also, the control and state variables can be approximated as functions of time. Convergence of the algorithms is proved. Numerical results are given for several test examples to demonstrate the applicability and efficiency of the method.
Directory of Open Access Journals (Sweden)
Luigi Vecchione
2015-07-01
Full Text Available One of the most important issues in biomass biocatalytic gasification is the correct prediction of gasification products, with particular attention to the Topping Atmosphere Residues (TARs. In this work, performedwithin the European 7FP UNIfHY project, we develops and validate experimentally a model which is able of predicting the outputs,including TARs, of a steam-fluidized bed biomass gasifier. Pine wood was chosen as biomass feedstock: the products obtained in pyrolysis tests are the relevant model input. Hydrodynamics and chemical properties of the reacting system are considered: the hydrodynamic approach is based on the two phase theory of fluidization, meanwhile the chemical model is based on the kinetic equations for the heterogeneous and homogenous reactions. The derived differentials equations for the gasifier at steady state were implemented MATLAB. Solution was consecutively carried out using the Boubaker Polynomials Expansion Scheme by varying steam/biomass ratio (0.5-1 and operating temperature (750-850°C.The comparison between models and experimental results showed that the model is able of predicting gas mole fractions and production rate including most of the representative TARs compounds
Fridjine, S.; Amlouk, M.
In this study, we define a synthetic parameter: optothermal expansivity as a quantitative guide to evaluating and optimizing both the thermal and the optical performance of PV-T functional materials. The definition of this parameter, ψAB (Amlouk-Boubaker parameter), takes into account the thermal diffusivity and the optical effective absorptivity of the material. The values of this parameter, which seems to be a characteristic one, correspond to the total volume that contains a fixed amount of heat per unit time (m3 s-1) and can be considered as a 3D velocity of the transmitted heat inside the material. As the PV-T combined devices need to have simultaneous optical and thermal efficiency, we try to investigate some recently proposed materials (β-SnS2, In2S3, ZnS1-xSex|0 ≤xabacus.
Energy Technology Data Exchange (ETDEWEB)
Boubaker, K.; Amlouk, M. [Unite de Physique de Dispositifs a Semiconducteurs-UPDS-Faculte des Sciences de Tunis, Campus Universitaire, 2092 Tunis (Tunisia)
2010-10-15
This study is a prelude to the definition of a new synthetic parameter inserted in a 2D abacus. This parameter: the Amlouk-Boubaker optothermal expansivity <{psi}{sub AB}>, is defined, for a given PV-T material, as a thermal diffusivity-to-optical effective absorptivity ratio. This parameter's unit evokes a heat flow velocity inside the material. Consequently, the parameter {psi}{sub AB} could be combined with the already known bandgap energy E{sub g}, in order to establish a 2D abacus. A sketched scheme of the 2D abacus is proposed as a guide for investigation and evaluation of PV-T candidate materials like metal oxides, amorphous silicon, zinc-doped binary compounds, and hydrogenated amorphous carbon. Using this abacus, designers will be able to compare solar energy-related materials on the basis of conjoint optical and thermal efficiency. (author)
ASYMPTOTIC EXPANSIONS OF ZEROS FOR KRAWTCHOUK POLYNOMIALS WITH ERROR BOUNDS
Institute of Scientific and Technical Information of China (English)
ZHU Xiao-feng; LI Xiu-chun
2006-01-01
Krawtchouk polynomials are frequently applied in modern physics. Based on the results which were educed by Li and Wong, the asymptotic expansions of Krawtchouk polynomials are improved by using Airy function, and uniform asymptotic expansions are got. Furthermore, the asymptotic expansions of the zeros for Krawtchouk polynomials are again deduced by using the property of the zeros of Airy function, and their corresponding error bounds are discussed. The obtained results give the asymptotic property of Krawtchouk polynomials with their zeros, which are better than the results educed by Li and Wong.
Asymptotic expansions of Feynman integrals of exponentials with polynomial exponent
Kravtseva, A. K.; Smolyanov, O. G.; Shavgulidze, E. T.
2016-10-01
In the paper, an asymptotic expansion of path integrals of functionals having exponential form with polynomials in the exponent is constructed. The definition of the path integral in the sense of analytic continuation is considered.
Dubey, B; Zhao, T G; Jonsson, M; Rahmanov, H
2010-05-01
In this study, an analytical method is introduced for the identification of predator-prey populations time-dependent evolution in a Lotka-Volterra predator-prey model which takes into account the concept of accelerated-predator-satiety. Oppositely to most of the predator-prey problem models, the actual model does not suppose that the predation is strictly proportional to the prey density. In reference to some recent experimental results and particularly to the conclusions of May (1973) about predators which are 'never not hungry', an accelerated satiety function is matched with the initial conventional equations. Solutions are plotted and compared to some relevant ones. The obtained trends are in good agreement with many standard Lotka-Volterra solutions except for the asymptotic behaviour.
A GLOBALLY UNIFORM ASYMPTOTIC EXPANSION OF THE HERMITE POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
Shi Wei
2008-01-01
In this article, the author extends the validity of a uniform asymptotic ex-pansion of the Hermite polynomials HN(√2n+1α) to include all positive values of a.His method makes use of the rational functions introduced by Olde Daalhuis and Temme (SIAM J. Math. Anal., (1994), 25: 304-321). A new estimate for the remainder is given.
Global Monte Carlo Simulation with High Order Polynomial Expansions
Energy Technology Data Exchange (ETDEWEB)
William R. Martin; James Paul Holloway; Kaushik Banerjee; Jesse Cheatham; Jeremy Conlin
2007-12-13
The functional expansion technique (FET) was recently developed for Monte Carlo simulation. The basic idea of the FET is to expand a Monte Carlo tally in terms of a high order expansion, the coefficients of which can be estimated via the usual random walk process in a conventional Monte Carlo code. If the expansion basis is chosen carefully, the lowest order coefficient is simply the conventional histogram tally, corresponding to a flat mode. This research project studied the applicability of using the FET to estimate the fission source, from which fission sites can be sampled for the next generation. The idea is that individual fission sites contribute to expansion modes that may span the geometry being considered, possibly increasing the communication across a loosely coupled system and thereby improving convergence over the conventional fission bank approach used in most production Monte Carlo codes. The project examined a number of basis functions, including global Legendre polynomials as well as “local” piecewise polynomials such as finite element hat functions and higher order versions. The global FET showed an improvement in convergence over the conventional fission bank approach. The local FET methods showed some advantages versus global polynomials in handling geometries with discontinuous material properties. The conventional finite element hat functions had the disadvantage that the expansion coefficients could not be estimated directly but had to be obtained by solving a linear system whose matrix elements were estimated. An alternative fission matrix-based response matrix algorithm was formulated. Studies were made of two alternative applications of the FET, one based on the kernel density estimator and one based on Arnoldi’s method of minimized iterations. Preliminary results for both methods indicate improvements in fission source convergence. These developments indicate that the FET has promise for speeding up Monte Carlo fission source
Enhancing sparsity of Hermite polynomial expansions by iterative rotations
Energy Technology Data Exchange (ETDEWEB)
Yang, Xiu; Lei, Huan; Baker, Nathan A.; Lin, Guang
2016-02-01
Compressive sensing has become a powerful addition to uncertainty quantification in recent years. This paper identifies new bases for random variables through linear mappings such that the representation of the quantity of interest is more sparse with new basis functions associated with the new random variables. This sparsity increases both the efficiency and accuracy of the compressive sensing-based uncertainty quantification method. Specifically, we consider rotation- based linear mappings which are determined iteratively for Hermite polynomial expansions. We demonstrate the effectiveness of the new method with applications in solving stochastic partial differential equations and high-dimensional (O(100)) problems.
Uncertainty Quantification for Polynomial Systems via Bernstein Expansions
Crespo, Luis G.; Kenny, Sean P.; Giesy, Daniel P.
2012-01-01
This paper presents a unifying framework to uncertainty quantification for systems having polynomial response metrics that depend on both aleatory and epistemic uncertainties. The approach proposed, which is based on the Bernstein expansions of polynomials, enables bounding the range of moments and failure probabilities of response metrics as well as finding supersets of the extreme epistemic realizations where the limits of such ranges occur. These bounds and supersets, whose analytical structure renders them free of approximation error, can be made arbitrarily tight with additional computational effort. Furthermore, this framework enables determining the importance of particular uncertain parameters according to the extent to which they affect the first two moments of response metrics and failure probabilities. This analysis enables determining the parameters that should be considered uncertain as well as those that can be assumed to be constants without incurring significant error. The analytical nature of the approach eliminates the numerical error that characterizes the sampling-based techniques commonly used to propagate aleatory uncertainties as well as the possibility of under predicting the range of the statistic of interest that may result from searching for the best- and worstcase epistemic values via nonlinear optimization or sampling.
Polynomial Chaos Expansion Approach to Interest Rate Models
Directory of Open Access Journals (Sweden)
Luca Di Persio
2015-01-01
Full Text Available The Polynomial Chaos Expansion (PCE technique allows us to recover a finite second-order random variable exploiting suitable linear combinations of orthogonal polynomials which are functions of a given stochastic quantity ξ, hence acting as a kind of random basis. The PCE methodology has been developed as a mathematically rigorous Uncertainty Quantification (UQ method which aims at providing reliable numerical estimates for some uncertain physical quantities defining the dynamic of certain engineering models and their related simulations. In the present paper, we use the PCE approach in order to analyze some equity and interest rate models. In particular, we take into consideration those models which are based on, for example, the Geometric Brownian Motion, the Vasicek model, and the CIR model. We present theoretical as well as related concrete numerical approximation results considering, without loss of generality, the one-dimensional case. We also provide both an efficiency study and an accuracy study of our approach by comparing its outputs with the ones obtained adopting the Monte Carlo approach, both in its standard and its enhanced version.
Direct discriminant locality preserving projection with Hammerstein polynomial expansion.
Chen, Xi; Zhang, Jiashu; Li, Defang
2012-12-01
Discriminant locality preserving projection (DLPP) is a linear approach that encodes discriminant information into the objective of locality preserving projection and improves its classification ability. To enhance the nonlinear description ability of DLPP, we can optimize the objective function of DLPP in reproducing kernel Hilbert space to form a kernel-based discriminant locality preserving projection (KDLPP). However, KDLPP suffers the following problems: 1) larger computational burden; 2) no explicit mapping functions in KDLPP, which results in more computational burden when projecting a new sample into the low-dimensional subspace; and 3) KDLPP cannot obtain optimal discriminant vectors, which exceedingly optimize the objective of DLPP. To overcome the weaknesses of KDLPP, in this paper, a direct discriminant locality preserving projection with Hammerstein polynomial expansion (HPDDLPP) is proposed. The proposed HPDDLPP directly implements the objective of DLPP in high-dimensional second-order Hammerstein polynomial space without matrix inverse, which extracts the optimal discriminant vectors for DLPP without larger computational burden. Compared with some other related classical methods, experimental results for face and palmprint recognition problems indicate the effectiveness of the proposed HPDDLPP.
Luo, Qiu-Ming
2009-12-01
We investigate Fourier expansions for the Apostol-Bernoulli and Apostol-Euler polynomials using the Lipschitz summation formula and obtain their integral representations. We give some explicit formulas at rational arguments for these polynomials in terms of the Hurwitz zeta function. We also derive the integral representations for the classical Bernoulli and Euler polynomials and related known results.
Expansions of one density via polynomials orthogonal with respect to the other
Szabłowski, Paweł J
2010-01-01
We expand Chebyshev polynomials and some of its linear combination in linear combinations of q-Hermite, Rogers and Al Salam-Chihara polynomials and vice versa. We use these expansions to obtain expansions of the some densities, including q-Normal and some related to it, in infinite series of orthogonal polynomials allowing deeper analysis, discovering new properties. On the way we find an easy proof of expansion of of Poisson-Mehler kernels for q-Hermite polynomials and also its inverse. We also formulate simple rule relating one set of orthogonal polynomials to the other given the properties of the ratio of the respective densities of measures orthogonalizing these polynomials sets.
Uncertainty Quantification for Airfoil Icing using Polynomial Chaos Expansions
DeGennaro, Anthony M; Martinelli, Luigi
2014-01-01
The formation and accretion of ice on the leading edge of a wing can be detrimental to airplane performance. Complicating this reality is the fact that even a small amount of uncertainty in the shape of the accreted ice may result in a large amount of uncertainty in aerodynamic performance metrics (e.g., stall angle of attack). The main focus of this work concerns using the techniques of Polynomial Chaos Expansions (PCE) to quantify icing uncertainty much more quickly than traditional methods (e.g., Monte Carlo). First, we present a brief survey of the literature concerning the physics of wing icing, with the intention of giving a certain amount of intuition for the physical process. Next, we give a brief overview of the background theory of PCE. Finally, we compare the results of Monte Carlo simulations to PCE-based uncertainty quantification for several different airfoil icing scenarios. The results are in good agreement and confirm that PCE methods are much more efficient for the canonical airfoil icing un...
On Continued Fraction Expansion of Real Roots of Polynomial Systems
DEFF Research Database (Denmark)
Mantzaflaris, Angelos; Mourrain, Bernard; Tsigaridas, Elias
2011-01-01
We elaborate on a correspondence between the coefficients of a multivariate polynomial represented in the Bernstein basis and in a tensor-monomial basis, which leads to homography representations of polynomial functions that use only integer arithmetic (in contrast to the Bernstein basis) and are...
Koçak, H.; Dahong, Z.; Yildirim, A.
2011-05-01
In this study, a range-free method is proposed in order to determine the Antoine constants for a given material (salicylic acid). The advantage of this method is mainly yielding analytical expressions which fit different temperature ranges.
Konakli, Katerina; Sudret, Bruno
2015-01-01
The growing need for uncertainty analysis of complex computational models has led to an expanding use of meta-models across engineering and sciences. The efficiency of meta-modeling techniques relies on their ability to provide statistically-equivalent analytical representations based on relatively few evaluations of the original model. Polynomial chaos expansions (PCE) have proven a powerful tool for developing meta-models in a wide range of applications; the key idea thereof is to expand th...
Simulation of stochastic systems via polynomial chaos expansions and convex optimization
Fagiano, Lorenzo
2012-01-01
Polynomial Chaos Expansions represent a powerful tool to simulate stochastic models of dynamical systems. Yet, deriving the expansion's coefficients for complex systems might require a significant and non-trivial manipulation of the model, or the computation of large numbers of simulation runs, rendering the approach too time consuming and impracticable for applications with more than a handful of random variables. We introduce a novel computationally tractable technique for computing the coefficients of polynomial chaos expansions. The approach exploits a regularization technique with a particular choice of weighting matrices, which allow to take into account the specific features of Polynomial Chaos expansions. The method, completely based on convex optimization, can be applied to problems with a large number of random variables and uses a modest number of Monte Carlo simulations, while avoiding model manipulations. Additional information on the stochastic process, when available, can be also incorporated i...
Energy Technology Data Exchange (ETDEWEB)
Konakli, Katerina, E-mail: konakli@ibk.baug.ethz.ch; Sudret, Bruno
2016-09-15
The growing need for uncertainty analysis of complex computational models has led to an expanding use of meta-models across engineering and sciences. The efficiency of meta-modeling techniques relies on their ability to provide statistically-equivalent analytical representations based on relatively few evaluations of the original model. Polynomial chaos expansions (PCE) have proven a powerful tool for developing meta-models in a wide range of applications; the key idea thereof is to expand the model response onto a basis made of multivariate polynomials obtained as tensor products of appropriate univariate polynomials. The classical PCE approach nevertheless faces the “curse of dimensionality”, namely the exponential increase of the basis size with increasing input dimension. To address this limitation, the sparse PCE technique has been proposed, in which the expansion is carried out on only a few relevant basis terms that are automatically selected by a suitable algorithm. An alternative for developing meta-models with polynomial functions in high-dimensional problems is offered by the newly emerged low-rank approximations (LRA) approach. By exploiting the tensor–product structure of the multivariate basis, LRA can provide polynomial representations in highly compressed formats. Through extensive numerical investigations, we herein first shed light on issues relating to the construction of canonical LRA with a particular greedy algorithm involving a sequential updating of the polynomial coefficients along separate dimensions. Specifically, we examine the selection of optimal rank, stopping criteria in the updating of the polynomial coefficients and error estimation. In the sequel, we confront canonical LRA to sparse PCE in structural-mechanics and heat-conduction applications based on finite-element solutions. Canonical LRA exhibit smaller errors than sparse PCE in cases when the number of available model evaluations is small with respect to the input
Konakli, Katerina; Sudret, Bruno
2016-09-01
The growing need for uncertainty analysis of complex computational models has led to an expanding use of meta-models across engineering and sciences. The efficiency of meta-modeling techniques relies on their ability to provide statistically-equivalent analytical representations based on relatively few evaluations of the original model. Polynomial chaos expansions (PCE) have proven a powerful tool for developing meta-models in a wide range of applications; the key idea thereof is to expand the model response onto a basis made of multivariate polynomials obtained as tensor products of appropriate univariate polynomials. The classical PCE approach nevertheless faces the "curse of dimensionality", namely the exponential increase of the basis size with increasing input dimension. To address this limitation, the sparse PCE technique has been proposed, in which the expansion is carried out on only a few relevant basis terms that are automatically selected by a suitable algorithm. An alternative for developing meta-models with polynomial functions in high-dimensional problems is offered by the newly emerged low-rank approximations (LRA) approach. By exploiting the tensor-product structure of the multivariate basis, LRA can provide polynomial representations in highly compressed formats. Through extensive numerical investigations, we herein first shed light on issues relating to the construction of canonical LRA with a particular greedy algorithm involving a sequential updating of the polynomial coefficients along separate dimensions. Specifically, we examine the selection of optimal rank, stopping criteria in the updating of the polynomial coefficients and error estimation. In the sequel, we confront canonical LRA to sparse PCE in structural-mechanics and heat-conduction applications based on finite-element solutions. Canonical LRA exhibit smaller errors than sparse PCE in cases when the number of available model evaluations is small with respect to the input dimension, a
Jacquelin, E.; Adhikari, S.; Sinou, J.-J.; Friswell, M. I.
2015-11-01
Polynomial chaos solution for the frequency response of linear non-proportionally damped dynamic systems has been considered. It has been observed that for lightly damped systems the convergence of the solution can be very poor in the vicinity of the deterministic resonance frequencies. To address this, Aitken's transformation and its generalizations are suggested. The proposed approach is successfully applied to the sequences defined by the first two moments of the responses, and this process significantly accelerates the polynomial chaos convergence. In particular, a 2-dof system with respectively 1 and 2 parameter uncertainties has been studied. The first two moments of the frequency response were calculated by Monte Carlo simulation, polynomial chaos expansion and Aitken's transformation of the polynomial chaos expansion. Whereas 200 polynomials are required to have a good agreement with Monte Carlo results around the deterministic eigenfrequencies, less than 50 polynomials transformed by the Aitken's method are enough. This latter result is improved if a generalization of Aitken's method (recursive Aitken's transformation, Shank's transformation) is applied. With the proposed convergence acceleration, polynomial chaos may be reconsidered as an efficient method to estimate the first two moments of a random dynamic response.
Apostol, Tom M. (Editor)
1991-01-01
In this 'Project Mathematics! series, sponsored by California Institute for Technology (CalTech), the mathematical concept of polynomials in rectangular coordinate (x, y) systems are explored. sing film footage of real life applications and computer animation sequences, the history of, the application of, and the different linear coordinate systems for quadratic, cubic, intersecting, and higher degree of polynomials are discussed.
Computation of Higher-Order Moments of Generalized Polynomial Chaos Expansions
Faverjon, Béatrice
2016-01-01
Because of the high complexity of steady-state or transient fluid flow solvers, non-intrusive uncertainty quantification techniques have been developed in aerodynamic simulations in order to compute the output quantities of interest that are required to evaluate the objective function of an optimization process, for example. The latter is commonly expressed in terms of moments of the quantities of interest, such as the mean, standard deviation, or even higher-order moments (skewness, kurtosis...). Polynomial surrogate models based on homogeneous chaos expansions have often been implemented in this respect. The original approach of uncertainty quantification using such polynomial expansions is however intrusive. It is based on a Galerkin-type projection formulation of the model equations to derive the governing equations for the polynomial expansion coefficients of the output quantities of interest. Both the intrusive and non-intrusive approaches call for the computation of third-order, even fourth-order momen...
Compressive sampling of polynomial chaos expansions: Convergence analysis and sampling strategies
Energy Technology Data Exchange (ETDEWEB)
Hampton, Jerrad; Doostan, Alireza, E-mail: alireza.doostan@colorado.edu
2015-01-01
Sampling orthogonal polynomial bases via Monte Carlo is of interest for uncertainty quantification of models with random inputs, using Polynomial Chaos (PC) expansions. It is known that bounding a probabilistic parameter, referred to as coherence, yields a bound on the number of samples necessary to identify coefficients in a sparse PC expansion via solution to an ℓ{sub 1}-minimization problem. Utilizing results for orthogonal polynomials, we bound the coherence parameter for polynomials of Hermite and Legendre type under their respective natural sampling distribution. In both polynomial bases we identify an importance sampling distribution which yields a bound with weaker dependence on the order of the approximation. For more general orthonormal bases, we propose the coherence-optimal sampling: a Markov Chain Monte Carlo sampling, which directly uses the basis functions under consideration to achieve a statistical optimality among all sampling schemes with identical support. We demonstrate these different sampling strategies numerically in both high-order and high-dimensional, manufactured PC expansions. In addition, the quality of each sampling method is compared in the identification of solutions to two differential equations, one with a high-dimensional random input and the other with a high-order PC expansion. In both cases, the coherence-optimal sampling scheme leads to similar or considerably improved accuracy.
Zernike expansion of derivatives and Laplacians of the Zernike circle polynomials.
Janssen, A J E M
2014-07-01
The partial derivatives and Laplacians of the Zernike circle polynomials occur in various places in the literature on computational optics. In a number of cases, the expansion of these derivatives and Laplacians in the circle polynomials are required. For the first-order partial derivatives, analytic results are scattered in the literature. Results start as early as 1942 in Nijboer's thesis and continue until present day, with some emphasis on recursive computation schemes. A brief historic account of these results is given in the present paper. By choosing the unnormalized version of the circle polynomials, with exponential rather than trigonometric azimuthal dependence, and by a proper combination of the two partial derivatives, a concise form of the expressions emerges. This form is appropriate for the formulation and solution of a model wavefront sensing problem of reconstructing a wavefront on the level of its expansion coefficients from (measurements of the expansion coefficients of) the partial derivatives. It turns out that the least-squares estimation problem arising here decouples per azimuthal order m, and per m the generalized inverse solution assumes a concise analytic form so that singular value decompositions are avoided. The preferred version of the circle polynomials, with proper combination of the partial derivatives, also leads to a concise analytic result for the Zernike expansion of the Laplacian of the circle polynomials. From these expansions, the properties of the Laplacian as a mapping from the space of circle polynomials of maximal degree N, as required in the study of the Neumann problem associated with the transport-of-intensity equation, can be read off within a single glance. Furthermore, the inverse of the Laplacian on this space is shown to have a concise analytic form.
Stress Response of Offshore Structures by Equivalent Polynomial Expansion Techniques
DEFF Research Database (Denmark)
Sigurdsson, Gudfinnur; Nielsen, Søren R.K.
This paper concerns an investigation of the effects of nonlinearity of drag loading on offshore structures excited by 2D wave fields, where the nonlinear term in the Morison equation is replaced by an equivalent cubic expansion. The equivalent cubic expansion coefficients for the equivalent drag...... model are obtained using the least mean square procedure. Numerical results are given. The displacement response and stress response processes obtained using the above loading model are compared with simulation results and those obtained from equivalent linearization of the drag term....
Energy Technology Data Exchange (ETDEWEB)
Lüchow, Arne, E-mail: luechow@rwth-aachen.de [Institute of Physical Chemistry, RWTH Aachen University, 52056 Aachen (Germany); Jülich Aachen Research Alliance (JARA-HPC), 52056 Aachen (Germany); Sturm, Alexander; Schulte, Christoph; Haghighi Mood, Kaveh [Institute of Physical Chemistry, RWTH Aachen University, 52056 Aachen (Germany)
2015-02-28
Jastrow correlation factors play an important role in quantum Monte Carlo calculations. Together with an orbital based antisymmetric function, they allow the construction of highly accurate correlation wave functions. In this paper, a generic expansion of the Jastrow correlation function in terms of polynomials that satisfy both the electron exchange symmetry constraint and the cusp conditions is presented. In particular, an expansion of the three-body electron-electron-nucleus contribution in terms of cuspless homogeneous symmetric polynomials is proposed. The polynomials can be expressed in fairly arbitrary scaling function allowing a generic implementation of the Jastrow factor. It is demonstrated with a few examples that the new Jastrow factor achieves 85%–90% of the total correlation energy in a variational quantum Monte Carlo calculation and more than 90% of the diffusion Monte Carlo correlation energy.
Indian Academy of Sciences (India)
A K Chattopadhyay; C V S Rao
2003-07-01
Here we describe the superiority of Bessel function as base function for radial expansion over Zernicke polynomial in the tomographic reconstruction technique. The causes for the superiority have been described in detail. The superiority has been shown both with simulated data for Kadomtsev’s model for saw-tooth oscillation and real experimental x-ray data from W7-AS Stellarator.
Limitations of polynomial chaos expansions in the Bayesian solution of inverse problems
Energy Technology Data Exchange (ETDEWEB)
Lu, Fei [Lawrence Berkeley National Laboratory (United States); Department of Mathematics, University of California, Berkeley (United States); Morzfeld, Matthias, E-mail: mmo@math.lbl.gov [Lawrence Berkeley National Laboratory (United States); Department of Mathematics, University of California, Berkeley (United States); Tu, Xuemin [Department of Mathematics, University of Kansas (United States); Chorin, Alexandre J. [Lawrence Berkeley National Laboratory (United States); Department of Mathematics, University of California, Berkeley (United States)
2015-02-01
Polynomial chaos expansions are used to reduce the computational cost in the Bayesian solutions of inverse problems by creating a surrogate posterior that can be evaluated inexpensively. We show, by analysis and example, that when the data contain significant information beyond what is assumed in the prior, the surrogate posterior can be very different from the posterior, and the resulting estimates become inaccurate. One can improve the accuracy by adaptively increasing the order of the polynomial chaos, but the cost may increase too fast for this to be cost effective compared to Monte Carlo sampling without a surrogate posterior.
Siminovitch, David; Untidt, Thomas; Nielsen, Niels Chr
2004-01-01
Our recent exact effective Hamiltonian theory (EEHT) for exact analysis of nuclear magnetic resonance (NMR) experiments relied on a novel entanglement of unitary exponential operators via finite expansion of the logarithmic mapping function. In the present study, we introduce simple alternant quotient expressions for the coefficients of the polynomial matrix expansion of these entangled operators. These expressions facilitate an extension of our previous closed solution to the Baker-Campbell-Hausdorff problem for SU(N) systems from Nfunction. The general applicability of these expressions is demonstrated by several examples with relevance for NMR spectroscopy. The specific form of the alternant quotients is also used to demonstrate the fundamentally important equivalence of Sylvester's theorem (also known as the spectral theorem) and the EEHT expansion.
Institute of Scientific and Technical Information of China (English)
E. H. Doha; S. I. El-Soubhy
2001-01-01
The formula of expressing the coefficients of an expansion ofultraspherical polynomials that has been integrated an arbitrary number of times in terms of the coefficients of the original expansion is stated in a more compact form and proved in a simpler way than the formula of Phillips and Karageorghis (1990). A new formula is proved for the q times integration of ultraspherical polynomials, of which the Chebyshev polynomials of the first and second kinds and Legendre polynomials are important special cases. An application of these formulae for solving ordinary differential equations with varying coefficients is discussed.CLC Number：O17 Document ID：AAuthor Resume：E. H. Doha，e-mail: eiddoha@frcu, eun. eg References：[1]Canuto,C. ,Spectral Methods in Fluid Dynamics,Springer,Belrin,1988.[2]Doha,E.H.,An Accurate Solution of Parabolic Equations by Expansion in Ultraspherical Polynomials,Comput. Math. Appl. ,19(1990),75-88.[3]Doha,E. H. ,The Coefficients of Differentiated Expansions and Derivatives of Ultraspherical Polynomials,Comput. Math. Appl.,21(1991),115-122.[4]Doha,E.H. ,The Chebyshev Coefficients of General order Derivatives of an Infinitely Differen-tiable Function in Two or Three Variables,Ann. Univ. Sci. Budapest. Sect. Comput. ,13(1992),83-91.[5]Doha,E. H.,On the Cefficients of Differentiable Expansions of Double and Triple Legendre Polynomials,Ann Univ. Sci. Budapest. Sect. Comput. ,15(1995),25-35.[6]Doha,E.H. ,The Ultraspherical Coefficients of the Moments of a General-Order Derivatives of an Infinitely Differentiable Function,J. Comput. Math. ,89(1998),53-72.[7]Doha,E.H. ,The Coefficients of Differentiated Expansions of Double and Triple Ultraspherical Polynomials,Annales Univ. Sci. Budapest.,Sect. Comp.,19(200),57-73.[8]Doha,E.H. and Al-Kholi,F. M. R. ,An Efficient Double Legerdre Spectral Method for Parabolic and Elliptic Partial Differential Equations,Intern. J. Computer. Math. (toAppear).[9]Fox,L. and Parker,I.B. ,Chebyshev Polynomials in
Piecewise polynomial chaos expansion with an application to brake squeal of a linear brake system
Sarrouy, E.; Dessombz, O.; Sinou, J.-J.
2013-02-01
This paper proposes numerical developments based on polynomial chaos (PC) expansions to process stochastic eigenvalue problems efficiently. These developments are applied to the problem of linear stability calculations for a simplified brake system: the stability of a finite element model of a brake is investigated when its friction coefficient or the contact stiffness are modeled as random parameters. Getting rid of the statistical point of view of the PC method but keeping the principle of a polynomial decomposition of eigenvalues and eigenvectors, the stochastic space is decomposed into several elements to realize a low degree piecewise polynomial approximation of these quantities. An approach relying on continuation principles is compared to the classical dichotomy method to build the partition. Moreover, a criterion for testing accuracy of the decomposition over each cell of the partition without requiring evaluation of exact eigenmodes is proposed and implemented. Several random distributions are tested, including a uniform-like law for description of friction coefficient variation. Results are compared to Monte Carlo simulations so as to determine the method accuracy and efficiency. Some general rules relative to the influence of the friction coefficient or the contact stiffness are also inferred from these calculations.
Efficient linear precoding for massive MIMO systems using truncated polynomial expansion
Müller, Axel
2014-06-01
Massive multiple-input multiple-output (MIMO) techniques have been proposed as a solution to satisfy many requirements of next generation cellular systems. One downside of massive MIMO is the increased complexity of computing the precoding, especially since the relatively \\'antenna-efficient\\' regularized zero-forcing (RZF) is preferred to simple maximum ratio transmission. We develop in this paper a new class of precoders for single-cell massive MIMO systems. It is based on truncated polynomial expansion (TPE) and mimics the advantages of RZF, while offering reduced and scalable computational complexity that can be implemented in a convenient parallel fashion. Using random matrix theory we provide a closed-form expression of the signal-to-interference-and-noise ratio under TPE precoding and compare it to previous works on RZF. Furthermore, the sum rate maximizing polynomial coefficients in TPE precoding are calculated. By simulation, we find that to maintain a fixed peruser rate loss as compared to RZF, the polynomial degree does not need to scale with the system, but it should be increased with the quality of the channel knowledge and signal-to-noise ratio. © 2014 IEEE.
A reduced polynomial chaos expansion method for the stochastic ﬁnite element analysis
Indian Academy of Sciences (India)
B Pascual; S Adhikari
2012-06-01
The stochastic ﬁnite element analysis of elliptic type partial differential equations is considered. A reduced method of the spectral stochastic ﬁnite element method using polynomial chaos is proposed. The method is based on the spectral decomposition of the deterministic system matrix. The reduction is achieved by retaining only the dominant eigenvalues and eigenvectors. The response of the reduced system is expanded as a series of Hermite polynomials, and a Galerkin error minimization approach is applied to obtain the deterministic coefﬁcients of the expansion. The moments and probability density function of the solution are obtained by a process similar to the classical spectral stochastic ﬁnite element method. The method is illustrated using three carefully selected numerical examples, namely, bending of a stochastic beam, ﬂow through porous media with stochastic permeability and transverse bending of a plate with stochastic properties. The results obtained from the proposed method are compared with classical polynomial chaos and direct Monte Carlo simulation results.
Fast and accurate sensitivity analysis of IMPT treatment plans using Polynomial Chaos Expansion
Perkó, Zoltán; van der Voort, Sebastian R.; van de Water, Steven; Hartman, Charlotte M. H.; Hoogeman, Mischa; Lathouwers, Danny
2016-06-01
The highly conformal planned dose distribution achievable in intensity modulated proton therapy (IMPT) can severely be compromised by uncertainties in patient setup and proton range. While several robust optimization approaches have been presented to address this issue, appropriate methods to accurately estimate the robustness of treatment plans are still lacking. To fill this gap we present Polynomial Chaos Expansion (PCE) techniques which are easily applicable and create a meta-model of the dose engine by approximating the dose in every voxel with multidimensional polynomials. This Polynomial Chaos (PC) model can be built in an automated fashion relatively cheaply and subsequently it can be used to perform comprehensive robustness analysis. We adapted PC to provide among others the expected dose, the dose variance, accurate probability distribution of dose-volume histogram (DVH) metrics (e.g. minimum tumor or maximum organ dose), exact bandwidths of DVHs, and to separate the effects of random and systematic errors. We present the outcome of our verification experiments based on 6 head-and-neck (HN) patients, and exemplify the usefulness of PCE by comparing a robust and a non-robust treatment plan for a selected HN case. The results suggest that PCE is highly valuable for both research and clinical applications.
Linear precoding based on polynomial expansion: reducing complexity in massive MIMO
Mueller, Axel
2016-02-29
Massive multiple-input multiple-output (MIMO) techniques have the potential to bring tremendous improvements in spectral efficiency to future communication systems. Counterintuitively, the practical issues of having uncertain channel knowledge, high propagation losses, and implementing optimal non-linear precoding are solved more or less automatically by enlarging system dimensions. However, the computational precoding complexity grows with the system dimensions. For example, the close-to-optimal and relatively “antenna-efficient” regularized zero-forcing (RZF) precoding is very complicated to implement in practice, since it requires fast inversions of large matrices in every coherence period. Motivated by the high performance of RZF, we propose to replace the matrix inversion and multiplication by a truncated polynomial expansion (TPE), thereby obtaining the new TPE precoding scheme which is more suitable for real-time hardware implementation and significantly reduces the delay to the first transmitted symbol. The degree of the matrix polynomial can be adapted to the available hardware resources and enables smooth transition between simple maximum ratio transmission and more advanced RZF. By deriving new random matrix results, we obtain a deterministic expression for the asymptotic signal-to-interference-and-noise ratio (SINR) achieved by TPE precoding in massive MIMO systems. Furthermore, we provide a closed-form expression for the polynomial coefficients that maximizes this SINR. To maintain a fixed per-user rate loss as compared to RZF, the polynomial degree does not need to scale with the system, but it should be increased with the quality of the channel knowledge and the signal-to-noise ratio.
Karajan, N; Otto, D; Oladyshkin, S; Ehlers, W
2014-10-01
A possibility to simulate the mechanical behaviour of the human spine is given by modelling the stiffer structures, i.e. the vertebrae, as a discrete multi-body system (MBS), whereas the softer connecting tissue, i.e. the softer intervertebral discs (IVD), is represented in a continuum-mechanical sense using the finite-element method (FEM). From a modelling point of view, the mechanical behaviour of the IVD can be included into the MBS in two different ways. They can either be computed online in a so-called co-simulation of a MBS and a FEM or offline in a pre-computation step, where a representation of the discrete mechanical response of the IVD needs to be defined in terms of the applied degrees of freedom (DOF) of the MBS. For both methods, an appropriate homogenisation step needs to be applied to obtain the discrete mechanical response of the IVD, i.e. the resulting forces and moments. The goal of this paper was to present an efficient method to approximate the mechanical response of an IVD in an offline computation. In a previous paper (Karajan et al. in Biomech Model Mechanobiol 12(3):453-466, 2012), it was proven that a cubic polynomial for the homogenised forces and moments of the FE model is a suitable choice to approximate the purely elastic response as a coupled function of the DOF of the MBS. In this contribution, the polynomial chaos expansion (PCE) is applied to generate these high-dimensional polynomials. Following this, the main challenge is to determine suitable deformation states of the IVD for pre-computation, such that the polynomials can be constructed with high accuracy and low numerical cost. For the sake of a simple verification, the coupling method and the PCE are applied to the same simplified motion segment of the spine as was used in the previous paper, i.e. two cylindrical vertebrae and a cylindrical IVD in between. In a next step, the loading rates are included as variables in the polynomial response functions to account for a more
A robust and efficient stepwise regression method for building sparse polynomial chaos expansions
Abraham, Simon; Raisee, Mehrdad; Ghorbaniasl, Ghader; Contino, Francesco; Lacor, Chris
2017-03-01
Polynomial Chaos (PC) expansions are widely used in various engineering fields for quantifying uncertainties arising from uncertain parameters. The computational cost of classical PC solution schemes is unaffordable as the number of deterministic simulations to be calculated grows dramatically with the number of stochastic dimension. This considerably restricts the practical use of PC at the industrial level. A common approach to address such problems is to make use of sparse PC expansions. This paper presents a non-intrusive regression-based method for building sparse PC expansions. The most important PC contributions are detected sequentially through an automatic search procedure. The variable selection criterion is based on efficient tools relevant to probabilistic method. Two benchmark analytical functions are used to validate the proposed algorithm. The computational efficiency of the method is then illustrated by a more realistic CFD application, consisting of the non-deterministic flow around a transonic airfoil subject to geometrical uncertainties. To assess the performance of the developed methodology, a detailed comparison is made with the well established LAR-based selection technique. The results show that the developed sparse regression technique is able to identify the most significant PC contributions describing the problem. Moreover, the most important stochastic features are captured at a reduced computational cost compared to the LAR method. The results also demonstrate the superior robustness of the method by repeating the analyses using random experimental designs.
Identification of nonlinear vibrating structures by polynomial expansion in the z-domain
Fasana, Alessandro; Garibaldi, Luigi; Marchesiello, Stefano
2017-02-01
A new method in the frequency domain for the identification of nonlinear vibrating structures is described, by adopting the perspective of nonlinearities as internal feedback forces. The technique is based on a polynomial expansion representation of the frequency response function of the underlying linear system, relying on a z-domain formulation. A least squares approach is adopted to take into account the information of all the frequency response functions but, when large data sets are used, the solution of the resulting system of algebraic linear equations can be a difficult task. A procedure to drastically reduce the matrix dimensions and consequently the computational cost - which largely depends on the number of spectral lines - is adopted, leading to a compact and well conditioned problem. The robustness and numerical performances of the method are demonstrated by its implementation on simulated data from single and two degree of freedom systems with typical nonlinear characteristics.
Polynomial expansion of the precoder for power minimization in large-scale MIMO systems
Sifaou, Houssem
2016-07-26
This work focuses on the downlink of a single-cell large-scale MIMO system in which the base station equipped with M antennas serves K single-antenna users. In particular, we are interested in reducing the implementation complexity of the optimal linear precoder (OLP) that minimizes the total power consumption while ensuring target user rates. As most precoding schemes, a major difficulty towards the implementation of OLP is that it requires fast inversions of large matrices at every new channel realizations. To overcome this issue, we aim at designing a linear precoding scheme providing the same performance of OLP but with lower complexity. This is achieved by applying the truncated polynomial expansion (TPE) concept on a per-user basis. To get a further leap in complexity reduction and allow for closed-form expressions of the per-user weighting coefficients, we resort to the asymptotic regime in which M and K grow large with a bounded ratio. Numerical results are used to show that the proposed TPE precoding scheme achieves the same performance of OLP with a significantly lower implementation complexity. © 2016 IEEE.
Kersaudy, Pierric; Sudret, Bruno; Varsier, Nadège; Picon, Odile; Wiart, Joe
2015-04-01
In numerical dosimetry, the recent advances in high performance computing led to a strong reduction of the required computational time to assess the specific absorption rate (SAR) characterizing the human exposure to electromagnetic waves. However, this procedure remains time-consuming and a single simulation can request several hours. As a consequence, the influence of uncertain input parameters on the SAR cannot be analyzed using crude Monte Carlo simulation. The solution presented here to perform such an analysis is surrogate modeling. This paper proposes a novel approach to build such a surrogate model from a design of experiments. Considering a sparse representation of the polynomial chaos expansions using least-angle regression as a selection algorithm to retain the most influential polynomials, this paper proposes to use the selected polynomials as regression functions for the universal Kriging model. The leave-one-out cross validation is used to select the optimal number of polynomials in the deterministic part of the Kriging model. The proposed approach, called LARS-Kriging-PC modeling, is applied to three benchmark examples and then to a full-scale metamodeling problem involving the exposure of a numerical fetus model to a femtocell device. The performances of the LARS-Kriging-PC are compared to an ordinary Kriging model and to a classical sparse polynomial chaos expansion. The LARS-Kriging-PC appears to have better performances than the two other approaches. A significant accuracy improvement is observed compared to the ordinary Kriging or to the sparse polynomial chaos depending on the studied case. This approach seems to be an optimal solution between the two other classical approaches. A global sensitivity analysis is finally performed on the LARS-Kriging-PC model of the fetus exposure problem.
Energy Technology Data Exchange (ETDEWEB)
Kersaudy, Pierric, E-mail: pierric.kersaudy@orange.com [Orange Labs, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France); Whist Lab, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France); ESYCOM, Université Paris-Est Marne-la-Vallée, 5 boulevard Descartes, 77700 Marne-la-Vallée (France); Sudret, Bruno [ETH Zürich, Chair of Risk, Safety and Uncertainty Quantification, Stefano-Franscini-Platz 5, 8093 Zürich (Switzerland); Varsier, Nadège [Orange Labs, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France); Whist Lab, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France); Picon, Odile [ESYCOM, Université Paris-Est Marne-la-Vallée, 5 boulevard Descartes, 77700 Marne-la-Vallée (France); Wiart, Joe [Orange Labs, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France); Whist Lab, 38 avenue du Général Leclerc, 92130 Issy-les-Moulineaux (France)
2015-04-01
In numerical dosimetry, the recent advances in high performance computing led to a strong reduction of the required computational time to assess the specific absorption rate (SAR) characterizing the human exposure to electromagnetic waves. However, this procedure remains time-consuming and a single simulation can request several hours. As a consequence, the influence of uncertain input parameters on the SAR cannot be analyzed using crude Monte Carlo simulation. The solution presented here to perform such an analysis is surrogate modeling. This paper proposes a novel approach to build such a surrogate model from a design of experiments. Considering a sparse representation of the polynomial chaos expansions using least-angle regression as a selection algorithm to retain the most influential polynomials, this paper proposes to use the selected polynomials as regression functions for the universal Kriging model. The leave-one-out cross validation is used to select the optimal number of polynomials in the deterministic part of the Kriging model. The proposed approach, called LARS-Kriging-PC modeling, is applied to three benchmark examples and then to a full-scale metamodeling problem involving the exposure of a numerical fetus model to a femtocell device. The performances of the LARS-Kriging-PC are compared to an ordinary Kriging model and to a classical sparse polynomial chaos expansion. The LARS-Kriging-PC appears to have better performances than the two other approaches. A significant accuracy improvement is observed compared to the ordinary Kriging or to the sparse polynomial chaos depending on the studied case. This approach seems to be an optimal solution between the two other classical approaches. A global sensitivity analysis is finally performed on the LARS-Kriging-PC model of the fetus exposure problem.
Dolgov, Sergey
2015-11-03
We apply the tensor train (TT) decomposition to construct the tensor product polynomial chaos expansion (PCE) of a random field, to solve the stochastic elliptic diffusion PDE with the stochastic Galerkin discretization, and to compute some quantities of interest (mean, variance, and exceedance probabilities). We assume that the random diffusion coefficient is given as a smooth transformation of a Gaussian random field. In this case, the PCE is delivered by a complicated formula, which lacks an analytic TT representation. To construct its TT approximation numerically, we develop the new block TT cross algorithm, a method that computes the whole TT decomposition from a few evaluations of the PCE formula. The new method is conceptually similar to the adaptive cross approximation in the TT format but is more efficient when several tensors must be stored in the same TT representation, which is the case for the PCE. In addition, we demonstrate how to assemble the stochastic Galerkin matrix and to compute the solution of the elliptic equation and its postprocessing, staying in the TT format. We compare our technique with the traditional sparse polynomial chaos and the Monte Carlo approaches. In the tensor product polynomial chaos, the polynomial degree is bounded for each random variable independently. This provides higher accuracy than the sparse polynomial set or the Monte Carlo method, but the cardinality of the tensor product set grows exponentially with the number of random variables. However, when the PCE coefficients are implicitly approximated in the TT format, the computations with the full tensor product polynomial set become possible. In the numerical experiments, we confirm that the new methodology is competitive in a wide range of parameters, especially where high accuracy and high polynomial degrees are required.
Dolgov, S.
2015-03-11
We apply the Tensor Train (TT) decomposition to construct the tensor product Polynomial Chaos Expansion (PCE) of a random field, to solve the stochastic elliptic diffusion PDE with the stochastic Galerkin discretization, and to compute some quantities of interest (mean, variance, exceedance probabilities). We assume that the random diffusion coefficient is given as a smooth transformation of a Gaussian random field. In this case, the PCE is delivered by a complicated formula, which lacks an analytic TT representation. To construct its TT approximation numerically, we develop the new block TT cross algorithm, a method that computes the whole TT decomposition from a few evaluations of the PCE formula. The new method is conceptually similar to the adaptive cross approximation in the TT format, but is more efficient when several tensors must be stored in the same TT representation, which is the case for the PCE. Besides, we demonstrate how to assemble the stochastic Galerkin matrix and to compute the solution of the elliptic equation and its post-processing, staying in the TT format. We compare our technique with the traditional sparse polynomial chaos and the Monte Carlo approaches. In the tensor product polynomial chaos, the polynomial degree is bounded for each random variable independently. This provides higher accuracy than the sparse polynomial set or the Monte Carlo method, but the cardinality of the tensor product set grows exponentially with the number of random variables. However, when the PCE coefficients are implicitly approximated in the TT format, the computations with the full tensor product polynomial set become possible. In the numerical experiments, we confirm that the new methodology is competitive in a wide range of parameters, especially where high accuracy and high polynomial degrees are required.
Directory of Open Access Journals (Sweden)
Mahmoud Paripour
2014-08-01
Full Text Available In this paper, the Bernstein polynomials are used to approximatethe solutions of linear integral equations with multiple time lags (IEMTL through expansion methods (collocation method, partition method, Galerkin method. The method is discussed in detail and illustrated by solving some numerical examples. Comparison between the exact and approximated results obtained from these methods is carried out
Ciriello, V.; Di Federico, V.; Riva, M.; Cadini, F.; De Sanctis, J.; Zio, E.; Guadagnini, A.
2012-04-01
We perform a Global Sensitivity Analysis (GSA) of a transport model used to compute the peak radionuclide concentration at a given control location in a randomly heterogeneous aquifer, following a release from a near surface repository of radioactive waste and subsequent contaminant migration within the host porous medium. We illustrate how uncertainty stemming from incomplete characterization of (a) the correlation scale of the variogram of hydraulic conductivity, (b) the partition coefficient associated with sorption of the migrating radionuclide, and (c) the effective dispersivity at the scale of interest propagates to the first two (ensemble) moments of the peak solute concentration detected at a target location within a two-dimensional randomly heterogeneous hydraulic conductivity field. We treat the uncertain system parameters as independent random variables and perform a variance-based GSA within a numerical Monte Carlo framework. Groundwater flow and transport are solved by randomly sampling the space of the uncertain parameters for an ensemble of generated hydraulic conductivity realizations. The Sobol indices are adopted as sensitivity measures. These are calculated by employing a Polynomial Chaos Expansion (PCE) technique. The PCE-based representation of the response surface of the adopted transport model is then adopted as a surrogate model of the transport process to reduce the computational burden associated with a standard Monte Carlo solution of the original governing equations. This methodology allows identifying the relative influence of the selected uncertain parameters on the target (ensemble) moments of peak concentrations. Our results suggest that the ensemble mean of peak concentration is strongly influenced by the partition coefficient and the longitudinal dispersivity for the scenario analyzed. On the other hand, the hydraulic conductivity correlation scale plays an important role in the variance of the calculated peak concentration values
A Posteriori Error Analysis of Stochastic Differential Equations Using Polynomial Chaos Expansions
Butler, T.
2011-01-01
We develop computable a posteriori error estimates for linear functionals of a solution to a general nonlinear stochastic differential equation with random model/source parameters. These error estimates are based on a variational analysis applied to stochastic Galerkin methods for forward and adjoint problems. The result is a representation for the error estimate as a polynomial in the random model/source parameter. The advantage of this method is that we use polynomial chaos representations for the forward and adjoint systems to cheaply produce error estimates by simple evaluation of a polynomial. By comparison, the typical method of producing such estimates requires repeated forward/adjoint solves for each new choice of random parameter. We present numerical examples showing that there is excellent agreement between these methods. © 2011 Society for Industrial and Applied Mathematics.
Mason, JC
2002-01-01
Chebyshev polynomials crop up in virtually every area of numerical analysis, and they hold particular importance in recent advances in subjects such as orthogonal polynomials, polynomial approximation, numerical integration, and spectral methods. Yet no book dedicated to Chebyshev polynomials has been published since 1990, and even that work focused primarily on the theoretical aspects. A broad, up-to-date treatment is long overdue.Providing highly readable exposition on the subject''s state of the art, Chebyshev Polynomials is just such a treatment. It includes rigorous yet down-to-earth coverage of the theory along with an in-depth look at the properties of all four kinds of Chebyshev polynomials-properties that lead to a range of results in areas such as approximation, series expansions, interpolation, quadrature, and integral equations. Problems in each chapter, ranging in difficulty from elementary to quite advanced, reinforce the concepts and methods presented.Far from being an esoteric subject, Chebysh...
Multiplication of a Schubert polynomial by a Stanley symmetric polynomial
Assaf, Sami
2017-01-01
We prove, combinatorially, that the product of a Schubert polynomial by a Stanley symmetric polynomial is a truncated Schubert polynomial. Using Monk's rule, we derive a nonnegative combinatorial formula for the Schubert polynomial expansion of a truncated Schubert polynomial. Combining these results, we give a nonnegative combinatorial rule for the product of a Schubert and a Schur polynomial in the Schubert basis.
Alvarez, G.; Şen, C.; Furukawa, N.; Motome, Y.; Dagotto, E.
2005-05-01
A software library is presented for the polynomial expansion method (PEM) of the density of states (DOS) introduced in [Y. Motome, N. Furukawa, J. Phys. Soc. Japan 68 (1999) 3853; N. Furukawa, Y. Motome, H. Nakata, Comput. Phys. Comm. 142 (2001) 410]. The library provides all necessary functions for the use of the PEM and its truncated version (TPEM) in a model independent way. The PEM/TPEM replaces the exact diagonalization of the one electron sector in models for fermions coupled to classical fields. The computational cost of the algorithm is O(N)—with N the number of lattice sites—for the TPEM [N. Furukawa, Y. Motome, J. Phys. Soc. Japan 73 (2004) 1482] which should be contrasted with the computational cost of the diagonalization technique that scales as O(N). The method is applied for the first time to a double exchange model with finite Hund coupling and also to diluted spin-fermion models. Program summaryTitle of library:TPEM Catalogue identifier: ADVK Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVK Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland No. of lines in distributed program, including test data, etc.: 1707 No. of bytes in distributed program, including test data, etc.: 13 644 Distribution format:tar.gz Operating system:Linux, UNIX Number of files:4 plus 1 test program Programming language used:C Computer:PC Nature of the physical problem:The study of correlated electrons coupled to classical fields appears in the treatment of many materials of much current interest in condensed matter theory, e.g., manganites, diluted magnetic semiconductors and high temperature superconductors among others. Method of solution: Typically an exact diagonalization of the electronic sector is performed in this type of models for each configuration of classical fields, which are integrated using a classical Monte Carlo algorithm. A polynomial expansion of the density of states is able to replace the exact
Mahata, Avik; Mukhopadhyay, Tanmoy; Adhikari, Sondipon
2016-03-01
Nano-twinned structures are mechanically stronger, ductile and stable than its non-twinned form. We have investigated the effect of varying twin spacing and twin boundary width (TBW) on the yield strength of the nano-twinned copper in a probabilistic framework. An efficient surrogate modelling approach based on polynomial chaos expansion has been proposed for the analysis. Effectively utilising 15 sets of expensive molecular dynamics simulations, thousands of outputs have been obtained corresponding to different sets of twin spacing and twin width using virtual experiments based on the surrogates. One of the major outcomes of this work is that there exists an optimal combination of twin boundary spacing and twin width until which the strength can be increased and after that critical point the nanowires weaken. This study also reveals that the yield strength of nano-twinned copper is more sensitive to TBW than twin spacing. Such robust inferences have been possible to be drawn only because of applying the surrogate modelling approach, which makes it feasible to obtain results corresponding to 40 000 combinations of different twin boundary spacing and twin width in a computationally efficient framework.
Freud, Géza
1971-01-01
Orthogonal Polynomials contains an up-to-date survey of the general theory of orthogonal polynomials. It deals with the problem of polynomials and reveals that the sequence of these polynomials forms an orthogonal system with respect to a non-negative m-distribution defined on the real numerical axis. Comprised of five chapters, the book begins with the fundamental properties of orthogonal polynomials. After discussing the momentum problem, it then explains the quadrature procedure, the convergence theory, and G. Szegő's theory. This book is useful for those who intend to use it as referenc
Energy Technology Data Exchange (ETDEWEB)
Blackett, S.A. [Univ. of Auckland (New Zealand). Dept of Engineering Science
1996-02-01
Numerical analysis is an important part of Engineering. Frequently relationships are not adequately understood, or too complicated to be represented by theoretical formulae. Instead, empirical approximations based on observed relationships can be used for simple fast and accurate evaluations. Historically, storage of data has been a large constraint on approximately methods. So the challenge is to find a sufficiently accurate representation of data which is valid over as large a range as possible while requiring the storage of only a few numerical values. Polynomials, popular as approximation functions because of their simplicity, can be used to represent simple data. Equation 1.1 shows a simple 3rd order polynomial approximation. However, just increasing the order and number of terms included in a polynomial approximation does not improve the overall result. Although the function may fit exactly to observed data, between these points it is likely that the approximation is increasingly less smooth and probably inadequate. An alternative to adding further terms to the approximation is to make the approximation rational. Equation 1.2 shows a rational polynomial, 3rd order in the numerator and denominator. A rational polynomial approximation allows poles and this can greatly enhance an approximation. In Sections 2 and 3 two different methods for fitting rational polynomials to a given data set are detailed. In Section 4, consideration is given to different rational polynomials used on adjacent regions. Section 5 shows the performance of the rational polynomial algorithms. Conclusions are presented in Section 6.
Dobbs, David E.
2010-01-01
This note develops and implements the theory of polynomial asymptotes to (graphs of) rational functions, as a generalization of the classical topics of horizontal asymptotes and oblique/slant asymptotes. Applications are given to hyperbolic asymptotes. Prerequisites include the division algorithm for polynomials with coefficients in the field of…
Befriending Askey-Wilson polynomials
Szabłowski, Paweł J
2011-01-01
Although our main interest is with the Askey-Wilson (AW) polynomials we recall and review four other families of the so-called Askey-Wilson scheme of polynomials. We do this for completeness as well as for better exposition of AW properties. Our main results concentrate on the complex parameters case, revealing new fascinating symmetries between the variables and some of the parameters. In particular we express Askey-Wilson polynomials as linear combinations of Al-Salam--Chihara (ASC) polynomials which together with the obtained earlier expansion of the Askey-Wilson density forms complete generalization of the situation met in the case of Al-Salam--Chihara and q-Hermite polynomials and the Poisson-Mehler expansion formula. As a by-product we get useful identities involving ASC polynomials. Finally by certain re-scaling of variables and parameters we arrive to AW polynomials and AW densities that have clear probabilistic interpretation. We recall some known and present some believed to be unknown identities an...
Institute of Scientific and Technical Information of China (English)
王雷
2008-01-01
<正>Polynomial functions are among the sim- plest expressions in algebra.They are easy to evaluate:only addition and repeated multipli- cation are required.Because of this,they are often used to approximate other more compli-
Haglund's conjecture on 3-column Macdonald polynomials
Blasiak, Jonah
2014-01-01
We prove a positive combinatorial formula for the Schur expansion of LLT polynomials indexed by a 3-tuple of skew shapes. This verifies a conjecture of Haglund. The proof requires expressing a noncommutative Schur function as a positive sum of monomials in Lam's algebra of ribbon Schur operators. Combining this result with the expression of Haglund, Haiman, and Loehr for transformed Macdonald polynomials in terms of LLT polynomials then yields a positive combinatorial rule for transformed Mac...
Koornwinder, T.H.
2012-01-01
Askey-Wilson polynomial refers to a four-parameter family of q-hypergeometric orthogonal polynomials which contains all families of classical orthogonal polynomials (in the wide sense) as special or limit cases.
Polynomial-Chaos-based Kriging
Schöbi, R; Sudret, B.; Wiart, J.
2015-01-01
International audience; Computer simulation has become the standard tool in many engineering fields for designing and optimizing systems, as well as for assessing their reliability. Optimization and uncertainty quantification problems typically require a large number of runs of the computational model at hand, which may not be feasible with high-fidelity models directly. Thus surrogate models (a.k.a metamodels) have been increasingly investigated in the last decade. Polynomial Chaos Expansion...
Reddy, A Satyanarayana
2011-01-01
A graph $X$ is said to be a pattern polynomial graph if its adjacency algebra is a coherent algebra. In this study we will find a necessary and sufficient condition for a graph to be a pattern polynomial graph. Some of the properties of the graphs which are polynomials in the pattern polynomial graph have been studied. We also identify known graph classes which are pattern polynomial graphs.
New classes of test polynomials of polynomial algebras
Institute of Scientific and Technical Information of China (English)
冯克勤; 余解台
1999-01-01
A polynomial p in a polynomial algebra over a field is called a test polynomial if any endomorphism of the polynomial algebra that fixes p is an automorphism. some classes of new test polynomials recognizing nonlinear automorphisms of polynomial algebras are given. In the odd prime characteristic case, test polynomials recognizing non-semisimple automorphisms are also constructed.
ON ABEL-GONTSCHAROFF-GOULD'S POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
He Tianxiao; Leetsch C. Hsu; Peter J. S. Shiue
2003-01-01
In this paper a connective study of Gould's annihilation coefficients and Abel-Gontscharoff polynomials is presented. It is shown that Gould's annihilation coefficients and Abel-Gontscharoff polynomials are actually equivalent to each other under certain linear substitutions for the variables. Moreover, a pair of related expansion formulas involving Gontscharoff's remainder and a new form of it are demonstrated, and also illustrated with several examples.
Polynomial Subtraction Method for Disconnected Quark Loops
Liu, Quan; Morgan, Ron
2014-01-01
The polynomial subtraction method, a new numerical approach for reducing the noise variance of Lattice QCD disconnected matrix elements calculation, is introduced in this paper. We use the MinRes polynomial expansion of the QCD matrix as the approximation to the matrix inverse and get a significant reduction in the variance calculation. We compare our results with that of the perturbative subtraction and find that the new strategy yields a faster decrease in variance which increases with quark mass.
Improved polynomial remainder sequences for Ore polynomials.
Jaroschek, Maximilian
2013-11-01
Polynomial remainder sequences contain the intermediate results of the Euclidean algorithm when applied to (non-)commutative polynomials. The running time of the algorithm is dependent on the size of the coefficients of the remainders. Different ways have been studied to make these as small as possible. The subresultant sequence of two polynomials is a polynomial remainder sequence in which the size of the coefficients is optimal in the generic case, but when taking the input from applications, the coefficients are often larger than necessary. We generalize two improvements of the subresultant sequence to Ore polynomials and derive a new bound for the minimal coefficient size. Our approach also yields a new proof for the results in the commutative case, providing a new point of view on the origin of the extraneous factors of the coefficients.
Factoring Polynomials and Fibonacci.
Schwartzman, Steven
1986-01-01
Discusses the factoring of polynomials and Fibonacci numbers, offering several challenges teachers can give students. For example, they can give students a polynomial containing large numbers and challenge them to factor it. (JN)
Palindromic random trigonometric polynomials
Conrey, J. Brian; Farmer, David W.; Imamoglu, Özlem
2008-01-01
We show that if a real trigonometric polynomial has few real roots, then the trigonometric polynomial obtained by writing the coefficients in reverse order must have many real roots. This is used to show that a class of random trigonometric polynomials has, on average, many real roots. In the case that the coefficients of a real trigonometric polynomial are independently and identically distributed, but with no other assumptions on the distribution, the expected fraction of real zeros is at l...
Polynomial Datapaths Optimization
Parta, Hojat
2014-01-01
The research presented focuses on optimization of polynomials using algebraic manipulations at the high level and digital arithmetic techniques at the implementation level. Previous methods lacked any algebraic understanding of the polynomials or only exposed limited potential. We have treated the polynomial optimization problem in abstract algebra allowing us algebraic freedom to transform polynomials. Unlike previous attempts where only a set of limited benchmarks have been used, we have fo...
Approximation for Transient of Nonlinear Circuits Using RHPM and BPES Methods
Directory of Open Access Journals (Sweden)
H. Vazquez-Leal
2013-01-01
Full Text Available The microelectronics area constantly demands better and improved circuit simulation tools. Therefore, in this paper, rational homotopy perturbation method and Boubaker Polynomials Expansion Scheme are applied to a differential equation from a nonlinear circuit. Comparing the results obtained by both techniques revealed that they are effective and convenient.
Blankertz, Raoul
2011-01-01
This diploma thesis is concerned with functional decomposition $f = g \\circ h$ of polynomials. First an algorithm is described which computes decompositions in polynomial time. This algorithm was originally proposed by Zippel (1991). A bound for the number of minimal collisions is derived. Finally a proof of a conjecture in von zur Gathen, Giesbrecht & Ziegler (2010) is given, which states a classification for a special class of decomposable polynomials.
A Bivariate Analogue to the Composed Product of Polynomials
Institute of Scientific and Technical Information of China (English)
Donald Mills; Kent M. Neuerburg
2003-01-01
The concept of a composed product for univariate polynomials has been explored extensively by Brawley, Brown, Carlitz, Gao,Mills, et al. Starting with these fundamental ideas andutilizing fractional power series representation(in particular, the Puiseux expansion) of bivariate polynomials, we generalize the univariate results. We define a bivariate composed sum,composed multiplication,and composed product (based on function composition). Further, we investigate the algebraic structure of certain classes of bivariate polynomials under these operations. We also generalize a result of Brawley and Carlitz concerningthe decomposition of polynomials into irreducibles.
Polynomial Fibonacci-Hessenberg matrices
Energy Technology Data Exchange (ETDEWEB)
Esmaeili, Morteza [Dept. of Mathematical Sciences, Isfahan University of Technology, 84156-83111 Isfahan (Iran, Islamic Republic of)], E-mail: emorteza@cc.iut.ac.ir; Esmaeili, Mostafa [Dept. of Electrical and Computer Engineering, Isfahan University of Technology, 84156-83111 Isfahan (Iran, Islamic Republic of)
2009-09-15
A Fibonacci-Hessenberg matrix with Fibonacci polynomial determinant is referred to as a polynomial Fibonacci-Hessenberg matrix. Several classes of polynomial Fibonacci-Hessenberg matrices are introduced. The notion of two-dimensional Fibonacci polynomial array is introduced and three classes of polynomial Fibonacci-Hessenberg matrices satisfying this property are given.
Polynomial Graphs and Symmetry
Goehle, Geoff; Kobayashi, Mitsuo
2013-01-01
Most quadratic functions are not even, but every parabola has symmetry with respect to some vertical line. Similarly, every cubic has rotational symmetry with respect to some point, though most cubics are not odd. We show that every polynomial has at most one point of symmetry and give conditions under which the polynomial has rotational or…
Nonnegativity of uncertain polynomials
Directory of Open Access Journals (Sweden)
iljak Dragoslav D.
1998-01-01
Full Text Available The purpose of this paper is to derive tests for robust nonnegativity of scalar and matrix polynomials, which are algebraic, recursive, and can be completed in finite number of steps. Polytopic families of polynomials are considered with various characterizations of parameter uncertainty including affine, multilinear, and polynomic structures. The zero exclusion condition for polynomial positivity is also proposed for general parameter dependencies. By reformulating the robust stability problem of complex polynomials as positivity of real polynomials, we obtain new sufficient conditions for robust stability involving multilinear structures, which can be tested using only real arithmetic. The obtained results are applied to robust matrix factorization, strict positive realness, and absolute stability of multivariable systems involving parameter dependent transfer function matrices.
Yu, Jiun-Hung
2012-01-01
Polynomial remainder codes are a large class of codes derived from the Chinese remainder theorem that includes Reed-Solomon codes as a special case. In this paper, we revisit these codes and study them more carefully than in previous work. We explicitly allow the code symbols to be polynomials of different degrees, which leads to two different notions of weight and distance. Algebraic decoding is studied in detail. If the moduli are not irreducible, the notion of an error locator polynomial is replaced by an error factor polynomial. We then obtain a collection of gcd-based decoding algorithms, some of which are not quite standard even when specialized to Reed-Solomon codes.
Chemical Reaction Networks for Computing Polynomials.
Salehi, Sayed Ahmad; Parhi, Keshab K; Riedel, Marc D
2017-01-20
Chemical reaction networks (CRNs) provide a fundamental model in the study of molecular systems. Widely used as formalism for the analysis of chemical and biochemical systems, CRNs have received renewed attention as a model for molecular computation. This paper demonstrates that, with a new encoding, CRNs can compute any set of polynomial functions subject only to the limitation that these functions must map the unit interval to itself. These polynomials can be expressed as linear combinations of Bernstein basis polynomials with positive coefficients less than or equal to 1. In the proposed encoding approach, each variable is represented using two molecular types: a type-0 and a type-1. The value is the ratio of the concentration of type-1 molecules to the sum of the concentrations of type-0 and type-1 molecules. The proposed encoding naturally exploits the expansion of a power-form polynomial into a Bernstein polynomial. Molecular encoders for converting any input in a standard representation to the fractional representation as well as decoders for converting the computed output from the fractional to a standard representation are presented. The method is illustrated first for generic CRNs; then chemical reactions designed for an example are mapped to DNA strand-displacement reactions.
Algorithms for Testing Monomials in Multivariate Polynomials
Chen, Zhixiang; Liu, Yang; Schweller, Robert
2010-01-01
This paper is our second step towards developing a theory of testing monomials in multivariate polynomials. The central question is to ask whether a polynomial represented by an arithmetic circuit has some types of monomials in its sum-product expansion. The complexity aspects of this problem and its variants have been investigated in our first paper by Chen and Fu (2010), laying a foundation for further study. In this paper, we present two pairs of algorithms. First, we prove that there is a randomized $O^*(p^k)$ time algorithm for testing $p$-monomials in an $n$-variate polynomial of degree $k$ represented by an arithmetic circuit, while a deterministic $O^*(6.4^k + p^k)$ time algorithm is devised when the circuit is a formula, here $p$ is a given prime number. Second, we present a deterministic $O^*(2^k)$ time algorithm for testing multilinear monomials in $\\Pi_m\\Sigma_2\\Pi_t\\times \\Pi_k\\Pi_3$ polynomials, while a randomized $O^*(1.5^k)$ algorithm is given for these polynomials. The first algorithm extends...
Inverse polynomial reconstruction method in DCT domain
Dadkhahi, Hamid; Gotchev, Atanas; Egiazarian, Karen
2012-12-01
The discrete cosine transform (DCT) offers superior energy compaction properties for a large class of functions and has been employed as a standard tool in many signal and image processing applications. However, it suffers from spurious behavior in the vicinity of edge discontinuities in piecewise smooth signals. To leverage the sparse representation provided by the DCT, in this article, we derive a framework for the inverse polynomial reconstruction in the DCT expansion. It yields the expansion of a piecewise smooth signal in terms of polynomial coefficients, obtained from the DCT representation of the same signal. Taking advantage of this framework, we show that it is feasible to recover piecewise smooth signals from a relatively small number of DCT coefficients with high accuracy. Furthermore, automatic methods based on minimum description length principle and cross-validation are devised to select the polynomial orders, as a requirement of the inverse polynomial reconstruction method in practical applications. The developed framework can considerably enhance the performance of the DCT in sparse representation of piecewise smooth signals. Numerical results show that denoising and image approximation algorithms based on the proposed framework indicate significant improvements over wavelet counterparts for this class of signals.
Fuzzy Morphological Polynomial Image Representation
Directory of Open Access Journals (Sweden)
Chin-Pan Huang
2010-01-01
Full Text Available A novel signal representation using fuzzy mathematical morphology is developed. We take advantage of the optimum fuzzy fitting and the efficient implementation of morphological operators to extract geometric information from signals. The new representation provides results analogous to those given by the polynomial transform. Geometrical decomposition of a signal is achieved by windowing and applying sequentially fuzzy morphological opening with structuring functions. The resulting representation is made to resemble an orthogonal expansion by constraining the results of opening to equate adapted structuring functions. Properties of the geometric decomposition are considered and used to calculate the adaptation parameters. Our procedure provides an efficient and flexible representation which can be efficiently implemented in parallel. The application of the representation is illustrated in data compression and fractal dimension estimation temporal signals and images.
Asymptotic analysis of the Nörlund and Stirling polynomials
Directory of Open Access Journals (Sweden)
Mark Daniel Ward
2012-04-01
Full Text Available We provide a full asymptotic analysis of the N{\\"o}rlund polynomials and Stirling polynomials. We give a general asymptotic expansion---to any desired degree of accuracy---when the parameter is not an integer. We use singularity analysis, Hankel contours, and transfer theory. This investigation was motivated by a need for such a complete asymptotic description, with parameter 1/2, during this author's recent solution of Wilf's 3rd (previously Unsolved Problem.
JACOBI POLYNOMIALS USED TO INVERT THE LAPLACE TRANSFORM
Institute of Scientific and Technical Information of China (English)
A.Al-Shuaibi; F.Al-Rawjih
2004-01-01
Given the Laplace transform F(s) of a function f(t), we develop a new algorithm to find an approximation to f(t) by the use of the classical Jacobi polynomials. The main contribution of our work is the development of a new and very effective method to determine the coefficients in the finite series expansion that approximation f(t) in terms of Jacobi polynomials. Some numerical examples are illustrated.
Additive and polynomial representations
Krantz, David H; Suppes, Patrick
1971-01-01
Additive and Polynomial Representations deals with major representation theorems in which the qualitative structure is reflected as some polynomial function of one or more numerical functions defined on the basic entities. Examples are additive expressions of a single measure (such as the probability of disjoint events being the sum of their probabilities), and additive expressions of two measures (such as the logarithm of momentum being the sum of log mass and log velocity terms). The book describes the three basic procedures of fundamental measurement as the mathematical pivot, as the utiliz
STABILITY OF SWITCHED POLYNOMIAL SYSTEMS
Institute of Scientific and Technical Information of China (English)
Zhiqiang LI; Yupeng QIAO; Hongsheng QI; Daizhan CHENG
2008-01-01
This paper investigates the stability of (switched) polynomial systems. Using semi-tensor product of matrices, the paper develops two tools for testing the stability of a (switched) polynomial system. One is to convert a product of multi-variable polynomials into a canonical form, and the other is an easily verifiable sufficient condition to justify whether a multi-variable polynomial is positive definite. Using these two tools, the authors construct a polynomial function as a candidate Lyapunov function and via testing its derivative the authors provide some sufficient conditions for the global stability of polynomial systems.
On Generalized Bell Polynomials
Directory of Open Access Journals (Sweden)
Roberto B. Corcino
2011-01-01
Full Text Available It is shown that the sequence of the generalized Bell polynomials Sn(x is convex under some restrictions of the parameters involved. A kind of recurrence relation for Sn(x is established, and some numbers related to the generalized Bell numbers and their properties are investigated.
Tricubic polynomial interpolation.
Birkhoff, G
1971-06-01
A new triangular "finite element" is described; it involves the 12-parameter family of all quartic polynomial functions that are "tricubic" in that their variation is cubic along any parallel to any side of the triangle. An interpolation scheme is described that approximates quite accurately any smooth function on any triangulated domain by a continuously differentiable function, tricubic on each triangular element.
Calculators and Polynomial Evaluation.
Weaver, J. F.
The intent of this paper is to suggest and illustrate how electronic hand-held calculators, especially non-programmable ones with limited data-storage capacity, can be used to advantage by students in one particular aspect of work with polynomial functions. The basic mathematical background upon which calculator application is built is summarized.…
Hetyei, Gábor
2010-01-01
We introduce the short toric polynomial associated to a graded Eulerian poset. This polynomial contains the same information as the two toric polynomials introduced by Stanley, but allows different algebraic manipulations. The intertwined recurrence defining Stanley's toric polynomials may be replaced by a single recurrence, in which the degree of the discarded terms is independent of the rank. A short toric variant of the formula by Bayer and Ehrenborg, expressing the toric $h$-vector in terms of the $cd$-index, may be stated in a rank-independent form, and it may be shown using weighted lattice path enumeration and the reflection principle. We use our techniques to derive a formula expressing the toric $h$-vector of a dual simplicial Eulerian poset in terms of its $f$-vector. This formula implies Gessel's formula for the toric $h$-vector of a cube, and may be used to prove that the nonnegativity of the toric $h$-vector of a simple polytope is a consequence of the Generalized Lower Bound Theorem holding for ...
Computing the Alexander Polynomial Numerically
DEFF Research Database (Denmark)
Hansen, Mikael Sonne
2006-01-01
Explains how to construct the Alexander Matrix and how this can be used to compute the Alexander polynomial numerically.......Explains how to construct the Alexander Matrix and how this can be used to compute the Alexander polynomial numerically....
ON PROPERTIES OF DIFFERENCE POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
Chen Zongxuan; Huang Zhibo; Zheng Xiumin
2011-01-01
We study the value distribution of difference polynomials of meromorphic functions, and extend classical theorems of Tumura-Clunie type to difference polynomials. We also consider the value distribution of f(z)f(z+c).
Chromatic polynomials for simplicial complexes
DEFF Research Database (Denmark)
Møller, Jesper Michael; Nord, Gesche
2016-01-01
In this note we consider s s -chromatic polynomials for finite simplicial complexes. When s=1 s=1 , the 1 1 -chromatic polynomial is just the usual graph chromatic polynomial of the 1 1 -skeleton. In general, the s s -chromatic polynomial depends on the s s -skeleton and its value at r r is the n...
R.J. Stroeker (Roel)
2002-01-01
textabstractA Q-derived polynomial is a univariate polynomial, defined over the rationals, with the property that its zeros, and those of all its derivatives are rational numbers. There is a conjecture that says that Q-derived polynomials of degree 4 with distinct roots for themselves and all their
R.J. Stroeker (Roel)
2006-01-01
textabstractA Q-derived polynomial is a univariate polynomial, defined over the rationals, with the property that its zeros, and those of all its derivatives are rational numbers. There is a conjecture that says that Q-derived polynomials of degree 4 with distinct roots for themselves and all their
Determinants and Polynomial Root Structure
De Pillis, L. G.
2005-01-01
A little known property of determinants is developed in a manner accessible to beginning undergraduates in linear algebra. Using the language of matrix theory, a classical result by Sylvester that describes when two polynomials have a common root is recaptured. Among results concerning the structure of polynomial roots, polynomials with pairs of…
Kuipers, J.
2012-06-01
New features of the symbolic algebra package Form 4 are discussed. Most importantly, these features include polynomial factorization and polynomial gcd computation. Examples of their use are shown. One of them is an exact version of Mincer which gives answers in terms of rational polynomials and 5 master integrals.
Interpolation and Polynomial Curve Fitting
Yang, Yajun; Gordon, Sheldon P.
2014-01-01
Two points determine a line. Three noncollinear points determine a quadratic function. Four points that do not lie on a lower-degree polynomial curve determine a cubic function. In general, n + 1 points uniquely determine a polynomial of degree n, presuming that they do not fall onto a polynomial of lower degree. The process of finding such a…
Parabolic refined invariants and Macdonald polynomials
Chuang, Wu-yen; Donagi, Ron; Pantev, Tony
2013-01-01
A string theoretic derivation is given for the conjecture of Hausel, Letellier, and Rodriguez-Villegas on the cohomology of character varieties with marked points. Their formula is identified with a refined BPS expansion in the stable pair theory of a local root stack, generalizing previous work of the first two authors in collaboration with G. Pan. Haiman's geometric construction for Macdonald polynomials is shown to emerge naturally in this context via geometric engineering. In particular this yields a new conjectural relation between Macdonald polynomials and refined local orbifold curve counting invariants. The string theoretic approach also leads to a new spectral cover construction for parabolic Higgs bundles in terms of holomorphic symplectic orbifolds.
Quantum Hurwitz numbers and Macdonald polynomials
Harnad, J
2015-01-01
Parametric families in the centre ${\\bf Z}({\\bf C}[S_n])$ of the group algebra of the symmetric group are obtained by identifying the indeterminates in the generating function for Macdonald polynomials as commuting Jucys-Murphy elements. Their eigenvalues provide coefficients in the double Schur function expansion of 2D Toda $\\tau$-functions of hypergeometric type. Expressing these in the basis of products of power sum symmetric functions, the coefficients may be interpreted geometrically as parametric families of quantum Hurwitz numbers, enumerating weighted branched coverings of the Riemann sphere. Combinatorially, they give quantum weighted sums over paths in the Cayley graph of $S_n$ generated by transpositions. Dual pairs of bases for the algebra of symmetric functions with respect to the scalar product in which the Macdonald polynomials are orthogonal provide both the geometrical and combinatorial significance of these quantum weighted enumerative invariants.
Polynomial interpolation methods for viscous flow calculations
Rubin, S. G.; Khosla, P. K.
1977-01-01
Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.
Schemes for Deterministic Polynomial Factoring
Ivanyos, Gábor; Saxena, Nitin
2008-01-01
In this work we relate the deterministic complexity of factoring polynomials (over finite fields) to certain combinatorial objects we call m-schemes. We extend the known conditional deterministic subexponential time polynomial factoring algorithm for finite fields to get an underlying m-scheme. We demonstrate how the properties of m-schemes relate to improvements in the deterministic complexity of factoring polynomials over finite fields assuming the generalized Riemann Hypothesis (GRH). In particular, we give the first deterministic polynomial time algorithm (assuming GRH) to find a nontrivial factor of a polynomial of prime degree n where (n-1) is a smooth number.
Complex Polynomial Vector Fields
DEFF Research Database (Denmark)
The two branches of dynamical systems, continuous and discrete, correspond to the study of differential equations (vector fields) and iteration of mappings respectively. In holomorphic dynamics, the systems studied are restricted to those described by holomorphic (complex analytic) functions...... vector fields. Since the class of complex polynomial vector fields in the plane is natural to consider, it is remarkable that its study has only begun very recently. There are numerous fundamental questions that are still open, both in the general classification of these vector fields, the decomposition...... of parameter spaces into structurally stable domains, and a description of the bifurcations. For this reason, the talk will focus on these questions for complex polynomial vector fields....
Directory of Open Access Journals (Sweden)
Hjalmar Rosengren
2006-12-01
Full Text Available We study multivariable Christoffel-Darboux kernels, which may be viewed as reproducing kernels for antisymmetric orthogonal polynomials, and also as correlation functions for products of characteristic polynomials of random Hermitian matrices. Using their interpretation as reproducing kernels, we obtain simple proofs of Pfaffian and determinant formulas, as well as Schur polynomial expansions, for such kernels. In subsequent work, these results are applied in combinatorics (enumeration of marked shifted tableaux and number theory (representation of integers as sums of squares.
Szabłowski, Paweł J
2010-01-01
We give equivalent forms of Askey-Wilson (AW) polynomials expressing them with a help of Al-Salam-Chihara polynomials. After restricting parameters of AW polynomials to complex conjugate pairs we give probabilistic interpretation of AW weight function and expand it in the series similar to Poisson-Mehler expansion formula and give its probabilistic interpretation. On the way (by setting certain parameter q to to 0) we get some formulae useful in rapidly developing so called 'free probability'.
A Characterization of Polynomials
DEFF Research Database (Denmark)
Andersen, Kurt Munk
1996-01-01
Given the problem:which functions f(x) are characterized by a relation of the form:f[x1,x2,...,xn]=h(x1+x2+...+xn), where n>1 and h(x) is a given function? Here f[x1,x2,...,xn] denotes the divided difference on n points x1,x2,...,xn of the function f(x).The answer is: f(x) is a polynomial of degree...
Some discrete multiple orthogonal polynomials
Arvesú, J.; Coussement, J.; van Assche, W.
2003-04-01
In this paper, we extend the theory of discrete orthogonal polynomials (on a linear lattice) to polynomials satisfying orthogonality conditions with respect to r positive discrete measures. First we recall the known results of the classical orthogonal polynomials of Charlier, Meixner, Kravchuk and Hahn (T.S. Chihara, An Introduction to Orthogonal Polynomials, Gordon and Breach, New York, 1978; R. Koekoek and R.F. Swarttouw, Reports of the Faculty of Technical Mathematics and Informatics No. 98-17, Delft, 1998; A.F. Nikiforov et al., Classical Orthogonal Polynomials of a Discrete Variable, Springer, Berlin, 1991). These polynomials have a lowering and raising operator, which give rise to a Rodrigues formula, a second order difference equation, and an explicit expression from which the coefficients of the three-term recurrence relation can be obtained. Then we consider r positive discrete measures and define two types of multiple orthogonal polynomials. The continuous case (Jacobi, Laguerre, Hermite, etc.) was studied by Van Assche and Coussement (J. Comput. Appl. Math. 127 (2001) 317-347) and Aptekarev et al. (Multiple orthogonal polynomials for classical weights, manuscript). The families of multiple orthogonal polynomials (of type II) that we will study have a raising operator and hence a Rodrigues formula. This will give us an explicit formula for the polynomials. Finally, there also exists a recurrence relation of order r+1 for these multiple orthogonal polynomials of type II. We compute the coefficients of the recurrence relation explicitly when r=2.
Spreading lengths of Hermite polynomials
Sánchez-Moreno, P; Manzano, D; Yáñez, R; 10.1016/j.cam.2009.09.043
2009-01-01
The Renyi, Shannon and Fisher spreading lengths of the classical or hypergeometric orthogonal polynomials, which are quantifiers of their distribution all over the orthogonality interval, are defined and investigated. These information-theoretic measures of the associated Rakhmanov probability density, which are direct measures of the polynomial spreading in the sense of having the same units as the variable, share interesting properties: invariance under translations and reflections, linear scaling and vanishing in the limit that the variable tends towards a given definite value. The expressions of the Renyi and Fisher lengths for the Hermite polynomials are computed in terms of the polynomial degree. The combinatorial multivariable Bell polynomials, which are shown to characterize the finite power of an arbitrary polynomial, play a relevant role for the computation of these information-theoretic lengths. Indeed these polynomials allow us to design an error-free computing approach for the entropic moments (w...
Oblivious Polynomial Evaluation
Institute of Scientific and Technical Information of China (English)
Hong-Da Li; Dong-Yao Ji; Deng-Guo Feng; Bao Li
2004-01-01
The problem of two-party oblivious polynomial evaluation(OPE)is studied,where one party(Alice)has a polynomial P(x)and the other party(Bob)with an input x wants to learn P(x)in such an oblivious way that Bob obtains P(x)without learning any additional information about P except what is implied by P(x)and Alice does not know Bob's input x.The former OPE protocols are based on an intractability assumption except for OT protocols.In fact,evaluating P(x)is equivalent to computing the product of the coefficient vectors(a0,...,an)and(1,...,xn).Using this idea,an efficient scale product protocol of two vectors is proposed first and then two OPE protocols are presented which do not need any other cryptographic assumption except for OT protocol.Compared with the existing OPE protocol,another characteristic of the proposed protocols is the degree of the polynomial is private.Another OPE protocol works in case of existence of untrusted third party.
Transition matrices for symmetric and quasisymmetric Hall-Littlewood polynomials
Loehr, Nicholas A; Warrington, Gregory S
2012-01-01
We introduce explicit combinatorial interpretations for the coefficients in some of the transition matrices between skew Hall-Littlewood polynomials P_lambda/mu(t), Hivert's quasisymmetric Hall-Littlewood polynomials G_gamma(t), and Gessel's fundamental and monomial quasisymmetric functions, F_alpha and M_beta. More specifically, we provide the following: (1) An expansion of the P_lambda in terms of the G_gamma, (2) expansions of the F_alpha and the M_beta in terms of the G_gamma and (3) an expansion of the P_lambda/mu in terms of the F_alpha. The F_alpha expansion of the P_lambda/mu is facilitated by introducing the set of starred tableaux.
Complex Polynomial Vector Fields
DEFF Research Database (Denmark)
Dias, Kealey
vector fields. Since the class of complex polynomial vector fields in the plane is natural to consider, it is remarkable that its study has only begun very recently. There are numerous fundamental questions that are still open, both in the general classification of these vector fields, the decomposition...... or meromorphic (allowing poles as singularities) functions. There already exists a well-developed theory for iterative holomorphic dynamical systems, and successful relations found between iteration theory and flows of vector fields have been one of the main motivations for the recent interest in holomorphic...
Symmetric functions and Hall polynomials
MacDonald, Ian Grant
1998-01-01
This reissued classic text is the acclaimed second edition of Professor Ian Macdonald's groundbreaking monograph on symmetric functions and Hall polynomials. The first edition was published in 1979, before being significantly expanded into the present edition in 1995. This text is widely regarded as the best source of information on Hall polynomials and what have come to be known as Macdonald polynomials, central to a number of key developments in mathematics and mathematical physics in the 21st century Macdonald polynomials gave rise to the subject of double affine Hecke algebras (or Cherednik algebras) important in representation theory. String theorists use Macdonald polynomials to attack the so-called AGT conjectures. Macdonald polynomials have been recently used to construct knot invariants. They are also a central tool for a theory of integrable stochastic models that have found a number of applications in probability, such as random matrices, directed polymers in random media, driven lattice gases, and...
Polynomial Regression on Riemannian Manifolds
Hinkle, Jacob; Fletcher, P Thomas; Joshi, Sarang
2012-01-01
In this paper we develop the theory of parametric polynomial regression in Riemannian manifolds and Lie groups. We show application of Riemannian polynomial regression to shape analysis in Kendall shape space. Results are presented, showing the power of polynomial regression on the classic rat skull growth data of Bookstein as well as the analysis of the shape changes associated with aging of the corpus callosum from the OASIS Alzheimer's study.
Witt Rings and Permutation Polynomials
Institute of Scientific and Technical Information of China (English)
Qifan Zhang
2005-01-01
Let p be a prime number. In this paper, the author sets up a canonical correspondence between polynomial functions over Z/p2Z and 3-tuples of polynomial functions over Z/pZ. Based on this correspondence, he proves and reproves some fundamental results on permutation polynomials mod pl. The main new result is the characterization of strong orthogonal systems over Z/p1Z.
Superoscillations with arbitrary polynomial shape
Chremmos, Ioannis; Fikioris, George
2015-07-01
We present a method for constructing superoscillatory functions the superoscillatory part of which approximates a given polynomial with arbitrarily small error in a fixed interval. These functions are obtained as the product of the polynomial with a sufficiently flat, bandlimited envelope function whose Fourier transform has at least N-1 continuous derivatives and an Nth derivative of bounded variation, N being the order of the polynomial. Polynomials of arbitrarily high order can be approximated if the Fourier transform of the envelope is smooth, i.e. a bump function.
Digital terrain modeling with the Chebyshev polynomials
Florinsky, I V
2015-01-01
Mathematical problems of digital terrain analysis include interpolation of digital elevation models (DEMs), DEM generalization and denoising, and computation of morphometric variables by calculation of partial derivatives of elevation. Traditionally, these procedures are based on numerical treatments of two-variable discrete functions of elevation. We developed a spectral analytical method and algorithm based on high-order orthogonal expansions using the Chebyshev polynomials of the first kind with the subsequent Fejer summation. The method and algorithm are intended for DEM analytical treatment, such as, DEM global approximation, denoising, and generalization as well as computation of morphometric variables by analytical calculation of partial derivatives. To test the method and algorithm, we used a DEM of the Northern Andes including 230,880 points (the elevation matrix 480 $\\times$ 481). DEMs were reconstructed with 480, 240, 120, 60, and 30 expansion coefficients. The first and second partial derivatives ...
Derivations and identities for Kravchuk polynomials
Bedratyuk, Leonid
2012-01-01
We introduce the notion of Kravchuk derivations of the polynomial algebra. We prove that any element of the kernel of the derivation gives a polynomial identity satisfied by the Kravchuk polynomials. Also, we prove that any kernel element of the basic Weitzenb\\"ok derivations yields a polynomial identity satisfied by the Kravchuk polynomials. We describe the corresponding intertwining maps.
Factorization of colored knot polynomials at roots of unity
Kononov, Ya.; Morozov, A.
2015-07-01
HOMFLY polynomials are the Wilson-loop averages in Chern-Simons theory and depend on four variables: the closed line (knot) in 3d space-time, representation R of the gauge group SU (N) and exponentiated coupling constant q. From analysis of a big variety of different knots we conclude that at q, which is a 2m-th root of unity, q2m = 1, HOMFLY polynomials in symmetric representations [ r ] satisfy recursion identity: Hr+m =Hr ṡHm for any A =qN, which is a generalization of the property Hr = H1r for special polynomials at m = 1. We conjecture a further generalization to arbitrary representation R, which, however, is checked only for torus knots. Next, Kashaev polynomial, which arises from HR at q2 = e 2 πi / | R |, turns equal to the special polynomial with A substituted by A| R |, provided R is a single-hook representations (including arbitrary symmetric) - what provides a q - A dual to the similar property of Alexander polynomial. All this implies non-trivial relations for the coefficients of the differential expansions, which are believed to provide reasonable coordinates in the space of knots - existence of such universal relations means that these variables are still not unconstrained.
Factorization of colored knot polynomials at roots of unity
Directory of Open Access Journals (Sweden)
Ya. Kononov
2015-07-01
Full Text Available HOMFLY polynomials are the Wilson-loop averages in Chern–Simons theory and depend on four variables: the closed line (knot in 3d space–time, representation R of the gauge group SU(N and exponentiated coupling constant q. From analysis of a big variety of different knots we conclude that at q, which is a 2m-th root of unity, q2m=1, HOMFLY polynomials in symmetric representations [r] satisfy recursion identity: Hr+m=Hr⋅Hm for any A=qN, which is a generalization of the property Hr=H1r for special polynomials at m=1. We conjecture a further generalization to arbitrary representation R, which, however, is checked only for torus knots. Next, Kashaev polynomial, which arises from HR at q2=e2πi/|R|, turns equal to the special polynomial with A substituted by A|R|, provided R is a single-hook representations (including arbitrary symmetric – what provides a q−A dual to the similar property of Alexander polynomial. All this implies non-trivial relations for the coefficients of the differential expansions, which are believed to provide reasonable coordinates in the space of knots – existence of such universal relations means that these variables are still not unconstrained.
Some New Formulae for Genocchi Numbers and Polynomials Involving Bernoulli and Euler Polynomials
Directory of Open Access Journals (Sweden)
Serkan Araci
2014-01-01
Full Text Available We give some new formulae for product of two Genocchi polynomials including Euler polynomials and Bernoulli polynomials. Moreover, we derive some applications for Genocchi polynomials to study a matrix formulation.
Parallel Construction of Irreducible Polynomials
DEFF Research Database (Denmark)
Frandsen, Gudmund Skovbjerg
Let arithmetic pseudo-NC^k denote the problems that can be solved by log space uniform arithmetic circuits over the finite prime field GF(p) of depth O(log^k (n + p)) and size polynomial in (n + p). We show that the problem of constructing an irreducible polynomial of specified degree over GF(p) ...
Complex Polynomial Vector Fields
DEFF Research Database (Denmark)
Dias, Kealey
The two branches of dynamical systems, continuous and discrete, correspond to the study of differential equations (vector fields) and iteration of mappings respectively. In holomorphic dynamics, the systems studied are restricted to those described by holomorphic (complex analytic) functions...... or meromorphic (allowing poles as singularities) functions. There already exists a well-developed theory for iterative holomorphic dynamical systems, and successful relations found between iteration theory and flows of vector fields have been one of the main motivations for the recent interest in holomorphic...... vector fields. Since the class of complex polynomial vector fields in the plane is natural to consider, it is remarkable that its study has only begun very recently. There are numerous fundamental questions that are still open, both in the general classification of these vector fields, the decomposition...
The number of polynomial solutions of polynomial Riccati equations
Gasull, Armengol; Torregrosa, Joan; Zhang, Xiang
2016-11-01
Consider real or complex polynomial Riccati differential equations a (x) y ˙ =b0 (x) +b1 (x) y +b2 (x)y2 with all the involved functions being polynomials of degree at most η. We prove that the maximum number of polynomial solutions is η + 1 (resp. 2) when η ≥ 1 (resp. η = 0) and that these bounds are sharp. For real trigonometric polynomial Riccati differential equations with all the functions being trigonometric polynomials of degree at most η ≥ 1 we prove a similar result. In this case, the maximum number of trigonometric polynomial solutions is 2η (resp. 3) when η ≥ 2 (resp. η = 1) and, again, these bounds are sharp. Although the proof of both results has the same starting point, the classical result that asserts that the cross ratio of four different solutions of a Riccati differential equation is constant, the trigonometric case is much more involved. The main reason is that the ring of trigonometric polynomials is not a unique factorization domain.
Prime power polynomial maps over finite fields
Berson, Joost
2012-01-01
We consider polynomial maps described by so-called prime power polynomials. These polynomials are defined using a fixed power of a prime number, say q. Considering invertible polynomial maps of this type over a characteristic zero field, we will only obtain (up to permutation of the variables) triangular maps, which are the most basic examples of polynomial automorphisms. However, over the finite field F_q automorphisms of this type have (in general) an entirely different structure. Namely, we will show that the prime power polynomial maps over F_q are in one-to-one correspondence with matrices having coefficients in a univariate polynomial ring over F_q. Furthermore, composition of polynomial maps translates to matrix multiplication, implying that invertible prime power polynomial maps correspond to invertible matrices. This alternate description of the prime power polynomial automorphism subgroup leads to the solution of many famous conjectures for this kind of polynomials and polynomial maps.
Hadamard Factorization of Stable Polynomials
Loredo-Villalobos, Carlos Arturo; Aguirre-Hernández, Baltazar
2011-11-01
The stable (Hurwitz) polynomials are important in the study of differential equations systems and control theory (see [7] and [19]). A property of these polynomials is related to Hadamard product. Consider two polynomials p,q ∈ R[x]:p(x) = anxn+an-1xn-1+...+a1x+a0q(x) = bmx m+bm-1xm-1+...+b1x+b0the Hadamard product (p × q) is defined as (p×q)(x) = akbkxk+ak-1bk-1xk-1+...+a1b1x+a0b0where k = min(m,n). Some results (see [16]) shows that if p,q ∈R[x] are stable polynomials then (p×q) is stable, also, i.e. the Hadamard product is closed; however, the reciprocal is not always true, that is, not all stable polynomial has a factorization into two stable polynomials the same degree n, if n> 4 (see [15]).In this work we will give some conditions to Hadamard factorization existence for stable polynomials.
Automatic differentiation for Fourier series and the radii polynomial approach
Lessard, Jean-Philippe; Mireles James, J. D.; Ransford, Julian
2016-11-01
In this work we develop a computer-assisted technique for proving existence of periodic solutions of nonlinear differential equations with non-polynomial nonlinearities. We exploit ideas from the theory of automatic differentiation in order to formulate an augmented polynomial system. We compute a numerical Fourier expansion of the periodic orbit for the augmented system, and prove the existence of a true solution nearby using an a-posteriori validation scheme (the radii polynomial approach). The problems considered here are given in terms of locally analytic vector fields (i.e. the field is analytic in a neighborhood of the periodic orbit) hence the computer-assisted proofs are formulated in a Banach space of sequences satisfying a geometric decay condition. In order to illustrate the use and utility of these ideas we implement a number of computer-assisted existence proofs for periodic orbits of the Planar Circular Restricted Three-Body Problem (PCRTBP).
Institute of Scientific and Technical Information of China (English)
GU Chuan-qing; PAN Bao-zhen; WU Bei-bei
2006-01-01
To solve Fredholm integral equations of the second kind, a generalized linear functional is introduced and a new function-valued Padé-type approximation is defined.By means of the power series expansion of the solution, this method can construct an approximate solution to solve the given integral equation. On the basis of the orthogonal polynomials, two useful determinant expressions of the numerator polynomial and the denominator polynomial for padé-type approximation are explicitly given.
Polynomial Regressions and Nonsense Inference
Directory of Open Access Journals (Sweden)
Daniel Ventosa-Santaulària
2013-11-01
Full Text Available Polynomial specifications are widely used, not only in applied economics, but also in epidemiology, physics, political analysis and psychology, just to mention a few examples. In many cases, the data employed to estimate such specifications are time series that may exhibit stochastic nonstationary behavior. We extend Phillips’ results (Phillips, P. Understanding spurious regressions in econometrics. J. Econom. 1986, 33, 311–340. by proving that an inference drawn from polynomial specifications, under stochastic nonstationarity, is misleading unless the variables cointegrate. We use a generalized polynomial specification as a vehicle to study its asymptotic and finite-sample properties. Our results, therefore, lead to a call to be cautious whenever practitioners estimate polynomial regressions.
Locally tame plane polynomial automorphisms
Berson, Joost; Furter, Jean-Philippe; Maubach, Stefan
2010-01-01
For automorphisms of a polynomial ring in two variables over a domain R, we show that local tameness implies global tameness provided that every 2-generated invertible R-module is free. We give many examples illustrating this property.
SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos
Ahlfeld, R.; Belkouchi, B.; Montomoli, F.
2016-09-01
A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrix is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10
Stochastic Estimation via Polynomial Chaos
2015-10-01
TΨ is a vector with P+1 elements. With these dimensions, (29) is solvable by standard numerical linear algebra techniques. The specific matrix...initial conditions for partial differential equations. Here, the elementary theory of the polynomial chaos is presented followed by the details of a...the elementary theory of the polynomial chaos is presented followed by the details of a number of example calculations where the statistical mean and
Properties of Leach-Flessas-Gorringe polynomials
Pursey, D. L.
1990-09-01
A generating function is obtained for the polynomials recently introduced by Leach, Flessas, and Gorringe [J. Math. Phys. 30, 406 (1989)], and is then used to relate the Leach-Flessas-Gorringe (or LFG) polynomials to Hermite polynomials. The generating function is also used to express a number of integrals involving the LFG polynomials as finite sums of parabolic cylinder functions.
Birth-death processes and associated polynomials
Doorn, van Erik A.
2003-01-01
We consider birth-death processes on the nonnegative integers and the corresponding sequences of orthogonal polynomials called birth-death polynomials. The sequence of associated polynomials linked with a sequence of birth-death polynomials and its orthogonalizing measure can be used in the analysis
On the Hermite-Apostol-Genocchi Polynomials
Kurt, Veli; Kurt, Burak
2011-09-01
In this study, we introduce and investigate the Hermite-Apostol-Genocchi polynomials by means of a suitable generating function. We establish several interesting properties of these general polynomials. Also, we prove two theorems between 2-dimensional Hermite polynomials and Hermite-Apostol-Genocchi polynomials.
Uniqueness and Zeros of -Shift Difference Polynomials
Indian Academy of Sciences (India)
Kai Liu; Xin-Ling Liu; Ting-Bin Cao
2011-08-01
In this paper, we consider the zero distributions of -shift difference polynomials of meromorphic functions with zero order, and obtain two theorems that extend the classical Hayman results on the zeros of differential polynomials to -shift difference polynomials. We also investigate the uniqueness problem of -shift difference polynomials that share a common value.
Sparse polynomial surrogates for aerodynamic computations with random inputs
Savin, Eric; Peter, Jacques
2015-01-01
This paper deals with some of the methodologies used to construct polynomial surrogate models based on generalized polynomial chaos (gPC) expansions for applications to uncertainty quantification (UQ) in aerodynamic computations. A core ingredient in gPC expansions is the choice of a dedicated sampling strategy, so as to define the most significant scenarios to be considered for the construction of such metamodels. A desirable feature of the proposed rules shall be their ability to handle several random inputs simultaneously. Methods to identify the relative "importance" of those variables or uncertain data shall be ideally considered as well. The present work is more particularly dedicated to the development of sampling strategies based on sparsity principles. Sparse multi-dimensional cubature rules based on general one-dimensional Gauss-Jacobi-type quadratures are first addressed. These sets are non nested, but they are well adapted to the probability density functions with compact support for the random in...
Multi-particle dynamical systems and polynomials
Demina, Maria V.; Kudryashov, Nikolai A.
2016-05-01
Polynomial dynamical systems describing interacting particles in the plane are studied. A method replacing integration of a polynomial multi-particle dynamical system by finding polynomial solutions of partial differential equations is introduced. The method enables one to integrate a wide class of polynomial multi-particle dynamical systems. The general solutions of certain dynamical systems related to linear second-order partial differential equations are found. As a by-product of our results, new families of orthogonal polynomials are derived.
DIFFERENTIAL CROSS SECTION ANALYSIS IN KAON PHOTOPRODUCTION USING ASSOCIATED LEGENDRE POLYNOMIALS
Energy Technology Data Exchange (ETDEWEB)
P. T. P. HUTAURUK, D. G. IRELAND, G. ROSNER
2009-04-01
Angular distributions of differential cross sections from the latest CLAS data sets,6 for the reaction γ + p→K+ + Λ have been analyzed using associated Legendre polynomials. This analysis is based upon theoretical calculations in Ref. 1 where all sixteen observables in kaon photoproduction can be classified into four Legendre classes. Each observable can be described by an expansion of associated Legendre polynomial functions. One of the questions to be addressed is how many associated Legendre polynomials are required to describe the data. In this preliminary analysis, we used data models with different numbers of associated Legendre polynomials. We then compared these models by calculating posterior probabilities of the models. We found that the CLAS data set needs no more than four associated Legendre polynomials to describe the differential cross section data. In addition, we also show the extracted coefficients of the best model.
Jacobi-Sobolev Orthogonal Polynomials: Asymptotics for N-Coherence of Measures
Directory of Open Access Journals (Sweden)
Marcellán Francisco
2011-01-01
Full Text Available Let us introduce the Sobolev-type inner product , where and , , with and for all A Mehler-Heine-type formula and the inner strong asymptotics on as well as some estimates for the polynomials orthogonal with respect to the above Sobolev inner product are obtained. Necessary conditions for the norm convergence of Fourier expansions in terms of such Sobolev orthogonal polynomials are given.
Range Image Flow using High-Order Polynomial Expansion
2013-09-01
give a special thanks to Dr. Steve Hobbs for his help with the high-order tensor calculations. MATLAB ® is a registered...that using multiple spatial scales and past information improve the final flow estimation, as we would expect. Also, we will port the MATLAB R...taken column- wise and diagonalized, and f is the range image data, taken column-wise. The values of these weights for a Velodyne R© and Odetic lidar
Asymptotic estimates for Apostol-Bernoulli and Apostol-Euler polynomials
Navas, Luis M; Varona, Juan L
2011-01-01
We analyze the asymptotic behavior of the Apostol-Bernoulli polynomials $\\mathcal{B}_{n}(x;\\lambda)$ in detail. The starting point is their Fourier series on $[0,1]$ which, it is shown, remains valid as an asymptotic expansion over compact subsets of the complex plane. This is used to determine explicit estimates on the constants in the approximation, and also to analyze oscillatory phenomena which arise in certain cases. These results are transferred to the Apostol-Euler polynomials $\\mathcal{E}_{n}(x;\\lambda)$ via a simple relation linking them to the Apostol-Bernoulli polynomials.
Plain Polynomial Arithmetic on GPU
Anisul Haque, Sardar; Moreno Maza, Marc
2012-10-01
As for serial code on CPUs, parallel code on GPUs for dense polynomial arithmetic relies on a combination of asymptotically fast and plain algorithms. Those are employed for data of large and small size, respectively. Parallelizing both types of algorithms is required in order to achieve peak performances. In this paper, we show that the plain dense polynomial multiplication can be efficiently parallelized on GPUs. Remarkably, it outperforms (highly optimized) FFT-based multiplication up to degree 212 while on CPU the same threshold is usually at 26. We also report on a GPU implementation of the Euclidean Algorithm which is both work-efficient and runs in linear time for input polynomials up to degree 218 thus showing the performance of the GCD algorithm based on systolic arrays.
Derivations and identities for Fibonacci and Lucas polynomials
Bedratyuk, Leonid
2012-01-01
We introduce the notion of Fibonacci and Lucas derivations of the polynomial algebras and prove that any element of kernel of the derivations defines a polynomial identity for the Fibonacci and Lucas polynomials. Also, we prove that any polynomial identity for Appel polynomial yields a polynomial identity for the Fibonacci and Lucas polynomials and describe the corresponding intertwining maps.
Tree modules and counting polynomials
Kinser, Ryan
2011-01-01
We give a formula for counting tree modules for the quiver S_g with g loops and one vertex in terms of tree modules on its universal cover. This formula, along with work of Helleloid and Rodriguez-Villegas, is used to show that the number of d-dimensional tree modules for S_g is polynomial in g with the same degree and leading coefficient as the counting polynomial A_{S_g}(d, q) for absolutely indecomposables over F_q, evaluated at q=1.
Orthogonal Polynomials and their Applications
Dehesa, Jesús; Marcellan, Francisco; Francia, José; Vinuesa, Jaime
1988-01-01
The Segovia meeting set out to stimulate an intensive exchange of ideas between experts in the area of orthogonal polynomials and its applications, to present recent research results and to reinforce the scientific and human relations among the increasingly international community working in orthogonal polynomials. This volume contains original research papers as well as survey papers about fundamental questions in the field (Nevai, Rakhmanov & López) and its relationship with other fields such as group theory (Koornwinder), Padé approximation (Brezinski), differential equations (Krall, Littlejohn) and numerical methods (Rivlin).
Symbolic computation of Appell polynomials using Maple
Directory of Open Access Journals (Sweden)
H. Alkahby
2001-07-01
Full Text Available This work focuses on the symbolic computation of Appell polynomials using the computer algebra system Maple. After describing the traditional approach of constructing Appell polynomials, the paper examines the operator method of constructing the same Appell polynomials. The operator approach enables us to express the Appell polynomial as Bessel function whose coefficients are Euler and Bernuolli numbers. We have also constructed algorithms using Maple to compute Appell polynomials based on the methods we have described. The achievement is the construction of Appell polynomials for any function of bounded variation.
Two new triangles of $q$-integers via $q$-Eulerian polynomials of type $A$ and $B$
Han, Guoniu; Zeng, Jiang
2012-01-01
The classical Eulerian polynomials can be expanded in the basis $t^{k-1}(1+t)^{n+1-2k}$ ($1\\leq k\\leq\\lfloor (n+1)/2\\rfloor$) with positive integral coefficients. This formula implies both the symmetry and the unimodality of the Eulerian polynomials. In this paper, we prove a $q$-analogue of this expansion for Carlitz's $q$-Eulerian polynomials as well as a similar formula for Chow-Gessel's $q$-Eulerian polynomials of type $B$. We shall give some applications of these two formulae, which involve two new sequences of polynomials in the variable $q$ with positive integral coefficients. An open problem is to give a combinatorial interpretation for these polynomials.
Polynomial Regressions and Nonsense Inference
DEFF Research Database (Denmark)
Ventosa-Santaulària, Daniel; Rodríguez-Caballero, Carlos Vladimir
Polynomial specifications are widely used, not only in applied economics, but also in epidemiology, physics, political analysis, and psychology, just to mention a few examples. In many cases, the data employed to estimate such estimations are time series that may exhibit stochastic nonstationary ...
Uniform approximation by (quantum) polynomials
Drucker, A.; de Wolf, R.
2011-01-01
We show that quantum algorithms can be used to re-prove a classical theorem in approximation theory, Jackson's Theorem, which gives a nearly-optimal quantitative version of Weierstrass's Theorem on uniform approximation of continuous functions by polynomials. We provide two proofs, based respectivel
On Modular Counting with Polynomials
DEFF Research Database (Denmark)
Hansen, Kristoffer Arnsfelt
2006-01-01
For any integers m and l, where m has r sufficiently large (depending on l) factors, that are powers of r distinct primes, we give a construction of a (symmetric) polynomial over Z_m of degree O(\\sqrt n) that is a generalized representation (commonly also called weak representation) of the MODl...
Two polynomial division inequalities in
Directory of Open Access Journals (Sweden)
Goetgheluck P
1998-01-01
Full Text Available This paper is a first attempt to give numerical values for constants and , in classical estimates and where is an algebraic polynomial of degree at most and denotes the -metric on . The basic tools are Markov and Bernstein inequalities.
Global Polynomial Kernel Hazard Estimation
DEFF Research Database (Denmark)
Hiabu, Munir; Miranda, Maria Dolores Martínez; Nielsen, Jens Perch
2015-01-01
This paper introduces a new bias reducing method for kernel hazard estimation. The method is called global polynomial adjustment (GPA). It is a global correction which is applicable to any kernel hazard estimator. The estimator works well from a theoretical point of view as it asymptotically redu...
Polynomial J-spectral factorization
Kwakernaak, Huibert; Sebek, Michael
1994-01-01
Several algorithms are presented for the J-spectral factorization of a para-Hermitian polynomial matrix. The four algorithms that are discussed are based on diagonalization, successive factor extraction, interpolation, and the solution of an algebraic Riccati equation, respectively. The paper includ
Global Polynomial Kernel Hazard Estimation
DEFF Research Database (Denmark)
Hiabu, Munir; Miranda, Maria Dolores Martínez; Nielsen, Jens Perch;
2015-01-01
This paper introduces a new bias reducing method for kernel hazard estimation. The method is called global polynomial adjustment (GPA). It is a global correction which is applicable to any kernel hazard estimator. The estimator works well from a theoretical point of view as it asymptotically...
A New Generalisation of Macdonald Polynomials
Garbali, Alexandr; de Gier, Jan; Wheeler, Michael
2017-01-01
We introduce a new family of symmetric multivariate polynomials, whose coefficients are meromorphic functions of two parameters (q, t) and polynomial in a further two parameters (u, v). We evaluate these polynomials explicitly as a matrix product. At u = v = 0 they reduce to Macdonald polynomials, while at q = 0, u = v = s they recover a family of inhomogeneous symmetric functions originally introduced by Borodin.
A Summation Formula for Macdonald Polynomials
de Gier, Jan; Wheeler, Michael
2016-03-01
We derive an explicit sum formula for symmetric Macdonald polynomials. Our expression contains multiple sums over the symmetric group and uses the action of Hecke generators on the ring of polynomials. In the special cases {t = 1} and {q = 0}, we recover known expressions for the monomial symmetric and Hall-Littlewood polynomials, respectively. Other specializations of our formula give new expressions for the Jack and q-Whittaker polynomials.
An analysis on the inversion of polynomials
M. F. González-Cardel; R. Díaz-Uribe
2006-01-01
In this work the application and the intervals of validity of an inverse polynomial, according to the method proposed by Arfken [1] for the inversion of series, is analyzed. It is shown that, for the inverse polynomial there exists a restricted domain whose longitude depends on the magnitude of the acceptable error when the inverse polynomial is used to approximate the inverse function of the original polynomial. A method for calculating the error of the approximation and its use in determini...
Application of Chebyshev Polynomial to simulated modeling
Institute of Scientific and Technical Information of China (English)
CHI Hai-hong; LI Dian-pu
2006-01-01
Chebyshev polynomial is widely used in many fields, and used usually as function approximation in numerical calculation. In this paper, Chebyshev polynomial expression of the propeller properties across four quadrants is given at first, then the expression of Chebyshev polynomial is transformed to ordinary polynomial for the need of simulation of propeller dynamics. On the basis of it,the dynamical models of propeller across four quadrants are given. The simulation results show the efficiency of mathematical model.
Positive trigonometric polynomials and signal processing applications
Dumitrescu, Bogdan
2007-01-01
Presents the results on positive trigonometric polynomials within a unitary framework; the theoretical results obtained partly from the general theory of real polynomials, partly from self-sustained developments. This book provides information on the theory of sum-of-squares trigonometric polynomials in two parts: theory and applications.
Lattice Platonic Solids and their Ehrhart polynomial
Ionascu, Eugen J
2011-01-01
First, we calculate the Ehrhart polynomial associated to an arbitrary cube with integer coordinates for its vertices. Then, we use this result to derive relationships between the Ehrhart polynomials for regular lattice tetrahedrons and those for regular lattice octahedrons. These relations allow one to reduce the calculation of these polynomials to only one coefficient.
General Eulerian Numbers and Eulerian Polynomials
Directory of Open Access Journals (Sweden)
Tingyao Xiong
2013-01-01
Full Text Available We will generalize the definitions of Eulerian numbers and Eulerian polynomials to general arithmetic progressions. Under the new definitions, we have been successful in extending several well-known properties of traditional Eulerian numbers and polynomials to the general Eulerian polynomials and numbers.
Energy Technology Data Exchange (ETDEWEB)
Vinet, Luc [Universite de Montreal, PO Box 6128, Station Centre-ville, Montreal QC H3C 3J7 (Canada); Zhedanov, Alexei [Donetsk Institute for Physics and Technology, Donetsk 83114 (Ukraine)
2009-10-30
We construct new families of elliptic solutions of the restricted Toda chain. The main tool is a special (so-called Stieltjes) ansatz for the moments of corresponding orthogonal polynomials. We show that the moments thus obtained are related to three types of Lame polynomials. The corresponding orthogonal polynomials can be considered as a generalization of the Stieltjes-Carlitz elliptic polynomials.
Frobenious-Euler Type Polynomials Related to Hermite-Bernoulli Polynomials
Kurt, Burak; Simsek, Yilmaz
2011-09-01
The aim of this paper is to define and investigate a new generating functions of the Frobenious-Euler polynomials and numbers. We establish some fundamental properties of these numbers and polynomials. We also derive relationship between these polynomials and Hermite-Apostol-Bernoulli polynomials and numbers. We also give some remarks and applications.
HOMFLY polynomials in representation [3, 1] for 3-strand braids
Mironov, A.; Morozov, A.; Morozov, An.; Sleptsov, A.
2016-09-01
This paper is a new step in the project of systematic description of colored knot polynomials started in [1]. In this paper, we managed to explicitly find the inclusive Racah matrix, i.e. the whole set of mixing matrices in channels R ⊗3 -→ Q with all possible Q, for R = [3 , 1]. The calculation is made possible by the use of a newly-developed efficient highest-weight method, still it remains tedious. The result allows one to evaluate and investigate [3 , 1]-colored polynomials for arbitrary 3-strand knots, and this confirms many previous conjectures on various factorizations, universality, and differential expansions. We consider in some detail the next-to-twist-knots three-strand family ( n, -1 | 1 , -1) and deduce its colored HOMFLY. Also confirmed and clarified is the eigenvalue hypothesis for the Racah matrices, which promises to provide a shortcut to generic formulas for arbitrary representations.
Knot polynomials in the first non-symmetric representation
Anokhina, A.; Mironov, A.; Morozov, A.; Morozov, And.
2014-05-01
We describe the explicit form and the hidden structure of the answer for the HOMFLY polynomial for the figure-8 and some other 3-strand knots in representation [21]. This is the first result for non-torus knots beyond (anti)symmetric representations, and its evaluation is far more complicated. We provide a whole variety of different arguments, allowing one to guess the answer for the figure-8 knot, which can be also partly used in more complicated situations. Finally we report the result of exact calculation for figure-8 and some other 3-strand knots based on the previously developed sophisticated technique of multi-strand calculations. We also discuss a formula for the superpolynomial in representation [21] for the figure-8 knot, which heavily relies on the conjectural form of superpolynomial expansion nearby the special polynomial point. Generalizations and details will be presented elsewhere.
Normal BGG solutions and polynomials
Cap, A; Hammerl, M
2012-01-01
First BGG operators are a large class of overdetermined linear differential operators intrinsically associated to a parabolic geometry on a manifold. The corresponding equations include those controlling infinitesimal automorphisms, higher symmetries, and many other widely studied PDE of geometric origin. The machinery of BGG sequences also singles out a subclass of solutions called normal solutions. These correspond to parallel tractor fields and hence to (certain) holonomy reductions of the canonical normal Cartan connection. Using the normal Cartan connection, we define a special class of local frames for any natural vector bundle associated to a parabolic geometry. We then prove that the coefficient functions of any normal solution of a first BGG operator with respect to such a frame are polynomials in the normal coordinates of the parabolic geometry. A bound on the degree of these polynomials in terms of representation theory data is derived. For geometries locally isomorphic to the homogeneous model of ...
BSDEs with polynomial growth generators
Directory of Open Access Journals (Sweden)
Philippe Briand
2000-01-01
Full Text Available In this paper, we give existence and uniqueness results for backward stochastic differential equations when the generator has a polynomial growth in the state variable. We deal with the case of a fixed terminal time, as well as the case of random terminal time. The need for this type of extension of the classical existence and uniqueness results comes from the desire to provide a probabilistic representation of the solutions of semilinear partial differential equations in the spirit of a nonlinear Feynman-Kac formula. Indeed, in many applications of interest, the nonlinearity is polynomial, e.g, the Allen-Cahn equation or the standard nonlinear heat and Schrödinger equations.
Leont'ev, V. K.
2015-11-01
A pseudo-Boolean function is an arbitrary mapping of the set of binary n-tuples to the real line. Such functions are a natural generalization of classical Boolean functions and find numerous applications in various applied studies. Specifically, the Fourier transform of a Boolean function is a pseudo-Boolean function. A number of facts associated with pseudo-Boolean polynomials are presented, and their applications to well-known discrete optimization problems are described.
Weak lensing tomography with orthogonal polynomials
Schaefer, Bjoern Malte
2011-01-01
The topic of this article is weak cosmic shear tomography where the line of sight-weighting is carried out with a set of specifically constructed orthogonal polynomials, dubbed TaRDiS (Tomography with orthogonAl Radial Distance polynomIal Systems). We investigate the properties of these polynomials and employ weak convergence spectra, which have been obtained by weighting with these polynomials, for the estimation of cosmological parameters. We quantify their power in constraining parameters in a Fisher-matrix technique and demonstrate how each polynomial projects out statistically independent information, and how the combination of multiple polynomials lifts degeneracies. The assumption of a reference cosmology is needed for the construction of the polynomials, and as a last point we investigate how errors in the construction with a wrong cosmological model propagate to misestimates in cosmological parameters. TaRDiS performs on a similar level as traditional tomographic methods and some key features of tomo...
Weak lensing tomography with orthogonal polynomials
Schäfer, Björn Malte; Heisenberg, Lavinia
2012-07-01
The topic of this paper is weak cosmic shear tomography where the line-of-sight weighting is carried out with a set of specifically constructed orthogonal polynomials, dubbed Tomography with Orthogonal Radial Distance Polynomial Systems (TaRDiS). We investigate the properties of these polynomials and employ weak convergence spectra, which have been obtained by weighting with these polynomials, for the estimation of cosmological parameters. We quantify their power in constraining parameters in a Fisher matrix technique and demonstrate how each polynomial projects out statistically independent information, and how the combination of multiple polynomials lifts degeneracies. The assumption of a reference cosmology is needed for the construction of the polynomials, and as a last point we investigate how errors in the construction with a wrong cosmological model propagate to misestimates in cosmological parameters. TaRDiS performs on a similar level as traditional tomographic methods and some key features of tomography are made easier to understand.
On Ternary Inclusion-Exclusion Polynomials
Bachman, Gennady
2010-01-01
Taking a combinatorial point of view on cyclotomic polynomials leads to a larger class of polynomials we shall call the inclusion-exclusion polynomials. This gives a more appropriate setting for certain types of questions about the coefficients of these polynomials. After establishing some basic properties of inclusion-exclusion polynomials we turn to a detailed study of the structure of ternary inclusion-exclusion polynomials. The latter subclass is exemplified by cyclotomic polynomials $\\Phi_{pqr}$, where $p
Stable piecewise polynomial vector fields
Directory of Open Access Journals (Sweden)
Claudio Pessoa
2012-09-01
Full Text Available Let $N={y>0}$ and $S={y<0}$ be the semi-planes of $mathbb{R}^2$ having as common boundary the line $D={y=0}$. Let $X$ and $Y$ be polynomial vector fields defined in $N$ and $S$, respectively, leading to a discontinuous piecewise polynomial vector field $Z=(X,Y$. This work pursues the stability and the transition analysis of solutions of $Z$ between $N$ and $S$, started by Filippov (1988 and Kozlova (1984 and reformulated by Sotomayor-Teixeira (1995 in terms of the regularization method. This method consists in analyzing a one parameter family of continuous vector fields $Z_{epsilon}$, defined by averaging $X$ and $Y$. This family approaches $Z$ when the parameter goes to zero. The results of Sotomayor-Teixeira and Sotomayor-Machado (2002 providing conditions on $(X,Y$ for the regularized vector fields to be structurally stable on planar compact connected regions are extended to discontinuous piecewise polynomial vector fields on $mathbb{R}^2$. Pertinent genericity results for vector fields satisfying the above stability conditions are also extended to the present case. A procedure for the study of discontinuous piecewise vector fields at infinity through a compactification is proposed here.
Projected Hartree Fock Theory as a Polynomial Similarity Transformation Theory of Single Excitations
Qiu, Yiheng; Scuseria, Gustavo E
2016-01-01
Spin-projected Hartree-Fock is introduced as a particle-hole excitation ansatz over a symmetry-adapted reference determinant. Remarkably, this expansion has an analytic expression that we were able to decipher. While the form of the polynomial expansion is universal, the excitation amplitudes need to be optimized. This is equivalent to the optimization of orbitals in the conventional projected Hartree-Fock framework of non-orthogonal determinants. Using the inverse of the particle-hole expansion, we similarity transform the Hamiltonian in a coupled-cluster style theory. The left eigenvector of the non-hermitian Hamiltonian is constructed in a similar particle-hole expansion fashion, and we show that to numerically reproduce variational projected Hartree-Fock results, one needs as many pair excitations in the bra as the number of strongly correlated entangled pairs in the system. This single-excitation polynomial similarity transformation theory is an alternative to our recently presented double excitation the...
New optical field operator expansion in number state representation
Institute of Scientific and Technical Information of China (English)
Yang Yang; Fan Hong-Yi
2013-01-01
By virtue of the Weyl ordering method,we find a new formalism of optical field operator expansion in number state representation.Miscellaneous optical fields' (coherent state,squeezed field,Wigner operator,etc.) new expansions are therefore exhibited.Some new generating functions of special polynomials are derived herewith.
Institute of Scientific and Technical Information of China (English)
无
2007-01-01
In this paper, we mainly study the relation of two cyclically reduced words w and w' on the condition they have the same trace polynomial (i.e., tr w= tr w' ). By defining an equivalence relation through such operators on words as inverse, cyclically left shift, and mirror, it is straightforward to get that w ～ w' implies tr w = tr w'. We show by a counter example that tr w = tr w' does not imply w ～ w'. And in two special cases, we prove that tr w = tr w' if and only if w ～ w'.
The stable computation of formal orthogonal polynomials
Beckermann, Bernhard
1996-12-01
For many applications - such as the look-ahead variants of the Lanczos algorithm - a sequence of formal (block-)orthogonal polynomials is required. Usually, one generates such a sequence by taking suitable polynomial combinations of a pair of basis polynomials. These basis polynomials are determined by a look-ahead generalization of the classical three term recurrence, where the polynomial coefficients are obtained by solving a small system of linear equations. In finite precision arithmetic, the numerical orthogonality of the polynomials depends on a good choice of the size of the small systems; this size is usually controlled by a heuristic argument such as the condition number of the small matrix of coefficients. However, quite often it happens that orthogonality gets lost.
Optimizing polynomials for floating-point implementation
De Dinechin, Florent
2008-01-01
The floating-point implementation of a function on an interval often reduces to polynomial approximation, the polynomial being typically provided by Remez algorithm. However, the floating-point evaluation of a Remez polynomial sometimes leads to catastrophic cancellations. This happens when some of the polynomial coefficients are very small in magnitude with respects to others. In this case, it is better to force these coefficients to zero, which also reduces the operation count. This technique, classically used for odd or even functions, may be generalized to a much larger class of functions. An algorithm is presented that forces to zero the smaller coefficients of the initial polynomial thanks to a modified Remez algorithm targeting an incomplete monomial basis. One advantage of this technique is that it is purely numerical, the function being used as a numerical black box. This algorithm is implemented within a larger polynomial implementation tool that is demonstrated on a range of examples, resulting in ...
A new Arnoldi approach for polynomial eigenproblems
Energy Technology Data Exchange (ETDEWEB)
Raeven, F.A.
1996-12-31
In this paper we introduce a new generalization of the method of Arnoldi for matrix polynomials. The new approach is compared with the approach of rewriting the polynomial problem into a linear eigenproblem and applying the standard method of Arnoldi to the linearised problem. The algorithm that can be applied directly to the polynomial eigenproblem turns out to be more efficient, both in storage and in computation.
Exceptional polynomials and SUSY quantum mechanics
Indian Academy of Sciences (India)
K V S Shiv Chaitanya; S Sree Ranjani; Prasanta K Panigrahi; R Radhakrishnan; V Srinivasan
2015-07-01
We show that for the quantum mechanical problem which admit classical Laguerre/Jacobi polynomials as solutions for the Schrödinger equations (SE), will also admit exceptional Laguerre/Jacobi polynomials as solutions having the same eigenvalues but with the ground state missing after a modification of the potential. Then, we claim that the existence of these exceptional polynomials leads to the presence of non-trivial supersymmetry.
On the verification of polynomial system solvers
Institute of Scientific and Technical Information of China (English)
Changbo CHEN; Marc MORENO MAZA; Wei PAN; Yuzhen XI
2008-01-01
We discuss the verification of mathematical software solving polynomial systems symbolically by way of triangular decomposition. Standard verification techniques are highly resource consuming and apply only to polynomial systems which are easy to solve. We exhibit a new approach which manipulates constructible sets represented by regular systems. We provide comparative benchmarks of different verification procedures applied to four solvers on a large set of well-known polynomial systems. Our experimental results illustrate the high effi-ciency of our new approach. In particular, we are able to verify triangular decompositions of polynomial systems which are not easy to solve.
Control to Facet for Polynomial Systems
DEFF Research Database (Denmark)
Sloth, Christoffer; Wisniewski, Rafael
2014-01-01
for the controller design are solved by searching for polynomials in Bernstein form. This allows the controller design problem to be formulated as a linear programming problem. Examples are provided that demonstrate the efficiency of the method for designing controls for polynomial systems.......This paper presents a solution to the control to facet problem for arbitrary polynomial vector fields defined on simplices. The novelty of the work is to use Bernstein coefficients of polynomials for determining certificates of positivity. Specifically, the constraints that are set up...
Multi-indexed (q)-Racah Polynomials
Odake, Satoru
2012-01-01
As the second stage of the project $multi-indexed orthogonal polynomials$, we present, in the framework of `discrete quantum mechanics' with real shifts in one dimension, the multi-indexed (q)-Racah polynomials. They are obtained from the (q)-Racah polynomials by multiple application of the discrete analogue of the Darboux transformations or the Crum-Krein-Adler deletion of `virtual state' vectors of type I and II, in a similar way to the multi-indexed Laguerre and Jacobi polynomials reported earlier. The virtual state vectors are the `solutions' of the matrix Schr\\"odinger equation with negative `eigenvalues', except for one of the two boundary points.
Polynomial chaotic inflation in supergravity revisited
Directory of Open Access Journals (Sweden)
Kazunori Nakayama
2014-10-01
Full Text Available We revisit a polynomial chaotic inflation model in supergravity which we proposed soon after the Planck first data release. Recently some issues have been raised in Ref. [12], concerning the validity of our polynomial chaotic inflation model. We study the inflaton dynamics in detail, and confirm that the inflaton potential is very well approximated by a polynomial potential for the parameters of our interest in any practical sense, and in particular, the spectral index and the tensor-to-scalar ratio can be estimated by single-field approximation. This justifies our analysis of the polynomial chaotic inflation in supergravity.
Asymptotics for a generalization of Hermite polynomials
Alfaro, M; Peña, A; Rezola, M L
2009-01-01
We consider a generalization of the classical Hermite polynomials by the addition of terms involving derivatives in the inner product. This type of generalization has been studied in the literature from the point of view of the algebraic properties. Thus, our aim is to study the asymptotics of this sequence of nonstandard orthogonal polynomials. In fact, we obtain Mehler--Heine type formulas for these polynomials and, as a consequence, we prove that there exists an acceleration of the convergence of the smallest positive zeros of these generalized Hermite polynomials towards the origin.
Relative risk regression models with inverse polynomials.
Ning, Yang; Woodward, Mark
2013-08-30
The proportional hazards model assumes that the log hazard ratio is a linear function of parameters. In the current paper, we model the log relative risk as an inverse polynomial, which is particularly suitable for modeling bounded and asymmetric functions. The parameters estimated by maximizing the partial likelihood are consistent and asymptotically normal. The advantages of the inverse polynomial model over the ordinary polynomial model and the fractional polynomial model for fitting various asymmetric log relative risk functions are shown by simulation. The utility of the method is further supported by analyzing two real data sets, addressing the specific question of the location of the minimum risk threshold.
Directory of Open Access Journals (Sweden)
Ryoo CS
2010-01-01
Full Text Available The purpose of this paper is to give some properties of several Bernstein type polynomials to represent the fermionic -adic integral on . From these properties, we derive some interesting identities on the Euler numbers and polynomials.
Algorithms in Solving Polynomial Inequalities
Directory of Open Access Journals (Sweden)
Christopher M. Cordero
2015-11-01
Full Text Available A new method to solve the solution set of polynomial inequalities was conducted. When −1 −2 >0 ℎ 1,2∈ ℝ 10 if n is even. Then, the solution set is ∈ ℝ ∈ −∞,1 ∪ ,+∞ ∪ ,+1 : }. Thus, when −1−2…−≥0, the solution is ∈ ℝ ∈−∞, 1∪ ,+∞∪, +1: }. If is odd, then the solution set is ∈ ℝ ∈ ,+∞ ∪ ,+1 : }. Thus, when −1 −2…−≥0, the solution set is ∈ ℝ ∈ ,+∞∪, +1: }. Let −1−2…−<0 if n is even. Then, the solution set is ∈ ℝ ∈ ,+1 ∶ }. Thus, when −1 −2…−≤0, then the solution set is ∈ ℝ ∈, +1: }. If is an odd, then the solution set is ∈ ℝ ∈ −∞,1 ∪ ,+1 : }. Thus, when −1 −2 … − ≤0, the solution set is ∈ ℝ ∈ −∞,1 ∪ ,+1 : }. This research provides a novel method in solving the solution set of polynomial inequalities, in addition to other existing methods.
Classification based polynomial image interpolation
Lenke, Sebastian; Schröder, Hartmut
2008-02-01
Due to the fast migration of high resolution displays for home and office environments there is a strong demand for high quality picture scaling. This is caused on the one hand by large picture sizes and on the other hand due to an enhanced visibility of picture artifacts on these displays [1]. There are many proposals for an enhanced spatial interpolation adaptively matched to picture contents like e.g. edges. The drawback of these approaches is the normally integer and often limited interpolation factor. In order to achieve rational factors there exist combinations of adaptive and non adaptive linear filters, but due to the non adaptive step the overall quality is notably limited. We present in this paper a content adaptive polyphase interpolation method which uses "offline" trained filter coefficients and an "online" linear filtering depending on a simple classification of the input situation. Furthermore we present a new approach to a content adaptive interpolation polynomial, which allows arbitrary polyphase interpolation factors at runtime and further improves the overall interpolation quality. The main goal of our new approach is to optimize interpolation quality by adapting higher order polynomials directly to the image content. In addition we derive filter constraints for enhanced picture quality. Furthermore we extend the classification based filtering to the temporal dimension in order to use it for an intermediate image interpolation.
Convergent series for lattice models with polynomial interactions
Ivanov, Aleksandr S.; Sazonov, Vasily K.
2017-01-01
The standard perturbative weak-coupling expansions in lattice models are asymptotic. The reason for this is hidden in the incorrect interchange of the summation and integration. However, substituting the Gaussian initial approximation of the perturbative expansions by a certain interacting model or regularizing original lattice integrals, one can construct desired convergent series. In this paper we develop methods, which are based on the joint and separate utilization of the regularization and new initial approximation. We prove, that the convergent series exist and can be expressed as re-summed standard perturbation theory for any model on the finite lattice with the polynomial interaction of even degree. We discuss properties of such series and study their applicability to practical computations on the example of the lattice ϕ4-model. We calculate expectation value using the convergent series, the comparison of the results with the Borel re-summation and Monte Carlo simulations shows a good agreement between all these methods.
Convergent series for lattice models with polynomial interactions
Ivanov, Aleksandr S
2016-01-01
The standard perturbative weak-coupling expansions in lattice models are asymptotic. The reason for this is hidden in the incorrect interchange of the summation and integration. However, substituting the Gaussian initial approximation of the perturbative expansions by a certain interacting model or regularizing original lattice integrals, one can construct desired convergent series. In this paper we develop methods, which are based on the joint and separate utilization of the regularization and new initial approximation. We prove, that the convergent series exist and can be expressed as the re-summed standard perturbation theory for any model on the finite lattice with the polynomial interaction of even degree. We discuss properties of such series and make them applicable to practical computations. The workability of the methods is demonstrated on the example of the lattice $\\phi^4$-model. We calculate the operator $\\langle\\phi_n^2\\rangle$ using the convergent series, the comparison of the results with the Bo...
ON FIRST INTEGRALS OF POLYNOMIAL AUTONOMOUS SYSTEMS
Institute of Scientific and Technical Information of China (English)
WANG Yuzhen; CHENG Daizhan; LI Chunwen
2002-01-01
Using Carleman linearization procedure, this paper investigates the problemof first integrals of polynomial autonomous systems and proposes a procedure to find thefirst integrals of polynomial family for the systems. A generalized eigenequation is obtainedand then the problem is reduced to the solvability of the eigenequation. The result is ageneralization of some known results.
Reliability polynomials crossing more than twice
Brown, J.I.; Koç, Y.; Kooij, R.E.
2011-01-01
In this paper we study all-terminal reliability polynomials of networks having the same number of nodes and the same number of links. First we show that the smallest possible size for a pair of networks that allows for two crossings of their reliability polynomials have seven nodes and fifteen edges
Sums of Powers of Fibonacci Polynomials
Indian Academy of Sciences (India)
Helmut Prodinger
2009-11-01
Using the explicit (Binet) formula for the Fibonacci polynomials, a summation formula for powers of Fibonacci polynomials is derived straightforwardly, which generalizes a recent result for squares that appeared in Proc. Ind. Acad. Sci. (Math. Sci.) 118 (2008) 27--41.
Notes on Schubert, Grothendieck and Key Polynomials
Kirillov, Anatol N.
2016-03-01
We introduce common generalization of (double) Schubert, Grothendieck, Demazure, dual and stable Grothendieck polynomials, and Di Francesco-Zinn-Justin polynomials. Our approach is based on the study of algebraic and combinatorial properties of the reduced rectangular plactic algebra and associated Cauchy kernels.
BOUNDS FOR THE ZEROS OF POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
W. M. Shah; A.Liman
2004-01-01
Let P(z) =n∑j=0 ajzj be a polynomial of degree n. In this paper we prove a more general result which interalia improves upon the bounds of a class of polynomials. We also prove a result which includes some extensions and generalizations of Enestrom-Kakeya theorem.
Several explicit formulae for Bernoulli polynomials
Komatsu, Takao; Pita Ruiz V., Claudio de J.
2016-01-01
We prove several explicit formulae for the $n$-th Bernoulli polynomial $B_{n}(x)$, in which $B_{n}(x)$ is equal to an affine combination of the polynomials $(x-1)^{n}$, $(x-2)^{n}$, $ldots$, $(x-k-1)^{n}$, where $k$ is any fixed positive integer greater or equal than $n$.
Elementary combinatorics of the HOMFLYPT polynomial
Chmutov, Sergei
2009-01-01
We explore Jaeger's state model for the HOMFLYPT polynomial. We reformulate this model in the language of Gauss diagrams and use it to obtain Gauss diagram formulas for a two-parameter family of Vassiliev invariants coming from the HOMFLYPT polynomial. These formulas are new already for invariants of degree 3.
A Note on Solvable Polynomial Algebras
Directory of Open Access Journals (Sweden)
Huishi Li
2014-03-01
Full Text Available In terms of their defining relations, solvable polynomial algebras introduced by Kandri-Rody and Weispfenning [J. Symbolic Comput., 9(1990] are characterized by employing Gr\\"obner bases of ideals in free algebras, thereby solvable polynomial algebras are completely determinable and constructible in a computational way.
Connections between the matching and chromatic polynomials
Directory of Open Access Journals (Sweden)
E. J. Farrell
1992-01-01
Full Text Available The main results established are (i a connection between the matching and chromatic polynomials and (ii a formula for the matching polynomial of a general complement of a subgraph of a graph. Some deductions on matching and chromatic equivalence and uniqueness are made.
On Polynomial Functions over Finite Commutative Rings
Institute of Scientific and Technical Information of China (English)
Jian Jun JIANG; Guo Hua PENG; Qi SUN; Qi Fan ZHANG
2006-01-01
Let R be an arbitrary finite commutative local ring. In this paper, we obtain a necessary and sufficient condition for a function over R to be a polynomial function. Before this paper, necessary and sufficient conditions for a function to be a polynomial function over some special finite commutative local rings were obtained.
A polynomial approach to nonlinear system controllability
Zheng, YF; Willems, JC; Zhang, CH
2001-01-01
This note uses a polynomial approach to present a necessary and sufficient condition for local controllability of single-input-single-output (SISO) nonlinear systems. The condition is presented in terms of common factors of a noncommutative polynomial expression. This result exposes controllability
The topology of Julia sets for polynomials
Institute of Scientific and Technical Information of China (English)
尹永成
2002-01-01
We prove that wandering components of the Julia set of a polynomial are singletons provided each critical point in a wandering Julia component is non-recurrent. This means a conjecture of Branner-Hubbard is true for this kind of polynomials.
Fractal Trigonometric Polynomials for Restricted Range Approximation
Chand, A. K. B.; Navascués, M. A.; Viswanathan, P.; Katiyar, S. K.
2016-05-01
One-sided approximation tackles the problem of approximation of a prescribed function by simple traditional functions such as polynomials or trigonometric functions that lie completely above or below it. In this paper, we use the concept of fractal interpolation function (FIF), precisely of fractal trigonometric polynomials, to construct one-sided uniform approximants for some classes of continuous functions.
New pole placement algorithm - Polynomial matrix approach
Shafai, B.; Keel, L. H.
1990-01-01
A simple and direct pole-placement algorithm is introduced for dynamical systems having a block companion matrix A. The algorithm utilizes well-established properties of matrix polynomials. Pole placement is achieved by appropriately assigning coefficient matrices of the corresponding matrix polynomial. This involves only matrix additions and multiplications without requiring matrix inversion. A numerical example is given for the purpose of illustration.
Distortion control of conjugacies between quadratic polynomials
Institute of Scientific and Technical Information of China (English)
无
2010-01-01
We use a new type of distortion control of univalent functions to give an alternative proof of Douady-Hubbard’s ray-landing theorem for quadratic Misiurewicz polynomials. The univalent maps arise from Thurston’s iterated algorithm on perturbation of such polynomials.
Uniqueness of meromorphic functions concerning differential polynomials
Institute of Scientific and Technical Information of China (English)
QIAO Lei
2007-01-01
Based on a unicity theorem for entire funcitions concerning differential polynomials proposed by M. L. Fang and W. Hong, we studied the uniqueness problem of two meromorphic functions whose differential polynomials share the same 1-point by proving two theorems and their related lemmas. The results extend and improve given by Fang and Hong's theorem.
Fostering Connections between Classes of Polynomial Functions.
Buck, Judy Curran
The typical path of instruction in high school algebra courses for the study of polynomial functions has been from linear functions, to quadratic functions, to polynomial functions of degree greater than two. This paper reports results of clinical interviews with an Algebra II student. The interviews were used to probe into the student's…
Colored HOMFLY polynomials can distinguish mutant knots
Nawata, Satoshi; Singh, Vivek Kumar
2015-01-01
We illustrate from the viewpoint of braiding operations on WZNW conformal blocks how colored HOMFLY polynomials with multiplicity structure can detect mutations. As an example, we explicitly evaluate the (2,1)-colored HOMFLY polynomials that distinguish a famous mutant pair, Kinoshita-Terasaka and Conway knot.
Indian Academy of Sciences (India)
V K Jain
2009-02-01
For a polynomial of degree , we have obtained an upper bound involving coefficients of the polynomial, for moduli of its zeros of smallest moduli, and then a refinement of the well-known Eneström–Kakeya theorem (under certain conditions).
Differential Krull dimension in differential polynomial extensions
Smirnov, Ilya
2011-01-01
We investigate the differential Krull dimension of differential polynomials over a differential ring. We prove a differential analogue of Jaffard's Special Chain Theorem and show that differential polynomial extensions of certain classes of differential rings have no anomaly of differential Krull dimension.
Polynomial weights and code constructions
DEFF Research Database (Denmark)
Massey, J; Costello, D; Justesen, Jørn
1973-01-01
polynomial included. This fundamental property is then used as the key to a variety of code constructions including 1) a simplified derivation of the binary Reed-Muller codes and, for any primepgreater than 2, a new extensive class ofp-ary "Reed-Muller codes," 2) a new class of "repeated-root" cyclic codes...... that are subcodes of the binary Reed-Muller codes and can be very simply instrumented, 3) a new class of constacyclic codes that are subcodes of thep-ary "Reed-Muller codes," 4) two new classes of binary convolutional codes with large "free distance" derived from known binary cyclic codes, 5) two new classes...... of long constraint length binary convolutional codes derived from2^r-ary Reed-Solomon codes, and 6) a new class ofq-ary "repeated-root" constacyclic codes with an algebraic decoding algorithm....
On the defect and stability of differential expansion
Kononov, Ya
2015-01-01
Empirical analysis of many colored knot polynomials, made possible by recent computational advances in Chern-Simons theory, reveals their stability: for any given negative N and any given knot the set of coefficients of the polynomial in r-th symmetric representation does not change with r, if it is large enough. This fact reflects the non-trivial and previously unknown properties of the differential expansion, and it turns out that from this point of view there are universality classes of knots, characterized by a single integer, which we call defect, and which is in fact related to the power of Alexander polynomial.
Morozov, Anton
2012-01-01
Like all other knot polynomials, the superpolynomials should be defined in arbitrary representation R of the gauge group in (refined) Chern-Simons theory. However, not a single example is yet known of a superpolynomial beyond symmetric or antisymmetric representations. We consider the expansion of the superpolynomial around the special polynomial in powers of (q-1) and (t-1) and suggest a simple formula for the first-order deviation, which is presumably valid for arbitrary representation. This formula can serve as a crucial lacking test of various formulas for non-trivial superpolynomials, which will appear in the literature in the near future.
A spin network generalization of the Jones Polynomial and Vassiliev invariants
Gambini, R; Pullin, J
1998-01-01
We apply the ideas of Alvarez and Labastida to the invariant of spin networks defined by Witten and Martin based on Chern-Simons theory. We show that it is possible to define ambient invariants of spin networks that (for the case of SU(2)) can be considered as extensions to spin networks of the Jones polynomial. Expansions of the coefficients of the polynomial yield primitive Vassiliev invariants. The resulting invariants are candidates for solutions of the Wheeler--DeWitt equations in the spin network representation of quantum gravity.
The Role of Orthogonal Polynomials in Tailoring Spherical Distributions to Kurtosis Requirements
Directory of Open Access Journals (Sweden)
Luca Bagnato
2016-08-01
Full Text Available This paper carries out an investigation of the orthogonal-polynomial approach to reshaping symmetric distributions to fit in with data requirements so as to cover the multivariate case. With this objective in mind, reference is made to the class of spherical distributions, given that they provide a natural multivariate generalization of univariate even densities. After showing how to tailor a spherical distribution via orthogonal polynomials to better comply with kurtosis requirements, we provide operational conditions for the positiveness of the resulting multivariate Gram–Charlier-like expansion, together with its kurtosis range. Finally, the approach proposed here is applied to some selected spherical distributions.
Sobolev orthogonal polynomials on a simplex
Aktas, Rabia
2011-01-01
The Jacobi polynomials on the simplex are orthogonal polynomials with respect to the weight function $W_\\bg(x) = x_1^{\\g_1} ... x_d^{\\g_d} (1- |x|)^{\\g_{d+1}}$ when all $\\g_i > -1$ and they are eigenfunctions of a second order partial differential operator $L_\\bg$. The singular cases that some, or all, $\\g_1,...,\\g_{d+1}$ are -1 are studied in this paper. Firstly a complete basis of polynomials that are eigenfunctions of $L_\\bg$ in each singular case is found. Secondly, these polynomials are shown to be orthogonal with respect to an inner product which is explicitly determined. This inner product involves derivatives of the functions, hence the name Sobolev orthogonal polynomials.
On Chebyshev polynomials and torus knots
Gavrilik, A M
2009-01-01
In this work we demonstrate that the q-numbers and their two-parameter generalization, the q,p-numbers, can be used to obtain some polynomial invariants for torus knots and links. First, we show that the q-numbers, which are closely connected with the Chebyshev polynomials, can also be related with the Alexander polynomials for the class T(s,2) of torus knots, s being an odd integer, and used for finding the corresponding skein relation. Then, we develop this procedure in order to obtain, with the help of q,p-numbers, the generalized two-variable Alexander polynomials, and prove their direct connection with the HOMFLY polynomials and the skein relation of the latter.
Baxter operator formalism for Macdonald polynomials
Gerasimov, Anton; Oblezin, Sergey
2012-01-01
We develop basic constructions of the Baxter operator formalism for the Macdonald polynomials. Precisely we construct a dual pair of mutually commuting Baxter operators such that the Macdonald polynomials are their common eigenfunctions. The dual pair of Baxter operators is closely related to the dual pair of recursive operators for Macdonald polynomials leading to various families of their integral representations. We also construct the Baxter operator formalism for the q-deformed Whittaker functions and the Jack polynomials obtained by degenerations of the Macdonald polynomials. This note provides a generalization of our previous results on the Baxter operator formalism for the Whittaker functions. It was demonstrated previously that Baxter operator formalism for the Whittaker functions has deep connections with representation theory. In particular the Baxter operators should be considered as elements of appropriate spherical Hecke algebras and their eigenvalues are identified with local Archimedean L-facto...
Orthogonal Polynomials from Hermitian Matrices II
Odake, Satoru
2016-01-01
This is the second part of the project `unified theory of classical orthogonal polynomials of a discrete variable derived from the eigenvalue problems of hermitian matrices.' In a previous paper, orthogonal polynomials having Jackson integral measures were not included, since such measures cannot be obtained from single infinite dimensional hermitian matrices. Here we show that Jackson integral measures for the polynomials of the big $q$-Jacobi family are the consequence of the recovery of self-adjointness of the unbounded Jacobi matrices governing the difference equations of these polynomials. The recovery of self-adjointness is achieved in an extended $\\ell^2$ Hilbert space on which a direct sum of two unbounded Jacobi matrices acts as a Hamiltonian or a difference Schr\\"odinger operator for an infinite dimensional eigenvalue problem. The polynomial appearing in the upper/lower end of Jackson integral constitutes the eigenvector of each of the two unbounded Jacobi matrix of the direct sum. We also point out...
Polynomials with Palindromic and Unimodal Coeﬃ cients
Institute of Scientific and Technical Information of China (English)
Hua SUN; Yi WANG; Hai Xia ZHANG
2015-01-01
Let f(q) = arqr +· · ·+asqs, with ar = 0 and as = 0, be a real polynomial. It is a palindromic polynomial of darga n if r+s = n and ar+i = as−i for all i. Polynomials of darga n form a linear subspace Pn(q) of R(q)n+1 of dimension ? n2 ?+1. We give transition matrices between two bases ?qj(1+q+· · ·+qn−2j)? , ?qj(1+q)n−2j? and the standard basis ?qj(1+qn−2j)? of Pn(q). We present some characterizations and sufficcient conditions for palindromic polynomials that can be expressed in terms of these two bases with nonnegative coefficients. We also point out the link between such polynomials and rank-generating functions of posets.
Jacob's ladders and new orthogonal systems generated by Jacobi polynomials
Moser, Jan
2010-01-01
Is is shown in this paper that there is a connection between the Riemann zeta-function $\\zf$ and the classical Jacobi's polynomials, i.e. the Legendre polynomials, Chebyshev polynomials of the first and the second kind,...
HIGHER ORDER MULTIVARIABLE NORLUND EULER-BERNOULLI POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
刘国栋
2002-01-01
The definitions of higher order multivariable Norlund Euler polynomials and Norlund Bernoulli polynomials are presented and some of their important properties are expounded. Some identities involving recurrence sequences and higher order multivariable Norlund Euler-Bernoulli polynomials are established.
Polynomial Interpolation in the Elliptic Curve Cryptosystem
Directory of Open Access Journals (Sweden)
Liew K. Jie
2011-01-01
Full Text Available Problem statement: In this research, we incorporate the polynomial interpolation method in the discrete logarithm problem based cryptosystem which is the elliptic curve cryptosystem. Approach: In this study, the polynomial interpolation method to be focused is the Lagrange polynomial interpolation which is the simplest polynomial interpolation method. This method will be incorporated in the encryption algorithm of the elliptic curve ElGamal cryptosystem. Results: The scheme modifies the elliptic curve ElGamal cryptosystem by adding few steps in the encryption algorithm. Two polynomials are constructed based on the encrypted points using Lagrange polynomial interpolation and encrypted for the second time using the proposed encryption method. We believe it is safe from the theoretical side as it still relies on the discrete logarithm problem of the elliptic curve. Conclusion/Recommendations: The modified scheme is expected to be more secure than the existing scheme as it offers double encryption techniques. On top of the existing encryption algorithm, we managed to encrypt one more time using the polynomial interpolation method. We also have provided detail examples based on the described algorithm.
Basis of symmetric polynomials for many-boson light-front wave functions.
Chabysheva, Sophia S; Hiller, John R
2014-12-01
We provide an algorithm for the construction of orthonormal multivariate polynomials that are symmetric with respect to the interchange of any two coordinates on the unit hypercube and are constrained to the hyperplane where the sum of the coordinates is one. These polynomials form a basis for the expansion of bosonic light-front momentum-space wave functions, as functions of longitudinal momentum, where momentum conservation guarantees that the fractions are on the interval [0,1] and sum to one. This generalizes earlier work on three-boson wave functions to wave functions for arbitrarily many identical bosons. A simple application in two-dimensional ϕ(4) theory illustrates the use of these polynomials.
Edgeworth expansions and rates of convergence for normalized sums: Chung's 1946 method revisited
2010-01-01
Abstract In this paper we revisit, correct and extend Chung?s 1946 method for deriving higher order Edgeworth expansions with respect to t-statistics and generalized self-normalized sums. Thereby we provide a set of formulas which allows the computation of the approximation of any order and specify the first four polynomials in the Edgeworth expansion the first two of which are well known. It turns out that knowledge of the first four polynomials is necessary and sufficient for cha...
On Chebyshev polynomials and torus knots
Gavrilik, A. M.; Pavlyuk, A. M.
2009-01-01
In this work we demonstrate that the q-numbers and their two-parameter generalization, the q,p-numbers, can be used to obtain some polynomial invariants for torus knots and links. First, we show that the q-numbers, which are closely connected with the Chebyshev polynomials, can also be related with the Alexander polynomials for the class T(s,2) of torus knots, s being an odd integer, and used for finding the corresponding skein relation. Then, we develop this procedure in order to obtain, wit...
Cycles are determined by their domination polynomials
Akbari, Saieed
2009-01-01
Let $G$ be a simple graph of order $n$. A dominating set of $G$ is a set $S$ of vertices of $G$ so that every vertex of $G$ is either in $S$ or adjacent to a vertex in $S$. The domination polynomial of $G$ is the polynomial $D(G,x)=\\sum_{i=1}^{n} d(G,i) x^{i}$, where $d(G,i)$ is the number of dominating sets of $G$ of size $i$. In this paper we show that cycles are determined by their domination polynomials.
Limits of zeros of polynomial sequences
Zhu, Xinyun; Grossman, George
2007-01-01
In the present paper we consider $F_k(x)=x^{k}-\\sum_{t=0}^{k-1}x^t,$ the characteristic polynomial of the $k$-th order Fibonacci sequence, the latter denoted $G(k,l).$ We determine the limits of the real roots of certain odd and even degree polynomials related to the derivatives and integrals of $F_k(x),$ that form infinite sequences of polynomials, of increasing degree. In particular, as $k \\to \\infty,$ the limiting values of the zeros are determined, for both odd and even cases. It is also ...
A bivariate chromatic polynomial for signed graphs
Beck, Matthias
2012-01-01
We study Dohmen--P\\"onitz--Tittmann's bivariate chromatic polynomial $c_\\Gamma(k,l)$ which counts all $(k+l)$-colorings of a graph $\\Gamma$ such that adjacent vertices get different colors if they are $\\le k$. Our first contribution is an extension of $c_\\Gamma(k,l)$ to signed graphs, for which we obtain an inclusion--exclusion formula and several special evaluations giving rise, e.g., to polynomials that encode balanced subgraphs. Our second goal is to derive combinatorial reciprocity theorems for $c_\\Gamma(k,l)$ and its signed-graph analogues, reminiscent of Stanley's reciprocity theorem linking chromatic polynomials to acyclic orientations.
A Polynomial Preconditioner for the CMRH Algorithm
Directory of Open Access Journals (Sweden)
Jiangzhou Lai
2011-01-01
Full Text Available Many large and sparse linear systems can be solved efficiently by restarted GMRES and CMRH methods Sadok 1999. The CMRH(m method is less expensive and requires slightly less storage than GMRES(m. But like GMRES, the restarted CMRH method may not converge. In order to remedy this defect, this paper presents a polynomial preconditioner for CMRH-based algorithm. Numerical experiments are given to show that the polynomial preconditioner is quite simple and easily constructed and the preconditioned CMRH(m with the polynomial preconditioner has better performance than CMRH(m.
On Calculation of Adomian Polynomials by MATLAB
Directory of Open Access Journals (Sweden)
Hossein ABOLGHASEMI
2011-01-01
Full Text Available Adomian Decomposition Method (ADM is an elegant technique to handle an extensive class of linear or nonlinear differential and integral equations. However, in case of nonlinear equations, ADM demands a special representation of each nonlinear term, namely, Adomian polynomials. The present paper introduces a novel MATLAB code which computes Adomian polynomials associated with several types of nonlinearities. The code exploits symbolic programming incorporated with a recently proposed alternative scheme to be straightforward and fast. For the sake of exemplification, Adomian polynomials of famous nonlinear operators, computed by the code, are given.
The Translated Dowling Polynomials and Numbers.
Mangontarum, Mahid M; Macodi-Ringia, Amila P; Abdulcarim, Normalah S
2014-01-01
More properties for the translated Whitney numbers of the second kind such as horizontal generating function, explicit formula, and exponential generating function are proposed. Using the translated Whitney numbers of the second kind, we will define the translated Dowling polynomials and numbers. Basic properties such as exponential generating functions and explicit formula for the translated Dowling polynomials and numbers are obtained. Convexity, integral representation, and other interesting identities are also investigated and presented. We show that the properties obtained are generalizations of some of the known results involving the classical Bell polynomials and numbers. Lastly, we established the Hankel transform of the translated Dowling numbers.
Polynomial threshold functions and Boolean threshold circuits
DEFF Research Database (Denmark)
Hansen, Kristoffer Arnsfelt; Podolskii, Vladimir V.
2013-01-01
We study the complexity of computing Boolean functions on general Boolean domains by polynomial threshold functions (PTFs). A typical example of a general Boolean domain is 12n . We are mainly interested in the length (the number of monomials) of PTFs, with their degree and weight being...... of secondary interest. We show that PTFs on general Boolean domains are tightly connected to depth two threshold circuits. Our main results in regard to this connection are: PTFs of polynomial length and polynomial degree compute exactly the functions computed by THRMAJ circuits. An exponential length lower...
Laurent polynomial moment problem: a case study
Pakovich, F; Zvonkin, A
2009-01-01
In recent years, the so-called polynomial moment problem, motivated by the classical Poincare center-focus problem, was thoroughly studied, and the answers to the main questions have been found. The study of a similar problem for rational functions is still at its very beginning. In this paper, we make certain progress in this direction; namely, we construct an example of a Laurent polynomial for which the solutions of the corresponding moment problem behave in a significantly more complicated way than it would be possible for a polynomial.
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.
More on rotations as spin matrix polynomials
Energy Technology Data Exchange (ETDEWEB)
Curtright, Thomas L. [Department of Physics, University of Miami, Coral Gables, Florida 33124-8046 (United States)
2015-09-15
Any nonsingular function of spin j matrices always reduces to a matrix polynomial of order 2j. The challenge is to find a convenient form for the coefficients of the matrix polynomial. The theory of biorthogonal systems is a useful framework to meet this challenge. Central factorial numbers play a key role in the theoretical development. Explicit polynomial coefficients for rotations expressed either as exponentials or as rational Cayley transforms are considered here. Structural features of the results are discussed and compared, and large j limits of the coefficients are examined.
Modular polynomials via isogeny volcanoes
Broker, Reinier; Sutherland, Andrew V
2010-01-01
We present a new algorithm to compute the classical modular polynomial Phi_n in the rings Z[X,Y] and (Z/mZ)[X,Y], for a prime n and any positive integer m. Our approach uses the graph of n-isogenies to efficiently compute Phi_n mod p for many primes p of a suitable form, and then applies the Chinese Remainder Theorem (CRT). Under the Generalized Riemann Hypothesis (GRH), we achieve an expected running time of O(n^3 (log n)^3 log log n), and compute Phi_n mod m using O(n^2 (log n)^2 + n^2 log m) space. We have used the new algorithm to compute Phi_n with n over 5000, and Phi_n mod m with n over 20000. We also consider several modular functions g for which Phi_n^g is smaller than Phi_n, allowing us to handle n over 60000.
Local Polynomial Estimation of Distribution Functions
Institute of Scientific and Technical Information of China (English)
LI Yong-hong; ZENG Xia
2007-01-01
Under the condition that the total distribution function is continuous and bounded on (-∞,∞), we constructed estimations for distribution and hazard functions with local polynomial method, and obtained the rate of strong convergence of the estimations.
Tutte Polynomial of Scale-Free Networks
Chen, Hanlin; Deng, Hanyuan
2016-05-01
The Tutte polynomial of a graph, or equivalently the q-state Potts model partition function, is a two-variable polynomial graph invariant of considerable importance in both statistical physics and combinatorics. The computation of this invariant for a graph is NP-hard in general. In this paper, we focus on two iteratively growing scale-free networks, which are ubiquitous in real-life systems. Based on their self-similar structures, we mainly obtain recursive formulas for the Tutte polynomials of two scale-free networks (lattices), one is fractal and "large world", while the other is non-fractal but possess the small-world property. Furthermore, we give some exact analytical expressions of the Tutte polynomial for several special points at ( x, y)-plane, such as, the number of spanning trees, the number of acyclic orientations, etc.
Generation of multivariate Hermite interpolating polynomials
Tavares, Santiago Alves
2005-01-01
Generation of Multivariate Hermite Interpolating Polynomials advances the study of approximate solutions to partial differential equations by presenting a novel approach that employs Hermite interpolating polynomials and bysupplying algorithms useful in applying this approach.Organized into three sections, the book begins with a thorough examination of constrained numbers, which form the basis for constructing interpolating polynomials. The author develops their geometric representation in coordinate systems in several dimensions and presents generating algorithms for each level number. He then discusses their applications in computing the derivative of the product of functions of several variables and in the construction of expression for n-dimensional natural numbers. Section II focuses on the construction of Hermite interpolating polynomials, from their characterizing properties and generating algorithms to a graphical analysis of their behavior. The final section of the book is dedicated to the applicatio...
Superconformal minimal models and admissible Jack polynomials
Blondeau-Fournier, Olivier; Ridout, David; Wood, Simon
2016-01-01
We give new proofs of the rationality of the N=1 superconformal minimal model vertex operator superalgebras and of the classification of their modules in both the Neveu-Schwarz and Ramond sectors. For this, we combine the standard free field realisation with the theory of Jack symmetric functions. A key role is played by Jack symmetric polynomials with a certain negative parameter that are labelled by admissible partitions. These polynomials are shown to describe free fermion correlators, suitably dressed by a symmetrising factor. The classification proofs concentrate on explicitly identifying Zhu's algebra and its twisted analogue. Interestingly, these identifications do not use an explicit expression for the non-trivial vacuum singular vector. While the latter is known to be expressible in terms of an Uglov symmetric polynomial or a linear combination of Jack superpolynomials, it turns out that standard Jack polynomials (and functions) suffice to prove the classification.
Hermite polynomials and quasi-classical asymptotics
Energy Technology Data Exchange (ETDEWEB)
Ali, S. Twareque, E-mail: twareque.ali@concordia.ca [Department of Mathematics and Statistics, Concordia University, Montréal, Québec H3G 1M8 (Canada); Engliš, Miroslav, E-mail: englis@math.cas.cz [Mathematics Institute, Silesian University in Opava, Na Rybníčku 1, 74601 Opava, Czech Republic and Mathematics Institute, Žitná 25, 11567 Prague 1 (Czech Republic)
2014-04-15
We study an unorthodox variant of the Berezin-Toeplitz type of quantization scheme, on a reproducing kernel Hilbert space generated by the real Hermite polynomials and work out the associated quasi-classical asymptotics.
Concentration for noncommutative polynomials in random matrices
2011-01-01
We present a concentration inequality for linear functionals of noncommutative polynomials in random matrices. Our hypotheses cover most standard ensembles, including Gaussian matrices, matrices with independent uniformly bounded entries and unitary or orthogonal matrices.
Thermodynamic characterization of networks using graph polynomials
Ye, Cheng; Peron, Thomas K DM; Silva, Filipi N; Rodrigues, Francisco A; Costa, Luciano da F; Torsello, Andrea; Hancock, Edwin R
2015-01-01
In this paper, we present a method for characterizing the evolution of time-varying complex networks by adopting a thermodynamic representation of network structure computed from a polynomial (or algebraic) characterization of graph structure. Commencing from a representation of graph structure based on a characteristic polynomial computed from the normalized Laplacian matrix, we show how the polynomial is linked to the Boltzmann partition function of a network. This allows us to compute a number of thermodynamic quantities for the network, including the average energy and entropy. Assuming that the system does not change volume, we can also compute the temperature, defined as the rate of change of entropy with energy. All three thermodynamic variables can be approximated using low-order Taylor series that can be computed using the traces of powers of the Laplacian matrix, avoiding explicit computation of the normalized Laplacian spectrum. These polynomial approximations allow a smoothed representation of the...
Characteristic Polynomials of Complex Random Matrix Models
Akemann, G
2003-01-01
We calculate the expectation value of an arbitrary product of characteristic polynomials of complex random matrices and their hermitian conjugates. Using the technique of orthogonal polynomials in the complex plane our result can be written in terms of a determinant containing these polynomials and their kernel. It generalizes the known expression for hermitian matrices and it also provides a generalization of the Christoffel formula to the complex plane. The derivation we present holds for complex matrix models with a general weight function at finite-N, where N is the size of the matrix. We give some explicit examples at finite-N for specific weight functions. The characteristic polynomials in the large-N limit at weak and strong non-hermiticity follow easily and they are universal in the weak limit. We also comment on the issue of the BMN large-N limit.
Twisted Polynomials and Forgery Attacks on GCM
DEFF Research Database (Denmark)
Abdelraheem, Mohamed Ahmed A. M. A.; Beelen, Peter; Bogdanov, Andrey;
2015-01-01
nonce misuse resistance, such as POET. The algebraic structure of polynomial hashing has given rise to security concerns: At CRYPTO 2008, Handschuh and Preneel describe key recovery attacks, and at FSE 2013, Procter and Cid provide a comprehensive framework for forgery attacks. Both approaches rely...... heavily on the ability to construct forgery polynomials having disjoint sets of roots, with many roots (“weak keys”) each. Constructing such polynomials beyond naïve approaches is crucial for these attacks, but still an open problem. In this paper, we comprehensively address this issue. We propose to use...... in an improved key recovery algorithm. As cryptanalytic applications of our twisted polynomials, we develop the first universal forgery attacks on GCM in the weak-key model that do not require nonce reuse. Moreover, we present universal weak-key forgeries for the nonce-misuse resistant AE scheme POET, which...
Limits of zeros of polynomial sequences
Zhu, Xinyun
2007-01-01
In the present paper we consider $F_k(x)=x^{k}-\\sum_{t=0}^{k-1}x^t,$ the characteristic polynomial of the $k$-th order Fibonacci sequence, the latter denoted $G(k,l).$ We determine the limits of the real roots of certain odd and even degree polynomials related to the derivatives and integrals of $F_k(x),$ that form infinite sequences of polynomials, of increasing degree. In particular, as $k \\to \\infty,$ the limiting values of the zeros are determined, for both odd and even cases. It is also shown, in both cases, that the convergence is monotone for sufficiently large degree. We give an upper bound for the modulus of the complex zeros of the polynomials for each sequence. This gives a general solution related to problems considered by Dubeau 1989, 1993, Miles 1960, Flores 1967, Miller 1971 and later by the second author in the present paper, and Narayan 1997.
Rational Convolution Roots of Isobaric Polynomials
Conci, Aura; Li, Huilan; MacHenry, Trueman
2014-01-01
In this paper, we exhibit two matrix representations of the rational roots of generalized Fibonacci polynomials (GFPs) under convolution product, in terms of determinants and permanents, respectively. The underlying root formulas for GFPs and for weighted isobaric polynomials (WIPs), which appeared in an earlier paper by MacHenry and Tudose, make use of two types of operators. These operators are derived from the generating functions for Stirling numbers of the first kind and second kind. Hen...
Blind Signature Scheme Based on Chebyshev Polynomials
Maheswara Rao Valluri
2011-01-01
A blind signature scheme is a cryptographic protocol to obtain a valid signature for a message from a signer such that signer’s view of the protocol can’t be linked to the resulting message signature pair. This paper presents blind signature scheme using Chebyshev polynomials. The security of the given scheme depends upon the intractability of the integer factorization problem and discrete logarithms ofChebyshev polynomials.
Positive maps, positive polynomials and entanglement witnesses
Skowronek, Lukasz
2009-01-01
We link the study of positive quantum maps, block positive operators, and entanglement witnesses with problems related to multivariate polynomials. For instance, we show how indecomposable block positive operators relate to biquadratic forms that are not sums of squares. Although the general problem of describing the set of positive maps remains open, in some particular cases we solve the corresponding polynomial inequalities and obtain explicit conditions for positivity.
Positive maps, positive polynomials and entanglement witnesses
Energy Technology Data Exchange (ETDEWEB)
Skowronek, Lukasz; Zyczkowski, Karol [Institute of Physics, Jagiellonian University, Krakow (Poland)], E-mail: lukasz.skowronek@uj.edu.pl, E-mail: karol@tatry.if.uj.edu.pl
2009-08-14
We link the study of positive quantum maps, block positive operators and entanglement witnesses with problems related to multivariate polynomials. For instance, we show how indecomposable block positive operators relate to biquadratic forms that are not sums of squares. Although the general problem of describing the set of positive maps remains open, in some particular cases we solve the corresponding polynomial inequalities and obtain explicit conditions for positivity.
Blind Signature Scheme Based on Chebyshev Polynomials
Directory of Open Access Journals (Sweden)
Maheswara Rao Valluri
2011-12-01
Full Text Available A blind signature scheme is a cryptographic protocol to obtain a valid signature for a message from a signer such that signer’s view of the protocol can’t be linked to the resulting message signature pair. This paper presents blind signature scheme using Chebyshev polynomials. The security of the given scheme depends upon the intractability of the integer factorization problem and discrete logarithms ofChebyshev polynomials.
Ferrers Matrices Characterized by the Rook Polynomials
Institute of Scientific and Technical Information of China (English)
MAHai-cheng; HUSheng-biao
2003-01-01
In this paper,we show that there exist precisely W(A) Ferrers matrices F(C1,C2,…,cm)such that the rook polynomials is equal to the rook polynomial of Ferrers matrix F(b1,b2,…,bm), where A={b1,b2-1,…,bm-m+1} is a repeated set,W(A) is weight of A.
Recursive Polynomial Remainder Sequence and its Subresultants
Terui, Akira
2008-01-01
We introduce concepts of "recursive polynomial remainder sequence (PRS)" and "recursive subresultant," along with investigation of their properties. A recursive PRS is defined as, if there exists the GCD (greatest common divisor) of initial polynomials, a sequence of PRSs calculated "recursively" for the GCD and its derivative until a constant is derived, and recursive subresultants are defined by determinants representing the coefficients in recursive PRS as functions of coefficients of init...
Subresultants in Recursive Polynomial Remainder Sequence
Terui, Akira
2008-01-01
We introduce concepts of "recursive polynomial remainder sequence (PRS)" and "recursive subresultant," and investigate their properties. In calculating PRS, if there exists the GCD (greatest common divisor) of initial polynomials, we calculate "recursively" with new PRS for the GCD and its derivative, until a constant is derived. We call such a PRS a recursive PRS. We define recursive subresultants to be determinants representing the coefficients in recursive PRS by coefficients of initial po...
Local fibred right adjoints are polynomial
DEFF Research Database (Denmark)
Kock, Anders; Kock, Joachim
2013-01-01
For any locally cartesian closed category E, we prove that a local fibred right adjoint between slices of E is given by a polynomial. The slices in question are taken in a well known fibred sense......For any locally cartesian closed category E, we prove that a local fibred right adjoint between slices of E is given by a polynomial. The slices in question are taken in a well known fibred sense...
Laguerre polynomials method in the valon model
Boroun, G R
2014-01-01
We used the Laguerre polynomials method for determination of the proton structure function in the valon model. We have examined the applicability of the valon model with respect to a very elegant method, where the structure of the proton is determined by expanding valon distributions and valon structure functions on Laguerre polynomials. We compared our results with the experimental data, GJR parameterization and DL model. Having checked, this method gives a good description for the proton structure function in valon model.
Online segmentation of time series based on polynomial least-squares approximations.
Fuchs, Erich; Gruber, Thiemo; Nitschke, Jiri; Sick, Bernhard
2010-12-01
The paper presents SwiftSeg, a novel technique for online time series segmentation and piecewise polynomial representation. The segmentation approach is based on a least-squares approximation of time series in sliding and/or growing time windows utilizing a basis of orthogonal polynomials. This allows the definition of fast update steps for the approximating polynomial, where the computational effort depends only on the degree of the approximating polynomial and not on the length of the time window. The coefficients of the orthogonal expansion of the approximating polynomial-obtained by means of the update steps-can be interpreted as optimal (in the least-squares sense) estimators for average, slope, curvature, change of curvature, etc., of the signal in the time window considered. These coefficients, as well as the approximation error, may be used in a very intuitive way to define segmentation criteria. The properties of SwiftSeg are evaluated by means of some artificial and real benchmark time series. It is compared to three different offline and online techniques to assess its accuracy and runtime. It is shown that SwiftSeg-which is suitable for many data streaming applications-offers high accuracy at very low computational costs.
Vector-Valued Jack Polynomials from Scratch
Directory of Open Access Journals (Sweden)
Jean-Gabriel Luque
2011-03-01
Full Text Available Vector-valued Jack polynomials associated to the symmetric group S_N are polynomials with multiplicities in an irreducible module of S_N and which are simultaneous eigenfunctions of the Cherednik-Dunkl operators with some additional properties concerning the leading monomial. These polynomials were introduced by Griffeth in the general setting of the complex reflections groups G(r,p,N and studied by one of the authors (C. Dunkl in the specialization r=p=1 (i.e. for the symmetric group. By adapting a construction due to Lascoux, we describe an algorithm allowing us to compute explicitly the Jack polynomials following a Yang-Baxter graph. We recover some properties already studied by C. Dunkl and restate them in terms of graphs together with additional new results. In particular, we investigate normalization, symmetrization and antisymmetrization, polynomials with minimal degree, restriction etc. We give also a shifted version of the construction and we discuss vanishing properties of the associated polynomials.
On λ-Bell polynomials associated with umbral calculus
Kim, T.; Kim, D. S.
2017-01-01
In this paper, we introduce some new λ-Bell polynomials and Bell polynomials of the second kind and investigate properties of these polynomials. Using our investigation, we derive some new identities for the two kinds of λ-Bell polynomials arising from umbral calculus.
A generalization of the dichromatic polynomial of a graph
1981-01-01
The Subgraph polynomial fo a graph pair (G,H), where H⫅G, is defined. By assigning particular weights to the variables, it is shown that this polynomial reduces to the dichromatic polynomial of G. This idea of a graph pair leads to a dual generalization of the dichromatic polynomial.
Interpolation on Real Algebraic Curves to Polynomial Data
Directory of Open Access Journals (Sweden)
Len Bos
2013-04-01
Full Text Available We discuss a polynomial interpolation problem where the data are of the form of a set of algebraic curves in R^2 on each of which is prescribed a polynomial. The object is then to construct a global bivariate polynomial that agrees with the given polynomials when restricted to the corresponding curves.
Interpolation Functions of -Extensions of Apostol's Type Euler Polynomials
Directory of Open Access Journals (Sweden)
Kim Young-Hee
2009-01-01
Full Text Available The main purpose of this paper is to present new -extensions of Apostol's type Euler polynomials using the fermionic -adic integral on . We define the - -Euler polynomials and obtain the interpolation functions and the Hurwitz type zeta functions of these polynomials. We define -extensions of Apostol type's Euler polynomials of higher order using the multivariate fermionic -adic integral on . We have the interpolation functions of these - -Euler polynomials. We also give -extensions of Apostol's type Euler polynomials of higher order and have the multiple Hurwitz type zeta functions of these - -Euler polynomials.
Colored knot amplitudes and Hall-Littlewood polynomials
Shakirov, Sh
2013-01-01
The amplitudes of refined Chern-Simons (CS) theory, colored by antisymmetric (or symmetric) representations, conjecturally generate the Lambda^r- (or S^r-) colored triply graded homology of (n,m) torus knots. This paper is devoted to the generalization of Rosso-Jones formula to refined amplitudes, that involves non-trivial Gamma-factors -- expansion coefficients in the Macdonald basis. We derive from refined CS theory a linear recursion w.r.t. transformations (n,m) -> (n, n+m) and (n,m) -> (m,-n) that fully determines these factors. Applying this recursion to (n,nk+1) torus knots colored by antisymmetric representations [1^r] we prove that their amplitudes are rectangular [n^r] Hall-Littlewood polynomials under k units of framing (a.k.a. the Garsia-Haiman nabla) operator. For symmetric representations [r], we find the dual -- q-Whittaker -- polynomials. These results confirm and give a colored extension of the observation of arXiv:1201.3339 that triply graded homology of many torus knots has a strikingly simp...
Notes on the Polynomial Identities in Random Overlap Structures
Sollich, Peter; Barra, Adriano
2012-04-01
In these notes we review first in some detail the concept of random overlap structure (ROSt) applied to fully connected and diluted spin glasses. We then sketch how to write down the general term of the expansion of the energy part from the Boltzmann ROSt (for the Sherrington-Kirkpatrick model) and the corresponding term from the RaMOSt, which is the diluted extension suitable for the Viana-Bray model. From the ROSt energy term, a set of polynomial identities (often known as Aizenman-Contucci or AC relations) is shown to hold rigorously at every order because of a recursive structure of these polynomials that we prove. We show also, however, that this set is smaller than the full set of AC identities that is already known. Furthermore, when investigating the RaMOSt energy for the diluted counterpart, at higher orders, combinations of such AC identities appear, ultimately suggesting a crucial role for the entropy in generating these constraints in spin glasses.
1970-10-01
sale: is disributici is unlimited = F’)RIWRD Seior Ignacio Soto, Rrecutive President, Instituto Mexicano del Cementc y Concreto , invited Mr. Bryant... Concreto , a.c., Kwidco, D. F., Mexico. Based on info.mation largely obtained from ACT Committee 223, Expansive ’ement. Concretes, ACI Journal, August 1Q70
Wavelet approach to accelerator problems. 1: Polynomial dynamics
Energy Technology Data Exchange (ETDEWEB)
Fedorova, A.; Zeitlin, M. [Russian Academy of Sciences, St. Petersburg (Russian Federation). Inst. of Problems of Mechanical Engineering; Parsa, Z. [Brookhaven National Lab., Upton, NY (United States). Dept. of Physics
1997-05-01
This is the first part of a series of talks in which the authors present applications of methods from wavelet analysis to polynomial approximations for a number of accelerator physics problems. In the general case they have the solution as a multiresolution expansion in the base of compactly supported wavelet basis. The solution is parameterized by solutions of two reduced algebraical problems, one is nonlinear and the second is some linear problem, which is obtained from one of the next wavelet constructions: Fast Wavelet Transform, Stationary Subdivision Schemes, the method of Connection Coefficients. In this paper the authors consider the problem of calculation of orbital motion in storage rings. The key point in the solution of this problem is the use of the methods of wavelet analysis, relatively novel set of mathematical methods, which gives one a possibility to work with well-localized bases in functional spaces and with the general type of operators (including pseudodifferential) in such bases.
Uncertainty Analysis via Failure Domain Characterization: Polynomial Requirement Functions
Crespo, Luis G.; Munoz, Cesar A.; Narkawicz, Anthony J.; Kenny, Sean P.; Giesy, Daniel P.
2011-01-01
This paper proposes an uncertainty analysis framework based on the characterization of the uncertain parameter space. This characterization enables the identification of worst-case uncertainty combinations and the approximation of the failure and safe domains with a high level of accuracy. Because these approximations are comprised of subsets of readily computable probability, they enable the calculation of arbitrarily tight upper and lower bounds to the failure probability. A Bernstein expansion approach is used to size hyper-rectangular subsets while a sum of squares programming approach is used to size quasi-ellipsoidal subsets. These methods are applicable to requirement functions whose functional dependency on the uncertainty is a known polynomial. Some of the most prominent features of the methodology are the substantial desensitization of the calculations from the uncertainty model assumed (i.e., the probability distribution describing the uncertainty) as well as the accommodation for changes in such a model with a practically insignificant amount of computational effort.
Explicit classes of permutation polynomials of F33m
Institute of Scientific and Technical Information of China (English)
无
2009-01-01
Permutation polynomials have been an interesting subject of study for a long time and have applications in many areas of mathematics and engineering. However, only a small number of specific classes of permutation polynomials are known so far. In this paper, six classes of linearized permutation polynomials and six classes of nonlinearized permutation polynomials over F33m are presented. These polynomials have simple shapes, and they are related to planar functions.
Explicit classes of permutation polynomials of F33m
Institute of Scientific and Technical Information of China (English)
DING CunSheng; XIANG Qing; YUAN Jin; YUAN PingZhi
2009-01-01
Permutation polynomials have been an interesting subject of study for a long time and have applications in many areas of mathematics and engineering. However, only a small number of specific classes of permutation polynomials are known so far. In this paper, six classes of linearized permutation polynomials and six classes of nonlinearized permutation polynomials over F33 are pre-sented. These polynomials have simple shapes, and they are related to planar functions.
Certain non-linear differential polynomials sharing a non zero polynomial
Directory of Open Access Journals (Sweden)
Majumder Sujoy
2015-10-01
functions sharing a nonzero polynomial and obtain two results which improves and generalizes the results due to L. Liu [Uniqueness of meromorphic functions and differential polynomials, Comput. Math. Appl., 56 (2008, 3236-3245.] and P. Sahoo [Uniqueness and weighted value sharing of meromorphic functions, Applied. Math. E-Notes., 11 (2011, 23-32.].
A new class of generalized polynomials associated with Hermite and Bernoulli polynomials
Directory of Open Access Journals (Sweden)
M. A. Pathan
2015-05-01
Full Text Available In this paper, we introduce a new class of generalized polynomials associated with the modified Milne-Thomson's polynomials Φ_{n}^{(α}(x,ν of degree n and order α introduced by Derre and Simsek.The concepts of Bernoulli numbers B_n, Bernoulli polynomials B_n(x, generalized Bernoulli numbers B_n(a,b, generalized Bernoulli polynomials B_n(x;a,b,c of Luo et al, Hermite-Bernoulli polynomials {_HB}_n(x,y of Dattoli et al and {_HB}_n^{(α} (x,y of Pathan are generalized to the one {_HB}_n^{(α}(x,y,a,b,c which is called the generalized polynomial depending on three positive real parameters. Numerous properties of these polynomials and some relationships between B_n, B_n(x, B_n(a,b, B_n(x;a,b,c and {}_HB_n^{(α}(x,y;a,b,c are established. Some implicit summation formulae and general symmetry identities are derived by using different analytical means and applying generating functions. These results extend some known summations and identities of generalized Bernoulli numbers and polynomials
Extending a Property of Cubic Polynomials to Higher-Degree Polynomials
Miller, David A.; Moseley, James
2012-01-01
In this paper, the authors examine a property that holds for all cubic polynomials given two zeros. This property is discovered after reviewing a variety of ways to determine the equation of a cubic polynomial given specific conditions through algebra and calculus. At the end of the article, they will connect the property to a very famous method…
Orthonormal Polynomial Approximation of Mineral Water Data with Errors in Both Variables
Bogdanova, Nina
2012-01-01
In this paper we introduce the data from mineral water probe with errors in both variables. For this case we apply our orthonormal polynomial expansion(OPEM) method to describe the data in the new error corridor. It receives the approximating curves and their derivatives including the errors by weighting approach. The numerical method and approximation results are presented and discussed. The special criteria are carried out for orthonormal and evaluated from it usual expansion. The numerical results are shown in tables and figures.
DEFF Research Database (Denmark)
Ribard, Nicolas; Wisniewski, Rafael; Sloth, Christoffer
2016-01-01
In the paper, we strive to develop an algorithm that simultaneously computes a polynomial control and a polynomial Lyapunov function. This ensures asymptotic stability of the designed feedback system. The above problem is translated to a certificate of positivity. To this end, we use the represen......In the paper, we strive to develop an algorithm that simultaneously computes a polynomial control and a polynomial Lyapunov function. This ensures asymptotic stability of the designed feedback system. The above problem is translated to a certificate of positivity. To this end, we use...... the representation of the given control system in Bernstein basis. Subsequently, the control synthesis problem is reduced to finite number of evaluations of a polynomial on vertices of cubes in the space of parameters representing admissible controls and Lyapunov functions....
Lattice harmonics expansion revisited
Kontrym-Sznajd, G.; Holas, A.
2017-04-01
The main subject of the work is to provide the most effective way of determining the expansion of some quantities into orthogonal polynomials, when these quantities are known only along some limited number of sampling directions. By comparing the commonly used Houston method with the method based on the orthogonality relation, some relationships, which define the applicability and correctness of these methods, are demonstrated. They are verified for various sets of sampling directions applicable for expanding quantities having the full symmetry of the Brillouin zone of cubic and non-cubic lattices. All results clearly show that the Houston method is always better than the orthogonality-relation one. For the cubic symmetry we present a few sets of special directions (SDs) showing how their construction and, next, a proper application depend on the choice of various sets of lattice harmonics. SDs are important mainly for experimentalists who want to reconstruct anisotropic quantities from their measurements, performed at a limited number of sampling directions.
Twisted Alexander polynomials of hyperbolic knots
Dunfield, Nathan M; Jackson, Nicholas
2011-01-01
We study a twisted Alexander polynomial naturally associated to a hyperbolic knot in an integer homology 3-sphere via a lift of the holonomy representation to SL(2, C). It is an unambiguous symmetric Laurent polynomial whose coefficients lie in the trace field of the knot. It contains information about genus, fibering, and chirality, and moreover is powerful enough to sometimes detect mutation. We calculated this invariant numerically for all 313,209 hyperbolic knots in S^3 with at most 15 crossings, and found that in all cases it gave a sharp bound on the genus of the knot and determined both fibering and chirality. We also study how such twisted Alexander polynomials vary as one moves around in an irreducible component X_0 of the SL(2, C)-character variety of the knot group. We show how to understand all of these polynomials at once in terms of a polynomial whose coefficients lie in the function field of X_0. We use this to help explain some of the patterns observed for knots in S^3, and explore a potential...
Leading log expansion of combinatorial Dyson Schwinger equations
Delage, Lucas
2016-01-01
We study combinatorial Dyson Schwinger equations, expressed in the Hopf algebra of words with a quasi shuffle product. We map them into an algebra of polynomials in one indeterminate L and show that the leading log expansion one obtains with such a mapping are simple power law like expression
Using Taylor Expansions to Prepare Students for Calculus
Lutzer, Carl V.
2011-01-01
We propose an alternative to the standard introduction to the derivative. Instead of using limits of difference quotients, students develop Taylor expansions of polynomials. This alternative allows students to develop many of the central ideas about the derivative at an intuitive level, using only skills and concepts from precalculus, and…
Uniform trigonometric polynomial B-spline curves
Institute of Scientific and Technical Information of China (English)
吕勇刚; 汪国昭; 杨勋年
2002-01-01
This paper presents a new kind of uniform spline curve, named trigonometric polynomialB-splines, over space Ω = span{sint, cost, tk-3,tk-4,…,t,1} of which k is an arbitrary integerlarger than or equal to 3. We show that trigonometric polynomial B-spline curves have many similarV properties to traditional B-splines. Based on the explicit representation of the curve we have also presented the subdivision formulae for this new kind of curve. Since the new spline can include both polynomial curves and trigonometric curves as special cases without rational form, it can be used as an efficient new model for geometric design in the fields of CAD/CAM.
Transversals of Complex Polynomial Vector Fields
DEFF Research Database (Denmark)
Dias, Kealey
by rotational constants. Transversals are a certain class of curves for such a family of vector fields that represent the bifurcation states for this family of vector fields. More specifically, transversals are curves that coincide with a homoclinic separatrix for some rotation of the vector field. Given......, an important step was proving that the transversals possessed a certain characteristic. Understanding transversals might be the key to proving other polynomial vector fields are generic, and they are important in understanding bifurcations of polynomial vector fields in general. We consider two important...... examples of rotated families to argue this. There will be discussed several open questions concerning the number of transversals that can appear for a certain degree d of a polynomial vector field, and furthermore how transversals are analyzed with respect to bifurcations around multiple equilibrium points....
Minimal residual method stronger than polynomial preconditioning
Energy Technology Data Exchange (ETDEWEB)
Faber, V.; Joubert, W.; Knill, E. [Los Alamos National Lab., NM (United States)] [and others
1994-12-31
Two popular methods for solving symmetric and nonsymmetric systems of equations are the minimal residual method, implemented by algorithms such as GMRES, and polynomial preconditioning methods. In this study results are given on the convergence rates of these methods for various classes of matrices. It is shown that for some matrices, such as normal matrices, the convergence rates for GMRES and for the optimal polynomial preconditioning are the same, and for other matrices such as the upper triangular Toeplitz matrices, it is at least assured that if one method converges then the other must converge. On the other hand, it is shown that matrices exist for which restarted GMRES always converges but any polynomial preconditioning of corresponding degree makes no progress toward the solution for some initial error. The implications of these results for these and other iterative methods are discussed.
Fast beampattern evaluation by polynomial rooting
Häcker, P.; Uhlich, S.; Yang, B.
2011-07-01
Current automotive radar systems measure the distance, the relative velocity and the direction of objects in their environment. This information enables the car to support the driver. The direction estimation capabilities of a sensor array depend on its beampattern. To find the array configuration leading to the best angle estimation by a global optimization algorithm, a huge amount of beampatterns have to be calculated to detect their maxima. In this paper, a novel algorithm is proposed to find all maxima of an array's beampattern fast and reliably, leading to accelerated array optimizations. The algorithm works for arrays having the sensors on a uniformly spaced grid. We use a general version of the gcd (greatest common divisor) function in order to write the problem as a polynomial. We differentiate and root the polynomial to get the extrema of the beampattern. In addition, we show a method to reduce the computational burden even more by decreasing the order of the polynomial.
A complete discrimination system for polynomials
Institute of Scientific and Technical Information of China (English)
杨路; 侯晓荣; 曾振柄
1996-01-01
Given a polynomial with symbolic/literal coefficients,a complete discrimination system is a set of explicit expressions in terms of the coefficients,which is sufficient for determining the numbers and multiplicities of the real and imaginary roots.Though it is of great significance,such a criterion for root-classification has never been given for polynomials with degrees greater than 4.The lack of efficient tools in this aspect extremely prevents computer implementations for Tarski’s and other methods in automated theorem proving.To remedy this defect,a generic algorithm is proposed to produce a complete discrimination system for a polynomial with any degrees.This result has extensive applications in various fields,and its efficiency was demonstrated by computer implementations.
Dominating Sets and Domination Polynomials of Paths
Directory of Open Access Journals (Sweden)
Saeid Alikhani
2009-01-01
Full Text Available Let G=(V,E be a simple graph. A set S⊆V is a dominating set of G, if every vertex in V\\S is adjacent to at least one vertex in S. Let 𝒫ni be the family of all dominating sets of a path Pn with cardinality i, and let d(Pn,j=|𝒫nj|. In this paper, we construct 𝒫ni, and obtain a recursive formula for d(Pn,i. Using this recursive formula, we consider the polynomial D(Pn,x=∑i=⌈n/3⌉nd(Pn,ixi, which we call domination polynomial of paths and obtain some properties of this polynomial.
Degroote, Matthias; Henderson, Thomas M.; Zhao, Jinmo; Dukelsky, Jorge; Scuseria, Gustavo E.
2016-03-01
We present a similarity transformation theory based on a polynomial form of a particle-hole pair excitation operator. In the weakly correlated limit, this polynomial becomes an exponential, leading to coupled cluster doubles. In the opposite strongly correlated limit, the polynomial becomes an extended Bessel expansion and yields the projected BCS wave function. In between, we interpolate using a single parameter. The effective Hamiltonian is non-Hermitian and this polynomial similarity transformation theory follows the philosophy of traditional coupled cluster, left projecting the transformed Hamiltonian onto subspaces of the Hilbert space in which the wave function variance is forced to be zero. Similarly, the interpolation parameter is obtained through minimizing the next residual in the projective hierarchy. We rationalize and demonstrate how and why coupled cluster doubles is ill suited to the strongly correlated limit, whereas the Bessel expansion remains well behaved. The model provides accurate wave functions with energy errors that in its best variant are smaller than 1% across all interaction strengths. The numerical cost is polynomial in system size and the theory can be straightforwardly applied to any realistic Hamiltonian.
Incomplete Bivariate Fibonacci and Lucas -Polynomials
Directory of Open Access Journals (Sweden)
Dursun Tasci
2012-01-01
Full Text Available We define the incomplete bivariate Fibonacci and Lucas -polynomials. In the case =1, =1, we obtain the incomplete Fibonacci and Lucas -numbers. If =2, =1, we have the incomplete Pell and Pell-Lucas -numbers. On choosing =1, =2, we get the incomplete generalized Jacobsthal number and besides for =1 the incomplete generalized Jacobsthal-Lucas numbers. In the case =1, =1, =1, we have the incomplete Fibonacci and Lucas numbers. If =1, =1, =1, =⌊(−1/(+1⌋, we obtain the Fibonacci and Lucas numbers. Also generating function and properties of the incomplete bivariate Fibonacci and Lucas -polynomials are given.
Spread polynomials, rotations and the butterfly effect
Goh, Shuxiang
2009-01-01
The spread between two lines in rational trigonometry replaces the concept of angle, allowing the complete specification of many geometrical and dynamical situations which have traditionally been viewed approximately. This paper investigates the case of powers of a rational spread rotation, and in particular, a curious periodicity in the prime power decomposition of the associated values of the spread polynomials, which are the analogs in rational trigonometry of the Chebyshev polynomials of the first kind. Rational trigonometry over finite fields plays a role, together with non-Euclidean geometries.
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.
Knot polynomial identities and quantum group coincidences
Morrison, Scott; Snyder, Noah
2010-01-01
We construct link invariants using the D_2n subfactor planar algebras, and use these to prove new identities relating certain specializations of colored Jones polynomials to specializations of other quantum knot polynomials. These identities can also be explained by coincidences between small modular categories involving the even parts of the D_2n planar algebras. We discuss the origins of these coincidences, explaining the role of SO level-rank duality, Kirby-Melvin symmetry, and properties of small Dynkin diagrams. One of these coincidences involves G_2 and does not appear to be related to level-rank duality.
Some Inequalities for the Derivative of Polynomials
Directory of Open Access Journals (Sweden)
Sunil Hans
2014-01-01
Full Text Available If pz=∑υ=0ncυzυ is a polynomial of degree n, having no zeros in z<1, then Aziz (1989 proved maxz=1p′z≤n/2Mα2+Mα+π21/2, where Mα=max1≤k≤npeiα+2kπ/n. In this paper, we consider a class of polynomial Pnμ of degree n, defined as pz=a0+∑υ=μnaυzυ and present certain generalizations of above inequality and some other well-known results.
Polynomial Kernelizations for $\\MINF_1$ and $\\MNP$
Kratsch, Stefan
2009-01-01
The relation of constant-factor approximability to fixed-parameter tractability and kernelization is a long-standing open question. We prove that two large classes of constant-factor approximable problems, namely $\\MINF_1$ and $\\MNP$, including the well-known subclass $\\MSNP$, admit polynomial kernelizations for their natural decision versions. This extends results of Cai and Chen (JCSS 1997), stating that the standard parameterizations of problems in $\\MSNP$ and $\\MINF_1$ are fixed-parameter tractable, and complements recent research on problems that do not admit polynomial kernelizations (Bodlaender et al. ICALP 2008).
The classification of polynomial basins of infinity
DeMarco, Laura
2011-01-01
We consider the problem of classifying the dynamics of complex polynomials $f: \\mathbb{C} \\to \\mathbb{C}$ restricted to their basins of infinity. We synthesize existing combinatorial tools --- tableaux, trees, and laminations --- into a new invariant of basin dynamics we call the pictograph. For polynomials with all critical points escaping to infinity, we obtain a complete description of the set of topological conjugacy classes. We give an algorithm for constructing abstract pictographs, and we provide an inductive algorithm for counting topological conjugacy classes with a given pictograph.
High degree interpolation polynomial in Newton form
Tal-Ezer, Hillel
1988-01-01
Polynomial interpolation is an essential subject in numerical analysis. Dealing with a real interval, it is well known that even if f(x) is an analytic function, interpolating at equally spaced points can diverge. On the other hand, interpolating at the zeroes of the corresponding Chebyshev polynomial will converge. Using the Newton formula, this result of convergence is true only on the theoretical level. It is shown that the algorithm which computes the divided differences is numerically stable only if: (1) the interpolating points are arranged in a different order, and (2) the size of the interval is 4.
The chromatic polynomial and list colorings
DEFF Research Database (Denmark)
Thomassen, Carsten
2009-01-01
We prove that, if a graph has a list of k available colors at every vertex, then the number of list-colorings is at least the chromatic polynomial evaluated at k when k is sufficiently large compared to the number of vertices of the graph.......We prove that, if a graph has a list of k available colors at every vertex, then the number of list-colorings is at least the chromatic polynomial evaluated at k when k is sufficiently large compared to the number of vertices of the graph....
Perturbations around the zeros of classical orthogonal polynomials
Sasaki, Ryu
2014-01-01
Starting from degree N solutions of a time dependent Schroedinger-like equation for classical orthogonal polynomials, a linear matrix equation describing perturbations around the N zeros of the polynomial is derived. The matrix has remarkable Diophantine properties. Its eigenvalues are independent of the zeros. The corresponding eigenvectors provide the representations of the lower degree (0,1,...,N-1) polynomials in terms of the zeros of the degree N polynomial. The results are valid universally for all the classical orthogonal polynomials, including the Askey scheme of hypergeometric orthogonal polynomials and its q-analogues.
Polynomial Representations for a Wavelet Model of Interest Rates
Directory of Open Access Journals (Sweden)
Dennis G. Llemit
2015-12-01
Full Text Available In this paper, we approximate a non – polynomial function which promises to be an essential tool in interest rates forecasting in the Philippines. We provide two numerical schemes in order to generate polynomial functions that approximate a new wavelet which is a modification of Morlet and Mexican Hat wavelets. The first is the Polynomial Least Squares method which approximates the underlying wavelet according to desired numerical errors. The second is the Chebyshev Polynomial approximation which generates the required function through a sequence of recursive and orthogonal polynomial functions. We seek to determine the lowest order polynomial representations of this wavelet corresponding to a set of error thresholds.
Institute of Scientific and Technical Information of China (English)
Tan Xiaogang; Wei Ping; Li Liping
2009-01-01
To detect higher order polynomial phase signals (HOPPSs), the smoothed-pseudo polynomial Wigner-Ville distribution (SP-PWVD), an improved version of the polynomial Wigner-Ville distribution (PWVD), is pre-sented using a separable kernel. By adjusting the lengths of the functions in the kernel, the balance between resolution retaining and interference suppressing can be adjusted conveniently. The proposed method with merits of interference terms reduction and noise suppression can provide time frequency representation of better readability and more accurate instantaneous frequency (IF) estimation with higher order SP-PWVD. The performance of the SP-PWVD is verified by computer simulations.
Borot, Gaëtan
2012-01-01
We propose a conjecture to compute the all-order asymptotic expansion of the colored Jones polynomial of the complement of a hyperbolic knot, J_N(q = exp(2u/N)) when N goes to infinity. Our conjecture claims that the asymptotic expansion of the colored Jones polynomial is a the formal wave function of an integrable system whose semiclassical spectral curve S would be the SL_2(C) character variety of the knot (the A-polynomial), and is formulated in the framework of the topological recursion. It takes as starting point the proposal made recently by Dijkgraaf, Fuji and Manabe (who kept only the perturbative part of the wave function, and found some discrepancies), but it also contains the non-perturbative parts, and solves the discrepancy problem. These non-perturbative corrections are derivatives of Theta functions associated to S, but the expansion is still in powers of 1/N due to the special properties of A-polynomials. We provide a detailed check for the figure-eight knot and the once-punctured torus bundle...
Convergent series for lattice models with polynomial interactions
Directory of Open Access Journals (Sweden)
Aleksandr S. Ivanov
2017-01-01
Full Text Available The standard perturbative weak-coupling expansions in lattice models are asymptotic. The reason for this is hidden in the incorrect interchange of the summation and integration. However, substituting the Gaussian initial approximation of the perturbative expansions by a certain interacting model or regularizing original lattice integrals, one can construct desired convergent series. In this paper we develop methods, which are based on the joint and separate utilization of the regularization and new initial approximation. We prove, that the convergent series exist and can be expressed as re-summed standard perturbation theory for any model on the finite lattice with the polynomial interaction of even degree. We discuss properties of such series and study their applicability to practical computations on the example of the lattice ϕ4-model. We calculate 〈ϕn2〉 expectation value using the convergent series, the comparison of the results with the Borel re-summation and Monte Carlo simulations shows a good agreement between all these methods.
Directory of Open Access Journals (Sweden)
Waleed M. Abd-Elhameed
2016-09-01
Full Text Available Herein, two numerical algorithms for solving some linear and nonlinear fractional-order differential equations are presented and analyzed. For this purpose, a novel operational matrix of fractional-order derivatives of Fibonacci polynomials was constructed and employed along with the application of the tau and collocation spectral methods. The convergence and error analysis of the suggested Fibonacci expansion were carefully investigated. Some numerical examples with comparisons are presented to ensure the efficiency, applicability and high accuracy of the proposed algorithms. Two accurate semi-analytic polynomial solutions for linear and nonlinear fractional differential equations are the result.
The 6 Vertex Model and Schubert Polynomials
Directory of Open Access Journals (Sweden)
Alain Lascoux
2007-02-01
Full Text Available We enumerate staircases with fixed left and right columns. These objects correspond to ice-configurations, or alternating sign matrices, with fixed top and bottom parts. The resulting partition functions are equal, up to a normalization factor, to some Schubert polynomials.
Dynamic system uncertainty propagation using polynomial chaos
Institute of Scientific and Technical Information of China (English)
Xiong Fenfen; Chen Shishi; Xiong Ying
2014-01-01
The classic polynomial chaos method (PCM), characterized as an intrusive methodology, has been applied to uncertainty propagation (UP) in many dynamic systems. However, the intrusive polynomial chaos method (IPCM) requires tedious modification of the governing equations, which might introduce errors and can be impractical. Alternative to IPCM, the non-intrusive polynomial chaos method (NIPCM) that avoids such modifications has been developed. In spite of the frequent application to dynamic problems, almost all the existing works about NIPCM for dynamic UP fail to elaborate the implementation process in a straightforward way, which is important to readers who are unfamiliar with the mathematics of the polynomial chaos theory. Meanwhile, very few works have compared NIPCM to IPCM in terms of their merits and applicability. Therefore, the mathematic procedure of dynamic UP via both methods considering parametric and initial condition uncertainties are comparatively discussed and studied in the present paper. Comparison of accuracy and efficiency in statistic moment estimation is made by applying the two methods to several dynamic UP problems. The relative merits of both approaches are discussed and summarized. The detailed description and insights gained with the two methods through this work are expected to be helpful to engineering designers in solving dynamic UP problems.
On the Schinzel Identity of Cyclotomic Polynomial
Institute of Scientific and Technical Information of China (English)
无
2000-01-01
@@For integer n>0, let n(x) denote the nth cyclotomic polynomial n(x)=tackrel{01 be an odd square-free number.Aurifeuille and Le Lasseur［1］ proved thatequationn(x)=An2(x)-(-1)n-12)nxBn2(x).equation
A note on Fibonacci-type polynomials
Amdeberhan, Tewodros
2008-01-01
We opt to study the convergence of maximal real roots of certain Fibonacci-type polynomials given by $G_n=x^kG_{n-1}+G_{n-2}$. The special cases $k=1$ and $k=2$ are found in [4] and [7], respectively.
Algebraic polynomial system solving and applications
Bleylevens, I.W.M.
2010-01-01
The problem of computing the solutions of a system of multivariate polynomial equations can be approached by the Stetter-Möller matrix method which casts the problem into a large eigenvalue problem. This Stetter-Möller matrix method forms the starting point for the development of computational proce
Quantum Hilbert matrices and orthogonal polynomials
DEFF Research Database (Denmark)
Andersen, Jørgen Ellegaard; Berg, Christian
2009-01-01
Using the notion of quantum integers associated with a complex number q≠0 , we define the quantum Hilbert matrix and various extensions. They are Hankel matrices corresponding to certain little q -Jacobi polynomials when |q|matrices...... of reciprocal Fibonacci numbers called Filbert matrices. We find a formula for the entries of the inverse quantum Hilbert matrix....
A recursive algorithm for Zernike polynomials
Davenport, J. W.
1982-01-01
The analysis of a function defined on a rotationally symmetric system, with either a circular or annular pupil is discussed. In order to numerically analyze such systems it is typical to expand the given function in terms of a class of orthogonal polynomials. Because of their particular properties, the Zernike polynomials are especially suited for numerical calculations. Developed is a recursive algorithm that can be used to generate the Zernike polynomials up to a given order. The algorithm is recursively defined over J where R(J,N) is the Zernike polynomial of degree N obtained by orthogonalizing the sequence R(J), R(J+2), ..., R(J+2N) over (epsilon, 1). The terms in the preceding row - the (J-1) row - up to the N+1 term is needed for generating the (J,N)th term. Thus, the algorith generates an upper left-triangular table. This algorithm was placed in the computer with the necessary support program also included.
The GCD property and irreduciable quadratic polynomials
Directory of Open Access Journals (Sweden)
Saroj Malik
1986-01-01
Full Text Available The proof of the following theorem is presented: If D is, respectively, a Krull domain, a Dedekind domain, or a Prüfer domain, then D is correspondingly a UFD, a PID, or a Bezout domain if and only if every irreducible quadratic polynomial in D[X] is a prime element.
The Tutte Polynomial of Some Matroids
Directory of Open Access Journals (Sweden)
Criel Merino
2012-01-01
graphs or matroids. In this work, we compile known formulas for the Tutte polynomial of some families of graphs and matroids. Also, we give brief explanations of the techniques that were used to find the formulas. Hopefully, this will be useful for researchers in Combinatorics and elsewhere.
Orthogonality Relations for Multivariate Krawtchouk Polynomials
Directory of Open Access Journals (Sweden)
Hiroshi Mizukawa
2011-02-01
Full Text Available The orthogonality relations of multivariate Krawtchouk polynomials are discussed. In case of two variables, the necessary and sufficient conditions of orthogonality is given by Grünbaum and Rahman in [SIGMA 6 (2010, 090, 12 pages]. In this study, a simple proof of the necessary and sufficient condition of orthogonality is given for a general case.
UNIQUENESS OF DIFFERENCE POLYNOMIALS OF MEROMORPHIC FUNCTIONS
Institute of Scientific and Technical Information of China (English)
刘永; 祁晓光
2014-01-01
In this article, we investigate the uniqueness problems of difference polynomials of meromorphic functions and obtain some results which can be viewed as discrete analogues of the results given by Shibazaki. Some examples are given to show the results in this article are best possible.
Dynamic system uncertainty propagation using polynomial chaos
Directory of Open Access Journals (Sweden)
Xiong Fenfen
2014-10-01
Full Text Available The classic polynomial chaos method (PCM, characterized as an intrusive methodology, has been applied to uncertainty propagation (UP in many dynamic systems. However, the intrusive polynomial chaos method (IPCM requires tedious modification of the governing equations, which might introduce errors and can be impractical. Alternative to IPCM, the non-intrusive polynomial chaos method (NIPCM that avoids such modifications has been developed. In spite of the frequent application to dynamic problems, almost all the existing works about NIPCM for dynamic UP fail to elaborate the implementation process in a straightforward way, which is important to readers who are unfamiliar with the mathematics of the polynomial chaos theory. Meanwhile, very few works have compared NIPCM to IPCM in terms of their merits and applicability. Therefore, the mathematic procedure of dynamic UP via both methods considering parametric and initial condition uncertainties are comparatively discussed and studied in the present paper. Comparison of accuracy and efficiency in statistic moment estimation is made by applying the two methods to several dynamic UP problems. The relative merits of both approaches are discussed and summarized. The detailed description and insights gained with the two methods through this work are expected to be helpful to engineering designers in solving dynamic UP problems.
Scalar Field Theories with Polynomial Shift Symmetries
Griffin, Tom; Horava, Petr; Yan, Ziqi
2014-01-01
We continue our study of naturalness in nonrelativistic QFTs of the Lifshitz type, focusing on scalar fields that can play the role of Nambu-Goldstone (NG) modes associated with spontaneous symmetry breaking. Such systems allow for an extension of the constant shift symmetry to a shift by a polynomial of degree $P$ in spatial coordinates. These "polynomial shift symmetries" in turn protect the technical naturalness of modes with a higher-order dispersion relation, and lead to a refinement of the proposed classification of infrared Gaussian fixed points available to describe NG modes in nonrelativistic theories. Generic interactions in such theories break the polynomial shift symmetry explicitly to the constant shift. It is thus natural to ask: Given a Gaussian fixed point with polynomial shift symmetry of degree $P$, what are the lowest-dimension operators that preserve this symmetry, and deform the theory into a self-interacting scalar field theory with the shift symmetry of degree $P$? To answer this (essen...
Ideals in Polynomial Near-rings
Institute of Scientific and Technical Information of China (English)
Mark Farag
2002-01-01
In this paper, we further explore the relationship between the ideals of N and those of N[x], where N is a zero-symmetric right near-ring with identity and N[x] is the polynomial near-ring introduced by Bagley in 1993.
Optimization of Cubic Polynomial Functions without Calculus
Taylor, Ronald D., Jr.; Hansen, Ryan
2008-01-01
In algebra and precalculus courses, students are often asked to find extreme values of polynomial functions in the context of solving an applied problem; but without the notion of derivative, something is lost. Either the functions are reduced to quadratics, since students know the formula for the vertex of a parabola, or solutions are…
Polynomial Vector Fields in One Complex Variable
DEFF Research Database (Denmark)
Branner, Bodil
In recent years Adrien Douady was interested in polynomial vector fields, both in relation to iteration theory and as a topic on their own. This talk is based on his work with Pierrette Sentenac, work of Xavier Buff and Tan Lei, and my own joint work with Kealey Dias....
Piecewise polynomial representations of genomic tracks.
Tarabichi, Maxime; Detours, Vincent; Konopka, Tomasz
2012-01-01
Genomic data from micro-array and sequencing projects consist of associations of measured values to chromosomal coordinates. These associations can be thought of as functions in one dimension and can thus be stored, analyzed, and interpreted as piecewise-polynomial curves. We present a general framework for building piecewise polynomial representations of genome-scale signals and illustrate some of its applications via examples. We show that piecewise constant segmentation, a typical step in copy-number analyses, can be carried out within this framework for both array and (DNA) sequencing data offering advantages over existing methods in each case. Higher-order polynomial curves can be used, for example, to detect trends and/or discontinuities in transcription levels from RNA-seq data. We give a concrete application of piecewise linear functions to diagnose and quantify alignment quality at exon borders (splice sites). Our software (source and object code) for building piecewise polynomial models is available at http://sourceforge.net/projects/locsmoc/.
Polynomial Asymptotes of the Second Kind
Dobbs, David E.
2011-01-01
This note uses the analytic notion of asymptotic functions to study when a function is asymptotic to a polynomial function. Along with associated existence and uniqueness results, this kind of asymptotic behaviour is related to the type of asymptote that was recently defined in a more geometric way. Applications are given to rational functions and…
Bernoulli Polynomials, Fourier Series and Zeta Numbers
DEFF Research Database (Denmark)
Scheufens, Ernst E
2013-01-01
Fourier series for Bernoulli polynomials are used to obtain information about values of the Riemann zeta function for integer arguments greater than one. If the argument is even we recover the well-known exact values, if the argument is odd we find integral representations and rapidly convergent...
Euler Polynomials, Fourier Series and Zeta Numbers
DEFF Research Database (Denmark)
Scheufens, Ernst E
2012-01-01
Fourier series for Euler polynomials is used to obtain information about values of the Riemann zeta function for integer arguments greater than one. If the argument is even we recover the well-known exact values, if the argument is odd we find integral representations and rapidly convergent series....
Polynomial Structure of Topological String Partition Functions
Zhou, Jie
2015-01-01
We review the polynomial structure of the topological string partition functions as solutions to the holomorphic anomaly equations. We also explain the connection between the ring of propagators defined from special K\\"ahler geometry and the ring of almost-holomorphic modular forms defined on modular curves.
On the optimal polynomial approximation of stochastic PDEs by galerkin and collocation methods
Beck, Joakim
2012-09-01
In this work we focus on the numerical approximation of the solution u of a linear elliptic PDE with stochastic coefficients. The problem is rewritten as a parametric PDE and the functional dependence of the solution on the parameters is approximated by multivariate polynomials. We first consider the stochastic Galerkin method, and rely on sharp estimates for the decay of the Fourier coefficients of the spectral expansion of u on an orthogonal polynomial basis to build a sequence of polynomial subspaces that features better convergence properties, in terms of error versus number of degrees of freedom, than standard choices such as Total Degree or Tensor Product subspaces. We consider then the Stochastic Collocation method, and use the previous estimates to introduce a new class of Sparse Grids, based on the idea of selecting a priori the most profitable hierarchical surpluses, that, again, features better convergence properties compared to standard Smolyak or tensor product grids. Numerical results show the effectiveness of the newly introduced polynomial spaces and sparse grids. © 2012 World Scientific Publishing Company.
Hong, X; Harris, C J
2000-01-01
This paper introduces a new neurofuzzy model construction algorithm for nonlinear dynamic systems based upon basis functions that are Bézier-Bernstein polynomial functions. This paper is generalized in that it copes with n-dimensional inputs by utilising an additive decomposition construction to overcome the curse of dimensionality associated with high n. This new construction algorithm also introduces univariate Bézier-Bernstein polynomial functions for the completeness of the generalized procedure. Like the B-spline expansion based neurofuzzy systems, Bézier-Bernstein polynomial function based neurofuzzy networks hold desirable properties such as nonnegativity of the basis functions, unity of support, and interpretability of basis function as fuzzy membership functions, moreover with the additional advantages of structural parsimony and Delaunay input space partition, essentially overcoming the curse of dimensionality associated with conventional fuzzy and RBF networks. This new modeling network is based on additive decomposition approach together with two separate basis function formation approaches for both univariate and bivariate Bézier-Bernstein polynomial functions used in model construction. The overall network weights are then learnt using conventional least squares methods. Numerical examples are included to demonstrate the effectiveness of this new data based modeling approach.
Polynomial Eigenvalue Solutions to Minimal Problems in Computer Vision.
Kukelova, Zuzana; Bujnak, Martin; Pajdla, Tomas
2012-07-01
We present a method for solving systems of polynomial equations appearing in computer vision. This method is based on polynomial eigenvalue solvers and is more straightforward and easier to implement than the state-of-the-art Gröbner basis method since eigenvalue problems are well studied, easy to understand, and efficient and robust algorithms for solving these problems are available. We provide a characterization of problems that can be efficiently solved as polynomial eigenvalue problems (PEPs) and present a resultant-based method for transforming a system of polynomial equations to a polynomial eigenvalue problem. We propose techniques that can be used to reduce the size of the computed polynomial eigenvalue problems. To show the applicability of the proposed polynomial eigenvalue method, we present the polynomial eigenvalue solutions to several important minimal relative pose problems.
On an Inequality Concerning the Polar Derivative of a Polynomial
Indian Academy of Sciences (India)
A Aziz; N A Rather
2007-08-01
In this paper, we present a correct proof of an -inequality concerning the polar derivative of a polynomial with restricted zeros. We also extend Zygmund’s inequality to the polar derivative of a polynomial.
d-Orthogonal Charlier Polynomials and the Weyl Algebra
Energy Technology Data Exchange (ETDEWEB)
Vinet, Luc [Centre de recherches mathematiques, Universite de Montreal, C.P. 6128, succ. Centre-ville, Montreal, Qc, H3C 3J7 (Canada); Zhedanov, Alexei, E-mail: luc.vinet@umontreal.ca, E-mail: zhedanov@yahoo.com [Department of Electronic and Kinetic Properties, Donetsk Institute for Physics and Technology R.Luxemburg str. 72, Donetsk, 83114 (Ukraine)
2011-03-01
It is shown that d-orthogonal Charlier polynomials arise as matrix elements of non unitary automorphisms of the Weyl algebra. The structural formulas that these polynomials obey are derived from this algebraic setting.
Combinatorial polynomials as moments, Hankel transforms and exponential Riordan arrays
Barry, Paul
2011-01-01
In the case of two combinatorial polynomials, we show that they can exhibited as moments of paramaterized families of orthogonal polynomials, and hence derive their Hankel transforms. Exponential Riordan arrays are the main vehicles used for this.
Irreducibility Results for Compositions of Polynomials in Several Variables
Indian Academy of Sciences (India)
Anca Iuliana Bonciocat; Alexandru Zaharescu
2005-05-01
We obtain explicit upper bounds for the number of irreducible factors for a class of compositions of polynomials in several variables over a given field. In particular, some irreducibility criteria are given for this class of compositions of polynomials.
Representations of Knot Groups and Twisted Alexander Polynomials
Institute of Scientific and Technical Information of China (English)
Xiao Song LIN
2001-01-01
We present a twisted version of the Alexander polynomial associated with a matrix representation of the knot group. Examples of two knots with the same Alexander module but differenttwisted Alexander polynomials are given.
Qiu, Yiheng; Henderson, Thomas M.; Scuseria, Gustavo E.
2016-09-01
Spin-projected Hartree-Fock is written as a particle-hole excitation ansatz over a symmetry-adapted reference determinant. Remarkably, this expansion has an analytic expression that we were able to decipher. While the form of the polynomial expansion is universal, the excitation amplitudes need to be optimized. This is equivalent to the optimization of orbitals in the conventional projected Hartree-Fock framework of non-orthogonal determinants. Using the inverse of the particle-hole expansion, we similarity transform the Hamiltonian in a coupled-cluster style theory. The left eigenvector of the non-Hermitian Hamiltonian is constructed in a similar particle-hole expansion fashion, and we show that to numerically reproduce variational projected Hartree-Fock results, one needs as many pair excitations in the bra as the number of strongly correlated entangled pairs in the system. This single-excitation polynomial similarity transformation theory is an alternative to our recently presented double excitation theory, but supports projected Hartree-Fock and coupled cluster simultaneously rather than interpolating between them.
Best polynomial degree reduction on q-lattices with applications to q-orthogonal polynomials
Ait-Haddou, Rachid
2015-06-07
We show that a weighted least squares approximation of q-Bézier coefficients provides the best polynomial degree reduction in the q-L
A 'missing' family of classical orthogonal polynomials
Energy Technology Data Exchange (ETDEWEB)
Vinet, Luc [Centre de Recherches Mathematiques, Universite de Montreal, PO Box 6128, Centre-ville Station, Montreal, Quebec, H3C 3J7 (Canada); Zhedanov, Alexei, E-mail: zhedanov@fti.dn.ua [Donetsk Institute for Physics and Technology, Donetsk 83114 (Ukraine)
2011-02-25
We study a family of 'classical' orthogonal polynomials which satisfy (apart from a three-term recurrence relation) an eigenvalue problem with a differential operator of Dunkl type. These polynomials can be obtained from the little q-Jacobi polynomials in the limit q = -1. We also show that these polynomials provide a nontrivial realization of the Askey-Wilson algebra for q = -1.
On the Lorentz degree of a product of polynomials
Ait-Haddou, Rachid
2015-01-01
In this note, we negatively answer two questions of T. Erdélyi (1991, 2010) on possible lower bounds on the Lorentz degree of product of two polynomials. We show that the correctness of one question for degree two polynomials is a direct consequence of a result of Barnard et al. (1991) on polynomials with nonnegative coefficients.
q-Extensions for the Apostol Type Polynomials
Directory of Open Access Journals (Sweden)
Nazim I. Mahmudov
2014-01-01
Full Text Available The aim of this work is to introduce an extension for q-standard notations. The q-Apostol type polynomials and study some of their properties. Besides, some relations between the mentioned polynomials and some other known polynomials are obtained.
Symmetry identities for 2-variable Apostol type and related polynomials
2015-01-01
In this article, certain symmetry identities for the 2-variable Apostol type polynomials are derived. By taking suitable values of parameters and indices, the symmetry identities for the special cases of the 2-variable Apostol type polynomials are established. Further, the symmetry identities for certain members belonging to the 2-variable Apostol type polynomials are also considered.
Approximating Exponential and Logarithmic Functions Using Polynomial Interpolation
Gordon, Sheldon P.; Yang, Yajun
2017-01-01
This article takes a closer look at the problem of approximating the exponential and logarithmic functions using polynomials. Either as an alternative to or a precursor to Taylor polynomial approximations at the precalculus level, interpolating polynomials are considered. A measure of error is given and the behaviour of the error function is…
The Roots of Adjoint Polynomial of the Graphs Contain Triangles
Institute of Scientific and Technical Information of China (English)
YECheng-fu
2004-01-01
We denote h(G,x) as the adjoint polynomial of graph G. In [5], Ma obtained the interpolation properties of the roots of adjoint polynomial of graphs containing triangles. By the properties, we prove the non-zero root of adjoint polynomial of Dn and Fn are single multiple.
Does the polynomial hierarchy collapse if onto functions are invertible?
H. Buhrman; L. Fortnow; M. Koucký; J.D. Rogers; N. Vereshchagin
2010-01-01
The class TFNP, defined by Megiddo and Papadimitriou, consists of multivalued functions with values that are polynomially verifiable and guaranteed to exist. Do we have evidence that such functions are hard, for example, if TFNP is computable in polynomial-time does this imply the polynomial-time hi
T. Kim; Choi, J.; Kim, Y. H.; C. S. Ryoo
2010-01-01
In this paper, we give a fermionic p-adic integral representions of Benstein polynomials associated with Euler numbers and polynomials. Finally, we give some interesting identities for the Euler numbers by using the properties of our integral represention.
Alvarez-Fernandez, Carlos; Manas, Manuel
2009-01-01
We consider the relation of the multi-component 2D Toda hierarchy with matrix orthogonal and biorthogonal polynomials. The multi-graded Hankel reduction of this hierarchy is considered and the corresponding generalized matrix orthogonal polynomials are studied. In particular for these polynomials we consider the recursion relations, and for rank one weights its relation with multiple orthogonal polynomials of mixed type with a type II normalization and the corresponding link with a Riemann--Hilbert problem.
Energy Technology Data Exchange (ETDEWEB)
Ceolin, C., E-mail: celina.ceolin@gmail.com [Universidade Federal de Santa Maria (UFSM), Frederico Westphalen, RS (Brazil). Centro de Educacao Superior Norte; Schramm, M.; Bodmann, B.E.J.; Vilhena, M.T., E-mail: celina.ceolin@gmail.com [Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, RS (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica
2015-07-01
Recently the stationary neutron diffusion equation in heterogeneous rectangular geometry was solved by the expansion of the scalar fluxes in polynomials in terms of the spatial variables (x; y), considering the two-group energy model. The focus of the present discussion consists in the study of an error analysis of the aforementioned solution. More specifically we show how the spatial subdomain segmentation is related to the degree of the polynomial and the Lipschitz constant. This relation allows to solve the 2-D neutron diffusion problem for second degree polynomials in each subdomain. This solution is exact at the knots where the Lipschitz cone is centered. Moreover, the solution has an analytical representation in each subdomain with supremum and infimum functions that shows the convergence of the solution. We illustrate the analysis with a selection of numerical case studies. (author)
HOMFLY polynomials in representation [3,1] for 3-strand braids
Mironov, A; Morozov, An; Sleptsov, A
2016-01-01
This paper is a new step in the project of systematic description of colored knot polynomials started in arXiv:1506.00339. In this paper, we managed to explicitly find the inclusive Racah matrix, i.e. the whole set of mixing matrices in channels R^3->Q with all possible Q, for R=[3,1]. The calculation is made possible by the use of a newly-developed efficient highest-weight method, still it remains tedious. The result allows one to evaluate and investigate [3,1]-colored polynomials for arbitrary 3-strand knots, and this confirms many previous conjectures on various factorizations, universality, and differential expansions. We consider in some detail the next-to-twist-knots three-strand family (n,-1|1,-1) and deduce its colored HOMFLY. Also confirmed and clarified is the eigenvalue hypothesis for the Racah matrices, which promises to provide a shortcut to generic formulas for arbitrary representations.
Bounding the Failure Probability Range of Polynomial Systems Subject to P-box Uncertainties
Crespo, Luis G.; Kenny, Sean P.; Giesy, Daniel P.
2012-01-01
This paper proposes a reliability analysis framework for systems subject to multiple design requirements that depend polynomially on the uncertainty. Uncertainty is prescribed by probability boxes, also known as p-boxes, whose distribution functions have free or fixed functional forms. An approach based on the Bernstein expansion of polynomials and optimization is proposed. In particular, we search for the elements of a multi-dimensional p-box that minimize (i.e., the best-case) and maximize (i.e., the worst-case) the probability of inner and outer bounding sets of the failure domain. This technique yields intervals that bound the range of failure probabilities. The offset between this bounding interval and the actual failure probability range can be made arbitrarily tight with additional computational effort.
Topological Expansion in the Complex Cubic Log-Gas Model: One-Cut Case
Bleher, Pavel; Deaño, Alfredo; Yattselev, Maxim
2016-09-01
We prove the topological expansion for the cubic log-gas partition function Z_N(t)= int _Γ \\cdots int _Γ prod _{1≤ jtopological expansion for log Z_N(t) in the one-cut phase region. The proof is based on the Riemann-Hilbert approach to semiclassical asymptotic expansions for the associated orthogonal polynomials and the theory of S-curves and quadratic differentials.
Modelling Trends in Ordered Correspondence Analysis Using Orthogonal Polynomials.
Lombardo, Rosaria; Beh, Eric J; Kroonenberg, Pieter M
2016-06-01
The core of the paper consists of the treatment of two special decompositions for correspondence analysis of two-way ordered contingency tables: the bivariate moment decomposition and the hybrid decomposition, both using orthogonal polynomials rather than the commonly used singular vectors. To this end, we will detail and explain the basic characteristics of a particular set of orthogonal polynomials, called Emerson polynomials. It is shown that such polynomials, when used as bases for the row and/or column spaces, can enhance the interpretations via linear, quadratic and higher-order moments of the ordered categories. To aid such interpretations, we propose a new type of graphical display-the polynomial biplot.
Discrete Fourier Analysis and Chebyshev Polynomials with G2 Group
Directory of Open Access Journals (Sweden)
Huiyuan Li
2012-10-01
Full Text Available The discrete Fourier analysis on the 30°-60°-90° triangle is deduced from the corresponding results on the regular hexagon by considering functions invariant under the group G2, which leads to the definition of four families generalized Chebyshev polynomials. The study of these polynomials leads to a Sturm-Liouville eigenvalue problem that contains two parameters, whose solutions are analogues of the Jacobi polynomials. Under a concept of m-degree and by introducing a new ordering among monomials, these polynomials are shown to share properties of the ordinary orthogonal polynomials. In particular, their common zeros generate cubature rules of Gauss type.
Polynomial threshold functions and Boolean threshold circuits
DEFF Research Database (Denmark)
Hansen, Kristoffer Arnsfelt; Podolskii, Vladimir V.
2013-01-01
of secondary interest. We show that PTFs on general Boolean domains are tightly connected to depth two threshold circuits. Our main results in regard to this connection are: PTFs of polynomial length and polynomial degree compute exactly the functions computed by THRMAJ circuits. An exponential length lower...... bound for PTFs that holds regardless of degree, thereby extending known lower bounds for THRMAJ circuits. We generalize two-party unbounded error communication complexity to the multi-party number-on-the-forehead setting, and show that communication lower bounds for 3-player protocols would yield size...... lower bounds for THRTHR circuits. We obtain several other results about PTFs. These include relationships between weight and degree of PTFs, and a degree lower bound for PTFs of constant length. We also consider a variant of PTFs over the max-plus algebra. We show that they are connected to PTFs over...
The Medusa Algorithm for Polynomial Matings
DEFF Research Database (Denmark)
Boyd, Suzanne Hruska; Henriksen, Christian
2012-01-01
The Medusa algorithm takes as input two postcritically finite quadratic polynomials and outputs the quadratic rational map which is the mating of the two polynomials (if it exists). Specifically, the output is a sequence of approximations for the parameters of the rational map, as well as an image...... of its Julia set. Whether these approximations converge is answered using Thurston's topological characterization of rational maps. This algorithm was designed by John Hamal Hubbard, and implemented in 1998 by Christian Henriksen and REU students David Farris and Kuon Ju Liu. In this paper we describe...... the algorithm and its implementation, discuss some output from the program (including many pictures) and related questions. Specifically, we include images and a discussion for some shared matings, Lattès examples, and tuning sequences of matings....
Characteristic polynomials in real Ginibre ensembles
Energy Technology Data Exchange (ETDEWEB)
Akemann, G; Phillips, M J [Department of Mathematical Sciences and BURSt Research Centre, Brunel University West London, UB8 3PH Uxbridge (United Kingdom); Sommers, H-J [Fachbereich Physik, Universitaet Duisburg-Essen, 47048 Duisburg (Germany)], E-mail: Gernot.Akemann@brunel.ac.uk, E-mail: Michael.Phillips@brunel.ac.uk, E-mail: H.J.Sommers@uni-due.de
2009-01-09
We calculate the average of two characteristic polynomials for the real Ginibre ensemble of asymmetric random matrices, and its chiral counterpart. Considered as quadratic forms they determine a skew-symmetric kernel from which all complex eigenvalue correlations can be derived. Our results are obtained in a very simple fashion without going to an eigenvalue representation, and are completely new in the chiral case. They hold for Gaussian ensembles which are partly symmetric, with kernels given in terms of Hermite and Laguerre polynomials respectively, depending on an asymmetry parameter. This allows us to interpolate between the maximally asymmetric real Ginibre and the Gaussian orthogonal ensemble, as well as their chiral counterparts. (fast track communication)
On Polynomial Sized MDP Succinct Policies
Liberatore, P
2011-01-01
Policies of Markov Decision Processes (MDPs) determine the next action to execute from the current state and, possibly, the history (the past states). When the number of states is large, succinct representations are often used to compactly represent both the MDPs and the policies in a reduced amount of space. In this paper, some problems related to the size of succinctly represented policies are analyzed. Namely, it is shown that some MDPs have policies that can only be represented in space super-polynomial in the size of the MDP, unless the polynomial hierarchy collapses. This fact motivates the study of the problem of deciding whether a given MDP has a policy of a given size and reward. Since some algorithms for MDPs work by finding a succinct representation of the value function, the problem of deciding the existence of a succinct representation of a value function of a given size and reward is also considered.
Eigenvalue conjecture and colored Alexander polynomials
Mironov, A
2016-01-01
We connect two important conjectures in the theory of knot polynomials. The first one is the property Al_R(q) = Al_{[1]}(q^{|R|}) for all single hook Young diagrams R, which is known to hold for all knots. The second conjecture claims that all the mixing matrices U_{i} in the relation {\\cal R}_i = U_i{\\cal R}_1U_i^{-1} between the i-th and the first generators {\\cal R}_i of the braid group are universally expressible through the eigenvalues of {\\cal R}_1. Since the above property of Alexander polynomials is very well tested, this relation provides a new support to the eigenvalue conjecture, especially for i>2, when its direct check by evaluation of the Racah matrices and their convolutions is technically difficult.
Tabulating knot polynomials for arborescent knots
Mironov, A; Morozov, An; Sleptsov, A; Ramadevi, P; Singh, Vivek Kumar
2016-01-01
Arborescent knots are the ones which can be represented in terms of double fat graphs or equivalently as tree Feynman diagrams. This is the class of knots for which the present knowledge is enough for lifting topological description to the level of effective analytical formulas. The paper describes the origin and structure of the new tables of colored knot polynomials, which will be posted at the dedicated site. Even if formal expressions are known in terms of modular transformation matrices, the computation in finite time requires additional ideas. We use the "family" approach, and apply it to arborescent knots in Rolfsen table by developing a Feynman diagram technique, associated with an auxiliary matrix model field theory. Gauge invariance in this theory helps to provide meaning to Racah matrices in the case of non-trivial multiplicities and explains the need for peculiar sign prescriptions in the calculation of [21]-colored HOMFLY polynomials.
Circular β ensembles, CMV representation, characteristic polynomials
Institute of Scientific and Technical Information of China (English)
SU ZhongGen
2009-01-01
In this note we first briefly review some recent progress in the study of the circular β ensemble on the unit circle, where 0 > 0 is a model parameter. In the special cases β = 1,2 and 4, this ensemble describes the joint probability density of eigenvalues of random orthogonal, unitary and sympletic matrices, respectively. For general β, Killip and Nenciu discovered a five-diagonal sparse matrix model, the CMV representation. This representation is new even in the case β = 2; and it has become a powerful tool for studying the circular β ensemble. We then give an elegant derivation for the moment identities of characteristic polynomials via the link with orthogonal polynomials on the unit circle.
Sparse DOA estimation with polynomial rooting
Xenaki, Angeliki; Gerstoft, Peter; Fernandez Grande, Efren
2015-01-01
Direction-of-arrival (DOA) estimation involves the localization of a few sources from a limited number of observations on an array of sensors. Thus, DOA estimation can be formulated as a sparse signal reconstruction problem and solved efficiently with compressive sensing (CS) to achieve highresolution imaging. Utilizing the dual optimal variables of the CS optimization problem, it is shown with Monte Carlo simulations that the DOAs are accurately reconstructed through polynomial rooting (Root...
Real meromorphic functions and linear differential polynomials
Institute of Scientific and Technical Information of China (English)
LANGLEY; J.; K.
2010-01-01
We determine all real meromorphic functions f in the plane such that f has finitely many zeros, the poles of f have bounded multiplicities, and f and F have finitely many non-real zeros, where F is a linear differential polynomial given by F = f (k) +Σk-1j=0ajf(j) , in which k≥2 and the coefficients aj are real numbers with a0≠0.
A Deterministic and Polynomial Modified Perceptron Algorithm
Directory of Open Access Journals (Sweden)
Olof Barr
2006-01-01
Full Text Available We construct a modified perceptron algorithm that is deterministic, polynomial and also as fast as previous known algorithms. The algorithm runs in time O(mn3lognlog(1/ρ, where m is the number of examples, n the number of dimensions and ρ is approximately the size of the margin. We also construct a non-deterministic modified perceptron algorithm running in timeO(mn2lognlog(1/ρ.
Time-reversal symmetry and random polynomials
Braun, D; Kus, M.; Zyczkowski, K.
1996-01-01
We analyze the density of roots of random polynomials where each complex coefficient is constructed of a random modulus and a fixed, deterministic phase. The density of roots is shown to possess a singular component only in the case for which the phases increase linearly with the index of coefficients. This means that, contrary to earlier belief, eigenvectors of a typical quantum chaotic system with some antiunitary symmetry will not display a clustering curve in the stellar representation. M...
Georeferencing CAMS data: Polynomial rectification and beyond
Yang, Xinghe
The Calibrated Airborne Multispectral Scanner (CAMS) is a sensor used in the commercial remote sensing program at NASA Stennis Space Center. In geographic applications of the CAMS data, accurate geometric rectification is essential for the analysis of the remotely sensed data and for the integration of the data into Geographic Information Systems (GIS). The commonly used rectification techniques such as the polynomial transformation and ortho rectification have been very successful in the field of remote sensing and GIS for most remote sensing data such as Landsat imagery, SPOT imagery and aerial photos. However, due to the geometric nature of the airborne line scanner which has high spatial frequency distortions, the polynomial model and the ortho rectification technique in current commercial software packages such as Erdas Imagine are not adequate for obtaining sufficient geometric accuracy. In this research, the geometric nature, especially the major distortions, of the CAMS data has been described. An analytical step-by-step geometric preprocessing has been utilized to deal with the potential high frequency distortions of the CAMS data. A generic sensor-independent photogrammetric model has been developed for the ortho-rectification of the CAMS data. Three generalized kernel classes and directional elliptical basis have been formulated into a rectification model of summation of multisurface functions, which is a significant extension to the traditional radial basis functions. The preprocessing mechanism has been fully incorporated into the polynomial, the triangle-based finite element analysis as well as the summation of multisurface functions. While the multisurface functions and the finite element analysis have the characteristics of localization, piecewise logic has been applied to the polynomial and photogrammetric methods, which can produce significant accuracy improvement over the global approach. A software module has been implemented with full
Normality and shared values concerning differential polynomials
Institute of Scientific and Technical Information of China (English)
无
2010-01-01
Let F be a family of functions meromorphic in a domain D, let P be a polynomial with either deg P≥3 or deg P = 2 and P having only one distinct zero, and let b be a finite nonzero complex number. If, each pair of functions f and g in F, P (f)f and P (g)g share b in D, then F is normal in D.
Completeness of the ring of polynomials
DEFF Research Database (Denmark)
Thorup, Anders
2015-01-01
Consider the polynomial ring R:=k[X1,…,Xn]R:=k[X1,…,Xn] in n≥2n≥2 variables over an uncountable field k. We prove that R is complete in its adic topology, that is, the translation invariant topology in which the non-zero ideals form a fundamental system of neighborhoods of 0. In addition we pro...
Piecewise polynomial solutions to linear inverse problems
DEFF Research Database (Denmark)
Hansen, Per Christian; Mosegaard, K.
1996-01-01
We have presented a new algorithm PP-TSVD that computes piecewise polynomial solutions to ill-posed problems, without a priori knowledge about the positions of the break points. In particular, we can compute piecewise constant functions that describe layered models. Such solutions are useful, e.g.......g., in seismological problems, and the algorithm can also be used as a preprocessor for other methods where break points/discontinuities must be incorporated explicitly....
On form factors and Macdonald polynomials
Lashkevich, Michael
2013-01-01
We are developing the algebraic construction for form factors of local operators in the sinh-Gordon theory proposed in [B.Feigin, M.Lashkeivch, 2008]. We show that the operators corresponding to the null vectors in this construction are given by the degenerate Macdonald polynomials with rectangular partitions and the parameters $t=-q$ on the unit circle. We obtain an integral representation for the null vectors and discuss its simple applications.
Moments, positive polynomials and their applications
Lasserre, Jean Bernard
2009-01-01
Many important applications in global optimization, algebra, probability and statistics, applied mathematics, control theory, financial mathematics, inverse problems, etc. can be modeled as a particular instance of the Generalized Moment Problem (GMP) . This book introduces a new general methodology to solve the GMP when its data are polynomials and basic semi-algebraic sets. This methodology combines semidefinite programming with recent results from real algebraic geometry to provide a hierarchy of semidefinite relaxations converging to the desired optimal value. Applied on appropriate cones,
Linear algebra for skew-polynomial matrices
Abramov, Sergei; Bronstein, Manuel
2002-01-01
We describe an algorithm for transforming skew-polynomial matrices over an Ore domain in row-reduced form, and show that this algorithm can be used to perform the standard calculations of linear algebra on such matrices (ranks, kernels, linear dependences, inhomogeneous solving). The main application of our algorithm is to desingularize recurrences and to compute the rational solutions of a large class of linear functional systems. It also turns out to be efficient when applied to ordinary co...
Polynomial approximation of functions in Sobolev spaces
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.
Some Orthogonal Polynomials in Four Variables
Directory of Open Access Journals (Sweden)
Charles F. Dunkl
2008-11-01
Full Text Available The symmetric group on 4 letters has the reflection group $D_3$ as an isomorphic image. This fact follows from the coincidence of the root systems $A_3$ and $D_3$. The isomorphism is used to construct an orthogonal basis of polynomials of 4 variables with 2 parameters. There is an associated quantum Calogero-Sutherland model of 4 identical particles on the line.
Algebraic polynomials and moments of stochastic integrals
Langovoy, Mikhail A
2011-01-01
We propose an algebraic method for proving estimates on moments of stochastic integrals. The method uses qualitative properties of roots of algebraic polynomials from certain general classes. As an application, we give a new proof of a variation of the Burkholder-Davis-Gundy inequality for the case of stochastic integrals with respect to real locally square integrable martingales. Further possible applications and extensions of the method are outlined.
A dynamic coefficient polynomial predistorter based on direct learning architecture
Institute of Scientific and Technical Information of China (English)
Li Bo; Ge Jianhua; Ai Bo
2008-01-01
A dynamic coefficient polynomial predistorter based on direct learning architecture is proposed. Compared to the existing polynomial predistorter, on the one hand, the proposed predistorter based on the direct learning architecture is more robust to initial conditions of the tap coefficients than that based on indirect learning architecture; on the other hand, by using two polynomial coefficient combinations, different polynomial coefficient combination can be selected when the input signal amplitude changes, which effectively decreases the estimate error. This paper introduces the direct learning architecture and gives the dynamic coefficient polynomial expression. A simplified nonlinear recursive least-squares (RLS) algorithm for polynomial coefficient estimation is also derived in detail. Computer simulations show that the proposed predistorter can attain 31dB, 28dB and 40dB spectrum suppression gain when our method is applied to the traveling wave tube amplifier (TWTA), solid state power amplifier (SSPA) and polynomial power amplifier (PA) model, respectively.
Characteristic polynomials of pseudo-Anosov maps
Birman, Joan; Kawamuro, Keiko
2010-01-01
We study the relationship between three different approaches to the action of a pseudo-Anosov mapping class $[F]$ on a surface: the original theorem of Thurston, its algorithmic proof by Bestvina-Handel, and related investigations of Penner-Harer. Bestvina and Handel represent $[F]$ as a suitably chosen homotopy equivalence $f: G\\to G$ of a finite graph, with an associated transition matrix $T$ whose largest eigenvalue is the dilatation of $[F]$. Extending a skew-symmetric form introduced by Penner and Harer to the setting of Bestvina and Handel, we show that the characteristic polynomial of $T$ is a monic and palindromic or anti-palindromic polynomial, possibly multiplied by a power of $x$. Moreover, it factors as a product of three polynomials. One of them reflects the action of $[F]$ on a certain symplectic space, the second one relates to the degeneracies of the skew-symmetric form, and the third one reflects the restriction of $f$ to the vertices of $G$. We give an application to the problem of deciding ...
Simplified Storm Surge Simulations Using Bernstein Polynomials
Beisiegel, Nicole; Behrens, Jörn
2016-04-01
Storm surge simulations are vital for forecasting, hazard assessment and eventually improving our understanding of Earth system processes. Discontinuous Galerkin (DG) methods have recently been explored in that context, because they are locally mass-conservative and in combination with suitable robust nodal filtering techniques (slope limiters) positivity-preserving and well-balanced for the still water state at rest. These filters manipulate interpolation point values in every time step in order to retain the desirable properties of the scheme. In particular, DG methods are able to represent prognostic variables such as the fluid height at high-order accuracy inside each element (triangle). For simulations that include wetting and drying, however, the high-order accuracy will destabilize the numerical model because point values on quadrature points may become negative during the computation if they do not coincide with interpolation points. This is why the model that we are presenting utilizes Bernstein polynomials as basis functions to model the wetting and drying. This has the advantage that negative pointvalues away from interpolation points are prevented, the model is stabilized and no additional time step restriction is introduced. Numerical tests show that the model is capable of simulating simplified storm surges. Furthermore, a comparison of model results with third-order Bernstein polynomials with results using traditional nodal Lagrange polynomials reveals an improvement in numerical convergence.
Deterministic Polynomial Factoring and Association Schemes
Arora, Manuel; Karpinski, Marek; Saxena, Nitin
2012-01-01
The problem of finding a nontrivial factor of a polynomial f(x) over a finite field F_q has many known efficient, but randomized, algorithms. The deterministic complexity of this problem is a famous open question even assuming the generalized Riemann hypothesis (GRH). In this work we improve the state of the art by focusing on prime degree polynomials; let n be the degree. If (n-1) has a `large' r-smooth divisor s, then we find a nontrivial factor of f(x) in deterministic poly(n^r,log q) time; assuming GRH and that s > sqrt{n/(2^r)}. Thus, for r = O(1) our algorithm is polynomial time. Further, for r > loglog n there are infinitely many prime degrees n for which our algorithm is applicable and better than the best known; assuming GRH. Our methods build on the algebraic-combinatorial framework of m-schemes initiated by Ivanyos, Karpinski and Saxena (ISSAC 2009). We show that the m-scheme on n points, implicitly appearing in our factoring algorithm, has an exceptional structure; leading us to the improved time ...
The bivariate Rogers Szegö polynomials
Chen, William Y. C.; Saad, Husam L.; Sun, Lisa H.
2007-06-01
We present an operator approach to deriving Mehler's formula and the Rogers formula for the bivariate Rogers-Szegö polynomials hn(x, y|q). The proof of Mehler's formula can be considered as a new approach to the nonsymmetric Poisson kernel formula for the continuous big q-Hermite polynomials Hn(x; a|q) due to Askey, Rahman and Suslov. Mehler's formula for hn(x, y|q) involves a 3phi2 sum and the Rogers formula involves a 2phi1 sum. The proofs of these results are based on parameter augmentation with respect to the q-exponential operator and the homogeneous q-shift operator in two variables. By extending recent results on the Rogers-Szegö polynomials hn(x|q) due to Hou, Lascoux and Mu, we obtain another Rogers-type formula for hn(x, y|q). Finally, we give a change of base formula for Hn(x; a|q) which can be used to evaluate some integrals by using the Askey-Wilson integral.
The bivariate Rogers-Szegoe polynomials
Energy Technology Data Exchange (ETDEWEB)
Chen, William Y C [Center for Combinatorics, LPMC, Nankai University, Tianjin 300071 (China); Saad, Husam L [Center for Combinatorics, LPMC, Nankai University, Tianjin 300071 (China); Sun, Lisa H [Center for Combinatorics, LPMC, Nankai University, Tianjin 300071 (China)
2007-06-08
We present an operator approach to deriving Mehler's formula and the Rogers formula for the bivariate Rogers-Szegoe polynomials h{sub n}(x, y vertical bar q). The proof of Mehler's formula can be considered as a new approach to the nonsymmetric Poisson kernel formula for the continuous big q-Hermite polynomials H{sub n}(x; a vertical bar q) due to Askey, Rahman and Suslov. Mehler's formula for h{sub n}(x, y vertical bar q) involves a {sub 3}{phi}{sub 2} sum and the Rogers formula involves a {sub 2}{phi}{sub 1} sum. The proofs of these results are based on parameter augmentation with respect to the q-exponential operator and the homogeneous q-shift operator in two variables. By extending recent results on the Rogers-Szegoe polynomials h{sub n}(x vertical bar q) due to Hou, Lascoux and Mu, we obtain another Rogers-type formula for h{sub n}(x, y vertical bar q). Finally, we give a change of base formula for H{sub n}(x; a vertical bar q) which can be used to evaluate some integrals by using the Askey-Wilson integral.
Nested Canalyzing, Unate Cascade, and Polynomial Functions.
Jarrah, Abdul Salam; Raposa, Blessilda; Laubenbacher, Reinhard
2007-09-15
This paper focuses on the study of certain classes of Boolean functions that have appeared in several different contexts. Nested canalyzing functions have been studied recently in the context of Boolean network models of gene regulatory networks. In the same context, polynomial functions over finite fields have been used to develop network inference methods for gene regulatory networks. Finally, unate cascade functions have been studied in the design of logic circuits and binary decision diagrams. This paper shows that the class of nested canalyzing functions is equal to that of unate cascade functions. Furthermore, it provides a description of nested canalyzing functions as a certain type of Boolean polynomial function. Using the polynomial framework one can show that the class of nested canalyzing functions, or, equivalently, the class of unate cascade functions, forms an algebraic variety which makes their analysis amenable to the use of techniques from algebraic geometry and computational algebra. As a corollary of the functional equivalence derived here, a formula in the literature for the number of unate cascade functions provides such a formula for the number of nested canalyzing functions.
Directory of Open Access Journals (Sweden)
Michael Basin
2011-04-01
Full Text Available In this paper, the mean-square filtering problem for polynomial system states confused with white Poisson noises over polynomial observations is studied proceeding from the general expression for the stochastic Ito differentials of the mean-square estimate and the error variance. In contrast to the previously obtained results, the paper deals with the general case of nonlinear polynomial states and observations with white Poisson noises. As a result, the Ito differentials for the mean-square estimate and error variance corresponding to the stated filtering problem are first derived. The procedure for obtaining an approximate closed-form finite-dimensional system of the filtering equations for any polynomial state over observations with any polynomial drift is then established. In the example, the obtained closed-form filter is applied to solve the third order sensor filtering problem for a quadratic state, assuming a conditionally Poisson initial condition for the extended third order state vector. The simulation results show that the designed filter yields a reliable and rapidly converging estimate.
Huang, Li
2016-11-01
Inspired by the recently proposed Legendre orthogonal polynomial representation for imaginary-time Green’s functions G(τ), we develop an alternate and superior representation for G(τ) and implement it in the hybridization expansion continuous-time quantum Monte Carlo impurity solver. This representation is based on the kernel polynomial method, which introduces some integral kernel functions to filter the numerical fluctuations caused by the explicit truncations of polynomial expansion series and can improve the computational precision significantly. As an illustration of the new representation, we re-examine the imaginary-time Green’s functions of the single-band Hubbard model in the framework of dynamical mean-field theory. The calculated results suggest that with carefully chosen integral kernel functions, whether the system is metallic or insulating, the Gibbs oscillations found in the previous Legendre orthogonal polynomial representation have been vastly suppressed and remarkable corrections to the measured Green’s functions have been obtained. Project supported by the National Natural Science Foundation of China (Grant No. 11504340).
Efficient modeling of photonic crystals with local Hermite polynomials
Boucher, C. R.; Li, Zehao; Albrecht, J. D.; Ram-Mohan, L. R.
2014-04-01
Developing compact algorithms for accurate electrodynamic calculations with minimal computational cost is an active area of research given the increasing complexity in the design of electromagnetic composite structures such as photonic crystals, metamaterials, optical interconnects, and on-chip routing. We show that electric and magnetic (EM) fields can be calculated using scalar Hermite interpolation polynomials as the numerical basis functions without having to invoke edge-based vector finite elements to suppress spurious solutions or to satisfy boundary conditions. This approach offers several fundamental advantages as evidenced through band structure solutions for periodic systems and through waveguide analysis. Compared with reciprocal space (plane wave expansion) methods for periodic systems, advantages are shown in computational costs, the ability to capture spatial complexity in the dielectric distributions, the demonstration of numerical convergence with scaling, and variational eigenfunctions free of numerical artifacts that arise from mixed-order real space basis sets or the inherent aberrations from transforming reciprocal space solutions of finite expansions. The photonic band structure of a simple crystal is used as a benchmark comparison and the ability to capture the effects of spatially complex dielectric distributions is treated using a complex pattern with highly irregular features that would stress spatial transform limits. This general method is applicable to a broad class of physical systems, e.g., to semiconducting lasers which require simultaneous modeling of transitions in quantum wells or dots together with EM cavity calculations, to modeling plasmonic structures in the presence of EM field emissions, and to on-chip propagation within monolithic integrated circuits.
Directory of Open Access Journals (Sweden)
Mario Faliva
2007-10-01
Full Text Available In this paper the issue of the inversion of a matrix polynomial about a unit root is tackled by resorting to Laurent expansion: The principal-part matrix coefficients associated with a simple and a second order pole are properly characterized and closed-form expressions are derived by virtue of a recent result on partitioned inversion (Faliva and Zoia, 2002. This eventually sheds more light on the analytical foundations of unit-root econometrics which in turn paves the way to an elegant unified representation theorem for (cointegrated processes up to the second order.
Energy Technology Data Exchange (ETDEWEB)
FEDOROVA,A.; ZEITLIN,M.; PARSA,Z.
2000-03-31
In this paper the authors present applications of methods from wavelet analysis to polynomial approximations for a number of accelerator physics problems. According to a variational approach in the general case they have the solution as a multiresolution (multiscales) expansion on the base of compactly supported wavelet basis. They give an extension of their results to the cases of periodic orbital particle motion and arbitrary variable coefficients. Then they consider more flexible variational method which is based on a biorthogonal wavelet approach. Also they consider a different variational approach, which is applied to each scale.
Time-reversal symmetry and random polynomials
Braun, D; Zyczkowski, K
1996-01-01
We analyze the density of roots of random polynomials where each complex coefficient is constructed of a random modulus and a fixed, deterministic phase. The density of roots is shown to possess a singular component only in the case for which the phases increase linearly with the index of coefficients. This means that, contrary to earlier belief, eigenvectors of a typical quantum chaotic system with some antiunitary symmetry will not display a clustering curve in the stellar representation. Moreover, a class of time-reverse invariant quantum systems is shown, for which spectra display fluctuations characteristic of orthogonal ensemble, while eigenvectors confer to predictions of unitary ensemble.
Conditional Density Approximations with Mixtures of Polynomials
DEFF Research Database (Denmark)
Varando, Gherardo; López-Cruz, Pedro L.; Nielsen, Thomas Dyhre
2015-01-01
Mixtures of polynomials (MoPs) are a non-parametric density estimation technique especially designed for hybrid Bayesian networks with continuous and discrete variables. Algorithms to learn one- and multi-dimensional (marginal) MoPs from data have recently been proposed. In this paper we introduce...... is found. We illustrate and study the methods using data sampled from known parametric distributions, and we demonstrate their applicability by learning models based on real neuroscience data. Finally, we compare the performance of the proposed methods with an approach for learning mixtures of truncated...
A Matricial Algorithm for Polynomial Refinement
King, Emily J
2011-01-01
In order to have a multiresolution analysis, the scaling function must be refinable. That is, it must be the linear combination of 2-dilation, $\\mathbb{Z}$-translates of itself. Refinable functions used in connection with wavelets are typically compactly supported. In 2002, David Larson posed the question, "Are all polynomials (of a single variable) finitely refinable?" That summer the author proved that the answer indeed was true using basic linear algebra. The result was presented in a number of talks but had not been typed up until now. The purpose of this short note is to record that particular proof.
Polynomial identities for ternary intermolecular recombination
Bremner, Murray R
2010-01-01
The operation of binary intermolecular recombination, originating in the theory of DNA computing, permits a natural generalization to n-ary operations which perform simultaneous recombination of n molecules. In the case n = 3, we use computer algebra to determine the polynomial identities of degree <= 9 satisfied by this trilinear nonassociative operation. Our approach requires computing a basis for the nullspace of a large integer matrix, and for this we compare two methods: (i) the row canonical form, and (ii) the Hermite normal form with lattice basis reduction. In the conclusion, we formulate some conjectures for the general case of n-ary intermolecular recombination.
Pade interpolation by F-polynomials and transfinite diameter
Coman, Dan
2011-01-01
We define $F$-polynomials as linear combinations of dilations by some frequencies of an entire function $F$. In this paper we use Pade interpolation of holomorphic functions in the unit disk by $F$-polynomials to obtain explicitly approximating $F$-polynomials with sharp estimates on their coefficients. We show that when frequencies lie in a compact set $K\\subset\\mathbb C$ then optimal choices for the frequencies of interpolating polynomials are similar to Fekete points. Moreover, the minimal norms of the interpolating operators form a sequence whose rate of growth is determined by the transfinite diameter of $K$. In case of the Laplace transforms of measures on $K$, we show that the coefficients of interpolating polynomials stay bounded provided that the frequencies are Fekete points. Finally, we give a sufficient condition for measures on the unit circle which ensures that the sums of the absolute values of the coefficients of interpolating polynomials stay bounded.
Guts of surfaces and the colored Jones polynomial
Futer, David; Purcell, Jessica S
2011-01-01
This work derives direct and concrete relations between colored Jones polynomials and the topology of incompressible spanning surfaces in knot and link complements. Under mild diagrammatic hypotheses that arise naturally in the study of knot polynomial invariants (A-adequacy), we prove that the growth of the degree of the colored Jones polynomials is a boundary slope of an essential surface in the knot complement. We show that certain coefficients of the polynomial measure how far this surface is from being a fiber in the knot complement; in particular, the surface is a fiber if and only if a particular coefficient vanishes. Our results also yield concrete relations between hyperbolic geometry and colored Jones polynomials: for certain families of links, coefficients of the polynomials determine the hyperbolic volume to within a factor of 4. Our approach is to generalize the checkerboard decompositions of alternating knots. Under mild diagrammatic hypotheses (A-adequacy), we show that the checkerboard knot su...
Computing the Tutte Polynomial in Vertex-Exponential Time
DEFF Research Database (Denmark)
Björklund, Andreas; Husfeldt, Thore; Kaski, Petteri;
2008-01-01
The deletion–contraction algorithm is perhapsthe most popular method for computing a host of fundamental graph invariants such as the chromatic, flow, and reliability polynomials in graph theory, the Jones polynomial of an alternating link in knot theory, and the partition functions of the models...... algorithm that computes the Tutte polynomial—and hence, all the aforementioned invariants and more—of an arbitrary graph in time within a polynomial factor of the number of connected vertex sets. The algorithm actually evaluates a multivariate generalization of the Tutte polynomial by making use...... of an identity due to Fortuin and Kasteleyn. We also provide a polynomial-space variant of the algorithm and give an analogous result for Chung and Graham's cover polynomial....
On the Tutte-Krushkal-Renardy polynomial for cell complexes
Bajo, Carlos; Chmutov, Sergei
2012-01-01
Recently V.Krushkal and D.Renardy generalized the Tutte polynomial from graphs to cell complexes. We show that evaluating this polynomial at the origin gives the number of cellular spanning trees in the sense of A.Duval, C.Klivans, and J.Martin. Moreover, after a slight modification, the Tutte-Krushkal-Renardy polynomial evaluated at the origin gives a weighted count of cellular spanning trees, and therefore its free term can be calculated by the cellular matrix-tree theorem of Duval et al. In the case of cell decomposition of a sphere, this modified polynomial satisfies the same duality identity as before. We find that evaluating the Tutte-Krushkal-Renardy along a certain line is the Bott polynomial. Finally we prove skein relations for the Tutte-Krushkal-Renardy polynomial.
Fock expansion of multimode pure Gaussian states
Energy Technology Data Exchange (ETDEWEB)
Cariolaro, Gianfranco; Pierobon, Gianfranco, E-mail: gianfranco.pierobon@unipd.it [Università di Padova, Padova (Italy)
2015-12-15
The Fock expansion of multimode pure Gaussian states is derived starting from their representation as displaced and squeezed multimode vacuum states. The approach is new and appears to be simpler and more general than previous ones starting from the phase-space representation given by the characteristic or Wigner function. Fock expansion is performed in terms of easily evaluable two-variable Hermite–Kampé de Fériet polynomials. A relatively simple and compact expression for the joint statistical distribution of the photon numbers in the different modes is obtained. In particular, this result enables one to give a simple characterization of separable and entangled states, as shown for two-mode and three-mode Gaussian states.
Multiple Twisted -Euler Numbers and Polynomials Associated with -Adic -Integrals
Directory of Open Access Journals (Sweden)
Jang Lee-Chae
2008-01-01
Full Text Available By using -adic -integrals on , we define multiple twisted -Euler numbers and polynomials. We also find Witt's type formula for multiple twisted -Euler numbers and discuss some characterizations of multiple twisted -Euler Zeta functions. In particular, we construct multiple twisted Barnes' type -Euler polynomials and multiple twisted Barnes' type -Euler Zeta functions. Finally, we define multiple twisted Dirichlet's type -Euler numbers and polynomials, and give Witt's type formula for them.
Representations of the Schroedinger group and matrix orthogonal polynomials
Energy Technology Data Exchange (ETDEWEB)
Vinet, Luc [Centre de recherches mathematiques, Universite de Montreal, CP 6128, succ. Centre-ville, Montreal, QC H3C 3J7 (Canada); Zhedanov, Alexei, E-mail: luc.vinet@umontreal.ca, E-mail: zhedanov@fti.dn.ua [Donetsk Institute for Physics and Technology, Donetsk 83114 (Ukraine)
2011-09-02
The representations of the Schroedinger group in one space dimension are explicitly constructed in the basis of the harmonic oscillator states. These representations are seen to involve matrix orthogonal polynomials in a discrete variable that have Charlier and Meixner polynomials as building blocks. The underlying Lie-theoretic framework allows for a systematic derivation of the structural formulas (recurrence relations, difference equations, Rodrigues' formula, etc) that these matrix orthogonal polynomials satisfy. (paper)
THE λ-GR(O)BNER BASES UNDER POLYNOMIAL COMPOSITION
Institute of Scientific and Technical Information of China (English)
Jinwang LIU; Dongmei LI; Xiaosong CHEN
2007-01-01
Polynomial composition is the operation of replacing variables in a polynomial with other polynomials. λ-Gr(o)bner basis is an especial Gr(o)bner basis. The main problem in the paper is: when does composition commute with λ-Gr(o)bner basis computation? We shall answer better the above question. This has a natural application in the computation of λ-Gr(o)bner bases.
Chromatic polynomials of graphs from Kac-Moody algebras
Venkatesh, R.; Viswanath, Sankaran
2013-01-01
We give a new interpretation of the chromatic polynomial of a simple graph G in terms of the Kac-Moody Lie algebra with Dynkin diagram G. We show that the chromatic polynomial is essentially the q-Kostant partition function of this Lie algebra evaluated on the sum of the simple roots. Applying the Peterson recurrence formula for root multiplicities, we obtain a new realization of the chromatic polynomial as a weighted sum of paths in the bond lattice of G.
On an Inequality of Pual Turan Concerning Polynomials-II
Institute of Scientific and Technical Information of China (English)
Abdullah Mir∗
2015-01-01
Let P(z) be a polynomial of degree n and for any complex number α, let DαP(z)=nP(z)+(α−z)P′(z) denote the polar derivative of the polynomial P(z) with respect to α. In this paper, we obtain inequalities for the polar derivative of a poly-nomial having all zeros inside a circle. Our results shall generalize and sharpen some well-known results of Turan, Govil, Dewan et al. and others.
Darboux polynomials for Lotka-Volterra systems in three dimensions
Christodoulides, Yiannis T
2008-01-01
We consider Lotka-Volterra systems in three dimensions depending on three real parameters. By using elementary algebraic methods we classify the Darboux polynomials (also known as second integrals) for such systems for various values of the parameters, and give the explicit form of the corresponding cofactors. More precisely, we show that a Darboux polynomial of degree greater than one is reducible. In fact, it is a product of linear Darboux polynomials and first integrals.
Some Results on the Polar Derivative of a Polynomial
Institute of Scientific and Technical Information of China (English)
Abdullah Mir; Bilal Dar
2014-01-01
Let P(z) be a polynomial of degree n and for any complex number α, let DαP(z)=nP(z)+(α-z)P′(z) denote the polar derivative of P(z) with respect toα. In this paper, we obtain certain inequalities for the polar derivative of a polynomial with restricted zeros. Our results generalize and sharpen some well-known polynomial inequalities.
Univariate Niho Bent Functions from o-Polynomials
Budaghyan, Lilya; Kholosha, Alexander; Carlet, Claude; Helleseth, Tor
2014-01-01
In this paper, we discover that any univariate Niho bent function is a sum of functions having the form of Leander-Kholosha bent functions with extra coefficients of the power terms. This allows immediately, knowing the terms of an o-polynomial, to obtain the powers of the additive terms in the polynomial representing corresponding bent function. However, the coefficients are calculated ambiguously. The explicit form is given for the bent functions obtained from quadratic and cubic o-polynomi...
Local polynomial Whittle estimation covering non-stationary fractional processes
DEFF Research Database (Denmark)
Nielsen, Frank
This paper extends the local polynomial Whittle estimator of Andrews & Sun (2004) to fractionally integrated processes covering stationary and non-stationary regions. We utilize the notion of the extended discrete Fourier transform and periodogram to extend the local polynomial Whittle estimator ...... study illustrates the performance of the proposed estimator compared to the classical local Whittle estimator and the local polynomial Whittle estimator. The empirical justi.cation of the proposed estimator is shown through an analysis of credit spreads....
Multivariate polynomial interpolation on Lissajous-Chebyshev nodes
Dencker, Peter; Erb, Wolfgang
2015-01-01
In this article, we study multivariate polynomial interpolation and quadrature rules on non-tensor product node sets related to Lissajous curves and Chebyshev varieties. After classifying multivariate Lissajous curves and the interpolation nodes linked to these curves, we derive a discrete orthogonality structure on these node sets. Using this orthogonality structure, we obtain unique polynomial interpolation in appropriately defined spaces of multivariate Chebyshev polynomials. Our results g...
Conventional modeling of the multilayer perceptron using polynomial basis functions
Chen, Mu-Song; Manry, Michael T.
1993-01-01
A technique for modeling the multilayer perceptron (MLP) neural network, in which input and hidden units are represented by polynomial basis functions (PBFs), is presented. The MLP output is expressed as a linear combination of the PBFs and can therefore be expressed as a polynomial function of its inputs. Thus, the MLP is isomorphic to conventional polynomial discriminant classifiers or Volterra filters. The modeling technique was successfully applied to several trained MLP networks.
Return times of polynomials as meta-Fibonacci numbers
Emerson, Nathaniel D.
2004-01-01
We consider generalized closest return times of a complex polynomial of degree at least two. Most previous studies on this subject have focused on the properties of polynomials with particular return times, especially the Fibonacci numbers. We study the general form of these closest return times. The main result of this paper is that these closest return times are meta-Fibonacci numbers. This result applies to the return times of a principal nest of a polynomial. Furthermore, we show that an ...
Directory of Open Access Journals (Sweden)
Tsugio Fukuchi
2014-06-01
Full Text Available The finite difference method (FDM based on Cartesian coordinate systems can be applied to numerical analyses over any complex domain. A complex domain is usually taken to mean that the geometry of an immersed body in a fluid is complex; here, it means simply an analytical domain of arbitrary configuration. In such an approach, we do not need to treat the outer and inner boundaries differently in numerical calculations; both are treated in the same way. Using a method that adopts algebraic polynomial interpolations in the calculation around near-wall elements, all the calculations over irregular domains reduce to those over regular domains. Discretization of the space differential in the FDM is usually derived using the Taylor series expansion; however, if we use the polynomial interpolation systematically, exceptional advantages are gained in deriving high-order differences. In using the polynomial interpolations, we can numerically solve the Poisson equation freely over any complex domain. Only a particular type of partial differential equation, Poisson's equations, is treated; however, the arguments put forward have wider generality in numerical calculations using the FDM.
Polynomial analysis of canopy spectra and biochemical component content inversion
Institute of Scientific and Technical Information of China (English)
YAN Chunyan; LIU Qiang; NIU Zheng; WANG Jihua; HUANG Wenjiang; LIU Liangyun
2005-01-01
A polynomial expression model was developed in this paper to describe directional canopy spectra, and the decomposition of the polynomial expression was used as a tool for retrieving biochemical component content from canopy multi-angle spectra. First, the basic formula of the polynomial expression was introduced and the physical meaning of its terms and coefficients was discussed. Based on this analysis, a complete polynomial expression model and its decomposition method were given. By decomposing the canopy spectra simulated with SAILH model, it shows that the polynomial expression can not only fit well the canopy spectra, but also show the contribution of every order scattering to the whole reflectance. Taking the first scattering coefficients a10 and a01 for example, the test results show that the polynomial coefficients reflect very well the hot spot phenomenon and the effects of viewing angles, LAI and leaf inclination angle on canopy spectra. By coupling the polynomial expression with leaf model PROSPECT, a canopy biochemical component content inversion model was given. In the simulated test, the canopy multi-angle spectra were simulated by two different models, SAILH and 4-SCALE respectively, then the biochemical component content was retrieved by inverting the coupled polynomial expression + PROSPECT model. Results of the simulated test are promising, and when applying the algorithm to measured corn canopy multi-angle spectra, we also get relatively accurate chlorophyll content. It shows that the polynomial analysis provides a new method to get biochemical component content independent of any specific canopy model.
Pseudorandom Numbers and Hash Functions from Iterations of Multivariate Polynomials
Ostafe, Alina
2009-01-01
Dynamical systems generated by iterations of multivariate polynomials with slow degree growth have proved to admit good estimates of exponential sums along their orbits which in turn lead to rather stronger bounds on the discrepancy for pseudorandom vectors generated by these iterations. Here we add new arguments to our original approach and also extend some of our recent constructions and results to more general orbits of polynomial iterations which may involve distinct polynomials as well. Using this construction we design a new class of hash functions from iterations of polynomials and use our estimates to motivate their "mixing" properties.
Automatically Discovering Relaxed Lyapunov Functions for Polynomial Dynamical Systems
Liu, Jiang; Zhao, Hengjun
2011-01-01
The notion of Lyapunov function plays a key role in design and verification of dynamical systems, as well as hybrid and cyber-physical systems. In this paper, to analyze the asymptotic stability of a dynamical system, we generalize standard Lyapunov functions to relaxed Lyapunov functions (RLFs), by considering higher order Lie derivatives of certain functions along the system's vector field. Furthermore, we present a complete method to automatically discovering polynomial RLFs for polynomial dynamical systems (PDSs). Our method is complete in the sense that it is able to discover all polynomial RLFs by enumerating all polynomial templates for any PDS.
An Analytic Formula for the A_2 Jack Polynomials
Directory of Open Access Journals (Sweden)
Vladimir V. Mangazeev
2007-01-01
Full Text Available In this letter I shall review my joint results with Vadim Kuznetsov and Evgeny Sklyanin [Indag. Math. 14 (2003, 451-482] on separation of variables (SoV for the $A_n$ Jack polynomials. This approach originated from the work [RIMS Kokyuroku 919 (1995, 27-34] where the integral representations for the $A_2$ Jack polynomials was derived. Using special polynomial bases I shall obtain a more explicit expression for the $A_2$ Jack polynomials in terms of generalised hypergeometric functions.
Multi-indexed Jacobi polynomials and Maya diagrams
Takemura, Kouichi
2014-11-01
Multi-indexed Jacobi polynomials are defined by the Wronskian of four types of eigenfunctions of the Pöschl-Teller Hamiltonian. We give a correspondence between multi-indexed Jacobi polynomials and pairs of Maya diagrams, and we show that any multi-indexed Jacobi polynomial is essentially equal to some multi-indexed Jacobi polynomial of two types of eigenfunction. As an application, we show a Wronskian-type formula of some special eigenstates of the deformed Pöschl-Teller Hamiltonian.
Casoratian Identities for the Wilson and Askey-Wilson Polynomials
Odake, Satoru
2013-01-01
Infinitely many Casoratian identities are derived for the Wilson and Askey-Wilson polynomials in parallel to the Wronskian identities for the Hermite, Laguerre and Jacobi polynomials, which were reported recently by the present authors. These identities form the basis of the equivalence between eigenstate adding and deleting Darboux transformations for solvable (discrete) quantum mechanical systems. Similar identities hold for various reduced form polynomials of the Wilson and Askey-Wilson polynomials, e.g. the continuous q-Jacobi, continuous (dual) (q-)Hahn, Meixner-Pollaczek, Al-Salam-Chihara, continuous (big) q-Hermite, etc.
Structure relations for monic orthogonal polynomials in two discrete variables
Rodal, J.; Area, I.; Godoy, E.
2008-04-01
In this paper, extensions of several relations linking differences of bivariate discrete orthogonal polynomials and polynomials themselves are given, by using an appropriate vector-matrix notation. Three-term recurrence relations are presented for the partial differences of the monic polynomial solutions of admissible second order partial difference equation of hypergeometric type. Structure relations, difference representations as well as lowering and raising operators are obtained. Finally, expressions for all matrix coefficients appearing in these finite-type relations are explicitly presented for a finite set of Hahn and Kravchuk orthogonal polynomials.
The distribution of zeros of general q-polynomials
Álvarez-Nodarse, R.; Buendía, E.; Dehesa, J. S.
1997-10-01
A general system of q-orthogonal polynomials is defined by means of its three-term recurrence relation. This system encompasses many of the known families of q-polynomials, among them the q-analogue of the classical orthogonal polynomials. The asymptotic density of zeros of the system is shown to be a simple and compact expression of the parameters which characterize the asymptotic behaviour of the coefficients of the recurrence relation. This result is applied to specific classes of polynomials known by the names q-Hahn, q-Kravchuk, q-Racah, q-Askey and Wilson, Al Salam - Carlitz and the celebrated little and big q-Jacobi.
Unimodularity of zeros of self-inversive polynomials
Lalin, Matilde N
2012-01-01
We generalise a necessary and sufficient condition given by Cohn for all the zeros of a self-inversive polynomial to be on the unit circle. Our theorem implies some sufficient conditions found by Lakatos, Losonczi and Schinzel. We apply our result to the study of a polynomial family closely related to Ramanujan polynomials, recently introduced by Gun, Murty and Rath, and studied by Murty, Smyth and Wang and Lal\\'in and Rogers. We prove that all polynomials in this family have their zeros on the unit circle, a result conjectured by Lal\\'in and Rogers on computational evidence.
Positive trigonometric polynomials and signal processing applications
Dumitrescu, Bogdan
2017-01-01
This revised edition is made up of two parts: theory and applications. Though many of the fundamental results are still valid and used, new and revised material is woven throughout the text. As with the original book, the theory of sum-of-squares trigonometric polynomials is presented unitarily based on the concept of Gram matrix (extended to Gram pair or Gram set). The programming environment has also evolved, and the books examples are changed accordingly. The applications section is organized as a collection of related problems that use systematically the theoretical results. All the problems are brought to a semi-definite programming form, ready to be solved with algorithms freely available, like those from the libraries SeDuMi, CVX and Pos3Poly. A new chapter discusses applications in super-resolution theory, where Bounded Real Lemma for trigonometric polynomials is an important tool. This revision is written to be more appealing and easier to use for new readers. < Features updated information on LMI...
Orthogonal Polynomials and $S$-curves
Rakhmanov, E A
2011-01-01
This paper is devoted to a study of $S$-curves, that is systems of curves in the complex plane whose equilibrium potential in a harmonic external field satisfies a special symmetry property ($S$-property). Such curves have many applications. In particular, they play a fundamental role in the theory of complex (non-hermitian) orthogonal polynomials. One of the main theorems on zero distribution of such polynomials asserts that the limit zero distribution is presented by an equilibrium measure of an $S$-curve associated with the problem if such a curve exists. These curves are also the starting point of the matrix Riemann-Hilbert approach to srtong asymptotics. Other approaches to the problem of strong asymptotics (differential equations, Riemann surfaces) are also related to $S$-curves or may be interpreted this way. Existence problem $S$-curve in a given class of curves in presence of a nontrivial external field presents certain challenge. We formulate and prove a version of existence theorem for the case whe...
On Factorization of Generalized Macdonald Polynomials
Kononov, Ya
2016-01-01
A remarkable feature of Schur functions -- the common eigenfunctions of cut-and-join operators from $W_\\infty$ -- is that they factorize at the peculiar two-parametric topological locus in the space of time-variables, what is known as the hook formula for quantum dimensions of representations of $U_q(SL_N)$ and plays a big role in various applications. This factorization survives at the level of Macdonald polynomials. We look for its further generalization to generalized Macdonald polynomials (GMP), associated in the same way with the toroidal Ding-Iohara-Miki algebras, which play the central role in modern studies in Seiberg-Witten-Nekrasov theory. In the simplest case of the first-coproduct eigenfunctions, where GMP depend on just two sets of time-variables, we discover a weak factorization -- on a one- (rather than four-) parametric slice of the topological locus, what is already a very non-trivial property, calling for proof and better understanding.
On factorization of generalized Macdonald polynomials
Energy Technology Data Exchange (ETDEWEB)
Kononov, Ya. [Landau Institute for Theoretical Physics, Chernogolovka (Russian Federation); HSE, Math Department, Moscow (Russian Federation); Morozov, A. [ITEP, Moscow (Russian Federation); Institute for Information Transmission Problems, Moscow (Russian Federation); National Research Nuclear University MEPhI, Moscow (Russian Federation)
2016-08-15
A remarkable feature of Schur functions - the common eigenfunctions of cut-and-join operators from W{sub ∞} - is that they factorize at the peculiar two-parametric topological locus in the space of time variables, which is known as the hook formula for quantum dimensions of representations of U{sub q}(SL{sub N}) and which plays a big role in various applications. This factorization survives at the level of Macdonald polynomials. We look for its further generalization to generalized Macdonald polynomials (GMPs), associated in the same way with the toroidal Ding-Iohara-Miki algebras, which play the central role in modern studies in Seiberg-Witten-Nekrasov theory. In the simplest case of the first-coproduct eigenfunctions, where GMP depend on just two sets of time variables, we discover a weak factorization - on a one- (rather than four-) parametric slice of the topological locus, which is already a very non-trivial property, calling for proof and better understanding. (orig.)
Tabulating knot polynomials for arborescent knots
Mironov, A.; Morozov, A.; Morozov, A.; Ramadevi, P.; Singh, Vivek Kumar; Sleptsov, A.
2017-02-01
Arborescent knots are those which can be represented in terms of double fat graphs or equivalently as tree Feynman diagrams. This is the class of knots for which the present knowledge is sufficient for lifting topological description to the level of effective analytical formulas. The paper describes the origin and structure of the new tables of colored knot polynomials, which will be posted at the dedicated site (http://knotebook.org). Even if formal expressions are known in terms of modular transformation matrices, the computation in finite time requires additional ideas. We use the ‘family’ approach, suggested in Mironov and Morozov (2015 Nucl. Phys. B 899 395–413), and apply it to arborescent knots in the Rolfsen table by developing a Feynman diagram technique, associated with an auxiliary matrix model field theory. Gauge invariance in this theory helps to provide meaning to Racah matrices in the case of non-trivial multiplicities and explains the need for peculiar sign prescriptions in the calculation of [21]-colored HOMFLY-PT polynomials.
On factorization of generalized Macdonald polynomials
Kononov, Ya.; Morozov, A.
2016-08-01
A remarkable feature of Schur functions—the common eigenfunctions of cut-and-join operators from W_∞ —is that they factorize at the peculiar two-parametric topological locus in the space of time variables, which is known as the hook formula for quantum dimensions of representations of U_q(SL_N) and which plays a big role in various applications. This factorization survives at the level of Macdonald polynomials. We look for its further generalization to generalized Macdonald polynomials (GMPs), associated in the same way with the toroidal Ding-Iohara-Miki algebras, which play the central role in modern studies in Seiberg-Witten-Nekrasov theory. In the simplest case of the first-coproduct eigenfunctions, where GMP depend on just two sets of time variables, we discover a weak factorization—on a one- (rather than four-) parametric slice of the topological locus, which is already a very non-trivial property, calling for proof and better understanding.
Seizure prediction using polynomial SVM classification.
Zisheng Zhang; Parhi, Keshab K
2015-08-01
This paper presents a novel patient-specific algorithm for prediction of seizures in epileptic patients with low hardware complexity and low power consumption. In the proposed approach, we first compute the spectrogram of the input fragmented EEG signals from a few electrodes. Each fragmented data clip is ten minutes in duration. Band powers, relative spectral powers and ratios of spectral powers are extracted as features. The features are then subjected to electrode selection and feature selection using classification and regression tree. The baseline experiment uses all features from selected electrodes and these features are then subjected to a radial basis function kernel support vector machine (RBF-SVM) classifier. The proposed method further selects a small number features from the selected electrodes and train a polynomial support vector machine (SVM) classifier with degree of 2 on these features. Prediction performances are compared between the baseline experiment and the proposed method. The algorithm is tested using intra-cranial EEG (iEEG) from the American Epilepsy Society Seizure Prediction Challenge database. The baseline experiment using a large number of features and RBF-SVM achieves a 100% sensitivity and an average AUC of 0.9985, while the proposed algorithm using only a small number of features and polynomial SVM with degree of 2 can achieve a sensitivity of 100.0%, an average area under curve (AUC) of 0.9795. For both experiments, only 10% of the available training data are used for training.
Solving Operator Equation Based on Expansion Approach
Directory of Open Access Journals (Sweden)
A. Aminataei
2014-01-01
Full Text Available To date, researchers usually use spectral and pseudospectral methods for only numerical approximation of ordinary and partial differential equations and also based on polynomial basis. But the principal importance of this paper is to develop the expansion approach based on general basis functions (in particular case polynomial basis for solving general operator equations, wherein the particular cases of our development are integral equations, ordinary differential equations, difference equations, partial differential equations, and fractional differential equations. In other words, this paper presents the expansion approach for solving general operator equations in the form Lu+Nu=g(x,x∈Γ, with respect to boundary condition Bu=λ, where L, N and B are linear, nonlinear, and boundary operators, respectively, related to a suitable Hilbert space, Γ is the domain of approximation, λ is an arbitrary constant, and g(x∈L2(Γ is an arbitrary function. Also the other importance of this paper is to introduce the general version of pseudospectral method based on general interpolation problem. Finally some experiments show the accuracy of our development and the error analysis is presented in L2(Γ norm.
Heller, J; Schmidt, C; van Rienen, U
2014-01-01
The electromagnetic properties of SRF cavities are mostly determined by their shape. Due to fabrication tolerances, tuning and limited resolution of measurement systems, the exact shape remains uncertain. In order to make assessments for the real life behaviour it is important to quantify how these geometrical uncertainties propagate through the mathematical system and influence certain electromagnetic properties, like the resonant frequencies of the structure’s eigenmodes. This can be done by using non-intrusive straightforward methods like Monte Carlo (MC) simulations. However, such simulations require a large number of deterministic problem solutions to obtain a sufficient accuracy. In order to avoid this scaling behaviour, the so-called generalized polynomial chaos (gPC) expansion is used. This technique allows for the relatively fast computation of uncertainty propagation for few uncertain parameters in the case of computationally expensive deterministic models. In this paper we use the gPC expansion t...
A Parallel Algorithm for Finding Roots of a Complex Polynomial
Institute of Scientific and Technical Information of China (English)
程锦松
1990-01-01
A distribution theory of the roots of a polynomial and a parallel algorithm for finding roots of a complex polynomial based on that theory are developed in this paper.With high parallelism,the algorithm is an improvement over the Wilf algorithm[3].
Weighted Approximation for Jackson-Matsuoka Polynomials on the Sphere
Directory of Open Access Journals (Sweden)
Guo Feng
2012-01-01
Full Text Available We consider the best approximation by Jackson-Matsuoka polynomials in the weighted Lp space on the unit sphere of Rd. Using the relation between K-functionals and modulus of smoothness on the sphere, we obtain the direct and inverse estimate of approximation by these polynomials for the h-spherical harmonics.
Pointwise Approximation Theorems for Combinations and Derivatives of Bernstein Polynomials
Institute of Scientific and Technical Information of China (English)
Lin Sen XIE
2005-01-01
We establish the pointwise approximation theorems for the combinations of Bernstein polynomials by the rth Ditzian-Totik modulus of smoothness ωγφ(f, t) where φ is an admissible step-weight function. An equivalence relation between the derivatives of these polynomials and the smoothness of functions is also obtained.
Polynomial coordinates and their behavior in higher dimensions
Berson, J.J.T.
2004-01-01
A coordinate is an element of a polynomial ring which is the first component of some automorphism of this ring. Understanding the structure of coordinates is still one of the major problems in the theory of polynomial automorphisms. It is already known, that in two variables over a field every coord
A second addition formula for continuous q-ultraspherical polynomials
Koornwinder, T.H.
2005-01-01
This paper provides the details of Remark 5.4 in the author's paper "Askey-Wilson polynomials as zonal spherical functions on the SU(2) quantum group", SIAM J. Math. Anal. 24 (1993), 795-813. In formula (5.9) of the 1993 paper a two-parameter class of Askey-Wilson polynomials was expanded as a finit
Mutations of Laurent Polynomials and Flat Families with Toric Fibers
Ilten, Nathan Owen
2012-01-01
We give a general criterion for two toric varieties to appear as fibers in a flat family over the projective line. We apply this to show that certain birational transformations mapping a Laurent polynomial to another Laurent polynomial correspond to deformations between the associated toric varieties.
Strict Positivstellens\\" atze for matrix polynomials with scalar constraints
Cimpric, Jaka
2010-01-01
We extend Krivine's strict positivstellensatz for usual (real multivariate) polynomials to symmetric matrix polynomials with scalar constraints. The proof is an elementary computation with Schur complements. Analogous extensions of Schm\\" udgen's and Putinar's strict positivstellensatz were recently proved by Hol and Scherer using methods from optimization theory.
Animating Nested Taylor Polynomials to Approximate a Function
Mazzone, Eric F.; Piper, Bruce R.
2010-01-01
The way that Taylor polynomials approximate functions can be demonstrated by moving the center point while keeping the degree fixed. These animations are particularly nice when the Taylor polynomials do not intersect and form a nested family. We prove a result that shows when this nesting occurs. The animations can be shown in class or…
POLYNOMIAL RECURRENCE FOR L（E）VY PROCESSES
Institute of Scientific and Technical Information of China (English)
ZHAO MINZHI; YING JIANGANG
2004-01-01
In this paper, the authors study the ω-transience and ω-recurrence for Lévy processes with any weight function ω, give a relation between ω-recurrence and the last exit times. As a special case, the polynomial recurrence and polynomial transience are also studied.
Approximation to Continuous Functions by a Kind of Interpolation Polynomials
Institute of Scientific and Technical Information of China (English)
Yuan Xue-gang; Wang De-hui
2001-01-01
In this paper, an interpolation polynomial operator Fn (f; l, x ) is constructed based on the zeros of a kind of Jacobi polynomials as the interpolation nodes. For any continuous function f(x)∈ Cb[1,1] (0≤b≤l) Fn(f; l,x) converges to f(x) uniformly, where l is an odd number.
Lm Extremal Polynomials Associated with Generalized Jacobi Weights
Institute of Scientific and Technical Information of China (English)
Ying-guang Shi
2003-01-01
Asymptotic estimations of the Christoffel type functions for Lm extremal polynomials with an even integer m associated with generalized Jacobi weights are established. Also, asymptotic behavior of the zeros of the Lm extremal polynomials and the Cotes numbers of the corresponding Turan quadrature formula is given.
On Period of the Sequence of Fibonacci Polynomials Modulo
Directory of Open Access Journals (Sweden)
İnci Gültekin
2013-01-01
Full Text Available It is shown that the sequence obtained by reducing modulo coefficient and exponent of each Fibonacci polynomials term is periodic. Also if is prime, then sequences of Fibonacci polynomial are compared with Wall numbers of Fibonacci sequences according to modulo . It is found that order of cyclic group generated with matrix is equal to the period of these sequences.
A polynomial f(R) inflation model
Huang, Qing-Guo
2014-01-01
Motivated by the ultraviolet complete quantum theory of gravity, for example string theory, we investigate a polynomial $f(R)$ inflation model in detail. We calculate the spectral index and tensor-to-scalar ratio in the $f(R)$ inflation model with the form of $f(R)=R+{R^2\\over 6M^2}+{\\lambda_n\\over 2n}{R^n\\over (3M^2)^{n-1}}$. If the dimensionless coupling $\\lambda_n$ is much smaller than one, there are two regions for achieving the slow-roll inflation if $\\lambda_n>0$ and only one region if $\\lambda_n<0$. Compared to Planck 2013, numerically we find that $|\\lambda_n|\\lesssim 10^{-2n+2.6}$ in the region I and $0<\\lambda_n\\lesssim 10^{-3n+5.6}$ in the region II.
Irreducible polynomials with prescribed sums of coefficients
Tuxanidy, Aleksandr; Wang, Qiang
2016-01-01
Let $q$ be a power of a prime, let $\\mathbb{F}_q$ be the finite field with $q$ elements and let $n \\geq 2$. For a polynomial $h(x) \\in \\mathbb{F}_q[x]$ of degree $n \\in \\mathbb{N}$ and a subset $W \\subseteq [0,n] := \\{0, 1, \\ldots, n\\}$, we define the sum-of-digits function $$S_W(h) = \\sum_{w \\in W}[x^{w}] h(x)$$ to be the sum of all the coefficients of $x^w$ in $h(x)$ with $w \\in W$. In the case when $q = 2$, we prove, except for a few genuine exceptions, that for any $c \\in \\mathbb{F}_2$ an...
The evolution of piecewise polynomial wave functions
Andrews, Mark
2017-01-01
For a non-relativistic particle, we consider the evolution of wave functions that consist of polynomial segments, usually joined smoothly together. These spline wave functions are compact (that is, they are initially zero outside a finite region), but they immediately extend over all available space as they evolve. The simplest splines are the square and triangular wave functions in one dimension, but very complicated splines have been used in physics. In general the evolution of such spline wave functions can be expressed in terms of antiderivatives of the propagator; in the case of a free particle or an oscillator, all the evolutions are expressed exactly in terms of Fresnel integrals. Some extensions of these methods to two and three dimensions are discussed.
Inequalities for a Polynomial and its Derivative
Indian Academy of Sciences (India)
V K Jain
2000-05-01
For an arbitrary entire function and any > 0, let $M(f, r):=\\max_{|z|=r}|f(z)|$. It is known that if is a polynomial of degree having no zeros in the open unit disc, and $m:=\\min_{|z|=1}|p(z)|$, then $$M(p',1)≤\\frac{n}{2}\\{M(p,1)-m\\},$$ $$M(p, R)≤\\left(\\frac{R^n+1}{2}\\right)M(p, 1)-m\\left(\\frac{R^n-1}{2}\\right), R>> 1.$$ It is also known that if has all its zeros in the closed unit disc, then $$M(p', 1)≥\\frac{n}{2}\\{M(p, 1)+m\\}.$$ The present paper contains certain generalizations of these inequalities
Study on the Grey Polynomial Geometric Programming
Institute of Scientific and Technical Information of China (English)
LUODang
2005-01-01
In the model of geometric programming, values of parameters cannot be gotten owing to data fluctuation and incompletion. But reasonable bounds of these parameters can be attained. This is to say, parameters of this model can be regarded as interval grey numbers. When the model contains grey numbers, it is hard for common programming method to solve them. By combining the common programming model with the grey system theory,and using some analysis strategies, a model of grey polynomial geometric programming, a model of 8 positioned geometric programming and their quasi-optimum solution or optimum solution are put forward. At the same time, we also developed an algorithm for the problem.This approach brings a new way for the application research of geometric programming. An example at the end of this paper shows the rationality and feasibility of the algorithm.
A polynomial algorithm for abstract maximum flow
Energy Technology Data Exchange (ETDEWEB)
McCormick, S.T. [Univ. of British Columbia, Vancouver, British Columbia (Canada)
1996-12-31
Ford and Fulkerson`s original 1956 max flow/min cut paper formulated max flow in terms of flows on paths, rather than the more familiar flows on arcs. In 1974 Hoffman pointed out that Ford and Fulkerson`s original proof was quite abstract, and applied to a wide range of max flow-like problems. In this abstract model we have capacitated elements, and linearly ordered subsets of elements called paths. When two paths share an element ({open_quote}cross{close_quote}), then there must be a path that is a subset of the first path up to the cross, and a subset of the second path after the cross. (Hoffman`s generalization of) Ford and Fulkerson`s proof showed that the max flow/min cut theorem still holds under this weak assumption. However, this proof is non-constructive. To get an algorithm, we assume that we have an oracle whose input is an arbitrary subset of elements, and whose output is either a path contained in that subset, or the statement that no such path exists. We then use complementary slackness to show how to augment any feasible set of path flows to a set with a strictly larger total flow value using a polynomial number of calls to the oracle. Then standard scaling techniques yield an overall polynomial algorithm for finding both a max flow and a min cut. Hoffman`s paper actually considers a sort of supermodular objective on the path flows, which allows him to include transportation problems and thus rain-cost flow in his frame-work. We also discuss extending our algorithm to this more general case.
Structured matrix based methods for approximate polynomial GCD
Boito, Paola
2011-01-01
Defining and computing a greatest common divisor of two polynomials with inexact coefficients is a classical problem in symbolic-numeric computation. The first part of this book reviews the main results that have been proposed so far in the literature. As usual with polynomial computations, the polynomial GCD problem can be expressed in matrix form: the second part of the book focuses on this point of view and analyses the structure of the relevant matrices, such as Toeplitz, Toepliz-block and displacement structures. New algorithms for the computation of approximate polynomial GCD are presented, along with extensive numerical tests. The use of matrix structure allows, in particular, to lower the asymptotic computational cost from cubic to quadratic order with respect to polynomial degree. .
Higher order branching of periodic orbits from polynomial isochrones
Directory of Open Access Journals (Sweden)
B. Toni
1999-09-01
Full Text Available We discuss the higher order local bifurcations of limit cycles from polynomial isochrones (linearizable centers when the linearizing transformation is explicitly known and yields a polynomial perturbation one-form. Using a method based on the relative cohomology decomposition of polynomial one-forms complemented with a step reduction process, we give an explicit formula for the overall upper bound of branch points of limit cycles in an arbitrary $n$ degree polynomial perturbation of the linear isochrone, and provide an algorithmic procedure to compute the upper bound at successive orders. We derive a complete analysis of the nonlinear cubic Hamiltonian isochrone and show that at most nine branch points of limit cycles can bifurcate in a cubic polynomial perturbation. Moreover, perturbations with exactly two, three, four, six, and nine local families of limit cycles may be constructed.
Asymptotic behaviour of zeros of exceptional Jacobi and Laguerre polynomials
Gómez-Ullate, David; Milson, Robert
2012-01-01
The location and asymptotic behaviour for large n of the zeros of exceptional Jacobi and Laguerre polynomials are discussed. The zeros of exceptional polynomials fall into two classes: the regular zeros, which lie in the interval of orthogonality and the exceptional zeros, which lie outside that interval. We show that the regular zeros have two interlacing properties: one is the natural interlacing between consecutive polynomials as a consequence of their Sturm-Liouville character, while the other one shows interlacing between the zeros of exceptional and classical polynomials. A generalization of the classical Heine-Mehler formula is provided for the exceptional polynomials, which allows to derive the asymptotic behaviour of their regular zeros. We also describe the location and the asymptotic behaviour of the exceptional zeros, which converge for large n to fixed values.
The lowest degree $0,1$-polynomial divisible by cyclotomic polynomial
Reddy, A Satyanarayana
2011-01-01
Let $n$ be an even positive integer with at most three distinct prime factors and let $\\ze_n$ be a primitive $n$-th root of unity. In this study, we made an attempt to find the lowest-degree $0,1$-polynomial $f(x) \\in \\Q[x]$ having at least three terms such that $f(\\ze_n)$ is a minimal vanishing sum of the $n$-th roots of unity.
Etessami, Kousha; Yannakakis, Mihalis
2012-01-01
We show that one can approximate the least fixed point solution for a multivariate system of monotone probabilistic max(min) polynomial equations, referred to as maxPPSs (and minPPSs, respectively), in time polynomial in both the encoding size of the system of equations and in log(1/epsilon), where epsilon > 0 is the desired additive error bound of the solution. (The model of computation is the standard Turing machine model.) We establish this result using a generalization of Newton's method which applies to maxPPSs and minPPSs, even though the underlying functions are only piecewise-differentiable. This generalizes our recent work which provided a P-time algorithm for purely probabilistic PPSs. These equations form the Bellman optimality equations for several important classes of infinite-state Markov Decision Processes (MDPs). Thus, as a corollary, we obtain the first polynomial time algorithms for computing to within arbitrary desired precision the optimal value vector for several classes of infinite-state...
Linear precoding based on polynomial expansion: reducing complexity in massive MIMO
Mueller, Axel; Kammoun, Abla; Björnson, Emil; Debbah, Merouane
2016-01-01
Massive multiple-input multiple-output (MIMO) techniques have the potential to bring tremendous improvements in spectral efficiency to future communication systems. Counterintuitively, the practical issues of having uncertain channel knowledge, high propagation losses, and implementing optimal non-linear precoding are solved more or less automatically by enlarging system dimensions. However, the computational precoding complexity grows with the system dimensions. For example, the close-to-opt...
Efficient linear precoding for massive MIMO systems using truncated polynomial expansion
Muller, Axel; Kammoun, Abla; Bjornson, Emil; Debbah, Merouane
2014-01-01
International audience; —Massive multiple-input multiple-output (MIMO) techniques have been proposed as a solution to satisfy many requirements of next generation cellular systems. One downside of massive MIMO is the increased complexity of computing the precoding, especially since the relatively "antenna-efficient" regu-larized zero-forcing (RZF) is preferred to simple maximum ratio transmission. We develop in this paper a new class of precoders for single-cell massive MIMO systems. It is ba...
Desai, A.; Witteveen, J.A.S.; Sarkar, S.
2013-01-01
The present study focuses on the uncertainty quantification of an aeroelastic instability system. This is a classical dynamical system often used to model the flow induced oscillation of flexible structures such as turbine blades. It is relevant as a preliminary fluid-structure interaction model, su
An Enhanced Asymptotic Expansion for the Stability of Nonlinear Elastic Structures
DEFF Research Database (Denmark)
Christensen, Claus Dencker; Byskov, Esben
2010-01-01
A new, enhanced asymptotic expansion applicable to stability of structures made of nonlinear elastic materials is established. The method utilizes “hyperbolic” terms instead of the conventional polynomial terms, covers full kinematic nonlinearity and is applied to nonlinear elastic Euler columns ...
Conformal Laplace superintegrable systems in 2D: polynomial invariant subspaces
Escobar-Ruiz, M. A.; Miller, Willard, Jr.
2016-07-01
2nd-order conformal superintegrable systems in n dimensions are Laplace equations on a manifold with an added scalar potential and 2n-1 independent 2nd order conformal symmetry operators. They encode all the information about Helmholtz (eigenvalue) superintegrable systems in an efficient manner: there is a 1-1 correspondence between Laplace superintegrable systems and Stäckel equivalence classes of Helmholtz superintegrable systems. In this paper we focus on superintegrable systems in two-dimensions, n = 2, where there are 44 Helmholtz systems, corresponding to 12 Laplace systems. For each Laplace equation we determine the possible two-variate polynomial subspaces that are invariant under the action of the Laplace operator, thus leading to families of polynomial eigenfunctions. We also study the behavior of the polynomial invariant subspaces under a Stäckel transform. The principal new results are the details of the polynomial variables and the conditions on parameters of the potential corresponding to polynomial solutions. The hidden gl 3-algebraic structure is exhibited for the exact and quasi-exact systems. For physically meaningful solutions, the orthogonality properties and normalizability of the polynomials are presented as well. Finally, for all Helmholtz superintegrable solvable systems we give a unified construction of one-dimensional (1D) and two-dimensional (2D) quasi-exactly solvable potentials possessing polynomial solutions, and a construction of new 2D PT-symmetric potentials is established.
Constructing general partial differential equations using polynomial and neural networks.
Zjavka, Ladislav; Pedrycz, Witold
2016-01-01
Sum fraction terms can approximate multi-variable functions on the basis of discrete observations, replacing a partial differential equation definition with polynomial elementary data relation descriptions. Artificial neural networks commonly transform the weighted sum of inputs to describe overall similarity relationships of trained and new testing input patterns. Differential polynomial neural networks form a new class of neural networks, which construct and solve an unknown general partial differential equation of a function of interest with selected substitution relative terms using non-linear multi-variable composite polynomials. The layers of the network generate simple and composite relative substitution terms whose convergent series combinations can describe partial dependent derivative changes of the input variables. This regression is based on trained generalized partial derivative data relations, decomposed into a multi-layer polynomial network structure. The sigmoidal function, commonly used as a nonlinear activation of artificial neurons, may transform some polynomial items together with the parameters with the aim to improve the polynomial derivative term series ability to approximate complicated periodic functions, as simple low order polynomials are not able to fully make up for the complete cycles. The similarity analysis facilitates substitutions for differential equations or can form dimensional units from data samples to describe real-world problems.
vs. a polynomial chaos-based MCMC
Siripatana, Adil
2014-08-01
Bayesian Inference of Manning\\'s n coefficient in a Storm Surge Model Framework: comparison between Kalman lter and polynomial based method Adil Siripatana Conventional coastal ocean models solve the shallow water equations, which describe the conservation of mass and momentum when the horizontal length scale is much greater than the vertical length scale. In this case vertical pressure gradients in the momentum equations are nearly hydrostatic. The outputs of coastal ocean models are thus sensitive to the bottom stress terms de ned through the formulation of Manning\\'s n coefficients. This thesis considers the Bayesian inference problem of the Manning\\'s n coefficient in the context of storm surge based on the coastal ocean ADCIRC model. In the first part of the thesis, we apply an ensemble-based Kalman filter, the singular evolutive interpolated Kalman (SEIK) filter to estimate both a constant Manning\\'s n coefficient and a 2-D parameterized Manning\\'s coefficient on one ideal and one of more realistic domain using observation system simulation experiments (OSSEs). We study the sensitivity of the system to the ensemble size. we also access the benefits from using an in ation factor on the filter performance. To study the limitation of the Guassian restricted assumption on the SEIK lter, 5 we also implemented in the second part of this thesis a Markov Chain Monte Carlo (MCMC) method based on a Generalized Polynomial chaos (gPc) approach for the estimation of the 1-D and 2-D Mannning\\'s n coe cient. The gPc is used to build a surrogate model that imitate the ADCIRC model in order to make the computational cost of implementing the MCMC with the ADCIRC model reasonable. We evaluate the performance of the MCMC-gPc approach and study its robustness to di erent OSSEs scenario. we also compare its estimates with those resulting from SEIK in term of parameter estimates and full distributions. we present a full analysis of the solution of these two methods, of the
Automorphisms of Algebras and Bochner's Property for Vector Orthogonal Polynomials
Horozov, Emil
2016-05-01
We construct new families of vector orthogonal polynomials that have the property to be eigenfunctions of some differential operator. They are extensions of the Hermite and Laguerre polynomial systems. A third family, whose first member has been found by Y. Ben Cheikh and K. Douak is also constructed. The ideas behind our approach lie in the studies of bispectral operators. We exploit automorphisms of associative algebras which transform elementary vector orthogonal polynomial systems which are eigenfunctions of a differential operator into other systems of this type.
Combinatorial theory of Macdonald polynomials I: Proof of Haglund's formula
Haglund, J.; Haiman, M.; Loehr, N
2005-01-01
Haglund recently proposed a combinatorial interpretation of the modified Macdonald polynomials H̃μ. We give a combinatorial proof of this conjecture, which establishes the existence and integrality of H̃μ. As corollaries, we obtain the cocharge formula of Lascoux and Schützenberger for Hall–Littlewood polynomials, a formula of Sahi and Knop for Jack's symmetric functions, a generalization of this result to the integral Macdonald polynomials Jμ, a formula for H̃μ in terms of Lascoux–Leclerc–Th...
Representations of non-negative polynomials via critical ideals
Hiep, Dang Tuan
2011-01-01
This paper studies the representations of a non-negative polynomial $f$ on a non-compact semi-algebraic set $K$ modulo its critical ideal. Under the assumptions that the semi-algebraic set $K$ is regular and $f$ satisfies the boundary Hessian conditions (BHC) at each zero of $f$ in $K$, we show that $f$ can be represented as a sum of squares (SOS) of real polynomials modulo its critical ideal if $f\\ge 0$ on $K$. In particular, we focus on the polynomial ring $\\mathbb R[x]$.
Polynomial system solving for decoding linear codes and algebraic cryptanalysis
2009-01-01
This thesis is devoted to applying symbolic methods to the problems of decoding linear codes and of algebraic cryptanalysis. The paradigm we employ here is as follows. We reformulate the initial problem in terms of systems of polynomial equations over a finite field. The solution(s) of such systems should yield a way to solve the initial problem. Our main tools for handling polynomials and polynomial systems in such a paradigm is the technique of Gröbner bases and normal form reductions. The ...
Generalized Freud's equation and level densities with polynomial potential
Boobna, Akshat; Ghosh, Saugata
2013-08-01
We study orthogonal polynomials with weight $\\exp[-NV(x)]$, where $V(x)=\\sum_{k=1}^{d}a_{2k}x^{2k}/2k$ is a polynomial of order 2d. We derive the generalised Freud's equations for $d=3$, 4 and 5 and using this obtain $R_{\\mu}=h_{\\mu}/h_{\\mu -1}$, where $h_{\\mu}$ is the normalization constant for the corresponding orthogonal polynomials. Moments of the density functions, expressed in terms of $R_{\\mu}$, are obtained using Freud's equation and using this, explicit results of level densities as $N\\rightarrow\\infty$ are derived.
Non-Hermitian oscillator Hamiltonians and multiple Charlier polynomials
Energy Technology Data Exchange (ETDEWEB)
Miki, Hiroshi, E-mail: miki@amp.i.kyoto-u.ac.jp [Department of Applied Mathematics and Physics, Graduate School of Informatics, Kyoto University, Sakyo-Ku, Kyoto 606 8501 (Japan); Vinet, Luc, E-mail: luc.vinet@umontreal.ca [Centre de recherches mathématiques, Université de Montréal, P.O. Box 6128, Centre-ville Station, Montréal (Québec), H3C 3J7 (Canada); Zhedanov, Alexei, E-mail: zhedanov@fti.dn.ua [Donetsk Institute for Physics and Technology, Donetsk 83 114 (Ukraine)
2011-12-05
A set of r non-Hermitian oscillator Hamiltonians in r dimensions is shown to be simultaneously diagonalizable. Their spectra are real and the common eigenstates are expressed in terms of multiple Charlier polynomials. An algebraic interpretation of these polynomials is thus achieved and the model is used to derive some of their properties. -- Highlights: ► A set of r non-Hermitian oscillator Hamiltonians in r dimensions is presented. ► Their spectra are real. ► The common eigenstates are expressed in terms of multiple Charlier polynomials (MCP). ► This “integrable” model allows to interpret structural formulas of the MCPs.
Chromatic Polynomials Of Some (m,l-Hyperwheels
Directory of Open Access Journals (Sweden)
Julian A. Allagan
2014-03-01
Full Text Available In this paper, using a standard method of computing the chromatic polynomial of hypergraphs, we introduce a new reduction theorem which allows us to find explicit formulae for the chromatic polynomials of some (complete non-uniform $(m,l-$hyperwheels and non-uniform $(m,l-$hyperfans. These hypergraphs, constructed through a ``join" graph operation, are some generalizations of the well-known wheel and fan graphs, respectively. Further, we revisit some results concerning these graphs and present their chromatic polynomials in a standard form that involves the Stirling numbers of the second kind.
Ratio Monotonicity of Polynomials Derived from Nondecreasing Sequences
Chen, William Y C; Zhou, Elaine L F
2010-01-01
The ratio monotonicity of a polynomial is a stronger property than log-concavity. Let P(x) be a polynomial with nonnegative and nondecreasing coefficients. We prove the ratio monotone property of P(x+1), which leads to the log-concavity of P(x+c) for any $c\\geq 1$ due to Llamas and Mart\\'{\\i}nez-Bernal. As a consequence, we obtain the ratio monotonicity of the Boros-Moll polynomials obtained by Chen and Xia without resorting to the recurrence relations of the coefficients.
Polynomials in Control Theory Parametrized by Their Roots
Directory of Open Access Journals (Sweden)
Baltazar Aguirre-Hernández
2012-01-01
Full Text Available The aim of this paper is to introduce the space of roots to study the topological properties of the spaces of polynomials. Instead of identifying a monic complex polynomial with the vector of its coefficients, we identify it with the set of its roots. Viète's map gives a homeomorphism between the space of roots and the space of coefficients and it gives an explicit formula to relate both spaces. Using this viewpoint we establish that the space of monic (Schur or Hurwitz aperiodic polynomials is contractible. Additionally we obtain a Boundary Theorem.
An eigenvalue problem for the associated Askey-Wilson polynomials
Bruder, Andrea; Suslov, Sergei K
2012-01-01
To derive an eigenvalue problem for the associated Askey-Wilson polynomials, we consider an auxiliary function in two variables which is related to the associated Askey-Wilson polynomials introduced by Ismail and Rahman. The Askey-Wilson operator, applied in each variable separately, maps this function to the ordinary Askey-Wilson polynomials with different sets of parameters. A third Askey-Wilson operator is found with the help of a computer algebra program which links the two, and an eigenvalue problem is stated.
q-Bernoulli numbers and q-Bernoulli polynomials revisited
Directory of Open Access Journals (Sweden)
Kim Taekyun
2011-01-01
Full Text Available Abstract This paper performs a further investigation on the q-Bernoulli numbers and q-Bernoulli polynomials given by Acikgöz et al. (Adv Differ Equ, Article ID 951764, 9, 2010, some incorrect properties are revised. It is point out that the generating function for the q-Bernoulli numbers and polynomials is unreasonable. By using the theorem of Kim (Kyushu J Math 48, 73-86, 1994 (see Equation 9, some new generating functions for the q-Bernoulli numbers and polynomials are shown. Mathematics Subject Classification (2000 11B68, 11S40, 11S80
Learning Read-constant Polynomials of Constant Degree modulo Composites
DEFF Research Database (Denmark)
Chattopadhyay, Arkadev; Gavaldá, Richard; Hansen, Kristoffer Arnsfelt;
2011-01-01
Boolean functions that have constant degree polynomial representation over a fixed finite ring form a natural and strict subclass of the complexity class \\textACC0ACC0. They are also precisely the functions computable efficiently by programs over fixed and finite nilpotent groups. This class...... is not known to be learnable in any reasonable learning model. In this paper, we provide a deterministic polynomial time algorithm for learning Boolean functions represented by polynomials of constant degree over arbitrary finite rings from membership queries, with the additional constraint that each variable...
Polynomial complexity algorithm for Max-Cut problem
Katkov, Mikhail
2010-01-01
The standard NP-complete max-cut problem is reformulated as a binary quadratic program xQx s.t x^2=1. This problem is further reformulated as global minimum of quartic polynomial (xQ'x - z)^2 + \\sum_i (x_i^2-1)^2+ \\alpha z^2, for some \\alpha. The global minimum is found by polynomial complexity semi-definite program. Numerical examples and code is provided. The resulting algorithm solves arbitrary max-cut problem in polynomial time, therefore P=NP.
Pamplona, Djenane C; Velloso, Raquel Q; Radwanski, Henrique N
2014-01-01
This article discusses skin expansion without considering cellular growth of the skin. An in vivo analysis was carried out that involved expansion at three different sites on one patient, allowing for the observation of the relaxation process. Those measurements were used to characterize the human skin of the thorax during the surgical process of skin expansion. A comparison between the in vivo results and the numerical finite elements model of the expansion was used to identify the material elastic parameters of the skin of the thorax of that patient. Delfino's constitutive equation was chosen to model the in vivo results. The skin is considered to be an isotropic, homogeneous, hyperelastic, and incompressible membrane. When the skin is extended, such as with expanders, the collagen fibers are also extended and cause stiffening in the skin, which results in increasing resistance to expansion or further stretching. We observed this phenomenon as an increase in the parameters as subsequent expansions continued. The number and shape of the skin expanders used in expansions were also studied, both mathematically and experimentally. The choice of the site where the expansion should be performed is discussed to enlighten problems that can lead to frustrated skin expansions. These results are very encouraging and provide insight into our understanding of the behavior of stretched skin by expansion. To our knowledge, this study has provided results that considerably improve our understanding of the behavior of human skin under expansion.
Lower bounds for polynomials using geometric programming
Ghasemi, Mehdi
2011-01-01
We make use of a result of Hurwitz and Reznick, and a consequence of this result due to Fidalgo and Kovacec, to determine a new sufficient condition for a polynomial $f\\in\\mathbb{R}[X_1,...,X_n]$ of even degree to be a sum of squares. This result generalizes a result of Lasserre and a result of Fidalgo and Kovacec, and it also generalizes the improvements of these results given in [6]. We apply this result to obtain a new lower bound $f_{gp}$ for $f$, and we explain how $f_{gp}$ can be computed using geometric programming. The lower bound $f_{gp}$ is generally not as good as the lower bound $f_{sos}$ introduced by Lasserre and Parrilo and Sturmfels, which is computed using semidefinite programming, but a run time comparison shows that, in practice, the computation of $f_{gp}$ is much faster. The computation is simplest when the highest degree term of $f$ has the form $\\sum_{i=1}^n a_iX_i^{2d}$, $a_i>0$, $i=1,...,n$. The lower bounds for $f$ established in [6] are obtained by evaluating the objective function ...
Integral Inequalities for Self-Reciprocal Polynomials
Indian Academy of Sciences (India)
Horst Alzer
2010-04-01
Let $n≥ 1$ be an integer and let $\\mathcal{P}_n$ be the class of polynomials of degree at most satisfying $z^nP(1/z)=P(z)$ for all $z\\in C$. Moreover, let be an integer with $1≤ r≤ n$. Then we have for all $P\\in\\mathcal{P}_n$: $$_n(r)\\int^{2}_0|P(e^{it})|^2dt≤\\int^{2}_0|P^{(r)}(e^{it})|^2dt≤_n(r)\\int^{2}_0|P(e^{it})|^2dt$$ with the best possible factors \\begin{equation*}_n(r)=\\begin{cases}\\prod^{r-1}_{j=0}\\left(\\frac{n}{2}-j\\right)^2, < \\text{if is even},\\\\ \\frac{1}{2}\\left[\\prod^{r-1}_{j=0}\\left(\\frac{n+1}{2}-j\\right)^2+\\prod^{r-1}_{j=0}\\left(\\frac{n-1}{2}-j\\right)^2\\right], < \\text{if is odd},\\end{cases}\\end{equation*} and \\begin{equation*}_n(r)=\\frac{1}{2}\\prod\\limits^{r-1}_{j=0}(n-j)^2.\\end{equation*} This refines and extends a result due to Aziz and Zargar (1997).
Polynomial Linear Programming with Gaussian Belief Propagation
Bickson, Danny; Shental, Ori; Dolev, Danny
2008-01-01
Interior-point methods are state-of-the-art algorithms for solving linear programming (LP) problems with polynomial complexity. Specifically, the Karmarkar algorithm typically solves LP problems in time O(n^{3.5}), where $n$ is the number of unknown variables. Karmarkar's celebrated algorithm is known to be an instance of the log-barrier method using the Newton iteration. The main computational overhead of this method is in inverting the Hessian matrix of the Newton iteration. In this contribution, we propose the application of the Gaussian belief propagation (GaBP) algorithm as part of an efficient and distributed LP solver that exploits the sparse and symmetric structure of the Hessian matrix and avoids the need for direct matrix inversion. This approach shifts the computation from realm of linear algebra to that of probabilistic inference on graphical models, thus applying GaBP as an efficient inference engine. Our construction is general and can be used for any interior-point algorithm which uses the Newt...
Polynomial super-gl(n) algebras
Jarvis, P D
2003-01-01
We introduce a class of finite dimensional nonlinear superalgebras $L = L_{\\bar{0}} + L_{\\bar{1}}$ providing gradings of $L_{\\bar{0}} = gl(n) \\simeq sl(n) + gl(1)$. Odd generators close by anticommutation on polynomials (of degree $>1$) in the $gl(n)$ generators. Specifically, we investigate `type I' super-$gl(n)$ algebras, having odd generators transforming in a single irreducible representation of $gl(n)$ together with its contragredient. Admissible structure constants are discussed in terms of available $gl(n)$ couplings, and various special cases and candidate superalgebras are identified and exemplified via concrete oscillator constructions. For the case of the $n$-dimensional defining representation, with odd generators $Q_{a}, \\bar{Q}{}^{b}$, and even generators ${E^{a}}_{b}$, $a,b = 1,...,n$, a three parameter family of quadratic super-$gl(n)$ algebras (deformations of $sl(n/1)$) is defined. In general, additional covariant Serre-type conditions are imposed, in order that the Jacobi identities be fulf...
Polynomial-time solutions to image segmentation
Energy Technology Data Exchange (ETDEWEB)
Asano, Tetsuo [Osaka Electro-Communication Univ., Neyagawa (Japan); Chen, D.Z. [Notre Dame, South Bend, IN (United States); Katoh, Naoki [Kobe Univ. of Commerce (Japan)
1996-12-31
Separating an object in an image from its background is a central problem (called segmentation) in pattern recognition and computer vision. In this paper, we study the complexity of the segmentation problem, assuming that the object forms a connected region in an intensity image. We show that the optimization problem of separating a connected region in an n-pixel grid is NP-hard under the interclass variance, a criterion that is used in discriminant analysis. More importantly, we consider the basic case in which the object is separated by two x-monotone curves (i.e., the object itself is x-monotone), and present polynomial-time algorithms for computing exact and approximate optimal segmentation. Our main algorithm for exact optimal segmentation by two x-monotone curves runs in O(n{sup 2}) time; this algorithm is based on several techniques such as a parametric optimization formulation, a hand-probing algorithm for the convex hull of an unknown point set, and dynamic programming using fast matrix searching. Our efficient approximation scheme obtains an {epsilon}-approximate solution in O({epsilon}{sup -1} n log L) time, where {epsilon} is any fixed constant with 1 > {epsilon} > 0, and L is the total sum of the absolute values of brightness levels of the image.
Polynomial Time Algorithms for Minimum Energy Scheduling
Baptiste, Philippe; Durr, Christoph
2009-01-01
The aim of power management policies is to reduce the amount of energy consumed by computer systems while maintaining satisfactory level of performance. One common method for saving energy is to simply suspend the system during the idle times. No energy is consumed in the suspend mode. However, the process of waking up the system itself requires a certain fixed amount of energy, and thus suspending the system is beneficial only if the idle time is long enough to compensate for this additional energy expenditure. In the specific problem studied in the paper, we have a set of jobs with release times and deadlines that need to be executed on a single processor. Preemptions are allowed. The processor requires energy L to be woken up and, when it is on, it uses one unit of energy per one unit of time. It has been an open problem whether a schedule minimizing the overall energy consumption can be computed in polynomial time. We solve this problem in positive, by providing an O(n^5)-time algorithm. In addition we pr...
Polynomial Method for PLL Controller Optimization
Directory of Open Access Journals (Sweden)
Tsung-Yu Chiou
2011-06-01
Full Text Available The Phase-Locked Loop (PLL is a key component of modern electronic communication and control systems. PLL is designed to extract signals from transmission channels. It plays an important role in systems where it is required to estimate the phase of a received signal, such as carrier tracking from global positioning system satellites. In order to robustly provide centimeter-level accuracy, it is crucial for the PLL to estimate the instantaneous phase of an incoming signal which is usually buried in random noise or some type of interference. This paper presents an approach that utilizes the recent development in the semi-definite programming and sum-of-squares field. A Lyapunov function will be searched as the certificate of the pull-in range of the PLL system. Moreover, a polynomial design procedure is proposed to further refine the controller parameters for system response away from the equilibrium point. Several simulation results as well as an experiment result are provided to show the effectiveness of this approach.
Constructing Polynomial Spectral Models for Stars
Rix, Hans-Walter; Conroy, Charlie; Hogg, David W
2016-01-01
Stellar spectra depend on the stellar parameters and on dozens of photospheric elemental abundances. Simultaneous fitting of these $\\mathcal{N}\\sim \\,$10-40 model labels to observed spectra has been deemed unfeasible, because the number of ab initio spectral model grid calculations scales exponentially with $\\mathcal{N}$. We suggest instead the construction of a polynomial spectral model (PSM) of order $\\mathcal{O}$ for the model flux at each wavelength. Building this approximation requires a minimum of only ${\\mathcal{N}+\\mathcal{O}\\choose\\mathcal{O}}$ calculations: e.g. a quadratic spectral model ($\\mathcal{O}=\\,$2), which can then fit $\\mathcal{N}=\\,$20 labels simultaneously, can be constructed from as few as 231 ab initio spectral model calculations; in practice, a somewhat larger number ($\\sim\\,$300-1000) of randomly chosen models lead to a better performing PSM. Such a PSM can be a good approximation to ab initio spectral models only over a limited portion of label space, which will vary case by case. Y...
Classification of complex polynomial vector fields in one complex variable
DEFF Research Database (Denmark)
Branner, Bodil; Dias, Kealey
2010-01-01
, the main result of the paper. This result is an extension and refinement of Douady et al. (Champs de vecteurs polynomiaux sur C. Unpublished manuscript) classification of the structurally stable polynomial vector fields. We further review some general concepts for completeness and show that vector fields......This paper classifies the global structure of monic and centred one-variable complex polynomial vector fields. The classification is achieved by means of combinatorial and analytic data. More specifically, given a polynomial vector field, we construct a combinatorial invariant, describing...... the topology, and a set of analytic invariants, describing the geometry. Conversely, given admissible combinatorial and analytic data sets, we show using surgery the existence of a unique monic and centred polynomial vector field realizing the given invariants. This is the content of the Structure Theorem...
Sparse Volterra and Polynomial Regression Models: Recoverability and Estimation
Kekatos, Vassilis
2011-01-01
Volterra and polynomial regression models play a major role in nonlinear system identification and inference tasks. Exciting applications ranging from neuroscience to genome-wide association analysis build on these models with the additional requirement of parsimony. This requirement has high interpretative value, but unfortunately cannot be met by least-squares based or kernel regression methods. To this end, compressed sampling (CS) approaches, already successful in linear regression settings, can offer a viable alternative. The viability of CS for sparse Volterra and polynomial models is the core theme of this work. A common sparse regression task is initially posed for the two models. Building on (weighted) Lasso-based schemes, an adaptive RLS-type algorithm is developed for sparse polynomial regressions. The identifiability of polynomial models is critically challenged by dimensionality. However, following the CS principle, when these models are sparse, they could be recovered by far fewer measurements. ...
Force prediction in cold rolling mills by polynomial methods
Directory of Open Access Journals (Sweden)
Nicu ROMAN
2007-12-01
Full Text Available A method for steel and aluminium strip thickness control is provided including a new technique for predictive rolling force estimation method by statistic model based on polynomial techniques.
Generalized Rayleigh and Jacobi Processes and Exceptional Orthogonal Polynomials
Chou, C.-I.; Ho, C.-L.
2013-09-01
We present four types of infinitely many exactly solvable Fokker-Planck equations, which are related to the newly discovered exceptional orthogonal polynomials. They represent the deformed versions of the Rayleigh process and the Jacobi process.
Log-concavity of the genus polynomials of Ringel Ladders
Directory of Open Access Journals (Sweden)
Jonathan L Gross
2015-10-01
Full Text Available A Ringel ladder can be formed by a self-bar-amalgamation operation on a symmetric ladder, that is, by joining the root vertices on its end-rungs. The present authors have previously derived criteria under which linear chains of copies of one or more graphs have log-concave genus polyno- mials. Herein we establish Ringel ladders as the first significant non-linear infinite family of graphs known to have log-concave genus polynomials. We construct an algebraic representation of self-bar-amalgamation as a matrix operation, to be applied to a vector representation of the partitioned genus distribution of a symmetric ladder. Analysis of the resulting genus polynomial involves the use of Chebyshev polynomials. This paper continues our quest to affirm the quarter-century-old conjecture that all graphs have log-concave genus polynomials.
Maximum likelihood polynomial regression for robust speech recognition
Institute of Scientific and Technical Information of China (English)
LU Yong; WU Zhenyang
2011-01-01
The linear hypothesis is the main disadvantage of maximum likelihood linear re- gression （MLLR）. This paper applies the polynomial regression method to model adaptation and establishes a nonlinear model adaptation algorithm using maximum likelihood polyno
Certified counting of roots of random univariate polynomials
Cleveland, Joseph; Hauenstein, Jonathan D; Haywood, Ian; Mehta, Dhagash; Morse, Anthony; Robol, Leonardo; Schlenk, Taylor
2014-01-01
A challenging problem in computational mathematics is to compute roots of a high-degree univariate random polynomial. We combine an efficient multiprecision implementation for solving high-degree random polynomials with two certification methods, namely Smale's $\\alpha$-theory and one based on Gerschgorin's theorem, for showing that a given numerical approximation is in the quadratic convergence region of Newton's method of some exact solution. With this combination, we can certifiably count the number of real roots of random polynomials. We quantify the difference between the two certification procedures and list the salient features of both of them. After benchmarking on random polynomials where the coefficients are drawn from the Gaussian distribution, we obtain novel experimental results for the Cauchy distribution case.
Kauffman polynomials of some links and invariants of 3-manifolds
Institute of Scientific and Technical Information of China (English)
李起升
2002-01-01
Kauffman bracket polynomials of the so-called generalized tree-like links are studied. An algorithm of Witten type invariants, which was defined by Blanchet and Habegger et al. of more general 3-manifolds is given.
Estimation in the polynomial errors-in-variables model
Institute of Scientific and Technical Information of China (English)
无
2002-01-01
Estimators are presented for the coefficients of the polynomial errors-in-variables (EV) model when replicated observations are taken at some experimental points. These estimators are shown to be strongly consistent under mild conditions.
Guts of surfaces and the colored Jones polynomial
Futer, David; Purcell, Jessica
2013-01-01
This monograph derives direct and concrete relations between colored Jones polynomials and the topology of incompressible spanning surfaces in knot and link complements. Under mild diagrammatic hypotheses, we prove that the growth of the degree of the colored Jones polynomials is a boundary slope of an essential surface in the knot complement. We show that certain coefficients of the polynomial measure how far this surface is from being a fiber for the knot; in particular, the surface is a fiber if and only if a particular coefficient vanishes. We also relate hyperbolic volume to colored Jones polynomials. Our method is to generalize the checkerboard decompositions of alternating knots. Under mild diagrammatic hypotheses, we show that these surfaces are essential, and obtain an ideal polyhedral decomposition of their complement. We use normal surface theory to relate the pieces of the JSJ decomposition of the complement to the combinatorics of certain surface spines (state graphs). Since state graphs have p...
Pure states, positive matrix polynomials and sums of hermitian squares
Klep, Igor
2009-01-01
Let M be an archimedean quadratic module of real t-by-t matrix polynomials in n variables, and let S be the set of all real n-tuples where each element of M is positive semidefinite. Our key finding is a natural bijection between the set of pure states of M and the cartesian product of S with the real projective (t-1)-space. This leads us to conceptual proofs of positivity certificates for matrix polynomials, including the recent seminal result of Hol and Scherer: If a symmetric matrix polynomial is positive definite on S, then it belongs to M. We also discuss what happens for non-symmetric matrix polynomials or in the absence of the archimedean assumption, and review some of the related classical results. The methods employed are both algebraic and functional analytic.
A complete algorithm to find exact minimal polynomial by approximations
Qin, Xiaolin; Chen, Jingwei; Zhang, Jingzhong
2010-01-01
We present a complete algorithm for finding an exact minimal polynomial from its approximate value by using an improved parameterized integer relation construction method. Our result is superior to the existence of error controlling on obtaining an exact rational number from its approximation. The algorithm is applicable for finding exact minimal polynomial of an algebraic number by its approximate root. This also enables us to provide an efficient method of converting the rational approximation representation to the minimal polynomial representation, and devise a simple algorithm to factor multivariate polynomials with rational coefficients. Compared with the subsistent methods, our method combines advantage of high efficiency in numerical computation, and exact, stable results in symbolic computation. we also discuss some applications to some transcendental numbers by approximations. Moreover, the Digits of our algorithm is far less than the LLL-lattice basis reduction technique in theory. In this paper, we...
Quantization of gauge fields, graph polynomials and graph homology
Energy Technology Data Exchange (ETDEWEB)
Kreimer, Dirk, E-mail: kreimer@physik.hu-berlin.de [Humboldt University, 10099 Berlin (Germany); Sars, Matthias [Humboldt University, 10099 Berlin (Germany); Suijlekom, Walter D. van [Radboud University Nijmegen, 6525 AJ Nijmegen (Netherlands)
2013-09-15
We review quantization of gauge fields using algebraic properties of 3-regular graphs. We derive the Feynman integrand at n loops for a non-abelian gauge theory quantized in a covariant gauge from scalar integrands for connected 3-regular graphs, obtained from the two Symanzik polynomials. The transition to the full gauge theory amplitude is obtained by the use of a third, new, graph polynomial, the corolla polynomial. This implies effectively a covariant quantization without ghosts, where all the relevant signs of the ghost sector are incorporated in a double complex furnished by the corolla polynomial–we call it cycle homology–and by graph homology. -- Highlights: •We derive gauge theory Feynman from scalar field theory with 3-valent vertices. •We clarify the role of graph homology and cycle homology. •We use parametric renormalization and the new corolla polynomial.
Symmetric polynomials in information theory: Entropy and subentropy
Energy Technology Data Exchange (ETDEWEB)
Jozsa, Richard; Mitchison, Graeme [DAMTP, Centre for Mathematical Sciences, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA (United Kingdom)
2015-06-15
Entropy and other fundamental quantities of information theory are customarily expressed and manipulated as functions of probabilities. Here we study the entropy H and subentropy Q as functions of the elementary symmetric polynomials in the probabilities and reveal a series of remarkable properties. Derivatives of all orders are shown to satisfy a complete monotonicity property. H and Q themselves become multivariate Bernstein functions and we derive the density functions of their Levy-Khintchine representations. We also show that H and Q are Pick functions in each symmetric polynomial variable separately. Furthermore, we see that H and the intrinsically quantum informational quantity Q become surprisingly closely related in functional form, suggesting a special significance for the symmetric polynomials in quantum information theory. Using the symmetric polynomials, we also derive a series of further properties of H and Q.
Fast computation of interlace polynomials on graphs of bounded treewidth
Bläser, Markus
2009-01-01
We consider the multivariate interlace polynomial introduced by Courcelle (2008), which generalizes several interlace polynomials defined by Arratia, Bollobas and Sorkin (2004) and by Aigner and van der Holst (2004). We present an algorithm to compute the multivariate interlace polynomial of a graph with n vertices given a tree decomposition of the graph of width k. Our algorithm uses 2^{4.5k^2+O(k)}*n arithmetic operations and can be efficiently implemented in parallel. It is tailor-made for the interlace polynomial and uses linear algebra arguments concerning adjacency matrices of graphs. The best previously known result (Courcelle 2008) employs a general logical framework and leads to an algorithm with running time f(k)*n^4, where f(k) is doubly exponential in k.
Cubic Polynomials with Rational Roots and Critical Points
Gupta, Shiv K.; Szymanski, Waclaw
2010-01-01
If you want your students to graph a cubic polynomial, it is best to give them one with rational roots and critical points. In this paper, we describe completely all such cubics and explain how to generate them.
Discrete Darboux transformation for discrete polynomials of hypergeometric type
Bangerezako, Gaspard
1998-03-01
The Darboux transformation, well known in second-order differential operator theory, is applied to the difference equations satisfied by the discrete hypergeometric polynomials (Charlier, Meixner-Kravchuk, Hahn).
An Elementary Proof of the Polynomial Matrix Spectral Factorization Theorem
Ephremidze, Lasha
2010-01-01
A very simple and short proof of the polynomial matrix spectral factorization theorem (on the unit circle as well as on the real line) is presented, which relies on elementary complex analysis and linear algebra.
An operator approach to the Al-Salam-Carlitz polynomials
Chen, William Y. C.; Saad, Husam L.; Sun, Lisa H.
2010-04-01
We present an operator approach to Rogers-type formulas and Mehler's formula for the Al-Salam-Carlitz polynomials Un(x,y,a;q). By using the q-exponential operator, we obtain a Rogers-type formula, which leads to a linearization formula. With the aid of a bivariate augmentation operator, we get a simple derivation of Mehler's formula due to Al-Salam and Carlitz ["Some orthogonal q-polynomials," Math. Nachr. 30, 47 (1965)]. By means of the Cauchy companion augmentation operator, we obtain an equivalent form of Mehler's formula. We also give several identities on the generating functions for products of the Al-Salam-Carlitz polynomials, which are extensions of the formulas for the Rogers-Szegö polynomials.
SOME EXTREMAL PROPERTIES OF THE INTEGRAL OF LEGENDRE POLYNOMIALS
Institute of Scientific and Technical Information of China (English)
史应光; 王子玉
2001-01-01
Some extremal properties of the integral of Legendre polynomials are given, which are of independent interest. Meanwhile they show that a conjecture of P. Erdos[1] is plausible and maybe provides some means to prove this conjecture.
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].
The Differential Dimension Polynomial for Characterizable Differential Ideals
Lange-Hegermann, Markus
2014-01-01
We generalize the differential dimension polynomial from prime differential ideals to characterizable differential ideals. Its computation is algorithmic, its degree and leading coefficient remain differential birational invariants, and it decides equality of characterizable differential ideals contained in each other.
Using Tutte polynomials to analyze the structure of the benzodiazepines
Cadavid Muñoz, Juan José
2014-05-01
Graph theory in general and Tutte polynomials in particular, are implemented for analyzing the chemical structure of the benzodiazepines. Similarity analysis are used with the Tutte polynomials for finding other molecules that are similar to the benzodiazepines and therefore that might show similar psycho-active actions for medical purpose, in order to evade the drawbacks associated to the benzodiazepines based medicine. For each type of benzodiazepines, Tutte polynomials are computed and some numeric characteristics are obtained, such as the number of spanning trees and the number of spanning forests. Computations are done using the computer algebra Maple's GraphTheory package. The obtained analytical results are of great importance in pharmaceutical engineering. As a future research line, the usage of the chemistry computational program named Spartan, will be used to extent and compare it with the obtained results from the Tutte polynomials of benzodiazepines.
Raising and Lowering Operators for Askey-Wilson Polynomials
Directory of Open Access Journals (Sweden)
Siddhartha Sahi
2007-01-01
Full Text Available In this paper we describe two pairs of raising/lowering operators for Askey-Wilson polynomials, which result from constructions involving very different techniques. The first technique is quite elementary, and depends only on the ''classical'' properties of these polynomials, viz. the q-difference equation and the three term recurrence. The second technique is less elementary, and involves the one-variable version of the double affine Hecke algebra.
Fast Parallel Computation of Polynomials Using Few Processors
DEFF Research Database (Denmark)
Valiant, Leslie G.; Skyum, Sven; Berkowitz, S.;
1983-01-01
It is shown that any multivariate polynomial of degree $d$ that can be computed sequentially in $C$ steps can be computed in parallel in $O((\\log d)(\\log C + \\log d))$ steps using only $(Cd)^{O(1)} $ processors.......It is shown that any multivariate polynomial of degree $d$ that can be computed sequentially in $C$ steps can be computed in parallel in $O((\\log d)(\\log C + \\log d))$ steps using only $(Cd)^{O(1)} $ processors....
Fast parallel computation of polynomials using few processors
DEFF Research Database (Denmark)
Valiant, Leslie; Skyum, Sven
1981-01-01
It is shown that any multivariate polynomial that can be computed sequentially in C steps and has degree d can be computed in parallel in 0((log d) (log C + log d)) steps using only (Cd)0(1) processors.......It is shown that any multivariate polynomial that can be computed sequentially in C steps and has degree d can be computed in parallel in 0((log d) (log C + log d)) steps using only (Cd)0(1) processors....
Nonstandard decision methods for the solvability of real polynomial equations
Institute of Scientific and Technical Information of China (English)
曾广兴
1999-01-01
For a multivariate polynomial equation with coefficients in a computable ordered field, two criteria of this equation having real solutions are given. Based on the criteria, decision methods for the existence of real zeros and the semidefiniteness of binary polynomials are provided. With the aid of computers, these methods are used to solve several examples. The technique is to extend the original field involved in the question to a computable non-Archimedean ordered field containing infinitesimal elements.
First extension groups of Verma modules and $R$-polynomials
Abe, Noriyuki
2010-01-01
We study the first extension groups between Verma modules. There was a conjecture which claims that the dimensions of the higher extension groups between Verma modules are the coefficients of $R$-polynomials defined by Kazhdan-Lusztig. This conjecture was known as the Gabber-Joseph conjecture (although Gebber and Joseph did not state.) However, Boe gives a counterexample to this conjecture. In this paper, we study how far are the dimensions of extension groups from the coefficients of $R$-polynomials.
Multivariate Permutation Polynomial Systems and Nonlinear Pseudorandom Number Generators
Ostafe, Alina
2009-01-01
In this paper we study a class of dynamical systems generated by iterations of multivariate permutation polynomial systems which lead to polynomial growth of the degrees of these iterations. Using these estimates and the same techniques studied previously for inversive generators, we bound exponential sums along the orbits of these dynamical systems and show that they admit much stronger estimates on average over all initial values than in the general case and thus can be of use for pseudorandom number generation.