An algebraic fractional order differentiator for a class of signals satisfying a linear differential equation

KAUST Repository

Liu, Da-Yan; Tian, Yang; Boutat, Driss; Laleg-Kirati, Taous-Meriem

2015-01-01

This paper aims at designing a digital fractional order differentiator for a class of signals satisfying a linear differential equation to estimate fractional derivatives with an arbitrary order in noisy case, where the input can be unknown or known with noises. Firstly, an integer order differentiator for the input is constructed using a truncated Jacobi orthogonal series expansion. Then, a new algebraic formula for the Riemann-Liouville derivative is derived, which is enlightened by the algebraic parametric method. Secondly, a digital fractional order differentiator is proposed using a numerical integration method in discrete noisy case. Then, the noise error contribution is analyzed, where an error bound useful for the selection of the design parameter is provided. Finally, numerical examples illustrate the accuracy and the robustness of the proposed fractional order differentiator.

An algebraic fractional order differentiator for a class of signals satisfying a linear differential equation

KAUST Repository

Liu, Da-Yan

2015-04-30

This paper aims at designing a digital fractional order differentiator for a class of signals satisfying a linear differential equation to estimate fractional derivatives with an arbitrary order in noisy case, where the input can be unknown or known with noises. Firstly, an integer order differentiator for the input is constructed using a truncated Jacobi orthogonal series expansion. Then, a new algebraic formula for the Riemann-Liouville derivative is derived, which is enlightened by the algebraic parametric method. Secondly, a digital fractional order differentiator is proposed using a numerical integration method in discrete noisy case. Then, the noise error contribution is analyzed, where an error bound useful for the selection of the design parameter is provided. Finally, numerical examples illustrate the accuracy and the robustness of the proposed fractional order differentiator.

Sensitivity analysis and design optimization through automatic differentiation

International Nuclear Information System (INIS)

Hovland, Paul D; Norris, Boyana; Strout, Michelle Mills; Bhowmick, Sanjukta; Utke, Jean

2005-01-01

Automatic differentiation is a technique for transforming a program or subprogram that computes a function, including arbitrarily complex simulation codes, into one that computes the derivatives of that function. We describe the implementation and application of automatic differentiation tools. We highlight recent advances in the combinatorial algorithms and compiler technology that underlie successful implementation of automatic differentiation tools. We discuss applications of automatic differentiation in design optimization and sensitivity analysis. We also describe ongoing research in the design of language-independent source transformation infrastructures for automatic differentiation algorithms

Automatic differentiation algorithms in model analysis

NARCIS (Netherlands)

Huiskes, M.J.

2002-01-01

Title: Automatic differentiation algorithms in model analysis
Author: M.J. Huiskes
Date: 19 March, 2002

In this thesis automatic differentiation algorithms and derivative-based methods

Design of high-order all-optical temporal differentiators based on multiple-phase-shifted fiber Bragg gratings.

Science.gov (United States)

Kulishov, Mykola; Azaña, José

2007-05-14

A simple and general approach for designing practical all-optical (all-fiber) arbitrary-order time differentiators is introduced here for the first time. Specifically, we demonstrate that the Nth time derivative of an input optical waveform can be obtained by reflection of this waveform in a single uniform fiber Bragg grating (FBG) incorporating N &pi-phase shifts properly located along its grating profile. The general design procedure of an arbitrary-order optical time differentiator based on a multiple-phase-shifted FBG is described and numerically demonstrated for up to fourth-order time differentiation. Our simulations show that the proposed approach can provide optical operation bandwidths in the tens-of-GHz regime using readily feasible FBG structures.

Low rank approach to computing first and higher order derivatives using automatic differentiation

International Nuclear Information System (INIS)

Reed, J. A.; Abdel-Khalik, H. S.; Utke, J.

2012-01-01

This manuscript outlines a new approach for increasing the efficiency of applying automatic differentiation (AD) to large scale computational models. By using the principles of the Efficient Subspace Method (ESM), low rank approximations of the derivatives for first and higher orders can be calculated using minimized computational resources. The output obtained from nuclear reactor calculations typically has a much smaller numerical rank compared to the number of inputs and outputs. This rank deficiency can be exploited to reduce the number of derivatives that need to be calculated using AD. The effective rank can be determined according to ESM by computing derivatives with AD at random inputs. Reduced or pseudo variables are then defined and new derivatives are calculated with respect to the pseudo variables. Two different AD packages are used: OpenAD and Rapsodia. OpenAD is used to determine the effective rank and the subspace that contains the derivatives. Rapsodia is then used to calculate derivatives with respect to the pseudo variables for the desired order. The overall approach is applied to two simple problems and to MATWS, a safety code for sodium cooled reactors. (authors)

The arbitrary order design code Tlie 1.0

International Nuclear Information System (INIS)

Zeijts, J. van; Neri, Filippo

1993-01-01

We describe the arbitrary order charged particle transfer map code TLIE. This code is a general 6D relativistic design code with a MAD compatible input language and among others implements user defined functions and subroutines and nested fitting and optimization. First we describe the mathematics and physics in the code. Aside from generating maps for all the standard accelerator elements we describe an efficient method for generating nonlinear transfer maps for realistic magnet models. We have implemented the method to arbitrary order in our accelerator design code for cylindrical current sheet magnets. We also have implemented a self-consistent space-charge approach as in CHARLIE. Subsequently we give a description of the input language and finally, we give several examples from productions run, such as cases with stacked multipoles with overlapping fringe fields. (Author)

The differential equation of an arbitrary reflecting surface

Science.gov (United States)

Melka, Richard F.; Berrettini, Vincent D.; Yousif, Hashim A.

2018-05-01

A differential equation describing the reflection of a light ray incident upon an arbitrary reflecting surface is obtained using the law of reflection. The derived equation is written in terms of a parameter and the value of this parameter determines the nature of the reflecting surface. Under various parametric constraints, the solution of the differential equation leads to the various conic surfaces but is not generally solvable. In addition, the dynamics of the light reflections from the conic surfaces are executed in the Mathematica software. Our derivation is the converse of the traditional approach and our analysis assumes a relation between the object distance and the image distance. This leads to the differential equation of the reflecting surface.

Post-convergence automatic differentiation of iterative schemes

International Nuclear Information System (INIS)

Azmy, Y.Y.

1997-01-01

A new approach for performing automatic differentiation (AD) of computer codes that embody an iterative procedure, based on differentiating a single additional iteration upon achieving convergence, is described and implemented. This post-convergence automatic differentiation (PAD) technique results in better accuracy of the computed derivatives, as it eliminates part of the derivatives convergence error, and a large reduction in execution time, especially when many iterations are required to achieve convergence. In addition, it provides a way to compute derivatives of the converged solution without having to repeat the entire iterative process every time new parameters are considered. These advantages are demonstrated and the PAD technique is validated via a set of three linear and nonlinear codes used to solve neutron transport and fluid flow problems. The PAD technique reduces the execution time over direct AD by a factor of up to 30 and improves the accuracy of the derivatives by up to two orders of magnitude. The PAD technique's biggest disadvantage lies in the necessity to compute the iterative map's Jacobian, which for large problems can be prohibitive. Methods are discussed to alleviate this difficulty

ON PARTIAL DIFFERENTIAL AND DIFFERENCE EQUATIONS WITH SYMMETRIES DEPENDING ON ARBITRARY FUNCTIONS

Directory of Open Access Journals (Sweden)

Giorgio Gubbiotti

2016-06-01

Full Text Available In this note we present some ideas on when Lie symmetries, both point and generalized, can depend on arbitrary functions. We show a few examples, both in partial differential and partial difference equations where this happens. Moreover we show that the infinitesimal generators of generalized symmetries depending on arbitrary functions, both for continuous and discrete equations, effectively play the role of master symmetries.

FORSIM-6, Automatic Solution of Coupled Differential Equation System

International Nuclear Information System (INIS)

Carver, M.B.; Stewart, D.G.; Blair, J.M.; Selander, W.N.

1983-01-01

1 - Description of problem or function: The FORSIM program is a versatile package which automates the solution of coupled differential equation systems. The independent variables are time, and up to three space coordinates, and the equations may be any mixture of partial and/or ordinary differential equations. The philosophy of the program is to provide a tool which will solve a system of differential equations for a user who has basic but unspecialized knowledge of numerical analysis and FORTRAN. The equations to be solved, together with the initial conditions and any special instructions, may be specified by the user in a single FORTRAN subroutine, although he may write a number of routines if this is more suitable. These are then loaded with the control routines, which perform the solution and any requested input and output. 2 - Method of solution: Partial differential equations are automatically converted into sets of coupled ordinary differential equations by variable order discretization in the spatial dimensions. These and other ordinary differential equations are integrated continuously in time using efficient variable order, variable step, error-controlled algorithms

Automatic Clustering Using FSDE-Forced Strategy Differential Evolution

Science.gov (United States)

Yasid, A.

2018-01-01

Clustering analysis is important in datamining for unsupervised data, cause no adequate prior knowledge. One of the important tasks is defining the number of clusters without user involvement that is known as automatic clustering. This study intends on acquiring cluster number automatically utilizing forced strategy differential evolution (AC-FSDE). Two mutation parameters, namely: constant parameter and variable parameter are employed to boost differential evolution performance. Four well-known benchmark datasets were used to evaluate the algorithm. Moreover, the result is compared with other state of the art automatic clustering methods. The experiment results evidence that AC-FSDE is better or competitive with other existing automatic clustering algorithm.

Operator overloading as an enabling technology for automatic differentiation

International Nuclear Information System (INIS)

Corliss, G.F.; Griewank, A.

1993-01-01

We present an example of the science that is enabled by object-oriented programming techniques. Scientific computation often needs derivatives for solving nonlinear systems such as those arising in many PDE algorithms, optimization, parameter identification, stiff ordinary differential equations, or sensitivity analysis. Automatic differentiation computes derivatives accurately and efficiently by applying the chain rule to each arithmetic operation or elementary function. Operator overloading enables the techniques of either the forward or the reverse mode of automatic differentiation to be applied to real-world scientific problems. We illustrate automatic differentiation with an example drawn from a model of unsaturated flow in a porous medium. The problem arises from planning for the long-term storage of radioactive waste

Planar real polynomial differential systems of degree n > 3 having a weak focus of high order

International Nuclear Information System (INIS)

Llibre, J.; Rabanal, R.

2008-06-01

We construct planar polynomial differential systems of even (respectively odd) degree n > 3, of the form linear plus a nonlinear homogeneous part of degree n having a weak focus of order n 2 -1 (respectively (n 2 -1)/2 ) at the origin. As far as we know this provides the highest order known until now for a weak focus of a polynomial differential system of arbitrary degree n. (author)

Differential equation for genus-two characters in arbitrary rational conformal field theories

International Nuclear Information System (INIS)

Mathur, S.D.; Sen, A.

1989-01-01

We develop a general method for deriving ordinary differential equations for the genus-two ''characters'' of an arbitrary rational conformal field theory using the hyperelliptic representation of the genus-two moduli space. We illustrate our method by explicitly deriving the character differential equations for k=1 SU(2), G 2 , and F 4 WZW models. Our method provides an intrinsic definition of conformal field theories on higher genus Riemann surfaces. (orig.)

Nonlinear differential equations for the wavefront surface at arbitrary Hartmann-plane distances.

Science.gov (United States)

Téllez-Quiñones, Alejandro; Malacara-Doblado, Daniel; Flores-Hernández, Ricardo; Gutiérrez-Hernández, David A; León-Rodríguez, Miguel

2016-03-20

In the Hartmann test, a wave aberration function W is estimated from the information of the spot diagram drawn in an observation plane. The distance from a reference plane to the observation plane, the Hartmann-plane distance, is typically chosen as z=f, where f is the radius of a reference sphere. The function W and the transversal aberrations {X,Y} calculated at the plane z=f are related by two well-known linear differential equations. Here, we propose two nonlinear differential equations to denote a more general relation between W and the transversal aberrations {U,V} calculated at any arbitrary Hartmann-plane distance z=r. We also show how to directly estimate the wavefront surface w from the information of {U,V}. The use of arbitrary r values could improve the reliability of the measurements of W, or w, when finding difficulties in adequate ray identification at z=f.

The orthogonal gradients method: A radial basis functions method for solving partial differential equations on arbitrary surfaces

KAUST Repository

Piret, Cé cile

2012-01-01

Much work has been done on reconstructing arbitrary surfaces using the radial basis function (RBF) method, but one can hardly find any work done on the use of RBFs to solve partial differential equations (PDEs) on arbitrary surfaces. In this paper

TMB: Automatic Differentiation and Laplace Approximation

Directory of Open Access Journals (Sweden)

Kasper Kristensen

2016-04-01

Full Text Available TMB is an open source R package that enables quick implementation of complex nonlinear random effects (latent variable models in a manner similar to the established AD Model Builder package (ADMB, http://admb-project.org/; Fournier et al. 2011. In addition, it offers easy access to parallel computations. The user defines the joint likelihood for the data and the random effects as a C++ template function, while all the other operations are done in R; e.g., reading in the data. The package evaluates and maximizes the Laplace approximation of the marginal likelihood where the random effects are automatically integrated out. This approximation, and its derivatives, are obtained using automatic differentiation (up to order three of the joint likelihood. The computations are designed to be fast for problems with many random effects (≈ 106 and parameters (≈ 103 . Computation times using ADMB and TMB are compared on a suite of examples ranging from simple models to large spatial models where the random effects are a Gaussian random field. Speedups ranging from 1.5 to about 100 are obtained with increasing gains for large problems. The package and examples are available at http://tmb-project.org/.

Multipole expansion of acoustical Bessel beams with arbitrary order and location.

Science.gov (United States)

Gong, Zhixiong; Marston, Philip L; Li, Wei; Chai, Yingbin

2017-06-01

An exact solution of expansion coefficients for a T-matrix method interacting with acoustic scattering of arbitrary order Bessel beams from an obstacle of arbitrary location is derived analytically. Because of the failure of the addition theorem for spherical harmonics for expansion coefficients of helicoidal Bessel beams, an addition theorem for cylindrical Bessel functions is introduced. Meanwhile, an analytical expression for the integral of products including Bessel and associated Legendre functions is applied to eliminate the integration over the polar angle. Note that this multipole expansion may also benefit other scattering methods and expansions of incident waves, for instance, partial-wave series solutions.

Photonic arbitrary waveform generator based on Taylor synthesis method

DEFF Research Database (Denmark)

Liao, Shasha; Ding, Yunhong; Dong, Jianji

2016-01-01

Arbitrary waveform generation has been widely used in optical communication, radar system and many other applications. We propose and experimentally demonstrate a silicon-on-insulator (SOI) on chip optical arbitrary waveform generator, which is based on Taylor synthesis method. In our scheme......, a Gaussian pulse is launched to some cascaded microrings to obtain first-, second- and third-order differentiations. By controlling amplitude and phase of the initial pulse and successive differentiations, we can realize an arbitrary waveform generator according to Taylor expansion. We obtain several typical...... waveforms such as square waveform, triangular waveform, flat-top waveform, sawtooth waveform, Gaussian waveform and so on. Unlike other schemes based on Fourier synthesis or frequency-to-time mapping, our scheme is based on Taylor synthesis method. Our scheme does not require any spectral disperser or large...

The orthogonal gradients method: A radial basis functions method for solving partial differential equations on arbitrary surfaces

KAUST Repository

Piret, Cécile

2012-05-01

Much work has been done on reconstructing arbitrary surfaces using the radial basis function (RBF) method, but one can hardly find any work done on the use of RBFs to solve partial differential equations (PDEs) on arbitrary surfaces. In this paper, we investigate methods to solve PDEs on arbitrary stationary surfaces embedded in . R3 using the RBF method. We present three RBF-based methods that easily discretize surface differential operators. We take advantage of the meshfree character of RBFs, which give us a high accuracy and the flexibility to represent the most complex geometries in any dimension. Two out of the three methods, which we call the orthogonal gradients (OGr) methods are the result of our work and are hereby presented for the first time. © 2012 Elsevier Inc.

Optimality Conditions in Differentiable Vector Optimization via Second-Order Tangent Sets

International Nuclear Information System (INIS)

Jimenez, Bienvenido; Novo, Vicente

2004-01-01

We provide second-order necessary and sufficient conditions for a point to be an efficient element of a set with respect to a cone in a normed space, so that there is only a small gap between necessary and sufficient conditions. To this aim, we use the common second-order tangent set and the asymptotic second-order cone utilized by Penot. As an application we establish second-order necessary conditions for a point to be a solution of a vector optimization problem with an arbitrary feasible set and a twice Frechet differentiable objective function between two normed spaces. We also establish second-order sufficient conditions when the initial space is finite-dimensional so that there is no gap with necessary conditions. Lagrange multiplier rules are also given

Tangent: Automatic Differentiation Using Source Code Transformation in Python

OpenAIRE

van Merriënboer, Bart; Wiltschko, Alexander B.; Moldovan, Dan

2017-01-01

Automatic differentiation (AD) is an essential primitive for machine learning programming systems. Tangent is a new library that performs AD using source code transformation (SCT) in Python. It takes numeric functions written in a syntactic subset of Python and NumPy as input, and generates new Python functions which calculate a derivative. This approach to automatic differentiation is different from existing packages popular in machine learning, such as TensorFlow and Autograd. Advantages ar...

Automatic differentiation tools in the dynamic simulation of chemical engineering processes

Directory of Open Access Journals (Sweden)

Castro M.C.

2000-01-01

Full Text Available Automatic Differentiation is a relatively recent technique developed for the differentiation of functions applicable directly to the source code to compute the function written in standard programming languages. That technique permits the automatization of the differentiation step, crucial for dynamic simulation and optimization of processes. The values for the derivatives obtained with AD are exact (to roundoff. The theoretical exactness of the AD comes from the fact that it uses the same rules of differentiation as in differential calculus, but these rules are applied to an algorithmic specification of the function rather than to a formula. The main purpose of this contribution is to discuss the impact of Automatic Differentiation in the field of dynamic simulation of chemical engineering processes. The influence of the differentiation technique on the behavior of the integration code, the performance of the generated code and the incorporation of AD tools in consistent initialization tools are discussed from the viewpoint of dynamic simulation of typical models in chemical engineering.

On existence of soliton solutions of arbitrary-order system of nonlinear Schrodinger equations

International Nuclear Information System (INIS)

Zhestkov, S.V.

2003-01-01

The soliton solutions are constructed for the system of arbitrary-order coupled nonlinear Schrodinger equations . The necessary and sufficient conditions of existence of these solutions are obtained. It is shown that the maximum number of solitons in nondegenerate case is 4L, where L is order of the system. (author)

The Mehler-Fock transform of general order and arbitrary index and its inversion

Directory of Open Access Journals (Sweden)

Cyril Nasim

1984-01-01

Full Text Available An integral transform involving the associated Legendre function of zero order, P−12+iτ(x, x∈[1,∞, as the kernel (considered as a function of τ, is called Mehler-Fock transform. Some generalizations, involving the function P−12+iτμ(x, where the order μ is an arbitrary complex number, including the case when μ=0,1,2,… have been known for some time. In this present note, we define a general Mehler-Fock transform involving, as the kernel, the Legendre function P−12+tμ(x, of general order μ and an arbitrary index −12+t, t=σ+iτ, −∞<τ<∞. Then we develop a symmetric inversion formulae for these transforms. Many well-known results are derived as special cases of this general form. These transforms are widely used for solving many axisymmetric potential problems.

Automatic differentiation for gradient-based optimization of radiatively heated microelectronics manufacturing equipment

Energy Technology Data Exchange (ETDEWEB)

Moen, C.D.; Spence, P.A.; Meza, J.C.; Plantenga, T.D.

1996-12-31

Automatic differentiation is applied to the optimal design of microelectronic manufacturing equipment. The performance of nonlinear, least-squares optimization methods is compared between numerical and analytical gradient approaches. The optimization calculations are performed by running large finite-element codes in an object-oriented optimization environment. The Adifor automatic differentiation tool is used to generate analytic derivatives for the finite-element codes. The performance results support previous observations that automatic differentiation becomes beneficial as the number of optimization parameters increases. The increase in speed, relative to numerical differences, has a limited value and results are reported for two different analysis codes.

Automatic Differentiation and its Program Realization

Czech Academy of Sciences Publication Activity Database

Hartman, J.; Lukšan, Ladislav; Zítko, J.

2009-01-01

Roč. 45, č. 5 (2009), s. 865-883 ISSN 0023-5954 R&D Projects: GA AV ČR IAA1030405 Institutional research plan: CEZ:AV0Z10300504 Keywords : automatic differentiation * modeling languages * systems of optimization Subject RIV: BA - General Mathematics Impact factor: 0.445, year: 2009 http://dml.cz/handle/10338.dmlcz/140037

Applications of automatic differentiation in topology optimization

DEFF Research Database (Denmark)

Nørgaard, Sebastian A.; Sagebaum, Max; Gauger, Nicolas R.

2017-01-01

The goal of this article is to demonstrate the applicability and to discuss the advantages and disadvantages of automatic differentiation in topology optimization. The technique makes it possible to wholly or partially automate the evaluation of derivatives for optimization problems and is demons...

Applications of automatic differentiation in topology optimization

DEFF Research Database (Denmark)

Nørgaard, Sebastian A.; Sagebaum, Max; Gauger, Nicolas R.

2017-01-01

and is demonstrated on two separate, previously published types of problems in topology optimization. Two separate software packages for automatic differentiation, CoDiPack and Tapenade are considered, and their performance and usability trade-offs are discussed and compared to a hand coded adjoint gradient...

New finite volume methods for approximating partial differential equations on arbitrary meshes

International Nuclear Information System (INIS)

Hermeline, F.

2008-12-01

This dissertation presents some new methods of finite volume type for approximating partial differential equations on arbitrary meshes. The main idea lies in solving twice the problem to be dealt with. One addresses the elliptic equations with variable (anisotropic, antisymmetric, discontinuous) coefficients, the parabolic linear or non linear equations (heat equation, radiative diffusion, magnetic diffusion with Hall effect), the wave type equations (Maxwell, acoustics), the elasticity and Stokes'equations. Numerous numerical experiments show the good behaviour of this type of method. (author)

Nonrelativistic equations of motion for particles with arbitrary spin

International Nuclear Information System (INIS)

Fushchich, V.I.; Nikitin, A.G.

1981-01-01

First- and second-order Galileo-invariant systems of differential equations which describe the motion of nonrelativistic particles of arbitrary spin are derived. The equations can be derived from a Lagrangian and describe the dipole, quadrupole, and spin-orbit interaction of the particles with an external field; these interactions have traditionally been regarded as purely relativistic effects. The problem of the motion of a nonrelativistic particle of arbitrary spin in a homogeneous magnetic field is solved exactly on the basis of the obtained equations. The generators of all classes of irreducible representations of the Galileo group are found

Fast and Simple Method for Evaluation of Polarization Correction to Propagation Constant of Arbitrary Order Guided Modes in Optical Fibers with Arbitrary Refractive Index Profile

Directory of Open Access Journals (Sweden)

Anton Bourdine

2015-01-01

Full Text Available This work presents fast and simple method for evaluation of polarization correction to scalar propagation constant of arbitrary order guided modes propagating over weakly guiding optical fibers. Proposed solution is based on earlier on developed modified Gaussian approximation extended for analysis of weakly guiding optical fibers with arbitrary refractive index profile in the core region bounded by single solid outer cladding. Some results are presented that illustrate the decreasing of computational error during the estimation of propagation constant when polarization corrections are taken into account. Analytical expressions for the first and second derivatives of polarization correction are derived and presented.

Some operational tools for solving fractional and higher integer order differential equations: A survey on their mutual relations

Science.gov (United States)

Kiryakova, Virginia S.

2012-11-01

The Laplace Transform (LT) serves as a basis of the Operational Calculus (OC), widely explored by engineers and applied scientists in solving mathematical models for their practical needs. This transform is closely related to the exponential and trigonometric functions (exp, cos, sin) and to the classical differentiation and integration operators, reducing them to simple algebraic operations. Thus, the classical LT and the OC give useful tool to handle differential equations and systems with constant coefficients. Several generalizations of the LT have been introduced to allow solving, in a similar way, of differential equations with variable coefficients and of higher integer orders, as well as of fractional (arbitrary non-integer) orders. Note that fractional order mathematical models are recently widely used to describe better various systems and phenomena of the real world. This paper surveys briefly some of our results on classes of such integral transforms, that can be obtained from the LT by means of "transmutations" which are operators of the generalized fractional calculus (GFC). On the list of these Laplace-type integral transforms, we consider the Borel-Dzrbashjan, Meijer, Krätzel, Obrechkoff, generalized Obrechkoff (multi-index Borel-Dzrbashjan) transforms, etc. All of them are G- and H-integral transforms of convolutional type, having as kernels Meijer's G- or Fox's H-functions. Besides, some special functions (also being G- and H-functions), among them - the generalized Bessel-type and Mittag-Leffler (M-L) type functions, are generating Gel'fond-Leontiev (G-L) operators of generalized differentiation and integration, which happen to be also operators of GFC. Our integral transforms have operational properties analogous to those of the LT - they do algebrize the G-L generalized integrations and differentiations, and thus can serve for solving wide classes of differential equations with variable coefficients of arbitrary, including non-integer order

Nonlocal symmetries of a class of scalar and coupled nonlinear ordinary differential equations of any order

International Nuclear Information System (INIS)

Pradeep, R Gladwin; Chandrasekar, V K; Senthilvelan, M; Lakshmanan, M

2011-01-01

In this paper, we devise a systematic procedure to obtain nonlocal symmetries of a class of scalar nonlinear ordinary differential equations (ODEs) of arbitrary order related to linear ODEs through nonlocal relations. The procedure makes use of the Lie point symmetries of the linear ODEs and the nonlocal connection to deduce the nonlocal symmetries of the corresponding nonlinear ODEs. Using these nonlocal symmetries, we obtain reduction transformations and reduced equations to specific examples. We find that the reduced equations can be explicitly integrated to deduce the general solutions for these cases. We also extend this procedure to coupled higher order nonlinear ODEs with specific reference to second-order nonlinear ODEs. (paper)

Application of Arbitrary-Order Hilbert Spectral Analysis to Passive Scalar Turbulence

International Nuclear Information System (INIS)

Huang, Y X; Lu, Z M; Liu, Y L; Schmitt, F G; Gagne, Y

2011-01-01

In previous work [Huang et al., PRE 82, 26319, 2010], we found that the passive scalar turbulence field maybe less intermittent than what we believed before. Here we apply the same method, namely arbitrary-order Hilbert spectral analysis, to a passive scalar (temperature) time series with a Taylor's microscale Reynolds number Re λ ≅ 3000. We find that with increasing Reynolds number, the discrepancy of scaling exponents between Hilbert ξ θ (q) and Kolmogorov-Obukhov-Corrsin (KOC) theory is increasing, and consequently the discrepancy between Hilbert and structure function could disappear at infinite Reynolds number.

High-Order Hyperbolic Residual-Distribution Schemes on Arbitrary Triangular Grids

Science.gov (United States)

Mazaheri, Alireza; Nishikawa, Hiroaki

2015-01-01

In this paper, we construct high-order hyperbolic residual-distribution schemes for general advection-diffusion problems on arbitrary triangular grids. We demonstrate that the second-order accuracy of the hyperbolic schemes can be greatly improved by requiring the scheme to preserve exact quadratic solutions. We also show that the improved second-order scheme can be easily extended to third-order by further requiring the exactness for cubic solutions. We construct these schemes based on the LDA and the SUPG methodology formulated in the framework of the residual-distribution method. For both second- and third-order-schemes, we construct a fully implicit solver by the exact residual Jacobian of the second-order scheme, and demonstrate rapid convergence of 10-15 iterations to reduce the residuals by 10 orders of magnitude. We demonstrate also that these schemes can be constructed based on a separate treatment of the advective and diffusive terms, which paves the way for the construction of hyperbolic residual-distribution schemes for the compressible Navier-Stokes equations. Numerical results show that these schemes produce exceptionally accurate and smooth solution gradients on highly skewed and anisotropic triangular grids, including curved boundary problems, using linear elements. We also present Fourier analysis performed on the constructed linear system and show that an under-relaxation parameter is needed for stabilization of Gauss-Seidel relaxation.

Proposal for arbitrary-order temporal integration of ultrafast optical signals using a single uniform-period fiber Bragg grating.

Science.gov (United States)

Asghari, Mohammad H; Azaña, José

2008-07-01

A simple and practical all-fiber design for implementing arbitrary-order temporal integration of ultrafast optical waveforms is proposed and numerically investigated. We demonstrate that an ultrafast photonics integrator of any desired integration order can be implemented using a uniform-period fiber Bragg grating (FBG) with a properly designed amplitude-only grating apodization profile. In particular, the grating coupling strength must vary according to the (N-1) power of the fiber distance for implementing an Nth-order photonics integrator (N=1,2,...). This approach requires the same level of practical difficulty for realizing any given integration order. The proposed integration devices operate over a limited time window, which is approximately fixed by the round-trip propagation time in the FBG. Ultrafast arbitrary-order all-optical integrators capable of accurate operation over nanosecond time windows can be implemented using readily feasible FBGs.

Exact Solutions for Certain Nonlinear Autonomous Ordinary Differential Equations of the Second Order and Families of Two-Dimensional Autonomous Systems

Directory of Open Access Journals (Sweden)

M. P. Markakis

2010-01-01

Full Text Available Certain nonlinear autonomous ordinary differential equations of the second order are reduced to Abel equations of the first kind ((Ab-1 equations. Based on the results of a previous work, concerning a closed-form solution of a general (Ab-1 equation, and introducing an arbitrary function, exact one-parameter families of solutions are derived for the original autonomous equations, for the most of which only first integrals (in closed or parametric form have been obtained so far. Two-dimensional autonomous systems of differential equations of the first order, equivalent to the considered herein autonomous forms, are constructed and solved by means of the developed analysis.

Multiple scattering of a zero-order Bessel beam with arbitrary incidence by an aggregate of uniaxial anisotropic spheres

International Nuclear Information System (INIS)

Li, Z.J.; Wu, Z.S.; Qu, T.; Shang, Q.C.; Bai, L.

2016-01-01

Based on the generalized multiparticle Mie theory, multiple scattering of an aggregate of uniaxial anisotropic spheres illuminated by a zero-order Bessel beam (ZOBB) with arbitrary propagation direction is investigated. The particle size and configuration are arbitrary. The arbitrary incident Bessel beam is expanded in terms of spherical vector wave functions (SVWFs). Utilizing the vector addition theorem of SVWFs, interactive and total scattering coefficients are derived through the continuous boundary conditions on which the interaction of the particles is considered. The accuracy of the theory and codes are verified by comparing results with those obtained for arbitrary plane wave incidence by CST simulation, and for ZOBB incidence by a numerical method. The effects of angle of incidence, pseudo-polarization angle, half-conical angle, beam center position, and permittivity tensor elements on the radar cross sections (RCSs) of several types of collective uniaxial anisotropic spheres, such as a linear chain, a 4×4×4 cube-shaped array, and other periodical structures consisting of massive spheres, are numerically analyzed. Selected results on the properties of typical particles such as TiO 2 , SiO 2 , or other particle lattices are calculated. This investigation could provide an effective test for further research on the scattering characteristics of an aggregate of anisotropic spheres by a high-order Bessel vortex beam. The results have important application in optical tweezers and particle manipulation. - Highlights: • Scattering of Bessel beam by an aggregate of uniaxial anisotropic spheres is studied. • The zero-order Bessel beam propagates and polarizes along arbitrary direction. • The accuracy of expansion coefficients, the scattering theory and codes is verified. • Effects of various parameters on scattering properties are numerically discussed. • Scattering properties of several type of periodical array are numerically analyzed.

Automatic Differentiation and Deep Learning

CERN Multimedia

CERN. Geneva

2018-01-01

Statistical learning has been getting more and more interest from the particle-physics community in recent times, with neural networks and gradient-based optimization being a focus. In this talk we shall discuss three things: automatic differention tools: tools to quickly build DAGs of computation that are fully differentiable. We shall focus on one such tool "PyTorch".  Easy deployment of trained neural networks into large systems with many constraints: for example, deploying a model at the reconstruction phase where the neural network has to be integrated into CERN's bulk data-processing C++-only environment Some recent models in deep learning for segmentation and generation that might be useful for particle physics problems.

LTSN solution of the adjoint neutron transport equation with arbitrary source for high order of quadrature in a homogeneous slab

International Nuclear Information System (INIS)

Goncalves, Glenio A.; Orengo, Gilberto; Vilhena, Marco Tullio M.B. de; Graca, Claudio O.

2002-01-01

In this work we present the LTS N solution of the adjoint transport equation for an arbitrary source, testing the aptness of this analytical solution for high order of quadrature in transport problems and comparing some preliminary results with the ANISN computations in a homogeneneous slab geometry. In order to do that we apply the new formulation for the LTS N method based on the invariance projection property, becoming possible to handle problems with arbitrary sources and demanding high order of quadrature or deep penetration. This new approach for the LTS N method is important both for direct and adjoint transport calculations and its development was inspired by the necessity of using generalized adjoint sources for important calculations. Although the mathematical convergence has been proved for an arbitrary source, when the quadrature order or deep penetration is required the LTS N method presents computational overflow even for simple sources (sin, cos, exp, polynomial). With the new formulation we eliminate this drawback and in this work we report the numerical simulations testing the new approach

Transformation between surface spherical harmonic expansion of arbitrary high degree and order and double Fourier series on sphere

Science.gov (United States)

Fukushima, Toshio

2018-02-01

In order to accelerate the spherical harmonic synthesis and/or analysis of arbitrary function on the unit sphere, we developed a pair of procedures to transform between a truncated spherical harmonic expansion and the corresponding two-dimensional Fourier series. First, we obtained an analytic expression of the sine/cosine series coefficient of the 4 π fully normalized associated Legendre function in terms of the rectangle values of the Wigner d function. Then, we elaborated the existing method to transform the coefficients of the surface spherical harmonic expansion to those of the double Fourier series so as to be capable with arbitrary high degree and order. Next, we created a new method to transform inversely a given double Fourier series to the corresponding surface spherical harmonic expansion. The key of the new method is a couple of new recurrence formulas to compute the inverse transformation coefficients: a decreasing-order, fixed-degree, and fixed-wavenumber three-term formula for general terms, and an increasing-degree-and-order and fixed-wavenumber two-term formula for diagonal terms. Meanwhile, the two seed values are analytically prepared. Both of the forward and inverse transformation procedures are confirmed to be sufficiently accurate and applicable to an extremely high degree/order/wavenumber as 2^{30} {≈ } 10^9. The developed procedures will be useful not only in the synthesis and analysis of the spherical harmonic expansion of arbitrary high degree and order, but also in the evaluation of the derivatives and integrals of the spherical harmonic expansion.

Schwarzian conditions for linear differential operators with selected differential Galois groups

International Nuclear Information System (INIS)

Abdelaziz, Y; Maillard, J-M

2017-01-01

We show that non-linear Schwarzian differential equations emerging from covariance symmetry conditions imposed on linear differential operators with hypergeometric function solutions can be generalized to arbitrary order linear differential operators with polynomial coefficients having selected differential Galois groups. For order three and order four linear differential operators we show that this pullback invariance up to conjugation eventually reduces to symmetric powers of an underlying order-two operator. We give, precisely, the conditions to have modular correspondences solutions for such Schwarzian differential equations, which was an open question in a previous paper. We analyze in detail a pullbacked hypergeometric example generalizing modular forms, that ushers a pullback invariance up to operator homomorphisms. We finally consider the more general problem of the equivalence of two different order-four linear differential Calabi–Yau operators up to pullbacks and conjugation, and clarify the cases where they have the same Yukawa couplings. (paper)

Schwarzian conditions for linear differential operators with selected differential Galois groups

Science.gov (United States)

Abdelaziz, Y.; Maillard, J.-M.

2017-11-01

We show that non-linear Schwarzian differential equations emerging from covariance symmetry conditions imposed on linear differential operators with hypergeometric function solutions can be generalized to arbitrary order linear differential operators with polynomial coefficients having selected differential Galois groups. For order three and order four linear differential operators we show that this pullback invariance up to conjugation eventually reduces to symmetric powers of an underlying order-two operator. We give, precisely, the conditions to have modular correspondences solutions for such Schwarzian differential equations, which was an open question in a previous paper. We analyze in detail a pullbacked hypergeometric example generalizing modular forms, that ushers a pullback invariance up to operator homomorphisms. We finally consider the more general problem of the equivalence of two different order-four linear differential Calabi-Yau operators up to pullbacks and conjugation, and clarify the cases where they have the same Yukawa couplings.

Finding Higher Order Differentials of MISTY1

Science.gov (United States)

Tsunoo, Yukiyasu; Saito, Teruo; Kawabata, Takeshi; Nakagawa, Hirokatsu

MISTY1 is a 64-bit block cipher that has provable security against differential and linear cryptanalysis. MISTY1 is one of the algorithms selected in the European NESSIE project, and it is recommended for Japanese e-Government ciphers by the CRYPTREC project. In this paper, we report on 12th order differentials in 3-round MISTY1 with FL functions and 44th order differentials in 4-round MISTY1 with FL functions both previously unknown. We also report that both data complexity and computational complexity of higher order differential attacks on 6-round MISTY1 with FL functions and 7-round MISTY1 with FL functions using the 46th order differential can be reduced to as much as 1/22 of the previous values by using multiple 44th order differentials simultaneously.

Carlson iterating rational approximation and performance analysis of fractional operator with arbitrary order

International Nuclear Information System (INIS)

He Qiu-Yan; Yuan Xiao; Yu Bo

2017-01-01

The performance analysis of the generalized Carlson iterating process, which can realize the rational approximation of fractional operator with arbitrary order, is presented in this paper. The reasons why the generalized Carlson iterating function possesses more excellent properties such as self-similarity and exponential symmetry are also explained. K-index, P-index, O-index, and complexity index are introduced to contribute to performance analysis. Considering nine different operational orders and choosing an appropriate rational initial impedance for a certain operational order, these rational approximation impedance functions calculated by the iterating function meet computational rationality, positive reality, and operational validity. Then they are capable of having the operational performance of fractional operators and being physical realization. The approximation performance of the impedance function to the ideal fractional operator and the circuit network complexity are also exhibited. (paper)

Automatic differentiation as a tool in engineering design

Science.gov (United States)

Barthelemy, Jean-Francois; Hall, Laura E.

1992-01-01

Automatic Differentiation (AD) is a tool that systematically implements the chain rule of differentiation to obtain the derivatives of functions calculated by computer programs. AD is assessed as a tool for engineering design. The forward and reverse modes of AD, their computing requirements, as well as approaches to implementing AD are discussed. The application of two different tools to two medium-size structural analysis problems to generate sensitivity information typically necessary in an optimization or design situation is also discussed. The observation is made that AD is to be preferred to finite differencing in most cases, as long as sufficient computer storage is available; in some instances, AD may be the alternative to consider in lieu of analytical sensitivity analysis.

Parallel computation of automatic differentiation applied to magnetic field calculations

International Nuclear Information System (INIS)

Hinkins, R.L.; Lawrence Berkeley Lab., CA

1994-09-01

The author presents a parallelization of an accelerator physics application to simulate magnetic field in three dimensions. The problem involves the evaluation of high order derivatives with respect to two variables of a multivariate function. Automatic differentiation software had been used with some success, but the computation time was prohibitive. The implementation runs on several platforms, including a network of workstations using PVM, a MasPar using MPFortran, and a CM-5 using CMFortran. A careful examination of the code led to several optimizations that improved its serial performance by a factor of 8.7. The parallelization produced further improvements, especially on the MasPar with a speedup factor of 620. As a result a problem that took six days on a SPARC 10/41 now runs in minutes on the MasPar, making it feasible for physicists at Lawrence Berkeley Laboratory to simulate larger magnets

Application of automatic change of interval to de Vogelaere's method of the solution of the differential equation y'' = f (x, y)

International Nuclear Information System (INIS)

Rogers, M.H.

1960-11-01

The paper gives an extension to de Vogelaere's method for the solution of systems of second order differential equations from which first derivatives are absent. The extension is a description of the way in which automatic change in step-length can be made to give a prescribed accuracy at each step. (author)

Skyrmion dynamics in single-hole Neel ordered doped two-dimensional antiferromagnets with arbitrary spin

International Nuclear Information System (INIS)

Moura, A.R.; Pereira, A.R.; Moura-Melo, W.A.; Pires, A.S.T.

2008-01-01

We develop an effective theory to study the skyrmion dynamics in the presence of a hole (removed spins from the lattice) in Neel ordered two-dimensional antiferromagnets with arbitrary spin value S. The general equation of motion for the 'mass center' of this structure is obtained. The frequency of small amplitude oscillations of pinned skyrmions around the defect center is calculated. It is proportional to the hole size and inversely proportional to the square of the skyrmion size

Numerically stable, scalable formulas for parallel and online computation of higher-order multivariate central moments with arbitrary weights

Energy Technology Data Exchange (ETDEWEB)

Pebay, Philippe [Sandia National Laboratories (SNL-CA), Livermore, CA (United States); Terriberry, Timothy B. [Xiph.Org Foundation, Arlington, VA (United States); Kolla, Hemanth [Sandia National Laboratories (SNL-CA), Livermore, CA (United States); Bennett, Janine [Sandia National Laboratories (SNL-CA), Livermore, CA (United States)

2016-03-29

Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Formulas such as these, are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearly the full representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.

Automatically ordering events and times in text

CERN Document Server

Derczynski, Leon R A

2017-01-01

The book offers a detailed guide to temporal ordering, exploring open problems in the field and providing solutions and extensive analysis. It addresses the challenge of automatically ordering events and times in text. Aided by TimeML, it also describes and presents concepts relating to time in easy-to-compute terms. Working out the order that events and times happen has proven difficult for computers, since the language used to discuss time can be vague and complex. Mapping out these concepts for a computational system, which does not have its own inherent idea of time, is, unsurprisingly, tough. Solving this problem enables powerful systems that can plan, reason about events, and construct stories of their own accord, as well as understand the complex narratives that humans express and comprehend so naturally. This book presents a theory and data-driven analysis of temporal ordering, leading to the identification of exactly what is difficult about the task. It then proposes and evaluates machine-learning so...

Automatic differentiation in geophysical inverse problems

Science.gov (United States)

Sambridge, M.; Rickwood, P.; Rawlinson, N.; Sommacal, S.

2007-07-01

Automatic differentiation (AD) is the technique whereby output variables of a computer code evaluating any complicated function (e.g. the solution to a differential equation) can be differentiated with respect to the input variables. Often AD tools take the form of source to source translators and produce computer code without the need for deriving and hand coding of explicit mathematical formulae by the user. The power of AD lies in the fact that it combines the generality of finite difference techniques and the accuracy and efficiency of analytical derivatives, while at the same time eliminating `human' coding errors. It also provides the possibility of accurate, efficient derivative calculation from complex `forward' codes where no analytical derivatives are possible and finite difference techniques are too cumbersome. AD is already having a major impact in areas such as optimization, meteorology and oceanography. Similarly it has considerable potential for use in non-linear inverse problems in geophysics where linearization is desirable, or for sensitivity analysis of large numerical simulation codes, for example, wave propagation and geodynamic modelling. At present, however, AD tools appear to be little used in the geosciences. Here we report on experiments using a state of the art AD tool to perform source to source code translation in a range of geoscience problems. These include calculating derivatives for Gibbs free energy minimization, seismic receiver function inversion, and seismic ray tracing. Issues of accuracy and efficiency are discussed.

AD Model Builder: using automatic differentiation for statistical inference of highly parameterized complex nonlinear models

DEFF Research Database (Denmark)

Fournier, David A.; Skaug, Hans J.; Ancheta, Johnoel

2011-01-01

Many criteria for statistical parameter estimation, such as maximum likelihood, are formulated as a nonlinear optimization problem.Automatic Differentiation Model Builder (ADMB) is a programming framework based on automatic differentiation, aimed at highly nonlinear models with a large number...... of such a feature is the generic implementation of Laplace approximation of high-dimensional integrals for use in latent variable models. We also review the literature in which ADMB has been used, and discuss future development of ADMB as an open source project. Overall, the main advantages ofADMB are flexibility...

On the asymptotic expansions of solutions of an nth order linear differential equation with power coefficients

International Nuclear Information System (INIS)

Paris, R.B.; Wood, A.D.

1984-11-01

The asymptotic expansions of solutions of a class of linear ordinary differential equations of arbitrary order n, containing a factor zsup(m) multiplying the lower order derivatives, are investigated for large values of z in the complex plane. Four classes of solutions are considered which exhibit the following behaviour as /z/ → infinity in certain sectors: (i) solutions whose behaviour is either exponentially large or algebraic (involving p ( < n) algebraic expansions), (ii) solutions which are exponentially small (iii) solutions with a single algebraic expansion and (iv) solutions which are even and odd functions of z whenever n+m is even. The asymptotic expansions of these solutions in a full neigbourhood of the point at infinity are obtained by means of the theory of the solutions in the case m=O developed in a previous paper

Study on time of flight property of electron optical systems by differential algebraic method

International Nuclear Information System (INIS)

Cheng Min; Tang Tiantong; Yao Zhenhua

2002-01-01

Differential algebraic method is a powerful and promising technique in computer numerical analysis. When applied to nonlinear dynamics systems, the arbitrary high-order transfer properties of the systems can be computed directly with high precision. In this paper, the principle of differential algebra is applied to study on the time of flight (TOF) property of electron optical systems and their arbitrary order TOF transfer properties can be numerically calculated out. As an example, TOF transfer properties of a uniform magnetic sector field analyzer have been studied by differential algebraic method. Relative errors of the first-order and second-order TOF transfer coefficients of the magnetic sector field analyzer are of the order 10 -11 or smaller compared with the analytic solutions. It is proved that differential algebraic TOF method is of high accuracy and very helpful for high-order TOF transfer property analysis of electron optical systems. (author)

An adaptive singular ES-FEM for mechanics problems with singular field of arbitrary order

OpenAIRE

Nguyen-Xuan, H.; Liu, G. R.; Bordas, Stéphane; Natarajan, S.; Rabczuk, T.

2013-01-01

This paper presents a singular edge-based smoothed finite element method (sES-FEM) for mechanics problems with singular stress fields of arbitrary order. The sES-FEM uses a basic mesh of three-noded linear triangular (T3) elements and a special layer of five-noded singular triangular elements (sT5) connected to the singular-point of the stress field. The sT5 element has an additional node on each of the two edges connected to the singular-point. It allows us to represent simple and efficient ...

Fractional order differentiation by integration with Jacobi polynomials

KAUST Repository

Liu, Dayan

2012-12-01

The differentiation by integration method with Jacobi polynomials was originally introduced by Mboup, Join and Fliess [22], [23]. This paper generalizes this method from the integer order to the fractional order for estimating the fractional order derivatives of noisy signals. The proposed fractional order differentiator is deduced from the Jacobi orthogonal polynomial filter and the Riemann-Liouville fractional order derivative definition. Exact and simple formula for this differentiator is given where an integral formula involving Jacobi polynomials and the noisy signal is used without complex mathematical deduction. Hence, it can be used both for continuous-time and discrete-time models. The comparison between our differentiator and the recently introduced digital fractional order Savitzky-Golay differentiator is given in numerical simulations so as to show its accuracy and robustness with respect to corrupting noises. © 2012 IEEE.

Fractional order differentiation by integration with Jacobi polynomials

KAUST Repository

Liu, Dayan; Gibaru, O.; Perruquetti, Wilfrid; Laleg-Kirati, Taous-Meriem

2012-01-01

The differentiation by integration method with Jacobi polynomials was originally introduced by Mboup, Join and Fliess [22], [23]. This paper generalizes this method from the integer order to the fractional order for estimating the fractional order derivatives of noisy signals. The proposed fractional order differentiator is deduced from the Jacobi orthogonal polynomial filter and the Riemann-Liouville fractional order derivative definition. Exact and simple formula for this differentiator is given where an integral formula involving Jacobi polynomials and the noisy signal is used without complex mathematical deduction. Hence, it can be used both for continuous-time and discrete-time models. The comparison between our differentiator and the recently introduced digital fractional order Savitzky-Golay differentiator is given in numerical simulations so as to show its accuracy and robustness with respect to corrupting noises. © 2012 IEEE.

Reducing the memory requirement in reverse mode automatic differentiation by solving TBR flow equations

International Nuclear Information System (INIS)

Naumann, U.

2002-01-01

The fast computation of gradients in reverse mode Automatic Differentiation (AD) requires the generation of adjoint versions of every statement in the original code. Due to the resulting reversal of the control flow certain intermediate values have to be made available in reverse order to compute the local partial derivatives. This can be achieved by storing these values or by recomputing them when they become required. In any case one is interested in minimizing the size of this set. Following an extensive introduction of the ''To-Be-Recorded'' (TBR) problem the authors present flow equations for propagating the TBR status of variables in the context of reverse mode AD of structured programs

Semiclassical approximation of the Wheeler-DeWitt equation: arbitrary orders and the question of unitarity

Science.gov (United States)

Kiefer, Claus; Wichmann, David

2018-06-01

We extend the Born-Oppenheimer type of approximation scheme for the Wheeler-DeWitt equation of canonical quantum gravity to arbitrary orders in the inverse Planck mass squared. We discuss in detail the origin of unitarity violation in this scheme and show that unitarity can be restored by an appropriate modification which requires back reaction from matter onto the gravitational sector. In our analysis, we heavily rely on the gauge aspects of the standard Born-Oppenheimer scheme in molecular physics.

Wavelet Methods for Solving Fractional Order Differential Equations

OpenAIRE

A. K. Gupta; S. Saha Ray

2014-01-01

Fractional calculus is a field of applied mathematics which deals with derivatives and integrals of arbitrary orders. The fractional calculus has gained considerable importance during the past decades mainly due to its application in diverse fields of science and engineering such as viscoelasticity, diffusion of biological population, signal processing, electromagnetism, fluid mechanics, electrochemistry, and many more. In this paper, we review different wavelet methods for solving both linea...

High Order Differential Frequency Hopping: Design and Analysis

Directory of Open Access Journals (Sweden)

Yong Li

2015-01-01

Full Text Available This paper considers spectrally efficient differential frequency hopping (DFH system design. Relying on time-frequency diversity over large spectrum and high speed frequency hopping, DFH systems are robust against hostile jamming interference. However, the spectral efficiency of conventional DFH systems is very low due to only using the frequency of each channel. To improve the system capacity, in this paper, we propose an innovative high order differential frequency hopping (HODFH scheme. Unlike in traditional DFH where the message is carried by the frequency relationship between the adjacent hops using one order differential coding, in HODFH, the message is carried by the frequency and phase relationship using two-order or higher order differential coding. As a result, system efficiency is increased significantly since the additional information transmission is achieved by the higher order differential coding at no extra cost on either bandwidth or power. Quantitative performance analysis on the proposed scheme demonstrates that transmission through the frequency and phase relationship using two-order or higher order differential coding essentially introduces another dimension to the signal space, and the corresponding coding gain can increase the system efficiency.

Adaptive finite element methods for differential equations

CERN Document Server

Bangerth, Wolfgang

2003-01-01

These Lecture Notes discuss concepts of `self-adaptivity' in the numerical solution of differential equations, with emphasis on Galerkin finite element methods. The key issues are a posteriori error estimation and it automatic mesh adaptation. Besides the traditional approach of energy-norm error control, a new duality-based technique, the Dual Weighted Residual method for goal-oriented error estimation, is discussed in detail. This method aims at economical computation of arbitrary quantities of physical interest by properly adapting the computational mesh. This is typically required in the design cycles of technical applications. For example, the drag coefficient of a body immersed in a viscous flow is computed, then it is minimized by varying certain control parameters, and finally the stability of the resulting flow is investigated by solving an eigenvalue problem. `Goal-oriented' adaptivity is designed to achieve these tasks with minimal cost. At the end of each chapter some exercises are posed in order ...

Field Method for Integrating the First Order Differential Equation

Institute of Scientific and Technical Information of China (English)

JIA Li-qun; ZHENG Shi-wang; ZHANG Yao-yu

2007-01-01

An important modern method in analytical mechanics for finding the integral, which is called the field-method, is used to research the solution of a differential equation of the first order. First, by introducing an intermediate variable, a more complicated differential equation of the first order can be expressed by two simple differential equations of the first order, then the field-method in analytical mechanics is introduced for solving the two differential equations of the first order. The conclusion shows that the field-method in analytical mechanics can be fully used to find the solutions of a differential equation of the first order, thus a new method for finding the solutions of the first order is provided.

The generation of arbitrary order, non-classical, Gauss-type quadrature for transport applications

International Nuclear Information System (INIS)

Spence, Peter J.

2015-01-01

A method is presented, based upon the Stieltjes method (1884), for the determination of non-classical Gauss-type quadrature rules, and the associated sets of abscissae and weights. The method is then used to generate a number of quadrature sets, to arbitrary order, which are primarily aimed at deterministic transport calculations. The quadrature rules and sets detailed include arbitrary order reproductions of those presented by Abu-Shumays in [4,8] (known as the QR sets, but labelled QRA here), in addition to a number of new rules and associated sets; these are generated in a similar way, and we label them the QRS quadrature sets. The method presented here shifts the inherent difficulty (encountered by Abu-Shumays) associated with solving the non-linear moment equations, particular to the required quadrature rule, to one of the determination of non-classical weight functions and the subsequent calculation of various associated inner products. Once a quadrature rule has been written in a standard form, with an associated weight function having been identified, the calculation of the required inner products is achieved using specific variable transformations, in addition to the use of rapid, highly accurate quadrature suited to this purpose. The associated non-classical Gauss quadrature sets can then be determined, and this can be done to any order very rapidly. In this paper, instead of listing weights and abscissae for the different quadrature sets detailed (of which there are a number), the MATLAB code written to generate them is included as Appendix D. The accuracy and efficacy (in a transport setting) of the quadrature sets presented is not tested in this paper (although the accuracy of the QRA quadrature sets has been studied in [12,13]), but comparisons to tabulated results listed in [8] are made. When comparisons are made with one of the azimuthal QRA sets detailed in [8], the inherent difficulty in the method of generation, used there, becomes apparent

The generation of arbitrary order, non-classical, Gauss-type quadrature for transport applications

Energy Technology Data Exchange (ETDEWEB)

Spence, Peter J., E-mail: peter.spence@awe.co.uk

2015-09-01

A method is presented, based upon the Stieltjes method (1884), for the determination of non-classical Gauss-type quadrature rules, and the associated sets of abscissae and weights. The method is then used to generate a number of quadrature sets, to arbitrary order, which are primarily aimed at deterministic transport calculations. The quadrature rules and sets detailed include arbitrary order reproductions of those presented by Abu-Shumays in [4,8] (known as the QR sets, but labelled QRA here), in addition to a number of new rules and associated sets; these are generated in a similar way, and we label them the QRS quadrature sets. The method presented here shifts the inherent difficulty (encountered by Abu-Shumays) associated with solving the non-linear moment equations, particular to the required quadrature rule, to one of the determination of non-classical weight functions and the subsequent calculation of various associated inner products. Once a quadrature rule has been written in a standard form, with an associated weight function having been identified, the calculation of the required inner products is achieved using specific variable transformations, in addition to the use of rapid, highly accurate quadrature suited to this purpose. The associated non-classical Gauss quadrature sets can then be determined, and this can be done to any order very rapidly. In this paper, instead of listing weights and abscissae for the different quadrature sets detailed (of which there are a number), the MATLAB code written to generate them is included as Appendix D. The accuracy and efficacy (in a transport setting) of the quadrature sets presented is not tested in this paper (although the accuracy of the QRA quadrature sets has been studied in [12,13]), but comparisons to tabulated results listed in [8] are made. When comparisons are made with one of the azimuthal QRA sets detailed in [8], the inherent difficulty in the method of generation, used there, becomes apparent

Fast RBF OGr for solving PDEs on arbitrary surfaces

Science.gov (United States)

Piret, Cécile; Dunn, Jarrett

2016-10-01

The Radial Basis Functions Orthogonal Gradients method (RBF-OGr) was introduced in [1] to discretize differential operators defined on arbitrary manifolds defined only by a point cloud. We take advantage of the meshfree character of RBFs, which give us a high accuracy and the flexibility to represent complex geometries in any spatial dimension. A large limitation of the RBF-OGr method was its large computational complexity, which greatly restricted the size of the point cloud. In this paper, we apply the RBF-Finite Difference (RBF-FD) technique to the RBF-OGr method for building sparse differentiation matrices discretizing continuous differential operators such as the Laplace-Beltrami operator. This method can be applied to solving PDEs on arbitrary surfaces embedded in ℛ3. We illustrate the accuracy of our new method by solving the heat equation on the unit sphere.

High-order asynchrony-tolerant finite difference schemes for partial differential equations

Science.gov (United States)

Aditya, Konduri; Donzis, Diego A.

2017-12-01

Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a method based on finite-difference schemes to solve partial differential equations in an asynchronous fashion - synchronization between PEs is relaxed at a mathematical level. While standard schemes can maintain their stability in the presence of asynchrony, their accuracy is drastically affected. In this work, we present a general methodology to derive asynchrony-tolerant (AT) finite difference schemes of arbitrary order of accuracy, which can maintain their accuracy when synchronizations are relaxed. We show that there are several choices available in selecting a stencil to derive these schemes and discuss their effect on numerical and computational performance. We provide a simple classification of schemes based on the stencil and derive schemes that are representative of different classes. Their numerical error is rigorously analyzed within a statistical framework to obtain the overall accuracy of the solution. Results from numerical experiments are used to validate the performance of the schemes.

Relativistic equation of the orbit of a particle in a arbitrary central force field

International Nuclear Information System (INIS)

Aaron, Francisc D.

2005-01-01

The equation of the orbit of a relativistic particle moving in an arbitrary central force field is derived. Straightforward generalizations of well-known first and second order differential equations are given. It is pointed out that the relativistic equation of the orbit has the same form as in the non-relativistic case, the only changes consisting in the appearance of additional terms proportional to 1/c 2 in both potential and total energies. (author)

On solutions of variable-order fractional differential equations

Directory of Open Access Journals (Sweden)

Ali Akgül

2017-01-01

solutions to fractional differential equations are compelling to get in real applications, due to the nonlocality and complexity of the fractional differential operators, especially for variable-order fractional differential equations. Therefore, it is significant to enhanced numerical methods for fractional differential equations. In this work, we consider variable-order fractional differential equations by reproducing kernel method. There has been much attention in the use of reproducing kernels for the solutions to many problems in the recent years. We give two examples to demonstrate how efficiently our theory can be implemented in practice.

Equation for disentangling time-ordered exponentials with arbitrary quadratic generators

International Nuclear Information System (INIS)

Budanov, V.G.

1987-01-01

In many quantum-mechanical constructions, it is necessary to disentangle an operator-valued time-ordered exponential with time-dependent generators quadratic in the creation and annihilation operators. By disentangling, one understands the finding of the matrix elements of the time-ordered exponential or, in a more general formulation. The solution of the problem can also be reduced to calculation of a matrix time-ordered exponential that solves the corresponding classical problem. However, in either case the evolution equations in their usual form do not enable one to take into account explicitly the symmetry of the system. In this paper the methods of Weyl analysis are used to find an ordinary differential equation on a matrix Lie algebra that is invariant with respect to the adjoint action of the dynamical symmetry group of a quadratic Hamiltonian and replaces the operator evolution equation for the Green's function

The arbitrary order mimetic finite difference method for a diffusion equation with a non-symmetric diffusion tensor

Science.gov (United States)

Gyrya, V.; Lipnikov, K.

2017-11-01

We present the arbitrary order mimetic finite difference (MFD) discretization for the diffusion equation with non-symmetric tensorial diffusion coefficient in a mixed formulation on general polygonal meshes. The diffusion tensor is assumed to be positive definite. The asymmetry of the diffusion tensor requires changes to the standard MFD construction. We present new approach for the construction that guarantees positive definiteness of the non-symmetric mass matrix in the space of discrete velocities. The numerically observed convergence rate for the scalar quantity matches the predicted one in the case of the lowest order mimetic scheme. For higher orders schemes, we observed super-convergence by one order for the scalar variable which is consistent with the previously published result for a symmetric diffusion tensor. The new scheme was also tested on a time-dependent problem modeling the Hall effect in the resistive magnetohydrodynamics.

Higgs Boson Production at Hadron Colliders: Differential Cross Section Through Next-to-Next-to-Leading Order

International Nuclear Information System (INIS)

Anastasiou, C

2004-01-01

The authors present a calculation of the fully differential cross section for Higgs boson production in the gluon fusion channel through next-to-next-to-leading order in perturbative QCD. They apply the method introduced in [1] to compute double real emission corrections. The calculation permits arbitrary cuts on the final state in the reaction hh → H + X. it can be easily extended to include decays of the Higgs boson into observable final states. In this Letter, they discuss the most important features of the calculation, and present some examples of physical applications that illustrate the range of observables that can be studied using the result. They compute the NNLO rapidity distribution of the Higgs boson, and also calculate the NNLO rapidity distribution with a veto on jet activity

SIVA/DIVA- INITIAL VALUE ORDINARY DIFFERENTIAL EQUATION SOLUTION VIA A VARIABLE ORDER ADAMS METHOD

Science.gov (United States)

Krogh, F. T.

1994-01-01

The SIVA/DIVA package is a collection of subroutines for the solution of ordinary differential equations. There are versions for single precision and double precision arithmetic. These solutions are applicable to stiff or nonstiff differential equations of first or second order. SIVA/DIVA requires fewer evaluations of derivatives than other variable order Adams predictor-corrector methods. There is an option for the direct integration of second order equations which can make integration of trajectory problems significantly more efficient. Other capabilities of SIVA/DIVA include: monitoring a user supplied function which can be separate from the derivative; dynamically controlling the step size; displaying or not displaying output at initial, final, and step size change points; saving the estimated local error; and reverse communication where subroutines return to the user for output or computation of derivatives instead of automatically performing calculations. The user must supply SIVA/DIVA with: 1) the number of equations; 2) initial values for the dependent and independent variables, integration stepsize, error tolerance, etc.; and 3) the driver program and operational parameters necessary for subroutine execution. SIVA/DIVA contains an extensive diagnostic message library should errors occur during execution. SIVA/DIVA is written in FORTRAN 77 for batch execution and is machine independent. It has a central memory requirement of approximately 120K of 8 bit bytes. This program was developed in 1983 and last updated in 1987.

Hadrons of arbitrary spin and heavy quark symmetry

International Nuclear Information System (INIS)

Hussain, F.; Thompson, G.; Koerner, J.G.

1993-11-01

We present a general construction of the spin content of the Bethe-Salpeter amplitudes (covariant wave functions) for heavy hadrons with arbitrary orbital excitations, using representations of l x O(3, 1). These wave functions incorporate the symmetries manifest in the heavy quark limit. In the baryonic sector we clearly differentiate between the Λ and Σ-type excited baryons. We then use the trace formalism to evaluate the weak transitions of ground state heavy hadrons to arbitrary excited heavy hadrons. The contributions of excited states to the Bjorken sum rule are also worked out in detail. (author). 21 refs

Solutions to an advanced functional partial differential equation of the pantograph type.

Science.gov (United States)

Zaidi, Ali A; Van Brunt, B; Wake, G C

2015-07-08

A model for cells structured by size undergoing growth and division leads to an initial boundary value problem that involves a first-order linear partial differential equation with a functional term. Here, size can be interpreted as DNA content or mass. It has been observed experimentally and shown analytically that solutions for arbitrary initial cell distributions are asymptotic as time goes to infinity to a certain solution called the steady size distribution. The full solution to the problem for arbitrary initial distributions, however, is elusive owing to the presence of the functional term and the paucity of solution techniques for such problems. In this paper, we derive a solution to the problem for arbitrary initial cell distributions. The method employed exploits the hyperbolic character of the underlying differential operator, and the advanced nature of the functional argument to reduce the problem to a sequence of simple Cauchy problems. The existence of solutions for arbitrary initial distributions is established along with uniqueness. The asymptotic relationship with the steady size distribution is established, and because the solution is known explicitly, higher-order terms in the asymptotics can be readily obtained.

Higher order differential calculus on SLq(N)

International Nuclear Information System (INIS)

Heckenberger, I.; Schueler, A.

1997-01-01

Let Γ be a bicovariant first order differential calculus on a Hopf algebra A. There are three possibilities to construct a differential N 0 -graded Hopf algebra Γcirconflex which contains Γ as its first order part. In all cases Γcirconflex is a quotient Γcirconflex = Γ x /J of the tensor algebra by some suitable ideal. We distinguish three possible choices u J, s J, and w J, where the first one generates the universal differential calculus (over Γ) and the last one is Woronowicz' external algebra. Let q be a transcendental complex number and let Γ be one of the N 2 -dimensional bicovariant first order differential calculi on the quantum group SL q (N). Then for N ≥ 3 the three ideals coincide. For Woronowicz' external algebra we calculate the dimensions of the spaces of left-invariant and bi-invariant k-forms. In this case each bi-invariant form is closed. In case of 4D ± calculi on SL q (2) the universal calculus is strictly larger than the other two calculi. In particular, the bi-invariant 1-form is not closed. (author)

Analysis of an Nth-order nonlinear differential-delay equation

Science.gov (United States)

Vallée, Réal; Marriott, Christopher

1989-01-01

The problem of a nonlinear dynamical system with delay and an overall response time which is distributed among N individual components is analyzed. Such a system can generally be modeled by an Nth-order nonlinear differential delay equation. A linear-stability analysis as well as a numerical simulation of that equation are performed and a comparison is made with the experimental results. Finally, a parallel is established between the first-order differential equation with delay and the Nth-order differential equation without delay.

Automatic differentiation for design sensitivity analysis of structural systems using multiple processors

Science.gov (United States)

Nguyen, Duc T.; Storaasli, Olaf O.; Qin, Jiangning; Qamar, Ramzi

1994-01-01

An automatic differentiation tool (ADIFOR) is incorporated into a finite element based structural analysis program for shape and non-shape design sensitivity analysis of structural systems. The entire analysis and sensitivity procedures are parallelized and vectorized for high performance computation. Small scale examples to verify the accuracy of the proposed program and a medium scale example to demonstrate the parallel vector performance on multiple CRAY C90 processors are included.

On the generalized Hankel±Clifford transformation of arbitrary order

Indian Academy of Sciences (India)

Е5Ж where C Y ЕxyЖИЕxyЖЕ З Жa2 J └ Е2. БББББ xy p. Ж is a generalized Bessel±Clifford function of first kind of order Е └ Ж, satisfying the differential equation xyHH └ Е З └ 1ЖyHЗ. Е x└1 З 1Жy И 0, which can be inverted by formulae. Proc. Indian Acad. Sci. (Math. Sci.), Vol. 110, No. 3, August 2000, pp. 293±304.

Preparation and tomographic reconstruction of an arbitrary single-photon path qubit

International Nuclear Information System (INIS)

Baek, So-Young; Kim, Yoon-Ho

2011-01-01

We report methods for preparation and tomographic reconstruction of an arbitrary single-photon path qubit. The arbitrary single-photon path qubit is prepared losslessly by passing the heralded single-photon state from spontaneous parametric down-conversion through variable beam splitter. Quantum state tomography of the single-photon path qubit is implemented by introducing path-projection measurements based on the first-order single-photon quantum interference. Using the state preparation and path-projection measurements methods for the single-photon path qubit, we demonstrate preparation and complete tomographic reconstruction of the single-photon path qubit with arbitrary purity. -- Highlights: → We report methods for preparation and tomographic reconstruction of an arbitrary single-photon path qubit. → We implement path-projection measurements based on the first-order single-photon quantum interference. → We demonstrate preparation and complete tomographic reconstruction of the single-photon path qubit with arbitrary purity.

Arbitrary-order Hilbert Spectral Analysis and Intermittency in Solar Wind Density Fluctuations

Science.gov (United States)

Carbone, Francesco; Sorriso-Valvo, Luca; Alberti, Tommaso; Lepreti, Fabio; Chen, Christopher H. K.; Němeček, Zdenek; Šafránková, Jana

2018-05-01

The properties of inertial- and kinetic-range solar wind turbulence have been investigated with the arbitrary-order Hilbert spectral analysis method, applied to high-resolution density measurements. Due to the small sample size and to the presence of strong nonstationary behavior and large-scale structures, the classical analysis in terms of structure functions may prove to be unsuccessful in detecting the power-law behavior in the inertial range, and may underestimate the scaling exponents. However, the Hilbert spectral method provides an optimal estimation of the scaling exponents, which have been found to be close to those for velocity fluctuations in fully developed hydrodynamic turbulence. At smaller scales, below the proton gyroscale, the system loses its intermittent multiscaling properties and converges to a monofractal process. The resulting scaling exponents, obtained at small scales, are in good agreement with those of classical fractional Brownian motion, indicating a long-term memory in the process, and the absence of correlations around the spectral-break scale. These results provide important constraints on models of kinetic-range turbulence in the solar wind.

Non-asymptotic fractional order differentiators via an algebraic parametric method

KAUST Repository

Liu, Dayan; Gibaru, O.; Perruquetti, Wilfrid

2012-01-01

Recently, Mboup, Join and Fliess [27], [28] introduced non-asymptotic integer order differentiators by using an algebraic parametric estimation method [7], [8]. In this paper, in order to obtain non-asymptotic fractional order differentiators we apply this algebraic parametric method to truncated expansions of fractional Taylor series based on the Jumarie's modified Riemann-Liouville derivative [14]. Exact and simple formulae for these differentiators are given where a sliding integration window of a noisy signal involving Jacobi polynomials is used without complex mathematical deduction. The efficiency and the stability with respect to corrupting noises of the proposed fractional order differentiators are shown in numerical simulations. © 2012 IEEE.

Non-asymptotic fractional order differentiators via an algebraic parametric method

KAUST Repository

Liu, Dayan

2012-08-01

Recently, Mboup, Join and Fliess [27], [28] introduced non-asymptotic integer order differentiators by using an algebraic parametric estimation method [7], [8]. In this paper, in order to obtain non-asymptotic fractional order differentiators we apply this algebraic parametric method to truncated expansions of fractional Taylor series based on the Jumarie\\'s modified Riemann-Liouville derivative [14]. Exact and simple formulae for these differentiators are given where a sliding integration window of a noisy signal involving Jacobi polynomials is used without complex mathematical deduction. The efficiency and the stability with respect to corrupting noises of the proposed fractional order differentiators are shown in numerical simulations. © 2012 IEEE.

Explicit formulation for natural frequencies of double-beam system with arbitrary boundary conditions

Energy Technology Data Exchange (ETDEWEB)

Mirzabeigy, Alborz; Madoliat, Reza [Iran University of Science and Technology, Narmak, Tehran (Iran, Islamic Republic of); Dabbagh, Vahid [University of Malaya, Kuala Lumpur (Malaysia)

2017-02-15

In this paper, free transverse vibration of two parallel beams connected through Winkler type elastic layer is investigated. Euler- Bernoulli beam hypothesis has been applied and it is assumed that boundary conditions of upper and lower beams are similar while arbitrary without any limitation even for non-ideal boundary conditions. Material properties and cross-section geometry of beams could be different from each other. The motion of the system is described by a homogeneous set of two partial differential equations, which is solved by using the classical Bernoulli-Fourier method. Explicit expressions are derived for the natural frequencies. In order to verify accuracy of results, the problem once again solved using modified Adomian decomposition method. Comparison between results indicates excellent accuracy of proposed formulation for any arbitrary boundary conditions. Derived explicit formulation is simplest method to determine natural frequencies of double-beam systems with high level of accuracy in comparison with other methods in literature.

Unambiguous formalism for higher order Lagrangian field theories

International Nuclear Information System (INIS)

Campos, Cedric M; De Leon, Manuel; De Diego, David MartIn; Vankerschaver, Joris

2009-01-01

The aim of this paper is to propose an unambiguous intrinsic formalism for higher order field theories which avoids the arbitrariness in the generalization of the conventional description of field theories, and implies the existence of different Cartan forms and Legendre transformations. We propose a differential-geometric setting for the dynamics of a higher order field theory, based on the Skinner and Rusk formalism for mechanics. This approach incorporates aspects of both the Lagrangian and the Hamiltonian description, since the field equations are formulated using the Lagrangian on a higher order jet bundle and the canonical multisymplectic form on its affine dual. As both of these objects are uniquely defined, the Skinner-Rusk approach has the advantage that it does not suffer from the arbitrariness in conventional descriptions. The result is that we obtain a unique and global intrinsic version of the Euler-Lagrange equations for higher order field theories. Several examples illustrate our construction.

Parameters and Fractional Differentiation Orders Estimation for Linear Continuous-Time Non-Commensurate Fractional Order Systems

KAUST Repository

Belkhatir, Zehor; Laleg-Kirati, Taous-Meriem

2017-01-01

This paper proposes a two-stage estimation algorithm to solve the problem of joint estimation of the parameters and the fractional differentiation orders of a linear continuous-time fractional system with non-commensurate orders. The proposed algorithm combines the modulating functions and the first-order Newton methods. Sufficient conditions ensuring the convergence of the method are provided. An error analysis in the discrete case is performed. Moreover, the method is extended to the joint estimation of smooth unknown input and fractional differentiation orders. The performance of the proposed approach is illustrated with different numerical examples. Furthermore, a potential application of the algorithm is proposed which consists in the estimation of the differentiation orders of a fractional neurovascular model along with the neural activity considered as input for this model.

Parameters and Fractional Differentiation Orders Estimation for Linear Continuous-Time Non-Commensurate Fractional Order Systems

KAUST Repository

Belkhatir, Zehor

2017-05-31

This paper proposes a two-stage estimation algorithm to solve the problem of joint estimation of the parameters and the fractional differentiation orders of a linear continuous-time fractional system with non-commensurate orders. The proposed algorithm combines the modulating functions and the first-order Newton methods. Sufficient conditions ensuring the convergence of the method are provided. An error analysis in the discrete case is performed. Moreover, the method is extended to the joint estimation of smooth unknown input and fractional differentiation orders. The performance of the proposed approach is illustrated with different numerical examples. Furthermore, a potential application of the algorithm is proposed which consists in the estimation of the differentiation orders of a fractional neurovascular model along with the neural activity considered as input for this model.

Differential equations

CERN Document Server

Barbu, Viorel

2016-01-01

This textbook is a comprehensive treatment of ordinary differential equations, concisely presenting basic and essential results in a rigorous manner. Including various examples from physics, mechanics, natural sciences, engineering and automatic theory, Differential Equations is a bridge between the abstract theory of differential equations and applied systems theory. Particular attention is given to the existence and uniqueness of the Cauchy problem, linear differential systems, stability theory and applications to first-order partial differential equations. Upper undergraduate students and researchers in applied mathematics and systems theory with a background in advanced calculus will find this book particularly useful. Supplementary topics are covered in an appendix enabling the book to be completely self-contained.

Fuchs indices and the first integrals of nonlinear differential equations

International Nuclear Information System (INIS)

Kudryashov, Nikolai A.

2005-01-01

New method of finding the first integrals of nonlinear differential equations in polynomial form is presented. Basic idea of our approach is to use the scaling of solution of nonlinear differential equation and to find the dimensions of arbitrary constants in the Laurent expansion of the general solution. These dimensions allows us to obtain the scalings of members for the first integrals of nonlinear differential equations. Taking the polynomials with unknown coefficients into account we present the algorithm of finding the first integrals of nonlinear differential equations in the polynomial form. Our method is applied to look for the first integrals of eight nonlinear ordinary differential equations of the fourth order. The general solution of one of the fourth order ordinary differential equations is given

Phase integral approximation for coupled ordinary differential equations of the Schroedinger type

International Nuclear Information System (INIS)

Skorupski, Andrzej A.

2008-01-01

Four generalizations of the phase integral approximation (PIA) to sets of ordinary differential equations of Schroedinger type [u j '' (x)+Σ k=1 N R jk (x)u k (x)=0, j=1,2,...,N] are described. The recurrence relations for higher order corrections are given in a form valid to arbitrary order and for the matrix R(x)[≡(R jk (x))] either Hermitian or non-Hermitian. For Hermitian and negative definite R(x) matrices, a Wronskian conserving PIA theory is formulated, which generalizes Fulling's current conserving theory pertinent to positive definite R(x) matrices. The idea of a modification of the PIA, which is well known for one equation [u '' (x)+R(x)u(x)=0], is generalized to sets. A simplification of Wronskian or current conserving theories is proposed which in each order eliminates one integration from the formulas for higher order corrections. If the PIA is generated by a nondegenerate eigenvalue of the R(x) matrix, the eliminated integration is the only one present. In that case, the simplified theory becomes fully algorithmic and is generalized to non-Hermitian R(x) matrices. The general theory is illustrated by a few examples automatically generated by using the author's program in MATHEMATICA published in e-print arXiv:0710.5406 [math-ph

Carrier introduction to moire pattern for automatic fringe-order distinguishing

International Nuclear Information System (INIS)

Fang, J.; Laermann, K.H.

1992-01-01

This paper presents an automatic procedure of pseudo-colour encoding of moire fringe orders. A carrier consisting of parallel fringes is introduced before the specimen deforms. The carrier pattern is captured by a camera and then stored in computer as a standard image. The space of the carrier fringes is distored by the strains on the specimen as it is loaded. On a certain condition, the orders of the frequency-modulated carrier still vary monotonically so that they can be easyly distinguished. Both the standard fringe-carrier and the frequency-modulated fringe pattern are transformed into two digital images, of which every fringe is encoded by one of the pseudo-colour codes corresponding to the monotonical fringe orders. At each pixel, the difference between the colour sequences of two images is calculated to obtain the fringe order of pure deformation. The moire pattern of the in-plane displacement is restored as a pseudo-colour image by whose colour-change the variation of the fringe orders is displayed. (orig.)

Fractional Order Differentiation by Integration and Error Analysis in Noisy Environment

KAUST Repository

Liu, Dayan

2015-03-31

The integer order differentiation by integration method based on the Jacobi orthogonal polynomials for noisy signals was originally introduced by Mboup, Join and Fliess. We propose to extend this method from the integer order to the fractional order to estimate the fractional order derivatives of noisy signals. Firstly, two fractional order differentiators are deduced from the Jacobi orthogonal polynomial filter, using the Riemann-Liouville and the Caputo fractional order derivative definitions respectively. Exact and simple formulae for these differentiators are given by integral expressions. Hence, they can be used for both continuous-time and discrete-time models in on-line or off-line applications. Secondly, some error bounds are provided for the corresponding estimation errors. These bounds allow to study the design parameters\\' influence. The noise error contribution due to a large class of stochastic processes is studied in discrete case. The latter shows that the differentiator based on the Caputo fractional order derivative can cope with a class of noises, whose mean value and variance functions are polynomial time-varying. Thanks to the design parameters analysis, the proposed fractional order differentiators are significantly improved by admitting a time-delay. Thirdly, in order to reduce the calculation time for on-line applications, a recursive algorithm is proposed. Finally, the proposed differentiator based on the Riemann-Liouville fractional order derivative is used to estimate the state of a fractional order system and numerical simulations illustrate the accuracy and the robustness with respect to corrupting noises.

Robust fractional order differentiators using generalized modulating functions method

KAUST Repository

Liu, Dayan

2015-02-01

This paper aims at designing a fractional order differentiator for a class of signals satisfying a linear differential equation with unknown parameters. A generalized modulating functions method is proposed first to estimate the unknown parameters, then to derive accurate integral formulae for the left-sided Riemann-Liouville fractional derivatives of the studied signal. Unlike the improper integral in the definition of the left-sided Riemann-Liouville fractional derivative, the integrals in the proposed formulae can be proper and be considered as a low-pass filter by choosing appropriate modulating functions. Hence, digital fractional order differentiators applicable for on-line applications are deduced using a numerical integration method in discrete noisy case. Moreover, some error analysis are given for noise error contributions due to a class of stochastic processes. Finally, numerical examples are given to show the accuracy and robustness of the proposed fractional order differentiators.

Robust fractional order differentiators using generalized modulating functions method

KAUST Repository

Liu, Dayan; Laleg-Kirati, Taous-Meriem

2015-01-01

This paper aims at designing a fractional order differentiator for a class of signals satisfying a linear differential equation with unknown parameters. A generalized modulating functions method is proposed first to estimate the unknown parameters, then to derive accurate integral formulae for the left-sided Riemann-Liouville fractional derivatives of the studied signal. Unlike the improper integral in the definition of the left-sided Riemann-Liouville fractional derivative, the integrals in the proposed formulae can be proper and be considered as a low-pass filter by choosing appropriate modulating functions. Hence, digital fractional order differentiators applicable for on-line applications are deduced using a numerical integration method in discrete noisy case. Moreover, some error analysis are given for noise error contributions due to a class of stochastic processes. Finally, numerical examples are given to show the accuracy and robustness of the proposed fractional order differentiators.

Asymptotic behavior of second-order impulsive differential equations

Directory of Open Access Journals (Sweden)

Haifeng Liu

2011-02-01

Full Text Available In this article, we study the asymptotic behavior of all solutions of 2-th order nonlinear delay differential equation with impulses. Our main tools are impulsive differential inequalities and the Riccati transformation. We illustrate the results by an example.

Harmonic arbitrary waveform generator

Science.gov (United States)

Roberts, Brock Franklin

2017-11-28

High frequency arbitrary waveforms have applications in radar, communications, medical imaging, therapy, electronic warfare, and charged particle acceleration and control. State of the art arbitrary waveform generators are limited in the frequency they can operate by the speed of the Digital to Analog converters that directly create their arbitrary waveforms. The architecture of the Harmonic Arbitrary Waveform Generator allows the phase and amplitude of the high frequency content of waveforms to be controlled without taxing the Digital to Analog converters that control them. The Harmonic Arbitrary Waveform Generator converts a high frequency input, into a precision, adjustable, high frequency arbitrary waveform.

On realization of nonlinear systems described by higher-order differential equations

NARCIS (Netherlands)

van der Schaft, Arjan

1987-01-01

We consider systems of smooth nonlinear differential and algebraic equations in which some of the variables are distinguished as “external variables.” The realization problem is to replace the higher-order implicit differential equations by first-order explicit differential equations and the

Solution of linear ordinary differential equations by means of the method of variation of arbitrary constants

DEFF Research Database (Denmark)

Mejlbro, Leif

1997-01-01

An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians.......An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians....

Digital PI-PD controller design for arbitrary order systems: Dominant pole placement approach.

Science.gov (United States)

Dincel, Emre; Söylemez, Mehmet Turan

2018-05-02

In this paper, a digital PI-PD controller design method is proposed for arbitrary order systems with or without time-delay to achieve desired transient response in the closed-loop via dominant pole placement approach. The digital PI-PD controller design problem is solved by converting the original problem to the digital PID controller design problem. Firstly, parametrization of the digital PID controllers which assign dominant poles to desired location is done. After that the subset of digital PID controller parameters in which the remaining poles are located away from the dominant pole pair is found via Chebyshev polynomials. The obtained PID controller parameters are then transformed into the PI-PD controller parameters by considering the closed-loop controller zero and the design is completed. Success of the proposed design method is firstly demonstrated on an example transfer function and compared with the well-known PID controller methods from the literature through simulations. After that the design method is implemented on the fan and plate laboratory system in a real environment. Copyright © 2018 ISA. Published by Elsevier Ltd. All rights reserved.

Target recognition of ladar range images using even-order Zernike moments.

Science.gov (United States)

Liu, Zheng-Jun; Li, Qi; Xia, Zhi-Wei; Wang, Qi

2012-11-01

Ladar range images have attracted considerable attention in automatic target recognition fields. In this paper, Zernike moments (ZMs) are applied to classify the target of the range image from an arbitrary azimuth angle. However, ZMs suffer from high computational costs. To improve the performance of target recognition based on small samples, even-order ZMs with serial-parallel backpropagation neural networks (BPNNs) are applied to recognize the target of the range image. It is found that the rotation invariance and classified performance of the even-order ZMs are both better than for odd-order moments and for moments compressed by principal component analysis. The experimental results demonstrate that combining the even-order ZMs with serial-parallel BPNNs can significantly improve the recognition rate for small samples.

Automatic validation of numerical solutions

DEFF Research Database (Denmark)

Stauning, Ole

1997-01-01

This thesis is concerned with ``Automatic Validation of Numerical Solutions''. The basic theory of interval analysis and self-validating methods is introduced. The mean value enclosure is applied to discrete mappings for obtaining narrow enclosures of the iterates when applying these mappings...... differential equations, but in this thesis, we describe how to use the methods for enclosing iterates of discrete mappings, and then later use them for discretizing solutions of ordinary differential equations. The theory of automatic differentiation is introduced, and three methods for obtaining derivatives...... are described: The forward, the backward, and the Taylor expansion methods. The three methods have been implemented in the C++ program packages FADBAD/TADIFF. Some examples showing how to use the three metho ds are presented. A feature of FADBAD/TADIFF not present in other automatic differentiation packages...

Parameter optimization of differential evolution algorithm for automatic playlist generation problem

Science.gov (United States)

Alamag, Kaye Melina Natividad B.; Addawe, Joel M.

2017-11-01

With the digitalization of music, the number of collection of music increased largely and there is a need to create lists of music that filter the collection according to user preferences, thus giving rise to the Automatic Playlist Generation Problem (APGP). Previous attempts to solve this problem include the use of search and optimization algorithms. If a music database is very large, the algorithm to be used must be able to search the lists thoroughly taking into account the quality of the playlist given a set of user constraints. In this paper we perform an evolutionary meta-heuristic optimization algorithm, Differential Evolution (DE) using different combination of parameter values and select the best performing set when used to solve four standard test functions. Performance of the proposed algorithm is then compared with normal Genetic Algorithm (GA) and a hybrid GA with Tabu Search. Numerical simulations are carried out to show better results from Differential Evolution approach with the optimized parameter values.

Modeling Ability Differentiation in the Second-Order Factor Model

Science.gov (United States)

Molenaar, Dylan; Dolan, Conor V.; van der Maas, Han L. J.

2011-01-01

In this article we present factor models to test for ability differentiation. Ability differentiation predicts that the size of IQ subtest correlations decreases as a function of the general intelligence factor. In the Schmid-Leiman decomposition of the second-order factor model, we model differentiation by introducing heteroscedastic residuals,…

On nonlinear differential equation with exact solutions having various pole orders

International Nuclear Information System (INIS)

Kudryashov, N.A.

2015-01-01

We consider a nonlinear ordinary differential equation having solutions with various movable pole order on the complex plane. We show that the pole order of exact solution is determined by values of parameters of the equation. Exact solutions in the form of the solitary waves for the second order nonlinear differential equation are found taking into account the method of the logistic function. Exact solutions of differential equations are discussed and analyzed

Series-, Parallel-, and Inter-Connection of Solid-State Arbitrary Fractional-Order Capacitors: Theoretical Study and Experimental Verification

KAUST Repository

Kartci, Aslihan

2018-02-26

In the paper, general analytical formulas are introduced for the determination of equivalent impedance, magnitude, and phase, i.e. order, for n arbitrary fractional-order capacitors (FoCs) connected in series, parallel, and their interconnection. The approach presented helps to evaluate these relevant quantities in the fractional domain since the order of each element has a significant effect on the impedance of each FoC and their equivalent capacitance cannot be considered. Three types of solid-state fractional-order passive capacitors of different orders, using ferroelectric polymer and reduced Graphene Oxide-percolated P(VDF-TrFE-CFE) composite structures, are fabricated and characterized. Using an impedance analyzer, the behavior of the devices was found to be stable in the frequency range 0.2MHz–20MHz, with a phase angle deviation of ±4 degrees. Multiple numerical and experimental case studies are given, in particular for two and three connected FoCs. The fundamental issues of the measurement units of the FoCs connected in series and parallel are derived. A MATLAB open access source code is given in Appendix sec:append for easy calculation of the equivalent FoC magnitude and phase. The experimental results are in good agreement with the theoretical assumptions.

Hojman's theorem of the third-order ordinary differential equation

International Nuclear Information System (INIS)

Hong-Sheng, Lü; Hong-Bin, Zhang; Shu-Long, Gu

2009-01-01

This paper extends Hojman's conservation law to the third-order differential equation. A new conserved quantity is constructed based on the Lie group of transformation generators of the equations of motion. The generators contain variations of the time and generalized coordinates. Two independent non-trivial conserved quantities of the third-order ordinary differential equation are obtained. A simple example is presented to illustrate the applications of the results. (general)

Lagrange-Noether method for solving second-order differential equations

Institute of Scientific and Technical Information of China (English)

Wu Hui-Bin; Wu Run-Heng

2009-01-01

The purpose of this paper is to provide a new method called the Lagrange-Noether method for solving second-order differential equations. The method is,firstly,to write the second-order differential equations completely or partially in the form of Lagrange equations,and secondly,to obtain the integrals of the equations by using the Noether theory of the Lagrange system. An example is given to illustrate the application of the result.

An inverse method for non linear ablative thermics with experimentation of automatic differentiation

Energy Technology Data Exchange (ETDEWEB)

Alestra, S [Simulation Information Technology and Systems Engineering, EADS IW Toulouse (France); Collinet, J [Re-entry Systems and Technologies, EADS ASTRIUM ST, Les Mureaux (France); Dubois, F [Professor of Applied Mathematics, Conservatoire National des Arts et Metiers Paris (France)], E-mail: stephane.alestra@eads.net, E-mail: jean.collinet@astrium.eads.net, E-mail: fdubois@cnam.fr

2008-11-01

Thermal Protection System is a key element for atmospheric re-entry missions of aerospace vehicles. The high level of heat fluxes encountered in such missions has a direct effect on mass balance of the heat shield. Consequently, the identification of heat fluxes is of great industrial interest but is in flight only available by indirect methods based on temperature measurements. This paper is concerned with inverse analyses of highly evolutive heat fluxes. An inverse problem is used to estimate transient surface heat fluxes (convection coefficient), for degradable thermal material (ablation and pyrolysis), by using time domain temperature measurements on thermal protection. The inverse problem is formulated as a minimization problem involving an objective functional, through an optimization loop. An optimal control formulation (Lagrangian, adjoint and gradient steepest descent method combined with quasi-Newton method computations) is then developed and applied, using Monopyro, a transient one-dimensional thermal model with one moving boundary (ablative surface) that has been developed since many years by ASTRIUM-ST. To compute numerically the adjoint and gradient quantities, for the inverse problem in heat convection coefficient, we have used both an analytical manual differentiation and an Automatic Differentiation (AD) engine tool, Tapenade, developed at INRIA Sophia-Antipolis by the TROPICS team. Several validation test cases, using synthetic temperature measurements are carried out, by applying the results of the inverse method with minimization algorithm. Accurate results of identification on high fluxes test cases, and good agreement for temperatures restitutions, are obtained, without and with ablation and pyrolysis, using bad fluxes initial guesses. First encouraging results with an automatic differentiation procedure are also presented in this paper.

Linear matrix differential equations of higher-order and applications

Directory of Open Access Journals (Sweden)

Mustapha Rachidi

2008-07-01

Full Text Available In this article, we study linear differential equations of higher-order whose coefficients are square matrices. The combinatorial method for computing the matrix powers and exponential is adopted. New formulas representing auxiliary results are obtained. This allows us to prove properties of a large class of linear matrix differential equations of higher-order, in particular results of Apostol and Kolodner are recovered. Also illustrative examples and applications are presented.

A comparison of different methods to implement higher order derivatives of density functionals

Energy Technology Data Exchange (ETDEWEB)

van Dam, Hubertus J.J. [Brookhaven National Lab. (BNL), Upton, NY (United States)

2016-05-18

Density functional theory is the dominant approach in electronic structure methods today. To calculate properties higher order derivatives of the density functionals are required. These derivatives might be implemented manually,by automatic differentiation, or by symbolic algebra programs. Different authors have cited different reasons for using the particular method of their choice. This paper presents work where all three approaches were used and the strengths and weaknesses of each approach are considered. It is found that all three methods produce code that is suffficiently performanted for practical applications, despite the fact that our symbolic algebra generated code and our automatic differentiation code still have scope for significant optimization. The automatic differentiation approach is the best option for producing readable and maintainable code.

Oscillation criteria for third order delay nonlinear differential equations

Directory of Open Access Journals (Sweden)

E. M. Elabbasy

2012-01-01

via comparison with some first differential equations whose oscillatory characters are known. Our results generalize and improve some known results for oscillation of third order nonlinear differential equations. Some examples are given to illustrate the main results.

On Fractional Order Hybrid Differential Equations

Directory of Open Access Journals (Sweden)

Mohamed A. E. Herzallah

2014-01-01

Full Text Available We develop the theory of fractional hybrid differential equations with linear and nonlinear perturbations involving the Caputo fractional derivative of order 0<α<1. Using some fixed point theorems we prove the existence of mild solutions for two types of hybrid equations. Examples are given to illustrate the obtained results.

Design of an ultrafast all-optical differentiator based on a fiber Bragg grating in transmission.

Science.gov (United States)

Preciado, Miguel A; Muriel, Miguel A

2008-11-01

We propose and analyze a first-order optical differentiator based on a fiber Bragg grating (FBG) in transmission. It is shown in the examples that a simple uniform-period FBG in a very strong coupling regime (maximum reflectivity very close to 100%) can perform close to ideal temporal differentiation of the complex envelope of an arbitrary-input optical signal.

Applications of automatic differentiation in computational fluid dynamics

Science.gov (United States)

Green, Lawrence L.; Carle, A.; Bischof, C.; Haigler, Kara J.; Newman, Perry A.

1994-01-01

Automatic differentiation (AD) is a powerful computational method that provides for computing exact sensitivity derivatives (SD) from existing computer programs for multidisciplinary design optimization (MDO) or in sensitivity analysis. A pre-compiler AD tool for FORTRAN programs called ADIFOR has been developed. The ADIFOR tool has been easily and quickly applied by NASA Langley researchers to assess the feasibility and computational impact of AD in MDO with several different FORTRAN programs. These include a state-of-the-art three dimensional multigrid Navier-Stokes flow solver for wings or aircraft configurations in transonic turbulent flow. With ADIFOR the user specifies sets of independent and dependent variables with an existing computer code. ADIFOR then traces the dependency path throughout the code, applies the chain rule to formulate derivative expressions, and generates new code to compute the required SD matrix. The resulting codes have been verified to compute exact non-geometric and geometric SD for a variety of cases. in less time than is required to compute the SD matrix using centered divided differences.

Hankel complementary integral transformations of arbitrary order

Directory of Open Access Journals (Sweden)

M. Linares Linares

1992-01-01

Full Text Available Four selfreciprocal integral transformations of Hankel type are defined through(ℋi,μf(y=Fi(y=∫0∞αi(xℊi,μ(xyf(xdx,   ℋi,μ−1=ℋi,μ,where i=1,2,3,4; μ≥0; α1(x=x1+2μ, ℊ1,μ(x=x−μJμ(x, Jμ(x being the Bessel function of the first kind of order μ; α2(x=x1−2μ, ℊ2,μ(x=(−1μx2μℊ1,μ(x; α3(x=x−1−2μ, ℊ3,μ(x=x1+2μℊ1,μ(x, and α4(x=x−1+2μ, ℊ4,μ(x=(−1μxℊ1,μ(x. The simultaneous use of transformations ℋ1,μ, and ℋ2,μ, (which are denoted by ℋμ allows us to solve many problems of Mathematical Physics involving the differential operator Δμ=D2+(1+2μx−1D, whereas the pair of transformations ℋ3,μ and ℋ4,μ, (which we express by ℋμ* permits us to tackle those problems containing its adjoint operator Δμ*=D2−(1+2μx−1D+(1+2μx−2, no matter what the real value of μ be. These transformations are also investigated in a space of generalized functions according to the mixed Parseval equation∫0∞f(xg(xdx=∫0∞(ℋμf(y(ℋμ*g(ydy,which is now valid for all real μ.

Second-order differential-delay equation to describe a hybrid bistable device

Science.gov (United States)

Vallee, R.; Dubois, P.; Cote, M.; Delisle, C.

1987-08-01

The problem of a dynamical system with delayed feedback, a hybrid bistable device, characterized by n response times and described by an nth-order differential-delay equation (DDE) is discussed. Starting from a linear-stability analysis of the DDE, the effects of the second-order differential terms on the position of the first bifurcation and on the frequency of the resulting self-oscillation are shown. The effects of the third-order differential terms on the first bifurcation are also considered. Experimental results are shown to support the linear analysis.

Third order differential equations with delay

Directory of Open Access Journals (Sweden)

Petr Liška

2015-05-01

Full Text Available In this paper, we study the oscillation and asymptotic properties of solutions of certain nonlinear third order differential equations with delay. In particular, we extend results of I. Mojsej (Nonlinear Analysis 68, 2008 and we improve conditions on the property B of N. Parhi and S. Padhi (Indian J. Pure Appl. Math., 33, 2002.

Multi-order, automatic dispersion compensation for 1.28 Terabaud signals

Science.gov (United States)

Paquot, Yvan; Schröder, Jochen; Van Erps, Jürgen; Vo, Trung D.; Pelusi, Mark D.; Madden, Steve; Luther-Davies, Barry; Eggleton, Benjamin J.

2012-06-01

Transmitting ultra-high symbol rate optical signals remains a challenge due to their high sensitivity to fluctuations of GVD and higher orders of dispersion in the transmission link. Being able to cancel the impairments due to those fluctuations is a key requirement to make transmission of ultrashort optical pulses practical. We demonstrate an automatic compensation scheme able to keep an Optical Time Division Multiplexed (OTDM) signal stable at a bandwidth of up to 1.28 Tbaud in spite of external perturbations. Our approach is based on monitoring the signal with a photonic-chip-based all-optical RF-spectrum analyzer. The measurement of a single parameter extracted from the RF-spectrum is used to drive a multidimensional optimization algorithm. We apply the method to the real time simultaneous compensation for 2nd, 3rd and 4th order dispersion using an LCOS spectral pulse shaper (SPS) as a tunable dispersion compensator.

First-order partial differential equations

CERN Document Server

Rhee, Hyun-Ku; Amundson, Neal R

2001-01-01

This first volume of a highly regarded two-volume text is fully usable on its own. After going over some of the preliminaries, the authors discuss mathematical models that yield first-order partial differential equations; motivations, classifications, and some methods of solution; linear and semilinear equations; chromatographic equations with finite rate expressions; homogeneous and nonhomogeneous quasilinear equations; formation and propagation of shocks; conservation equations, weak solutions, and shock layers; nonlinear equations; and variational problems. Exercises appear at the end of mo

Conventional hamiltonian for first-order differential systems

International Nuclear Information System (INIS)

Farias, J.R.

1984-01-01

Lagrangian systems corresponding to a given set of 2n first-order ordinary differential equations are singular ones. In despite this, it is shown that these systems can be put into a Hamiltonian form in the usual manner. (Author) [pt

Asymptotic behavior of solutions of linear multi-order fractional differential equation systems

OpenAIRE

Diethelm, Kai; Siegmund, Stefan; Tuan, H. T.

2017-01-01

In this paper, we investigate some aspects of the qualitative theory for multi-order fractional differential equation systems. First, we obtain a fundamental result on the existence and uniqueness for multi-order fractional differential equation systems. Next, a representation of solutions of homogeneous linear multi-order fractional differential equation systems in series form is provided. Finally, we give characteristics regarding the asymptotic behavior of solutions to some classes of line...

Processes of arbitrary order in quantum electrodynamics with a pair-creating external field

International Nuclear Information System (INIS)

Gitman, D.M.

1977-01-01

Dyson's perturbation theory analogue for quantum electrodynamical processes with arbitrary initial and final states in an external field creating pairs is discussed. The interaction with the field is taken into account exactly. The possibility of using Feynman diagrams, together with modified correspondence rules, for the representation of the above mentioned processes is demonstrated. (author)

Exact decoupling of the Dirac Hamiltonian. II. The generalized Douglas-Kroll-Hess transformation up to arbitrary order

International Nuclear Information System (INIS)

Reiher, Markus; Wolf, Alexander

2004-01-01

In order to achieve exact decoupling of the Dirac Hamiltonian within a unitary transformation scheme, we have discussed in part I of this series that either a purely numerical iterative technique (the Barysz-Sadlej-Snijders method) or a stepwise analytic approach (the Douglas-Kroll-Hess method) are possible. For the evaluation of Douglas-Kroll-Hess Hamiltonians up to a pre-defined order it was shown that a symbolic scheme has to be employed. In this work, an algorithm for this analytic derivation of Douglas-Kroll-Hess Hamiltonians up to any arbitrary order in the external potential is presented. We discuss how an estimate for the necessary order for exact decoupling (within machine precision) for a given system can be determined from the convergence behavior of the Douglas-Kroll-Hess expansion prior to a quantum chemical calculation. Once this maximum order has been accomplished, the spectrum of the positive-energy part of the decoupled Hamiltonian, e.g., for electronic bound states, cannot be distinguished from the corresponding part of the spectrum of the Dirac operator. An efficient scalar-relativistic implementation of the symbolic operations for the evaluation of the positive-energy part of the block-diagonal Hamiltonian is presented, and its accuracy is tested for ground-state energies of one-electron ions over the whole periodic table. Furthermore, the first many-electron calculations employing sixth up to fourteenth order DKH Hamiltonians are presented

Exact decoupling of the Dirac Hamiltonian. II. The generalized Douglas-Kroll-Hess transformation up to arbitrary order.

Science.gov (United States)

Reiher, Markus; Wolf, Alexander

2004-12-08

In order to achieve exact decoupling of the Dirac Hamiltonian within a unitary transformation scheme, we have discussed in part I of this series that either a purely numerical iterative technique (the Barysz-Sadlej-Snijders method) or a stepwise analytic approach (the Douglas-Kroll-Hess method) are possible. For the evaluation of Douglas-Kroll-Hess Hamiltonians up to a pre-defined order it was shown that a symbolic scheme has to be employed. In this work, an algorithm for this analytic derivation of Douglas-Kroll-Hess Hamiltonians up to any arbitrary order in the external potential is presented. We discuss how an estimate for the necessary order for exact decoupling (within machine precision) for a given system can be determined from the convergence behavior of the Douglas-Kroll-Hess expansion prior to a quantum chemical calculation. Once this maximum order has been accomplished, the spectrum of the positive-energy part of the decoupled Hamiltonian, e.g., for electronic bound states, cannot be distinguished from the corresponding part of the spectrum of the Dirac operator. An efficient scalar-relativistic implementation of the symbolic operations for the evaluation of the positive-energy part of the block-diagonal Hamiltonian is presented, and its accuracy is tested for ground-state energies of one-electron ions over the whole periodic table. Furthermore, the first many-electron calculations employing sixth up to fourteenth order DKH Hamiltonians are presented. (c) 2004 American Institute of Physics.

Higher Order Differential Attack on 6-Round MISTY1

Science.gov (United States)

Tsunoo, Yukiyasu; Saito, Teruo; Nakashima, Hiroki; Shigeri, Maki

MISTY1 is a 64-bit block cipher that has provable security against differential and linear cryptanalysis. MISTY1 is one of the algorithms selected in the European NESSIE project, and it has been recommended for Japanese e-Government ciphers by the CRYPTREC project. This paper reports a previously unknown higher order differential characteristic of 4-round MISTY1 with the FL functions. It also shows that a higher order differential attack that utilizes this newly discovered characteristic is successful against 6-round MISTY1 with the FL functions. This attack can recover a partial subkey with a data complexity of 253.7 and a computational complexity of 264.4, which is better than any previous cryptanalysis of MISTY1.

Automatic Differentiation in Quantum Chemistry with Applications to Fully Variational Hartree-Fock.

Science.gov (United States)

Tamayo-Mendoza, Teresa; Kreisbeck, Christoph; Lindh, Roland; Aspuru-Guzik, Alán

2018-05-23

Automatic differentiation (AD) is a powerful tool that allows calculating derivatives of implemented algorithms with respect to all of their parameters up to machine precision, without the need to explicitly add any additional functions. Thus, AD has great potential in quantum chemistry, where gradients are omnipresent but also difficult to obtain, and researchers typically spend a considerable amount of time finding suitable analytical forms when implementing derivatives. Here, we demonstrate that AD can be used to compute gradients with respect to any parameter throughout a complete quantum chemistry method. We present DiffiQult , a Hartree-Fock implementation, entirely differentiated with the use of AD tools. DiffiQult is a software package written in plain Python with minimal deviation from standard code which illustrates the capability of AD to save human effort and time in implementations of exact gradients in quantum chemistry. We leverage the obtained gradients to optimize the parameters of one-particle basis sets in the context of the floating Gaussian framework.

Aerodynamic design applying automatic differentiation and using robust variable fidelity optimization

Science.gov (United States)

Takemiya, Tetsushi

, and that (2) the AMF terminates optimization erroneously when the optimization problems have constraints. The first problem is due to inaccuracy in computing derivatives in the AMF, and the second problem is due to erroneous treatment of the trust region ratio, which sets the size of the domain for an optimization in the AMF. In order to solve the first problem of the AMF, automatic differentiation (AD) technique, which reads the codes of analysis models and automatically generates new derivative codes based on some mathematical rules, is applied. If derivatives are computed with the generated derivative code, they are analytical, and the required computational time is independent of the number of design variables, which is very advantageous for realistic aerospace engineering problems. However, if analysis models implement iterative computations such as computational fluid dynamics (CFD), which solves system partial differential equations iteratively, computing derivatives through the AD requires a massive memory size. The author solved this deficiency by modifying the AD approach and developing a more efficient implementation with CFD, and successfully applied the AD to general CFD software. In order to solve the second problem of the AMF, the governing equation of the trust region ratio, which is very strict against the violation of constraints, is modified so that it can accept the violation of constraints within some tolerance. By accepting violations of constraints during the optimization process, the AMF can continue optimization without terminating immaturely and eventually find the true optimum design point. With these modifications, the AMF is referred to as "Robust AMF," and it is applied to airfoil and wing aerodynamic design problems using Euler CFD software. The former problem has 21 design variables, and the latter 64. In both problems, derivatives computed with the proposed AD method are first compared with those computed with the finite

Pair production of arbitrary spin particles by electromagnetic fields

International Nuclear Information System (INIS)

Kruglov, S.I.

2006-01-01

The exact solutions of the wave equation for arbitrary spin particles in the field of the soliton-like electric impulse were obtained. The differential probability of pair production of particles by electromagnetic fields has been evaluated on the basis of the exact solutions. As a particular case, the particle pair production in the constant and uniform electric field were studied

The Cauchy problem for higher order abstract differential equations

CERN Document Server

Xiao, Ti-Jun

1998-01-01

This monograph is the first systematic exposition of the theory of the Cauchy problem for higher order abstract linear differential equations, which covers all the main aspects of the developed theory. The main results are complete with detailed proofs and established recently, containing the corresponding theorems for first and incomplete second order cases and therefore for operator semigroups and cosine functions. They will find applications in many fields. The special power of treating the higher order problems directly is demonstrated, as well as that of the vector-valued Laplace transforms in dealing with operator differential equations and operator families. The reader is expected to have a knowledge of complex and functional analysis.

Reduction of static field equation of Faddeev model to first order PDE

International Nuclear Information System (INIS)

Hirayama, Minoru; Shi Changguang

2007-01-01

A method to solve the static field equation of the Faddeev model is presented. For a special combination of the concerned field, we adopt a form which is compatible with the field equation and involves two arbitrary complex functions. As a result, the static field equation is reduced to a set of first order partial differential equations

Green's matrix for a second-order self-adjoint matrix differential operator

International Nuclear Information System (INIS)

Sisman, Tahsin Cagri; Tekin, Bayram

2010-01-01

A systematic construction of the Green's matrix for a second-order self-adjoint matrix differential operator from the linearly independent solutions of the corresponding homogeneous differential equation set is carried out. We follow the general approach of extracting the Green's matrix from the Green's matrix of the corresponding first-order system. This construction is required in the cases where the differential equation set cannot be turned to an algebraic equation set via transform techniques.

Lattice Boltzmann model for high-order nonlinear partial differential equations.

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂_{t}ϕ+∑_{k=1}^{m}α_{k}∂_{x}^{k}Π_{k}(ϕ)=0 (1≤k≤m≤6), α_{k} are constant coefficients, Π_{k}(ϕ) are some known differential functions of ϕ. As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K(n,n)-Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009)1672-179910.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009)PHYADX0378-437110.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Lattice Boltzmann model for high-order nonlinear partial differential equations

Science.gov (United States)

Chai, Zhenhua; He, Nanzhong; Guo, Zhaoli; Shi, Baochang

2018-01-01

In this paper, a general lattice Boltzmann (LB) model is proposed for the high-order nonlinear partial differential equation with the form ∂tϕ +∑k=1mαk∂xkΠk(ϕ ) =0 (1 ≤k ≤m ≤6 ), αk are constant coefficients, Πk(ϕ ) are some known differential functions of ϕ . As some special cases of the high-order nonlinear partial differential equation, the classical (m)KdV equation, KdV-Burgers equation, K (n ,n ) -Burgers equation, Kuramoto-Sivashinsky equation, and Kawahara equation can be solved by the present LB model. Compared to the available LB models, the most distinct characteristic of the present model is to introduce some suitable auxiliary moments such that the correct moments of equilibrium distribution function can be achieved. In addition, we also conducted a detailed Chapman-Enskog analysis, and found that the high-order nonlinear partial differential equation can be correctly recovered from the proposed LB model. Finally, a large number of simulations are performed, and it is found that the numerical results agree with the analytical solutions, and usually the present model is also more accurate than the existing LB models [H. Lai and C. Ma, Sci. China Ser. G 52, 1053 (2009), 10.1007/s11433-009-0149-3; H. Lai and C. Ma, Phys. A (Amsterdam) 388, 1405 (2009), 10.1016/j.physa.2009.01.005] for high-order nonlinear partial differential equations.

Arbitrary vehicle steering characteristics with changing ratio rack and pinion transmission

Directory of Open Access Journals (Sweden)

András G Bendefy

2015-12-01

Full Text Available In order to achieve arbitrary steering characteristics at vehicles, a steering mechanism was developed, in which changing ratio rack and pinion connections have been applied. In contrary to a regular steering mechanism where only a single rack is used, two racks were applied in order to make arbitrary characteristics possible. The turning wheel’s required motion functions had to be defined first, thereafter could we determine the changing ratio rack and pinion geometry which produces this motion. First, a simplified two-dimensional mockup was created in order to study the difficulties and possibilities of a real construction. Later, a fully functional assembly was designed and manufactured to make further experiments possible.

An Arbitrary First Order Theory Can Be Represented by a Program: A Theorem

Science.gov (United States)

Hosheleva, Olga

1997-01-01

How can we represent knowledge inside a computer? For formalized knowledge, classical logic seems to be the most adequate tool. Classical logic is behind all formalisms of classical mathematics, and behind many formalisms used in Artificial Intelligence. There is only one serious problem with classical logic: due to the famous Godel's theorem, classical logic is algorithmically undecidable; as a result, when the knowledge is represented in the form of logical statements, it is very difficult to check whether, based on this statement, a given query is true or not. To make knowledge representations more algorithmic, a special field of logic programming was invented. An important portion of logic programming is algorithmically decidable. To cover knowledge that cannot be represented in this portion, several extensions of the decidable fragments have been proposed. In the spirit of logic programming, these extensions are usually introduced in such a way that even if a general algorithm is not available, good heuristic methods exist. It is important to check whether the already proposed extensions are sufficient, or further extensions is necessary. In the present paper, we show that one particular extension, namely, logic programming with classical negation, introduced by M. Gelfond and V. Lifschitz, can represent (in some reasonable sense) an arbitrary first order logical theory.

Numerov iteration method for second order integral-differential equation

International Nuclear Information System (INIS)

Zeng Fanan; Zhang Jiaju; Zhao Xuan

1987-01-01

In this paper, Numerov iterative method for second order integral-differential equation and system of equations are constructed. Numerical examples show that this method is better than direct method (Gauss elimination method) in CPU time and memoy requireing. Therefore, this method is an efficient method for solving integral-differential equation in nuclear physics

Pseudospectral collocation methods for fourth order differential equations

Science.gov (United States)

Malek, Alaeddin; Phillips, Timothy N.

1994-01-01

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Quaternion based generalization of Chern–Simons theories in arbitrary dimensions

Directory of Open Access Journals (Sweden)

Alessandro D'Adda

2017-08-01

Full Text Available A generalization of Chern–Simons gauge theory is formulated in any dimension and arbitrary gauge group where gauge fields and gauge parameters are differential forms of any degree. The quaternion algebra structure of this formulation is shown to be equivalent to a three Z2-gradings structure, thus clarifying the quaternion role in the previous formulation.

Differential and difference equations a comparison of methods of solution

CERN Document Server

Maximon, Leonard C

2016-01-01

This book, intended for researchers and graduate students in physics, applied mathematics and engineering, presents a detailed comparison of the important methods of solution for linear differential and difference equations - variation of constants, reduction of order, Laplace transforms and generating functions - bringing out the similarities as well as the significant differences in the respective analyses. Equations of arbitrary order are studied, followed by a detailed analysis for equations of first and second order. Equations with polynomial coefficients are considered and explicit solutions for equations with linear coefficients are given, showing significant differences in the functional form of solutions of differential equations from those of difference equations. An alternative method of solution involving transformation of both the dependent and independent variables is given for both differential and difference equations. A comprehensive, detailed treatment of Green’s functions and the associat...

Convergence acceleration for time-independent first-order PDE using optimal PNB-approximations

Energy Technology Data Exchange (ETDEWEB)

Holmgren, S.; Branden, H. [Uppsala Univ. (Sweden)

1996-12-31

We consider solving time-independent (steady-state) flow problems in 2D or 3D governed by hyperbolic or {open_quotes}almost hyperbolic{close_quotes} systems of partial differential equations (PDE). Examples of such PDE are the Euler and the Navier-Stokes equations. The PDE is discretized using a finite difference or finite volume scheme with arbitrary order of accuracy. If the matrix B describes the discretized differential operator and u denotes the approximate solution, the discrete problem is given by a large system of equations.

Discussion of “Attitude of Physicians Towards Automatic Alerting in Computerized Physician Order Entry Systems”

DEFF Research Database (Denmark)

Bates, D. W.; Baysari, M. T.; Dugas, M.

2013-01-01

With these comments on the paper “Attitude of Physicians Towards Automatic Alerting in Computerized Physician Order Entry Systems”, written by Martin Jung and coauthors, with Dr. Elske Ammenwerth as senior author, the journal wants to stimulate a broad discussion on computerized physician order...... entry systems. An international group of experts have been invited by the editor of Methods to comment on this paper. Each of the invited commentaries forms one section of this paper....

First- and Second-Order Full-Differential in Edge Analysis of Images

Directory of Open Access Journals (Sweden)

Dong-Mei Pu

2014-01-01

mathematics. We propose and reformulate them with a uniform definition framework. Based on our observation and analysis with the difference, we propose an algorithm to detect the edge from image. Experiments on Corel5K and PASCAL VOC 2007 are done to show the difference between the first order and the second order. After comparison with Canny operator and the proposed first-order differential, the main result is that the second-order differential has the better performance in analysis of changes of the context of images with good selection of control parameter.

LOOP- SIMULATION OF THE AUTOMATIC FREQUENCY CONTROL SUBSYSTEM OF A DIFFERENTIAL MINIMUM SHIFT KEYING RECEIVER

Science.gov (United States)

Davarian, F.

1994-01-01

The LOOP computer program was written to simulate the Automatic Frequency Control (AFC) subsystem of a Differential Minimum Shift Keying (DMSK) receiver with a bit rate of 2400 baud. The AFC simulated by LOOP is a first order loop configuration with a first order R-C filter. NASA has been investigating the concept of mobile communications based on low-cost, low-power terminals linked via geostationary satellites. Studies have indicated that low bit rate transmission is suitable for this application, particularly from the frequency and power conservation point of view. A bit rate of 2400 BPS is attractive due to its applicability to the linear predictive coding of speech. Input to LOOP includes the following: 1) the initial frequency error; 2) the double-sided loop noise bandwidth; 3) the filter time constants; 4) the amount of intersymbol interference; and 5) the bit energy to noise spectral density. LOOP output includes: 1) the bit number and the frequency error of that bit; 2) the computed mean of the frequency error; and 3) the standard deviation of the frequency error. LOOP is written in MS SuperSoft FORTRAN 77 for interactive execution and has been implemented on an IBM PC operating under PC DOS with a memory requirement of approximately 40K of 8 bit bytes. This program was developed in 1986.

A compact, multichannel, and low noise arbitrary waveform generator.

Science.gov (United States)

Govorkov, S; Ivanov, B I; Il'ichev, E; Meyer, H-G

2014-05-01

A new type of high functionality, fast, compact, and easy programmable arbitrary waveform generator for low noise physical measurements is presented. The generator provides 7 fast differential waveform channels with a maximum bandwidth up to 200 MHz frequency. There are 6 fast pulse generators on the generator board with 78 ps time resolution in both duration and delay, 3 of them with amplitude control. The arbitrary waveform generator is additionally equipped with two auxiliary slow 16 bit analog-to-digital converters and four 16 bit digital-to-analog converters for low frequency applications. Electromagnetic shields are introduced to the power supply, digital, and analog compartments and with a proper filter design perform more than 110 dB digital noise isolation to the output signals. All the output channels of the board have 50 Ω SubMiniature version A termination. The generator board is suitable for use as a part of a high sensitive physical equipment, e.g., fast read out and manipulation of nuclear magnetic resonance or superconducting quantum systems and any other application, which requires electromagnetic interference free fast pulse and arbitrary waveform generation.

A compact, multichannel, and low noise arbitrary waveform generator

International Nuclear Information System (INIS)

Govorkov, S.; Ivanov, B. I.; Il'ichev, E.; Meyer, H.-G.

2014-01-01

A new type of high functionality, fast, compact, and easy programmable arbitrary waveform generator for low noise physical measurements is presented. The generator provides 7 fast differential waveform channels with a maximum bandwidth up to 200 MHz frequency. There are 6 fast pulse generators on the generator board with 78 ps time resolution in both duration and delay, 3 of them with amplitude control. The arbitrary waveform generator is additionally equipped with two auxiliary slow 16 bit analog-to-digital converters and four 16 bit digital-to-analog converters for low frequency applications. Electromagnetic shields are introduced to the power supply, digital, and analog compartments and with a proper filter design perform more than 110 dB digital noise isolation to the output signals. All the output channels of the board have 50 Ω SubMiniature version A termination. The generator board is suitable for use as a part of a high sensitive physical equipment, e.g., fast read out and manipulation of nuclear magnetic resonance or superconducting quantum systems and any other application, which requires electromagnetic interference free fast pulse and arbitrary waveform generation

The Oscillation of a Class of the Fractional-Order Delay Differential Equations

Directory of Open Access Journals (Sweden)

Qianli Lu

2014-01-01

Full Text Available Several oscillation results are proposed including necessary and sufficient conditions for the oscillation of fractional-order delay differential equations with constant coefficients, the sufficient or necessary and sufficient conditions for the oscillation of fractional-order delay differential equations by analysis method, and the sufficient or necessary and sufficient conditions for the oscillation of delay partial differential equation with three different boundary conditions. For this, α-exponential function which is a kind of functions that play the same role of the classical exponential functions of fractional-order derivatives is used.

Automatic generation of anatomic characteristics from cerebral aneurysm surface models.

Science.gov (United States)

Neugebauer, M; Lawonn, K; Beuing, O; Preim, B

2013-03-01

Computer-aided research on cerebral aneurysms often depends on a polygonal mesh representation of the vessel lumen. To support a differentiated, anatomy-aware analysis, it is necessary to derive anatomic descriptors from the surface model. We present an approach on automatic decomposition of the adjacent vessels into near- and far-vessel regions and computation of the axial plane. We also exemplarily present two applications of the geometric descriptors: automatic computation of a unique vessel order and automatic viewpoint selection. Approximation methods are employed to analyze vessel cross-sections and the vessel area profile along the centerline. The resulting transition zones between near- and far- vessel regions are used as input for an optimization process to compute the axial plane. The unique vessel order is defined via projection into the plane space of the axial plane. The viewing direction for the automatic viewpoint selection is derived from the normal vector of the axial plane. The approach was successfully applied to representative data sets exhibiting a broad variability with respect to the configuration of their adjacent vessels. A robustness analysis showed that the automatic decomposition is stable against noise. A survey with 4 medical experts showed a broad agreement with the automatically defined transition zones. Due to the general nature of the underlying algorithms, this approach is applicable to most of the likely aneurysm configurations in the cerebral vasculature. Additional geometric information obtained during automatic decomposition can support correction in case the automatic approach fails. The resulting descriptors can be used for various applications in the field of visualization, exploration and analysis of cerebral aneurysms.

Arbitrary temporal shape pulsed fiber laser based on SPGD algorithm

Science.gov (United States)

Jiang, Min; Su, Rongtao; Zhang, Pengfei; Zhou, Pu

2018-06-01

A novel adaptive pulse shaping method for a pulsed master oscillator power amplifier fiber laser to deliver an arbitrary pulse shape is demonstrated. Numerical simulation has been performed to validate the feasibility of the scheme and provide meaningful guidance for the design of the algorithm control parameters. In the proof-of-concept experiment, information on the temporal property of the laser is exchanged and evaluated through a local area network, and the laser adjusted the parameters of the seed laser according to the monitored output of the system automatically. Various pulse shapes, including a rectangular shape, ‘M’ shape, and elliptical shape are achieved through experimental iterations.

Ultrasound speckle reduction based on fractional order differentiation.

Science.gov (United States)

Shao, Dangguo; Zhou, Ting; Liu, Fan; Yi, Sanli; Xiang, Yan; Ma, Lei; Xiong, Xin; He, Jianfeng

2017-07-01

Ultrasound images show a granular pattern of noise known as speckle that diminishes their quality and results in difficulties in diagnosis. To preserve edges and features, this paper proposes a fractional differentiation-based image operator to reduce speckle in ultrasound. An image de-noising model based on fractional partial differential equations with balance relation between k (gradient modulus threshold that controls the conduction) and v (the order of fractional differentiation) was constructed by the effective combination of fractional calculus theory and a partial differential equation, and the numerical algorithm of it was achieved using a fractional differential mask operator. The proposed algorithm has better speckle reduction and structure preservation than the three existing methods [P-M model, the speckle reducing anisotropic diffusion (SRAD) technique, and the detail preserving anisotropic diffusion (DPAD) technique]. And it is significantly faster than bilateral filtering (BF) in producing virtually the same experimental results. Ultrasound phantom testing and in vivo imaging show that the proposed method can improve the quality of an ultrasound image in terms of tissue SNR, CNR, and FOM values.

Automatic Optimization for Large-Scale Real-Time Coastal Water Simulation

Directory of Open Access Journals (Sweden)

Shunli Wang

2016-01-01

Full Text Available We introduce an automatic optimization approach for the simulation of large-scale coastal water. To solve the singular problem of water waves obtained with the traditional model, a hybrid deep-shallow-water model is estimated by using an automatic coupling algorithm. It can handle arbitrary water depth and different underwater terrain. As a certain feature of coastal terrain, coastline is detected with the collision detection technology. Then, unnecessary water grid cells are simplified by the automatic simplification algorithm according to the depth. Finally, the model is calculated on Central Processing Unit (CPU and the simulation is implemented on Graphics Processing Unit (GPU. We show the effectiveness of our method with various results which achieve real-time rendering on consumer-level computer.

Dispersion in a thermal plasma including arbitrary degeneracy and quantum recoil.

Science.gov (United States)

Melrose, D B; Mushtaq, A

2010-11-01

The longitudinal response function for a thermal electron gas is calculated including two quantum effects exactly, degeneracy, and the quantum recoil. The Fermi-Dirac distribution is expanded in powers of a parameter that is small in the nondegenerate limit and the response function is evaluated in terms of the conventional plasma dispersion function to arbitrary order in this parameter. The infinite sum is performed in terms of polylogarithms in the long-wavelength and quasistatic limits, giving results that apply for arbitrary degeneracy. The results are applied to the dispersion relations for Langmuir waves and to screening, reproducing known results in the nondegenerate and completely degenerate limits, and generalizing them to arbitrary degeneracy.

OSCILLATION BEHAVIOR OF SOLUTIONS FOR EVEN ORDER NEUTRAL FUNCTIONAL DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

T.Candan

2006-01-01

Even order neutral functional differential equations are considered. Sufficient conditions for the oscillation behavior of solutions for this differential equation are presented. The new results are presented and some examples are also given.

Low Dimensional Vessiot-Guldberg-Lie Algebras of Second-Order Ordinary Differential Equations

Directory of Open Access Journals (Sweden)

Rutwig Campoamor-Stursberg

2016-03-01

Full Text Available A direct approach to non-linear second-order ordinary differential equations admitting a superposition principle is developed by means of Vessiot-Guldberg-Lie algebras of a dimension not exceeding three. This procedure allows us to describe generic types of second-order ordinary differential equations subjected to some constraints and admitting a given Lie algebra as Vessiot-Guldberg-Lie algebra. In particular, well-known types, such as the Milne-Pinney or Kummer-Schwarz equations, are recovered as special cases of this classification. The analogous problem for systems of second-order differential equations in the real plane is considered for a special case that enlarges the generalized Ermakov systems.

Higher order multi-term time-fractional partial differential equations involving Caputo-Fabrizio derivative

OpenAIRE

Erkinjon Karimov; Sardor Pirnafasov

2017-01-01

In this work we discuss higher order multi-term partial differential equation (PDE) with the Caputo-Fabrizio fractional derivative in time. Using method of separation of variables, we reduce fractional order partial differential equation to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Nonclassical Symmetries for Nonlinear Partial Differential Equations via Compatibility

International Nuclear Information System (INIS)

El-Sabbagh, Mostafa F.; Ahmad, Ali T.

2011-01-01

The determining equations for the nonclassical symmetry reductions of nonlinear partial differential equations with arbitrary order can be obtained by requiring the compatibility between the original equations and the invariant surface conditions. The (2+1)-dimensional shallow water wave equation, Boussinesq equation, and the dispersive wave equations in shallow water serve as examples illustrating how compatibility leads quickly and easily to the determining equations for their nonclassical symmetries. (general)

Differential geometry of group lattices

International Nuclear Information System (INIS)

Dimakis, Aristophanes; Mueller-Hoissen, Folkert

2003-01-01

In a series of publications we developed ''differential geometry'' on discrete sets based on concepts of noncommutative geometry. In particular, it turned out that first-order differential calculi (over the algebra of functions) on a discrete set are in bijective correspondence with digraph structures where the vertices are given by the elements of the set. A particular class of digraphs are Cayley graphs, also known as group lattices. They are determined by a discrete group G and a finite subset S. There is a distinguished subclass of ''bicovariant'' Cayley graphs with the property ad(S)S subset of S. We explore the properties of differential calculi which arise from Cayley graphs via the above correspondence. The first-order calculi extend to higher orders and then allow us to introduce further differential geometric structures. Furthermore, we explore the properties of ''discrete'' vector fields which describe deterministic flows on group lattices. A Lie derivative with respect to a discrete vector field and an inner product with forms is defined. The Lie-Cartan identity then holds on all forms for a certain subclass of discrete vector fields. We develop elements of gauge theory and construct an analog of the lattice gauge theory (Yang-Mills) action on an arbitrary group lattice. Also linear connections are considered and a simple geometric interpretation of the torsion is established. By taking a quotient with respect to some subgroup of the discrete group, generalized differential calculi associated with so-called Schreier diagrams are obtained

Higher order multi-term time-fractional partial differential equations involving Caputo-Fabrizio derivative

Directory of Open Access Journals (Sweden)

Erkinjon Karimov

2017-10-01

Full Text Available In this work we discuss higher order multi-term partial differential equation (PDE with the Caputo-Fabrizio fractional derivative in time. Using method of separation of variables, we reduce fractional order partial differential equation to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Automatic identification of inertial sensor placement on human body segments during walking

NARCIS (Netherlands)

Weenk, D.; van Beijnum, Bernhard J.F.; Baten, Christian T.M.; Hermens, Hermanus J.; Veltink, Petrus H.

2013-01-01

We present a novel method for the automatic identification of inertial sensors on human body segments during walking. This method allows the user to place (wireless) inertial sensors on arbitrary body segments. Next, the user walks for just a few seconds and the segment to which each sensor is

Arbitrary waveform modulated pulse EPR at 200 GHz

Science.gov (United States)

Kaminker, Ilia; Barnes, Ryan; Han, Songi

2017-06-01

We report here on the implementation of arbitrary waveform generation (AWG) capabilities at ∼200 GHz into an Electron Paramagnetic Resonance (EPR) and Dynamic Nuclear Polarization (DNP) instrument platform operating at 7 T. This is achieved with the integration of a 1 GHz, 2 channel, digital to analog converter (DAC) board that enables the generation of coherent arbitrary waveforms at Ku-band frequencies with 1 ns resolution into an existing architecture of a solid state amplifier multiplier chain (AMC). This allows for the generation of arbitrary phase- and amplitude-modulated waveforms at 200 GHz with >150 mW power. We find that the non-linearity of the AMC poses significant difficulties in generating amplitude-modulated pulses at 200 GHz. We demonstrate that in the power-limited regime of ω1 10 MHz) spin manipulation in incoherent (inversion), as well as coherent (echo formation) experiments. Highlights include the improvement by one order of magnitude in inversion bandwidth compared to that of conventional rectangular pulses, as well as a factor of two in improvement in the refocused echo intensity at 200 GHz.

Analysis of Caputo Impulsive Fractional Order Differential Equations with Applications

Directory of Open Access Journals (Sweden)

Lakshman Mahto

2013-01-01

Full Text Available We use Sadovskii's fixed point method to investigate the existence and uniqueness of solutions of Caputo impulsive fractional differential equations of order with one example of impulsive logistic model and few other examples as well. We also discuss Caputo impulsive fractional differential equations with finite delay. The results proven are new and compliment the existing one.

Dispersion in thermal plasma including arbitrary degeneracy and quantum recoil

International Nuclear Information System (INIS)

Mushtaq, A.; Melrose, D.B.

2012-01-01

The longitudinal response function for a thermal electron gas was calculated including two quantum effects exactly, degeneracy and the quantum recoil. The Fermi-Dirac distribution was expanded in powers of a parameter that is small in the non-degenerate limit and the response function was evaluated in terms of the conventional plasma dispersion function to arbitrary order in this parameter. The infinite sum was performed in terms of poly logarithms in the long-wavelength and quasi-static limits, giving results that apply for arbitrary degeneracy. The results were applied to the dispersion relations for Langmuir waves and to screening, reproducing known results in the non-degenerate and completely degenerate limits], and generalizing them to arbitrary degeneracy. The occupation number for the completely degenerate limit is shown. The importance of the results regarding to semiconductor plasmas were highlighted. (orig./A.B.)

Finding higher order Darboux polynomials for a family of rational first order ordinary differential equations

Science.gov (United States)

Avellar, J.; Claudino, A. L. G. C.; Duarte, L. G. S.; da Mota, L. A. C. P.

2015-10-01

For the Darbouxian methods we are studying here, in order to solve first order rational ordinary differential equations (1ODEs), the most costly (computationally) step is the finding of the needed Darboux polynomials. This can be so grave that it can render the whole approach unpractical. Hereby we introduce a simple heuristics to speed up this process for a class of 1ODEs.

A Variable Order Fractional Differential-Based Texture Enhancement Algorithm with Application in Medical Imaging.

Directory of Open Access Journals (Sweden)

Qiang Yu

Full Text Available Texture enhancement is one of the most important techniques in digital image processing and plays an essential role in medical imaging since textures discriminate information. Most image texture enhancement techniques use classical integral order differential mask operators or fractional differential mask operators using fixed fractional order. These masks can produce excessive enhancement of low spatial frequency content, insufficient enhancement of large spatial frequency content, and retention of high spatial frequency noise. To improve upon existing approaches of texture enhancement, we derive an improved Variable Order Fractional Centered Difference (VOFCD scheme which dynamically adjusts the fractional differential order instead of fixing it. The new VOFCD technique is based on the second order Riesz fractional differential operator using a Lagrange 3-point interpolation formula, for both grey scale and colour image enhancement. We then use this method to enhance photographs and a set of medical images related to patients with stroke and Parkinson's disease. The experiments show that our improved fractional differential mask has a higher signal to noise ratio value than the other fractional differential mask operators. Based on the corresponding quantitative analysis we conclude that the new method offers a superior texture enhancement over existing methods.

Numerical solution of second-order stochastic differential equations with Gaussian random parameters

Directory of Open Access Journals (Sweden)

Rahman Farnoosh

2014-07-01

Full Text Available In this paper, we present the numerical solution of ordinary differential equations (or SDEs, from each orderespecially second-order with time-varying and Gaussian random coefficients. We indicate a complete analysisfor second-order equations in specially case of scalar linear second-order equations (damped harmonicoscillators with additive or multiplicative noises. Making stochastic differential equations system from thisequation, it could be approximated or solved numerically by different numerical methods. In the case oflinear stochastic differential equations system by Computing fundamental matrix of this system, it could becalculated based on the exact solution of this system. Finally, this stochastic equation is solved by numericallymethod like E.M. and Milstein. Also its Asymptotic stability and statistical concepts like expectationand variance of solutions are discussed.

A Domain Specific Embedded Language in C++ for Automatic Differentiation, Projection, Integration and Variational Formulations

Directory of Open Access Journals (Sweden)

Christophe Prud'homme

2006-01-01

Full Text Available In this article, we present a domain specific embedded language in C++ that can be used in various contexts such as numerical projection onto a functional space, numerical integration, variational formulations and automatic differentiation. Albeit these tools operate in different ways, the language overcomes this difficulty by decoupling expression constructions from evaluation. The language is implemented using expression templates and meta-programming techniques and uses various Boost libraries. The language is exercised on a number of non-trivial examples and a benchmark presents the performance behavior on a few test problems.

ADM For Solving Linear Second-Order Fredholm Integro-Differential Equations

Science.gov (United States)

Karim, Mohd F.; Mohamad, Mahathir; Saifullah Rusiman, Mohd; Che-Him, Norziha; Roslan, Rozaini; Khalid, Kamil

2018-04-01

In this paper, we apply Adomian Decomposition Method (ADM) as numerically analyse linear second-order Fredholm Integro-differential Equations. The approximate solutions of the problems are calculated by Maple package. Some numerical examples have been considered to illustrate the ADM for solving this equation. The results are compared with the existing exact solution. Thus, the Adomian decomposition method can be the best alternative method for solving linear second-order Fredholm Integro-Differential equation. It converges to the exact solution quickly and in the same time reduces computational work for solving the equation. The result obtained by ADM shows the ability and efficiency for solving these equations.

A Four-Stage Fifth-Order Trigonometrically Fitted Semi-Implicit Hybrid Method for Solving Second-Order Delay Differential Equations

Directory of Open Access Journals (Sweden)

Sufia Zulfa Ahmad

2016-01-01

Full Text Available We derived a two-step, four-stage, and fifth-order semi-implicit hybrid method which can be used for solving special second-order ordinary differential equations. The method is then trigonometrically fitted so that it is suitable for solving problems which are oscillatory in nature. The methods are then used for solving oscillatory delay differential equations. Numerical results clearly show the efficiency of the new method when compared to the existing explicit and implicit methods in the scientific literature.

High-order quantum algorithm for solving linear differential equations

International Nuclear Information System (INIS)

Berry, Dominic W

2014-01-01

Linear differential equations are ubiquitous in science and engineering. Quantum computers can simulate quantum systems, which are described by a restricted type of linear differential equations. Here we extend quantum simulation algorithms to general inhomogeneous sparse linear differential equations, which describe many classical physical systems. We examine the use of high-order methods (where the error over a time step is a high power of the size of the time step) to improve the efficiency. These provide scaling close to Δt 2 in the evolution time Δt. As with other algorithms of this type, the solution is encoded in amplitudes of the quantum state, and it is possible to extract global features of the solution. (paper)

Attentional Bias for Pain and Sex, and Automatic Appraisals of Sexual Penetration: Differential Patterns in Dyspareunia vs Vaginismus?

Science.gov (United States)

Melles, Reinhilde J; Dewitte, Marieke D; Ter Kuile, Moniek M; Peters, Madelon M L; de Jong, Peter J

2016-08-01

Current information processing models propose that heightened attention bias for sex-related threats (eg, pain) and lowered automatic incentive processes ("wanting") may play an important role in the impairment of sexual arousal and the development of sexual dysfunctions such as genitopelvic pain/penetration disorder (GPPPD). Differential threat and incentive processing may also help explain the stronger persistence of coital avoidance in women with vaginismus compared to women with dyspareunia. As the first aim, we tested if women with GPPPD show (1) heightened attention for pain and sex, and (2) heightened threat and lower incentive associations with sexual penetration. Second, we examined whether the stronger persistence of coital avoidance in vaginismus vs dyspareunia might be explained by a stronger attentional bias or more dysfunctional automatic threat/incentive associations. Women with lifelong vaginismus (n = 37), dyspareunia (n = 29), and a no-symptoms comparison group (n = 51) completed a visual search task to assess attentional bias, and single target implicit-association tests to measure automatic sex-threat and sex-wanting associations. There were no group differences in attentional bias or automatic associations. Correlational analysis showed that slowed detection of sex stimuli and stronger automatic threat associations were related to lowered sexual arousal. The findings do not corroborate the view that attentional bias for pain or sex contributes to coital pain, or that differences in coital avoidance may be explained by differences in attentional bias or automatic threat/incentive associations. However, the correlational findings are consistent with the view that automatic threat associations and impaired attention for sex stimuli may interfere with the generation of sexual arousal. Copyright © 2016 International Society for Sexual Medicine. Published by Elsevier Inc. All rights reserved.

A NEW OSCILLATION CRITERION FOR FIRST ORDER NEUTRAL DELAY DIFFERENTIAL EQUATION

Institute of Scientific and Technical Information of China (English)

无

2006-01-01

In this paper, a new sufficient condition for the oscillation of all solutions of first order neutral delay differential equations is obtained. Secondly, the result can also be extended to a general neutral differential equation, and many known results in the literatures are improved.

Asymptotic behavior and stability of second order neutral delay differential equations

NARCIS (Netherlands)

Chen, G.L.; van Gaans, O.W.; Verduyn Lunel, Sjoerd

2014-01-01

We study the asymptotic behavior of a class of second order neutral delay differential equations by both a spectral projection method and an ordinary differential equation method approach. We discuss the relation of these two methods and illustrate some features using examples. Furthermore, a fixed

Seeley-Gilkey coefficients for fourth-order operators on Riemannian manifold

International Nuclear Information System (INIS)

Gusynin, V.P.

1990-01-01

The covariant pseudodifferential-operator method of Widom is developed for computing the coefficients in the heat kernel expansion. It allows one to calculate Seeley-Gilkey coefficients for both minimal and nonminimal differential operators acting on a vector bundle over a riemannian manifold. The coefficients for the fourth-order minimal operators in arbitrary dimensions of space are calculated. In contrast to the second-order operators the coefficients for the fourth-order (and higher) operators turn out to be essentially dependent on the space dimension. The algorithmic character of the method allows one to calculate the coefficients by computer using an analytical calculation system. The method also permits a simple generalization to manifolds with torsion and supermanifolds. (orig.)

On the Liouvillian solution of second-order linear differential equations and algebraic invariant curves

International Nuclear Information System (INIS)

Man, Yiu-Kwong

2010-01-01

In this communication, we present a method for computing the Liouvillian solution of second-order linear differential equations via algebraic invariant curves. The main idea is to integrate Kovacic's results on second-order linear differential equations with the Prelle-Singer method for computing first integrals of differential equations. Some examples on using this approach are provided. (fast track communication)

Double Hopf bifurcation in delay differential equations

Directory of Open Access Journals (Sweden)

Redouane Qesmi

2014-07-01

Full Text Available The paper addresses the computation of elements of double Hopf bifurcation for retarded functional differential equations (FDEs with parameters. We present an efficient method for computing, simultaneously, the coefficients of center manifolds and normal forms, in terms of the original FDEs, associated with the double Hopf singularity up to an arbitrary order. Finally, we apply our results to a nonlinear model with periodic delay. This shows the applicability of the methodology in the study of delay models arising in either natural or technological problems.

Self-consistent field theory simulations of polymers on arbitrary domains

Energy Technology Data Exchange (ETDEWEB)

Ouaknin, Gaddiel, E-mail: gaddielouaknin@umail.ucsb.edu [Department of Mechanical Engineering, University of California, Santa Barbara, CA 93106-5070 (United States); Laachi, Nabil; Delaney, Kris [Materials Research Laboratory, University of California, Santa Barbara, CA 93106-5080 (United States); Fredrickson, Glenn H. [Materials Research Laboratory, University of California, Santa Barbara, CA 93106-5080 (United States); Department of Chemical Engineering, University of California, Santa Barbara, CA 93106-5080 (United States); Department of Materials, University of California, Santa Barbara, CA 93106-5050 (United States); Gibou, Frederic [Department of Mechanical Engineering, University of California, Santa Barbara, CA 93106-5070 (United States); Department of Computer Science, University of California, Santa Barbara, CA 93106-5110 (United States)

2016-12-15

We introduce a framework for simulating the mesoscale self-assembly of block copolymers in arbitrary confined geometries subject to Neumann boundary conditions. We employ a hybrid finite difference/volume approach to discretize the mean-field equations on an irregular domain represented implicitly by a level-set function. The numerical treatment of the Neumann boundary conditions is sharp, i.e. it avoids an artificial smearing in the irregular domain boundary. This strategy enables the study of self-assembly in confined domains and enables the computation of physically meaningful quantities at the domain interface. In addition, we employ adaptive grids encoded with Quad-/Oc-trees in parallel to automatically refine the grid where the statistical fields vary rapidly as well as at the boundary of the confined domain. This approach results in a significant reduction in the number of degrees of freedom and makes the simulations in arbitrary domains using effective boundary conditions computationally efficient in terms of both speed and memory requirement. Finally, in the case of regular periodic domains, where pseudo-spectral approaches are superior to finite differences in terms of CPU time and accuracy, we use the adaptive strategy to store chain propagators, reducing the memory footprint without loss of accuracy in computed physical observables.

Asymptotically stable fourth-order accurate schemes for the diffusion equation on complex shapes

International Nuclear Information System (INIS)

Abarbanel, S.; Ditkowski, A.

1997-01-01

An algorithm which solves the multidimensional diffusion equation on complex shapes to fourth-order accuracy and is asymptotically stable in time is presented. This bounded-error result is achieved by constructing, on a rectangular grid, a differentiation matrix whose symmetric part is negative definite. The differentiation matrix accounts for the Dirichlet boundary condition by imposing penalty-like terms. Numerical examples in 2-D show that the method is effective even where standard schemes, stable by traditional definitions, fail. The ability of the paradigm to be applied to arbitrary geometric domains is an important feature of the algorithm. 5 refs., 14 figs

Robust automatic high resolution segmentation of SOFC anode porosity in 3D

DEFF Research Database (Denmark)

Jørgensen, Peter Stanley; Bowen, Jacob R.

2008-01-01

Routine use of 3D characterization of SOFCs by focused ion beam (FIB) serial sectioning is generally restricted by the time consuming task of manually delineating structures within each image slice. We apply advanced image analysis algorithms to automatically segment the porosity phase of an SOFC...... anode in 3D. The technique is based on numerical approximations to partial differential equations to evolve a 3D surface to the desired phase boundary. Vector fields derived from the experimentally acquired data are used as the driving force. The automatic segmentation compared to manual delineation...... reveals and good correspondence and the two approaches are quantitatively compared. It is concluded that the. automatic approach is more robust, more reproduceable and orders of magnitude quicker than manual segmentation of SOFC anode porosity for subsequent quantitative 3D analysis. Lastly...

Iterative oscillation results for second-order differential equations with advanced argument

Directory of Open Access Journals (Sweden)

Irena Jadlovska

2017-07-01

Full Text Available This article concerns the oscillation of solutions to a linear second-order differential equation with advanced argument. Sufficient oscillation conditions involving limit inferior are given which essentially improve known results. We base our technique on the iterative construction of solution estimates and some of the recent ideas developed for first-order advanced differential equations. We demonstrate the advantage of our results on Euler-type advanced equation. Using MATLAB software, a comparison of the effectiveness of newly obtained criteria as well as the necessary iteration length in particular cases are discussed.

Content-aware automatic cropping for consumer photos

Science.gov (United States)

Tang, Hao; Tretter, Daniel; Lin, Qian

2013-03-01

Consumer photos are typically authored once, but need to be retargeted for reuse in various situations. These include printing a photo on different size paper, changing the size and aspect ratio of an embedded photo to accommodate the dynamic content layout of web pages or documents, adapting a large photo for browsing on small displays such as mobile phone screens, and improving the aesthetic quality of a photo that was badly composed at the capture time. In this paper, we propose a novel, effective, and comprehensive content-aware automatic cropping (hereafter referred to as "autocrop") method for consumer photos to achieve the above purposes. Our autocrop method combines the state-of-the-art context-aware saliency detection algorithm, which aims to infer the likely intent of the photographer, and the "branch-and-bound" efficient subwindow search optimization technique, which seeks to locate the globally optimal cropping rectangle in a fast manner. Unlike most current autocrop methods, which can only crop a photo into an arbitrary rectangle, our autocrop method can automatically crop a photo into either a rectangle of arbitrary dimensions or a rectangle of the desired aspect ratio specified by the user. The aggressiveness of the cropping operation may be either automatically determined by the method or manually indicated by the user with ease. In addition, our autocrop method is extended to support the cropping of a photo into non-rectangular shapes such as polygons of any number of sides. It may also be potentially extended to return multiple cropping suggestions, which will enable the creation of new photos to enrich the original photo collections. Our experimental results show that the proposed autocrop method in this paper can generate high-quality crops for consumer photos of various types.

The exact equation of motion of a simple pendulum of arbitrary amplitude: a hypergeometric approach

International Nuclear Information System (INIS)

Qureshi, M I; Rafat, M; Azad, S Ismail

2010-01-01

The motion of a simple pendulum of arbitrary amplitude is usually treated by approximate methods. By using generalized hypergeometric functions, it is however possible to solve the problem exactly. In this paper, we provide the exact equation of motion of a simple pendulum of arbitrary amplitude. A new and exact expression for the time of swinging of a simple pendulum from the vertical position to an arbitrary angular position θ is given by equation (3.10). The time period of such a pendulum is also exactly expressible in terms of hypergeometric functions. The exact expressions thus obtained are used to plot the graphs that compare the exact time period T(θ 0 ) with the time period T(0) (based on simple harmonic approximation). We also compare the relative difference between T(0) and T(θ 0 ) found from the exact equation of motion with the usual perturbation theory estimate. The treatment is intended for graduate students, who have acquired some familiarity with the hypergeometric functions. This approach may also be profitably used by specialists who encounter during their investigations nonlinear differential equations similar in form to the pendulum equation. Such nonlinear differential equations could arise in diverse fields, such as acoustic vibrations, oscillations in small molecules, turbulence and electronic filters, among others.

Hopf bifurcation formula for first order differential-delay equations

Science.gov (United States)

Rand, Richard; Verdugo, Anael

2007-09-01

This work presents an explicit formula for determining the radius of a limit cycle which is born in a Hopf bifurcation in a class of first order constant coefficient differential-delay equations. The derivation is accomplished using Lindstedt's perturbation method.

Security Analysis of 7-Round MISTY1 against Higher Order Differential Attacks

Science.gov (United States)

Tsunoo, Yukiyasu; Saito, Teruo; Shigeri, Maki; Kawabata, Takeshi

MISTY1 is a 64-bit block cipher that has provable security against differential and linear cryptanalysis. MISTY1 is one of the algorithms selected in the European NESSIE project, and it has been recommended for Japanese e-Government ciphers by the CRYPTREC project. This paper shows that higher order differential attacks can be successful against 7-round versions of MISTY1 with FL functions. The attack on 7-round MISTY1 can recover a partial subkey with a data complexity of 254.1 and a computational complexity of 2120.8, which signifies the first successful attack on 7-round MISTY1 with no limitation such as a weak key. This paper also evaluates the complexity of this higher order differential attack on MISTY1 in which the key schedule is replaced by a pseudorandom function. It is shown that resistance to the higher order differential attack is not substantially improved even in 7-round MISTY1 in which the key schedule is replaced by a pseudorandom function.

Rational approximations to solutions of linear differential equations.

Science.gov (United States)

Chudnovsky, D V; Chudnovsky, G V

1983-08-01

Rational approximations of Padé and Padé type to solutions of differential equations are considered. One of the main results is a theorem stating that a simultaneous approximation to arbitrary solutions of linear differential equations over C(x) cannot be "better" than trivial ones implied by the Dirichlet box principle. This constitutes, in particular, the solution in the linear case of Kolchin's problem that the "Roth's theorem" holds for arbitrary solutions of algebraic differential equations. Complete effective proofs for several valuations are presented based on the Wronskian methods and graded subrings of Picard-Vessiot extensions.

High order Poisson Solver for unbounded flows

DEFF Research Database (Denmark)

Hejlesen, Mads Mølholm; Rasmussen, Johannes Tophøj; Chatelain, Philippe

2015-01-01

This paper presents a high order method for solving the unbounded Poisson equation on a regular mesh using a Green’s function solution. The high order convergence was achieved by formulating mollified integration kernels, that were derived from a filter regularisation of the solution field....... The method was implemented on a rectangular domain using fast Fourier transforms (FFT) to increase computational efficiency. The Poisson solver was extended to directly solve the derivatives of the solution. This is achieved either by including the differential operator in the integration kernel...... the equations of fluid mechanics as an example, but can be used in many physical problems to solve the Poisson equation on a rectangular unbounded domain. For the two-dimensional case we propose an infinitely smooth test function which allows for arbitrary high order convergence. Using Gaussian smoothing...

Seeley-Gilkey coefficients for the fourth-order operators on a Riemannian manifold

International Nuclear Information System (INIS)

Gusynin, V.P.

1989-01-01

A new covariant method for computing the coefficients in the heat kernel expansion is suggested. It allows one to calculate Seeley-Gilkey coefficients for both minimal and nonminimal differential operators acting on a vector bundle over a Riemannian manifold. The coefficients for the fourth-order minimal operators in arbitrary dimension of the space are calculated. In contrast to the second-order operators the coefficients for the fourth-order (and higher) operators turn out to be essentially dependent on the space dimension. The algorithmic character of the method suggested allows one to calculate coefficients by computer using the analytical calculation system. 19 refs.; 1 fig

From differential to difference equations for first order ODEs

Science.gov (United States)

Freed, Alan D.; Walker, Kevin P.

1991-01-01

When constructing an algorithm for the numerical integration of a differential equation, one should first convert the known ordinary differential equation (ODE) into an ordinary difference equation. Given this difference equation, one can develop an appropriate numerical algorithm. This technical note describes the derivation of two such ordinary difference equations applicable to a first order ODE. The implicit ordinary difference equation has the same asymptotic expansion as the ODE itself, whereas the explicit ordinary difference equation has an asymptotic that is similar in structure but different in value when compared with that of the ODE.

Nonlinear singular perturbation problems of arbitrary real orders

International Nuclear Information System (INIS)

Bijura, Angelina M.

2003-10-01

Higher order asymptotic solutions of singularly perturbed nonlinear fractional integral and derivatives of order 1/2 are investigated. It is particularly shown that whilst certain asymptotic expansions are applied successfully to linear equations and particular nonlinear problems, the standard formal asymptotic expansion is appropriate for the general class of nonlinear equations. This theory is then generalised to the general equation (of order β, 0 < β < 1). (author)

Riemann-Christoffel Tensor in Differential Geometry of Fractional Order Application to Fractal Space-Time

Science.gov (United States)

Jumarie, Guy

2013-04-01

By using fractional differences, one recently proposed an alternative to the formulation of fractional differential calculus, of which the main characteristics is a new fractional Taylor series and its companion Rolle's formula which apply to non-differentiable functions. The key is that now we have at hand a differential increment of fractional order which can be manipulated exactly like in the standard Leibniz differential calculus. Briefly the fractional derivative is the quotient of fractional increments. It has been proposed that this calculus can be used to construct a differential geometry on manifold of fractional order. The present paper, on the one hand, refines the framework, and on the other hand, contributes some new results related to arc length of fractional curves, area on fractional differentiable manifold, covariant fractal derivative, Riemann-Christoffel tensor of fractional order, fractional differential equations of fractional geodesic, strip modeling of fractal space time and its relation with Lorentz transformation. The relation with Nottale's fractal space-time theory then appears in quite a natural way.

Mathematic model analysis of Gaussian beam propagation through an arbitrary thickness random phase screen.

Science.gov (United States)

Tian, Yuzhen; Guo, Jin; Wang, Rui; Wang, Tingfeng

2011-09-12

In order to research the statistical properties of Gaussian beam propagation through an arbitrary thickness random phase screen for adaptive optics and laser communication application in the laboratory, we establish mathematic models of statistical quantities, which are based on the Rytov method and the thin phase screen model, involved in the propagation process. And the analytic results are developed for an arbitrary thickness phase screen based on the Kolmogorov power spectrum. The comparison between the arbitrary thickness phase screen and the thin phase screen shows that it is more suitable for our results to describe the generalized case, especially the scintillation index.

Hyers-Ulam stability for second-order linear differential equations with boundary conditions

Directory of Open Access Journals (Sweden)

Pasc Gavruta

2011-06-01

Full Text Available We prove the Hyers-Ulam stability of linear differential equations of second-order with boundary conditions or with initial conditions. That is, if y is an approximate solution of the differential equation $y''+ eta (x y = 0$ with $y(a = y(b =0$, then there exists an exact solution of the differential equation, near y.

Second Order Impulsive Retarded Differential Inclusions with Nonlocal Conditions

Directory of Open Access Journals (Sweden)

Hernán R. Henríquez

2014-01-01

Full Text Available In this work we establish some existence results for abstract second order Cauchy problems modeled by a retarded differential inclusion involving nonlocal and impulsive conditions. Our results are obtained by using fixed point theory for the measure of noncompactness.

Numerical solution of distributed order fractional differential equations

Science.gov (United States)

Katsikadelis, John T.

2014-02-01

In this paper a method for the numerical solution of distributed order FDEs (fractional differential equations) of a general form is presented. The method applies to both linear and nonlinear equations. The Caputo type fractional derivative is employed. The distributed order FDE is approximated with a multi-term FDE, which is then solved by adjusting appropriately the numerical method developed for multi-term FDEs by Katsikadelis. Several example equations are solved and the response of mechanical systems described by such equations is studied. The convergence and the accuracy of the method for linear and nonlinear equations are demonstrated through well corroborated numerical results.

Oscillation criteria for third order nonlinear delay differential equations with damping

Directory of Open Access Journals (Sweden)

Said R. Grace

2015-01-01

Full Text Available This note is concerned with the oscillation of third order nonlinear delay differential equations of the form \\[\\label{*} \\left( r_{2}(t\\left( r_{1}(ty^{\\prime}(t\\right^{\\prime}\\right^{\\prime}+p(ty^{\\prime}(t+q(tf(y(g(t=0.\\tag{\\(\\ast\\}\\] In the papers [A. Tiryaki, M. F. Aktas, Oscillation criteria of a certain class of third order nonlinear delay differential equations with damping, J. Math. Anal. Appl. 325 (2007, 54-68] and [M. F. Aktas, A. Tiryaki, A. Zafer, Oscillation criteria for third order nonlinear functional differential equations, Applied Math. Letters 23 (2010, 756-762], the authors established some sufficient conditions which insure that any solution of equation (\\(\\ast\\ oscillates or converges to zero, provided that the second order equation \\[\\left( r_{2}(tz^{\\prime }(t\\right^{\\prime}+\\left(p(t/r_{1}(t\\right z(t=0\\tag{\\(\\ast\\ast\\}\\] is nonoscillatory. Here, we shall improve and unify the results given in the above mentioned papers and present some new sufficient conditions which insure that any solution of equation (\\(\\ast\\ oscillates if equation (\\(\\ast\\ast\\ is nonoscillatory. We also establish results for the oscillation of equation (\\(\\ast\\ when equation (\\(\\ast\\ast\\ is oscillatory.

Finite difference method and algebraic polynomial interpolation for numerically solving Poisson's equation over arbitrary domains

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.

Third-order nonlinear differential operators preserving invariant subspaces of maximal dimension

International Nuclear Information System (INIS)

Qu Gai-Zhu; Zhang Shun-Li; Li Yao-Long

2014-01-01

In this paper, third-order nonlinear differential operators are studied. It is shown that they are quadratic forms when they preserve invariant subspaces of maximal dimension. A complete description of third-order quadratic operators with constant coefficients is obtained. One example is given to derive special solutions for evolution equations with third-order quadratic operators. (general)

On oscillation of second-order linear ordinary differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, A.; Šremr, Jiří

2011-01-01

Roč. 54, - (2011), s. 69-81 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z10190503 Keywords : linear second-order ordinary differential equation * Kamenev theorem * oscillation Subject RIV: BA - General Mathematics http://www.rmi.ge/jeomj/memoirs/vol54/abs54-4.htm

Growth of meromorphic solutions of higher-order linear differential equations

Directory of Open Access Journals (Sweden)

Wenjuan Chen

2009-01-01

Full Text Available In this paper, we investigate the higher-order linear differential equations with meromorphic coefficients. We improve and extend a result of M.S. Liu and C.L. Yuan, by using the estimates for the logarithmic derivative of a transcendental meromorphic function due to Gundersen, and the extended Winman-Valiron theory which proved by J. Wang and H.X. Yi. In addition, we also consider the nonhomogeneous linear differential equations.

On the entropy of random surfaces with arbitrary genus

International Nuclear Information System (INIS)

Kostov, I.K.; Krzywicki, A.

1987-01-01

We calculate the susceptibility critical exponent γ for Polyakov random surfaces with arbitrary genus, using the Liouville theory to one-loop order. Some rigorous results obtained for special dimensionalities in a discrete version of the model are also noted. In all cases γ grows linearly with the genus of the surface. (orig.)

Effective Hamiltonian for 2-dimensional arbitrary spin Ising model

International Nuclear Information System (INIS)

Sznajd, J.; Polska Akademia Nauk, Wroclaw. Inst. Niskich Temperatur i Badan Strukturalnych)

1983-08-01

The method of the reduction of the generalized arbitrary-spin 2-dimensional Ising model to spin-half Ising model is presented. The method is demonstrated in detail by calculating the effective interaction constants to the third order in cumulant expansion for the triangular spin-1 Ising model (the Blume-Emery-Griffiths model). (author)

Quantum electrodynamics with arbitrary charge on a noncommutative space

International Nuclear Information System (INIS)

Zhou Wanping; Long Zhengwen; Cai Shaohong

2009-01-01

Using the Seiberg-Witten map, we obtain a quantum electrodynamics on a noncommutative space, which has arbitrary charge and keep the gauge invariance to at the leading order in theta. The one-loop divergence and Compton scattering are reinvestigated. The noncommutative effects are larger than those in ordinary noncommutative quantum electrodynamics. (authors)

Reduction of observer order by differentiation, almost controllability subspace covers and minimal order PID-observers

NARCIS (Netherlands)

TRENTELMAN, HL

1984-01-01

This note generalizes the geometric theory around minimal and reduced order observers to the situation in which differentiation of certain components of the observed output is allowed. A geometric theory involving the notion of PID-observer is introduced, using the concept of almost complementary

All-optical temporal fractional order differentiator using an in-fiber ellipsoidal air-microcavity

Science.gov (United States)

Zhang, Lihong; Sun, Shuqian; Li, Ming; Zhu, Ninghua

2017-12-01

An all-optical temporal fractional order differentiator with ultrabroad bandwidth (~1.6 THz) and extremely simple fabrication is proposed and experimentally demonstrated based on an in-fiber ellipsoidal air-microcavity. The ellipsoidal air-microcavity is fabricated by splicing a single mode fiber (SMF) and a photonic crystal fiber (PCF) together using a simple arc-discharging technology. By changing the arc-discharging times, the propagation loss can be adjusted and then the differentiation order is tuned. A nearly Gaussian-like optical pulse with 3 dB bandwidth of 8 nm is launched into the differentiator and a 0.65 order differentiation of the input pulse is achieved with a processing error of 2.55%. Project supported by the the National Natural Science Foundation of China (Nos. 61522509, 61377002, 61535012), the National High-Tech Research & Development Program of China (No. SS2015AA011002), and the Beijing Natural Science Foundation (No. 4152052). Ming Li was supported in part by the Thousand Young Talent Program.

Oscillation of solutions of some higher order linear differential equations

Directory of Open Access Journals (Sweden)

Hong-Yan Xu

2009-11-01

Full Text Available In this paper, we deal with the order of growth and the hyper order of solutions of higher order linear differential equations $$f^{(k}+B_{k-1}f^{(k-1}+\\cdots+B_1f'+B_0f=F$$ where $B_j(z (j=0,1,\\ldots,k-1$ and $F$ are entire functions or polynomials. Some results are obtained which improve and extend previous results given by Z.-X. Chen, J. Wang, T.-B. Cao and C.-H. Li.

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer from the step size reduction faced by standard explicit methods. The family is based on the standard second order orthogonal Runge-Kutta-Chebyshev (ROCK2) methods for deterministic problems. The convergence, meansquare, and asymptotic stability properties of the methods are analyzed. Numerical experiments, including applications to nonlinear SDEs and parabolic stochastic partial differential equations are presented and confirm the theoretical results. © 2013 Society for Industrial and Applied Mathematics.

Analytic formulae for the Hartree-Fock order parameter at arbitrary p/q filling factors for the 2DEG in a magnetic field

International Nuclear Information System (INIS)

Cabo Monte Oca, A. de.

1994-07-01

Analytic expressions for order parameters are given for the previously introduced general class of Hartree Fock states at arbitrary filling factors ν=p/q for odd q values. The order parameters are expressed as sums of magnetic translations eigenvalues over the filled single electron states. Simple summation formulae for the band spectra in terms of the same eigenvalues are also presented. The energy per particle at ν=1/3 is calculated for various states differing in the way of filling of the 1/3 of the orbitals. The calculated energies are not competing with the usual CDW results. However the high degree of electron overlapping allows for the next corrections to modify this situation. The discussion suggests these Hartree-Fock Slater determinants as interesting alternatives for the Tao-Thouless parent states which may correct their anomalous symmetry and correlation functions properties. (author). 28 refs

ACCURATE ESTIMATES OF CHARACTERISTIC EXPONENTS FOR SECOND ORDER DIFFERENTIAL EQUATION

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper, a second order linear differential equation is considered, and an accurate estimate method of characteristic exponent for it is presented. Finally, we give some examples to verify the feasibility of our result.

Existence and Uniqueness of Solutions for Coupled Systems of Higher-Order Nonlinear Fractional Differential Equations

Directory of Open Access Journals (Sweden)

Ahmad Bashir

2010-01-01

Full Text Available We study an initial value problem for a coupled Caputo type nonlinear fractional differential system of higher order. As a first problem, the nonhomogeneous terms in the coupled fractional differential system depend on the fractional derivatives of lower orders only. Then the nonhomogeneous terms in the fractional differential system are allowed to depend on the unknown functions together with the fractional derivative of lower orders. Our method of analysis is based on the reduction of the given system to an equivalent system of integral equations. Applying the nonlinear alternative of Leray-Schauder, we prove the existence of solutions of the fractional differential system. The uniqueness of solutions of the fractional differential system is established by using the Banach contraction principle. An illustrative example is also presented.

Linear dynamic analysis of arbitrary thin shells modal superposition by using finite element method

International Nuclear Information System (INIS)

Goncalves Filho, O.J.A.

1978-11-01

The linear dynamic behaviour of arbitrary thin shells by the Finite Element Method is studied. Plane triangular elements with eighteen degrees of freedom each are used. The general equations of movement are obtained from the Hamilton Principle and solved by the Modal Superposition Method. The presence of a viscous type damping can be considered by means of percentages of the critical damping. An automatic computer program was developed to provide the vibratory properties and the dynamic response to several types of deterministic loadings, including temperature effects. The program was written in FORTRAN IV for the Burroughs B-6700 computer. (author)

Approximating second-order vector differential operators on distorted meshes in two space dimensions

International Nuclear Information System (INIS)

Hermeline, F.

2008-01-01

A new finite volume method is presented for approximating second-order vector differential operators in two space dimensions. This method allows distorted triangle or quadrilateral meshes to be used without the numerical results being too much altered. The matrices that need to be inverted are symmetric positive definite therefore, the most powerful linear solvers can be applied. The method has been tested on a few second-order vector partial differential equations coming from elasticity and fluids mechanics areas. These numerical experiments show that it is second-order accurate and locking-free. (authors)

Integrability of systems of two second-order ordinary differential equations admitting four-dimensional Lie algebras.

Science.gov (United States)

Gainetdinova, A A; Gazizov, R K

2017-01-01

We suggest an algorithm for integrating systems of two second-order ordinary differential equations with four symmetries. In particular, if the admitted transformation group has two second-order differential invariants, the corresponding system can be integrated by quadratures using invariant representation and the operator of invariant differentiation. Otherwise, the systems reduce to partially uncoupled forms and can also be integrated by quadratures.

Deterministic factor analysis: methods of integro-differentiation of non-integral order

Directory of Open Access Journals (Sweden)

Valentina V. Tarasova

2016-12-01

Full Text Available Objective to summarize the methods of deterministic factor economic analysis namely the differential calculus and the integral method. nbsp Methods mathematical methods for integrodifferentiation of nonintegral order the theory of derivatives and integrals of fractional nonintegral order. Results the basic concepts are formulated and the new methods are developed that take into account the memory and nonlocality effects in the quantitative description of the influence of individual factors on the change in the effective economic indicator. Two methods are proposed for integrodifferentiation of nonintegral order for the deterministic factor analysis of economic processes with memory and nonlocality. It is shown that the method of integrodifferentiation of nonintegral order can give more accurate results compared with standard methods method of differentiation using the first order derivatives and the integral method using the integration of the first order for a wide class of functions describing effective economic indicators. Scientific novelty the new methods of deterministic factor analysis are proposed the method of differential calculus of nonintegral order and the integral method of nonintegral order. Practical significance the basic concepts and formulas of the article can be used in scientific and analytical activity for factor analysis of economic processes. The proposed method for integrodifferentiation of nonintegral order extends the capabilities of the determined factorial economic analysis. The new quantitative method of deterministic factor analysis may become the beginning of quantitative studies of economic agents behavior with memory hereditarity and spatial nonlocality. The proposed methods of deterministic factor analysis can be used in the study of economic processes which follow the exponential law in which the indicators endogenous variables are power functions of the factors exogenous variables including the processes

Automatic low-order aberration correction based on geometrical optics for slab lasers.

Science.gov (United States)

Yu, Xin; Dong, Lizhi; Lai, Boheng; Yang, Ping; Liu, Yong; Kong, Qingfeng; Yang, Kangjian; Tang, Guomao; Xu, Bing

2017-02-20

In this paper, we present a method based on geometry optics to simultaneously correct low-order aberrations and reshape the beams of slab lasers. A coaxial optical system with three lenses is adapted. The positions of the three lenses are directly calculated based on the beam parameters detected by wavefront sensors. The initial sizes of the input beams are 1.8  mm×11  mm, and peak-to-valley (PV) values of the wavefront range up to several tens of microns. After automatic correction, the dimensions may reach nearly 22  mm×22  mm as expected, and PV values of the wavefront are less than 2 μm. The effectiveness and precision of this method are verified with experiments.

Higher-order meshing of implicit geometries, Part II: Approximations on manifolds

Science.gov (United States)

Fries, T. P.; Schöllhammer, D.

2017-11-01

A new concept for the higher-order accurate approximation of partial differential equations on manifolds is proposed where a surface mesh composed by higher-order elements is automatically generated based on level-set data. Thereby, it enables a completely automatic workflow from the geometric description to the numerical analysis without any user-intervention. A master level-set function defines the shape of the manifold through its zero-isosurface which is then restricted to a finite domain by additional level-set functions. It is ensured that the surface elements are sufficiently continuous and shape regular which is achieved by manipulating the background mesh. The numerical results show that optimal convergence rates are obtained with a moderate increase in the condition number compared to handcrafted surface meshes.

Data-driven uncertainty quantification using the arbitrary polynomial chaos expansion

International Nuclear Information System (INIS)

Oladyshkin, S.; Nowak, W.

2012-01-01

We discuss the arbitrary polynomial chaos (aPC), which has been subject of research in a few recent theoretical papers. Like all polynomial chaos expansion techniques, aPC approximates the dependence of simulation model output on model parameters by expansion in an orthogonal polynomial basis. The aPC generalizes chaos expansion techniques towards arbitrary distributions with arbitrary probability measures, which can be either discrete, continuous, or discretized continuous and can be specified either analytically (as probability density/cumulative distribution functions), numerically as histogram or as raw data sets. We show that the aPC at finite expansion order only demands the existence of a finite number of moments and does not require the complete knowledge or even existence of a probability density function. This avoids the necessity to assign parametric probability distributions that are not sufficiently supported by limited available data. Alternatively, it allows modellers to choose freely of technical constraints the shapes of their statistical assumptions. Our key idea is to align the complexity level and order of analysis with the reliability and detail level of statistical information on the input parameters. We provide conditions for existence and clarify the relation of the aPC to statistical moments of model parameters. We test the performance of the aPC with diverse statistical distributions and with raw data. In these exemplary test cases, we illustrate the convergence with increasing expansion order and, for the first time, with increasing reliability level of statistical input information. Our results indicate that the aPC shows an exponential convergence rate and converges faster than classical polynomial chaos expansion techniques.

High-order noise filtering in nontrivial quantum logic gates

CSIR Research Space (South Africa)

Green, T

2012-07-01

Full Text Available composed of arbitrary control sequences. We present a general method to calculate the ensemble-averaged entanglement fidelity to arbitrary order in terms of noise filter functions, and provide explicit expressions to fourth order in the noise strength...

Functional methods for arbitrary densities in curved spacetime

International Nuclear Information System (INIS)

Basler, M.

1993-01-01

This paper gives an introduction to the technique of functional differentiation and integration in curved spacetime, applied to examples from quantum field theory. Special attention is drawn on the choice of functional integral measure. Referring to a suggestion by Toms, fields are choosen as arbitrary scalar, spinorial or vectorial densities. The technique developed by Toms for a pure quadratic Lagrangian are extended to the calculation of the generating functional with external sources. Included are two examples of interacting theories, a self-interacting scalar field and a Yang-Mills theory. For these theories the complete set of Feynman graphs depending on the weight of variables is derived. (orig.)

High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new methods of weak order two, in particular, semi-implicit integrators well suited for stiff (meansquare stable) stochastic problems, and implicit integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of our methodology.

ON THE INSTABILITY OF SOLUTIONS TO A NONLINEAR VECTOR DIFFERENTIAL EQUATION OF FOURTH ORDER

Institute of Scientific and Technical Information of China (English)

无

2011-01-01

This paper presents a new result related to the instability of the zero solution to a nonlinear vector differential equation of fourth order.Our result includes and improves an instability result in the previous literature,which is related to the instability of the zero solution to a nonlinear scalar differential equation of fourth order.

Lie group classification of first-order delay ordinary differential equations

Science.gov (United States)

Dorodnitsyn, Vladimir A.; Kozlov, Roman; Meleshko, Sergey V.; Winternitz, Pavel

2018-05-01

A group classification of first-order delay ordinary differential equations (DODEs) accompanied by an equation for the delay parameter (delay relation) is presented. A subset of such systems (delay ordinary differential systems or DODSs), which consists of linear DODEs and solution-independent delay relations, have infinite-dimensional symmetry algebras—as do nonlinear ones that are linearizable by an invertible transformation of variables. Genuinely nonlinear DODSs have symmetry algebras of dimension n, . It is shown how exact analytical solutions of invariant DODSs can be obtained using symmetry reduction.

N-th order impulsive integro-differential equations in Banach spaces

Directory of Open Access Journals (Sweden)

Manfeng Hu

2004-03-01

Full Text Available We investigate the maximal and minimal solutions of initial value problem for N-th order nonlinear impulsive integro-differential equation in Banach space by establishing a comparison result and using the upper and lower solutions methods.

On nonnegative solutions of second order linear functional differential equations

Czech Academy of Sciences Publication Activity Database

Lomtatidze, Alexander; Vodstrčil, Petr

2004-01-01

Roč. 32, č. 1 (2004), s. 59-88 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z1019905 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics

Construction and accuracy of partial differential equation approximations to the chemical master equation.

Science.gov (United States)

Grima, Ramon

2011-11-01

The mesoscopic description of chemical kinetics, the chemical master equation, can be exactly solved in only a few simple cases. The analytical intractability stems from the discrete character of the equation, and hence considerable effort has been invested in the development of Fokker-Planck equations, second-order partial differential equation approximations to the master equation. We here consider two different types of higher-order partial differential approximations, one derived from the system-size expansion and the other from the Kramers-Moyal expansion, and derive the accuracy of their predictions for chemical reactive networks composed of arbitrary numbers of unimolecular and bimolecular reactions. In particular, we show that the partial differential equation approximation of order Q from the Kramers-Moyal expansion leads to estimates of the mean number of molecules accurate to order Ω(-(2Q-3)/2), of the variance of the fluctuations in the number of molecules accurate to order Ω(-(2Q-5)/2), and of skewness accurate to order Ω(-(Q-2)). We also show that for large Q, the accuracy in the estimates can be matched only by a partial differential equation approximation from the system-size expansion of approximate order 2Q. Hence, we conclude that partial differential approximations based on the Kramers-Moyal expansion generally lead to considerably more accurate estimates in the mean, variance, and skewness than approximations of the same order derived from the system-size expansion.

Oscillation of certain higher-order neutral partial functional differential equations.

Science.gov (United States)

Li, Wei Nian; Sheng, Weihong

2016-01-01

In this paper, we study the oscillation of certain higher-order neutral partial functional differential equations with the Robin boundary conditions. Some oscillation criteria are established. Two examples are given to illustrate the main results in the end of this paper.

Asymptotic integration of a linear fourth order differential equation of Poincaré type

Directory of Open Access Journals (Sweden)

Anibal Coronel

2015-11-01

Full Text Available This article deals with the asymptotic behavior of nonoscillatory solutions of fourth order linear differential equation where the coefficients are perturbations of constants. We define a change of variable and deduce that the new variable satisfies a third order nonlinear differential equation. We assume three hypotheses. The first hypothesis is related to the constant coefficients and set up that the characteristic polynomial associated with the fourth order linear equation has simple and real roots. The other two hypotheses are related to the behavior of the perturbation functions and establish asymptotic integral smallness conditions of the perturbations. Under these general hypotheses, we obtain four main results. The first two results are related to the application of a fixed point argument to prove that the nonlinear third order equation has a unique solution. The next result concerns with the asymptotic behavior of the solutions of the nonlinear third order equation. The fourth main theorem is introduced to establish the existence of a fundamental system of solutions and to precise the formulas for the asymptotic behavior of the linear fourth order differential equation. In addition, we present an example to show that the results introduced in this paper can be applied in situations where the assumptions of some classical theorems are not satisfied.

Unveiling Reality of the Mind: Cultural Arbitrary of Consumerism

Science.gov (United States)

Choi, Su-Jin

2012-01-01

This paper discusses the cultural arbitrary of consumerism by focusing on a personal realm. That is, I discuss what consumerism appeals to and how it flourishes in relation to our minds. I argue that we need to unveil reality of the mind, be aware of ourselves in relation to the perpetuation of consumerism, in order to critically intervene in the…

Differential operators in a Clifford analysis associated to differential equations with anti-monogenic right-hand sides

International Nuclear Information System (INIS)

Nguyen Thanh Van

2006-12-01

This paper deals with the initial value problem of the type φw / φt = L (t, x, w, φw / φx i ) (1) w(0, x) = φ(x) (2) where t is the time, L is a linear first order operator in a Clifford Analysis and φ is a generalized monogenic function. We give sufficient conditions on the coefficients of operator L under which L is associated to differential equations with anti-monogenic right-hand sides. For such operator L the initial problem (1),(2) is solvable for an arbitrary generalized monogenic initial function φ and the solution is also generalized monogenic for each t. (author)

Antiperiodic Boundary Value Problems for Second-Order Impulsive Ordinary Differential Equations

Directory of Open Access Journals (Sweden)

2009-02-01

Full Text Available We consider a second-order ordinary differential equation with antiperiodic boundary conditions and impulses. By using Schaefer's fixed-point theorem, some existence results are obtained.

Saturation behavior: a general relationship described by a simple second-order differential equation.

Science.gov (United States)

Kepner, Gordon R

2010-04-13

The numerous natural phenomena that exhibit saturation behavior, e.g., ligand binding and enzyme kinetics, have been approached, to date, via empirical and particular analyses. This paper presents a mechanism-free, and assumption-free, second-order differential equation, designed only to describe a typical relationship between the variables governing these phenomena. It develops a mathematical model for this relation, based solely on the analysis of the typical experimental data plot and its saturation characteristics. Its utility complements the traditional empirical approaches. For the general saturation curve, described in terms of its independent (x) and dependent (y) variables, a second-order differential equation is obtained that applies to any saturation phenomena. It shows that the driving factor for the basic saturation behavior is the probability of the interactive site being free, which is described quantitatively. Solving the equation relates the variables in terms of the two empirical constants common to all these phenomena, the initial slope of the data plot and the limiting value at saturation. A first-order differential equation for the slope emerged that led to the concept of the effective binding rate at the active site and its dependence on the calculable probability the interactive site is free. These results are illustrated using specific cases, including ligand binding and enzyme kinetics. This leads to a revised understanding of how to interpret the empirical constants, in terms of the variables pertinent to the phenomenon under study. The second-order differential equation revealed the basic underlying relations that describe these saturation phenomena, and the basic mathematical properties of the standard experimental data plot. It was shown how to integrate this differential equation, and define the common basic properties of these phenomena. The results regarding the importance of the slope and the new perspectives on the empirical

Arbitrary Inequality in Reputation Systems

Science.gov (United States)

Frey, Vincenz; van de Rijt, Arnout

2016-12-01

Trust is an essential condition for exchange. Large societies must substitute the trust traditionally provided through kinship and sanctions in small groups to make exchange possible. The rise of internet-supported reputation systems has been celebrated for providing trust at a global scale, enabling the massive volumes of transactions between distant strangers that are characteristic of modern human societies. Here we problematize an overlooked side-effect of reputation systems: Equally trustworthy individuals may realize highly unequal exchange volumes. We report the results of a laboratory experiment that shows emergent differentiation between ex ante equivalent individuals when information on performance in past exchanges is shared. This arbitrary inequality results from cumulative advantage in the reputation-building process: Random initial distinctions grow as parties of good repute are chosen over those lacking a reputation. We conjecture that reputation systems produce artificial concentration in a wide range of markets and leave superior but untried exchange alternatives unexploited.

Periodic solutions of certain third order nonlinear differential systems with delay

International Nuclear Information System (INIS)

Tejumola, H.O.; Afuwape, A.U.

1990-12-01

This paper investigates the existence of 2π-periodic solutions of systems of third-order nonlinear differential equations, with delay, under varied assumptions. The results obtained extend earlier works of Tejumola and generalize to third order systems those of Conti, Iannacci and Nkashama as well as DePascale and Iannacci and Iannacci and Nkashama. 16 refs

Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces

Directory of Open Access Journals (Sweden)

Yongjin Li

2013-08-01

Full Text Available We prove the Hyers-Ulam stability of linear second-order differential equations in complex Banach spaces. That is, if y is an approximate solution of the differential equation $y''+ alpha y'(t +eta y = 0$ or $y''+ alpha y'(t +eta y = f(t$, then there exists an exact solution of the differential equation near to y.

POSITIVE SOLUTIONS TO SEMI-LINEAR SECOND-ORDER ORDINARY DIFFERENTIAL EQUATIONS IN BANACH SPACE

Institute of Scientific and Technical Information of China (English)

无

2008-01-01

In this paper,we study the existence of positive periodic solution to some second- order semi-linear differential equation in Banach space.By the fixed point index theory, we prove that the semi-linear differential equation has two positive periodic solutions.

Global asymptotical ω-periodicity of a fractional-order non-autonomous neural networks.

Science.gov (United States)

Chen, Boshan; Chen, Jiejie

2015-08-01

We study the global asymptotic ω-periodicity for a fractional-order non-autonomous neural networks. Firstly, based on the Caputo fractional-order derivative it is shown that ω-periodic or autonomous fractional-order neural networks cannot generate exactly ω-periodic signals. Next, by using the contraction mapping principle we discuss the existence and uniqueness of S-asymptotically ω-periodic solution for a class of fractional-order non-autonomous neural networks. Then by using a fractional-order differential and integral inequality technique, we study global Mittag-Leffler stability and global asymptotical periodicity of the fractional-order non-autonomous neural networks, which shows that all paths of the networks, starting from arbitrary points and responding to persistent, nonconstant ω-periodic external inputs, asymptotically converge to the same nonconstant ω-periodic function that may be not a solution. Copyright © 2015 Elsevier Ltd. All rights reserved.

Spectroscopy of the Schwarzschild black hole at arbitrary frequencies.

Science.gov (United States)

Casals, Marc; Ottewill, Adrian

2012-09-14

Linear field perturbations of a black hole are described by the Green function of the wave equation that they obey. After Fourier decomposing the Green function, its two natural contributions are given by poles (quasinormal modes) and a largely unexplored branch cut in the complex frequency plane. We present new analytic methods for calculating the branch cut on a Schwarzschild black hole for arbitrary values of the frequency. The branch cut yields a power-law tail decay for late times in the response of a black hole to an initial perturbation. We determine explicitly the first three orders in the power-law and show that the branch cut also yields a new logarithmic behavior T(-2ℓ-5)lnT for late times. Before the tail sets in, the quasinormal modes dominate the black hole response. For electromagnetic perturbations, the quasinormal mode frequencies approach the branch cut at large overtone index n. We determine these frequencies up to n(-5/2) and, formally, to arbitrary order. Highly damped quasinormal modes are of particular interest in that they have been linked to quantum properties of black holes.

Einstein-Weyl spaces and third-order differential equations

Science.gov (United States)

Tod, K. P.

2000-08-01

The three-dimensional null-surface formalism of Tanimoto [M. Tanimoto, "On the null surface formalism," Report No. gr-qc/9703003 (1997)] and Forni et al. [Forni et al., "Null surfaces formation in 3D," J. Math Phys. (submitted)] are extended to describe Einstein-Weyl spaces, following Cartan [E. Cartan, "Les espaces généralisées et l'integration de certaines classes d'equations différentielles," C. R. Acad. Sci. 206, 1425-1429 (1938); "La geometria de las ecuaciones diferenciales de tercer order," Rev. Mat. Hispano-Am. 4, 1-31 (1941)]. In the resulting formalism, Einstein-Weyl spaces are obtained from a particular class of third-order differential equations. Some examples of the construction which include some new Einstein-Weyl spaces are given.

Free massless fermionic fields of arbitrary spin in d-dimensional anti-de Sitter space

Energy Technology Data Exchange (ETDEWEB)

Vasiliev, M A

1988-04-25

Free massless fermionic fields of arbitrary spins, corresponding to fully symmetric tensor-spinor irreducible representations of the flat little group SO(d-2), are described in d-dimensional anti-de Sitter space in terms of differential forms. Appropriate linearized higher-spin curvature 2-forms are found. Explicitly gauge invariant higher-spin actions are constructed in terms of these linearized curvatures.

Eight equation model for arbitrary shaped pipe conveying fluid

International Nuclear Information System (INIS)

Gale, J.; Tiselj, I.

2006-01-01

Linear eight-equation system for two-way coupling of single-phase fluid transient and arbitrary shaped one-dimensional pipeline movement is described and discussed. The governing phenomenon described with this system is also known as Fluid-Structure Interaction. Standard Skalak's four-equation model for axial coupling was improved with additional four Timoshenko's beam equations for description of flexural displacements and rotations. In addition to the conventional eight-equation system that enables coupling of straight sections, the applied mathematical model was improved for description of the arbitrary shaped pipeline located in two-dimensional plane. The applied model was solved with second-order accurate numerical method that is based on Godounov's characteristic upwind schemes. The model was successfully used for simulation of the rod impact induced transient and conventional instantaneous valve closure induced transient in the tank-pipe-valve system. (author)

Applying fractional order PID to design TCSC-based damping controller in coordination with automatic generation control of interconnected multi-source power system

Directory of Open Access Journals (Sweden)

Javad Morsali

2017-02-01

Full Text Available In this paper, fractional order proportional-integral-differential (FOPID controller is employed in the design of thyristor controlled series capacitor (TCSC-based damping controller in coordination with the secondary integral controller as automatic generation control (AGC loop. In doing so, the contribution of the TCSC in tie-line power exchange is extracted mathematically for small load disturbance. Adjustable parameters of the proposed FOPID-based TCSC damping controller and the AGC loop are optimized concurrently via an improved particle swarm optimization (IPSO algorithm which is reinforced by chaotic parameter and crossover operator to obtain a globally optimal solution. The powerful FOMCON toolbox is used along with MATLAB for handling fractional order modeling and control. An interconnected multi-source power system is simulated regarding the physical constraints of generation rate constraint (GRC nonlinearity and governor dead band (GDB effect. Simulation results using FOMCON toolbox demonstrate that the proposed FOPID-based TCSC damping controller achieves the greatest dynamic performance under different load perturbation patterns in comparison with phase lead-lag and classical PID-based TCSC damping controllers, all in coordination with the integral AGC. Moreover, sensitivity analyses are performed to show the robustness of the proposed controller under various uncertainty scenarios.

A Solution to the Fundamental Linear Fractional Order Differential Equation

Science.gov (United States)

Hartley, Tom T.; Lorenzo, Carl F.

1998-01-01

This paper provides a solution to the fundamental linear fractional order differential equation, namely, (sub c)d(sup q, sub t) + ax(t) = bu(t). The impulse response solution is shown to be a series, named the F-function, which generalizes the normal exponential function. The F-function provides the basis for a qth order "fractional pole". Complex plane behavior is elucidated and a simple example, the inductor terminated semi- infinite lossy line, is used to demonstrate the theory.

The Complementary Hankel Type Transformations Of Arbitrary Order

Directory of Open Access Journals (Sweden)

B.B. Waphare

2013-09-01

Full Text Available In this paper four self-reciprocal integral transformations of Hankel type are defined. The simultaneous use of trans-formations H1,α,β and H2,α,β (which are denoted by Hα,β allows us to solve many problems of Mathematical Physics involving the differential operator ∆α,β= D2+4αx−1D, whereas the pair of transformations H3,α,β and H4,α,β (which we express by Hα,β permits us to tackle those problems containing its adjoint operator, no matter what the real value of α − β be. These transformations are also investigated in a space of generalized functions according to the mixed Parseval equation.

Automatic Censoring CFAR Detector Based on Ordered Data Difference for Low-Flying Helicopter Safety

Directory of Open Access Journals (Sweden)

Wen Jiang

2016-07-01

Full Text Available Being equipped with a millimeter-wave radar allows a low-flying helicopter to sense the surroundings in real time, which significantly increases its safety. However, nonhomogeneous clutter environments, such as a multiple target situation and a clutter edge environment, can dramatically affect the radar signal detection performance. In order to improve the radar signal detection performance in nonhomogeneous clutter environments, this paper proposes a new automatic censored cell averaging CFAR detector. The proposed CFAR detector does not require any prior information about the background environment and uses the hypothesis test of the first-order difference (FOD result of ordered data to reject the unwanted samples in the reference window. After censoring the unwanted ranked cells, the remaining samples are combined to form an estimate of the background power level, thus getting better radar signal detection performance. The simulation results show that the FOD-CFAR detector provides low loss CFAR performance in a homogeneous environment and also performs robustly in nonhomogeneous environments. Furthermore, the measured results of a low-flying helicopter validate the basic performance of the proposed method.

Automatic Censoring CFAR Detector Based on Ordered Data Difference for Low-Flying Helicopter Safety

Science.gov (United States)

Jiang, Wen; Huang, Yulin; Yang, Jianyu

2016-01-01

Being equipped with a millimeter-wave radar allows a low-flying helicopter to sense the surroundings in real time, which significantly increases its safety. However, nonhomogeneous clutter environments, such as a multiple target situation and a clutter edge environment, can dramatically affect the radar signal detection performance. In order to improve the radar signal detection performance in nonhomogeneous clutter environments, this paper proposes a new automatic censored cell averaging CFAR detector. The proposed CFAR detector does not require any prior information about the background environment and uses the hypothesis test of the first-order difference (FOD) result of ordered data to reject the unwanted samples in the reference window. After censoring the unwanted ranked cells, the remaining samples are combined to form an estimate of the background power level, thus getting better radar signal detection performance. The simulation results show that the FOD-CFAR detector provides low loss CFAR performance in a homogeneous environment and also performs robustly in nonhomogeneous environments. Furthermore, the measured results of a low-flying helicopter validate the basic performance of the proposed method. PMID:27399714

Application of the Lie Symmetry Analysis for second-order fractional differential equations

Directory of Open Access Journals (Sweden)

Mousa Ilie

2017-12-01

Full Text Available Obtaining analytical or numerical solution of fractional differential equations is one of the troublesome and challenging issue among mathematicians and engineers, specifically in recent years. The purpose of this paper Lie Symmetry method is developed to solve second-order fractional differential equations, based on conformable fractional derivative. Some numerical examples are presented to illustrate the proposed approach.

Variations in the Solution of Linear First-Order Differential Equations. Classroom Notes

Science.gov (United States)

Seaman, Brian; Osler, Thomas J.

2004-01-01

A special project which can be given to students of ordinary differential equations is described in detail. Students create new differential equations by changing the dependent variable in the familiar linear first-order equation (dv/dx)+p(x)v=q(x) by means of a substitution v=f(y). The student then creates a table of the new equations and…

On Some Pursuit and Evasion Differential Game Problems for an Infinite Number of First-Order Differential Equations

Directory of Open Access Journals (Sweden)

Abbas Badakaya Ja'afaru

2012-01-01

Full Text Available We study pursuit and evasion differential game problems described by infinite number of first-order differential equations with function coefficients in Hilbert space l2. Problems involving integral, geometric, and mix constraints to the control functions of the players are considered. In each case, we give sufficient conditions for completion of pursuit and for which evasion is possible. Consequently, strategy of the pursuer and control function of the evader are constructed in an explicit form for every problem considered.

Symmetries of th-Order Approximate Stochastic Ordinary Differential Equations

OpenAIRE

Fredericks, E.; Mahomed, F. M.

2012-01-01

Symmetries of $n$ th-order approximate stochastic ordinary differential equations (SODEs) are studied. The determining equations of these SODEs are derived in an Itô calculus context. These determining equations are not stochastic in nature. SODEs are normally used to model nature (e.g., earthquakes) or for testing the safety and reliability of models in construction engineering when looking at the impact of random perturbations.

A characteristic based multiple balance approach for SN on arbitrary polygonal meshes

International Nuclear Information System (INIS)

Grove, R.E.; Pevey, R.E.

1995-01-01

The authors introduce a new approach for characteristic based S n transport on arbitrary polygonal meshes in XY geometry. They approximate a general surface as an arbitrary polygon and rotate to a coordinate system aligned with the direction of particle travel. They use exact moment balance equations on whole cells and subregions called slices and close the system by analytically solving the characteristic equation. The authors assume spatial functions for boundary conditions and cell sources and formulate analogous functions for outgoing edge and cell angular fluxes which exactly preserve spatial moments of the analytic solution. In principle, their approach provides the framework to extend characteristic methods formulated on rectangular grids to arbitrary polygonal meshes. The authors derive schemes based on step and linear spatial approximations. Their step characteristic scheme is mathematically equivalent to the Extended Step Characteristic (ESC) method but their approach and scheme differ in the geometry rotation and in the solution form. Their solutions are simple and permit edge-based transport sweep ordering

Plane waves and spherical means applied to partial differential equations

CERN Document Server

John, Fritz

2004-01-01

Elementary and self-contained, this heterogeneous collection of results on partial differential equations employs certain elementary identities for plane and spherical integrals of an arbitrary function, showing how a variety of results on fairly general differential equations follow from those identities. The first chapter deals with the decomposition of arbitrary functions into functions of the type of plane waves. Succeeding chapters introduce the first application of the Radon transformation and examine the solution of the initial value problem for homogeneous hyperbolic equations with con

Archimedeanization of ordered vector spaces

OpenAIRE

Emelyanov, Eduard Yu.

2014-01-01

In the case of an ordered vector space with an order unit, the Archimedeanization method has been developed recently by V.I Paulsen and M. Tomforde. We present a general version of the Archimedeanization which covers arbitrary ordered vector spaces.

Coherence measures in automatic time-migration velocity analysis

International Nuclear Information System (INIS)

Maciel, Jonathas S; Costa, Jessé C; Schleicher, Jörg

2012-01-01

Time-migration velocity analysis can be carried out automatically by evaluating the coherence of migrated seismic events in common-image gathers (CIGs). The performance of gradient methods for automatic time-migration velocity analysis depends on the coherence measures used as the objective function. We compare the results of four different coherence measures, being conventional semblance, differential semblance, an extended differential semblance using differences of more distant image traces and the product of the latter with conventional semblance. In our numerical experiments, the objective functions based on conventional semblance and on the product of conventional semblance with extended differential semblance provided the best velocity models, as evaluated by the flatness of the resulting CIGs. The method can be easily extended to anisotropic media. (paper)

Acoustic streaming: an arbitrary Lagrangian-Eulerian perspective.

Science.gov (United States)

Nama, Nitesh; Huang, Tony Jun; Costanzo, Francesco

2017-08-25

We analyse acoustic streaming flows using an arbitrary Lagrangian Eulerian (ALE) perspective. The formulation stems from an explicit separation of time scales resulting in two subproblems: a first-order problem, formulated in terms of the fluid displacement at the fast scale, and a second-order problem, formulated in terms of the Lagrangian flow velocity at the slow time scale. Following a rigorous time-averaging procedure, the second-order problem is shown to be intrinsically steady, and with exact boundary conditions at the oscillating walls. Also, as the second-order problem is solved directly for the Lagrangian velocity, the formulation does not need to employ the notion of Stokes drift, or any associated post-processing, thus facilitating a direct comparison with experiments. Because the first-order problem is formulated in terms of the displacement field, our formulation is directly applicable to more complex fluid-structure interaction problems in microacoustofluidic devices. After the formulation's exposition, we present numerical results that illustrate the advantages of the formulation with respect to current approaches.

Second-order sliding mode controller with model reference adaptation for automatic train operation

Science.gov (United States)

Ganesan, M.; Ezhilarasi, D.; Benni, Jijo

2017-11-01

In this paper, a new approach to model reference based adaptive second-order sliding mode control together with adaptive state feedback is presented to control the longitudinal dynamic motion of a high speed train for automatic train operation with the objective of minimal jerk travel by the passengers. The nonlinear dynamic model for the longitudinal motion of the train comprises of a locomotive and coach subsystems is constructed using multiple point-mass model by considering the forces acting on the vehicle. An adaptation scheme using Lyapunov criterion is derived to tune the controller gains by considering a linear, stable reference model that ensures the stability of the system in closed loop. The effectiveness of the controller tracking performance is tested under uncertain passenger load, coupler-draft gear parameters, propulsion resistance coefficients variations and environmental disturbances due to side wind and wet rail conditions. The results demonstrate improved tracking performance of the proposed control scheme with a least jerk under maximum parameter uncertainties when compared to constant gain second-order sliding mode control.

Joint estimation of the fractional differentiation orders and the unknown input for linear fractional non-commensurate system

KAUST Repository

Belkhatir, Zehor

2015-11-05

This paper deals with the joint estimation of the unknown input and the fractional differentiation orders of a linear fractional order system. A two-stage algorithm combining the modulating functions with a first-order Newton method is applied to solve this estimation problem. First, the modulating functions approach is used to estimate the unknown input for a given fractional differentiation orders. Then, the method is combined with a first-order Newton technique to identify the fractional orders jointly with the input. To show the efficiency of the proposed method, numerical examples illustrating the estimation of the neural activity, considered as input of a fractional model of the neurovascular coupling, along with the fractional differentiation orders are presented in both noise-free and noisy cases.

Arbitrary waveform generator to improve laser diode driver performance

Science.gov (United States)

Fulkerson, Jr, Edward Steven

2015-11-03

An arbitrary waveform generator modifies the input signal to a laser diode driver circuit in order to reduce the overshoot/undershoot and provide a "flat-top" signal to the laser diode driver circuit. The input signal is modified based on the original received signal and the feedback from the laser diode by measuring the actual current flowing in the laser diode after the original signal is applied to the laser diode.

On the maximal cut of Feynman integrals and the solution of their differential equations

Directory of Open Access Journals (Sweden)

Amedeo Primo

2017-03-01

Full Text Available The standard procedure for computing scalar multi-loop Feynman integrals consists in reducing them to a basis of so-called master integrals, derive differential equations in the external invariants satisfied by the latter and, finally, try to solve them as a Laurent series in ϵ=(4−d/2, where d are the space–time dimensions. The differential equations are, in general, coupled and can be solved using Euler's variation of constants, provided that a set of homogeneous solutions is known. Given an arbitrary differential equation of order higher than one, there exists no general method for finding its homogeneous solutions. In this paper we show that the maximal cut of the integrals under consideration provides one set of homogeneous solutions, simplifying substantially the solution of the differential equations.

Top-quark decay at next-to-next-to-leading order in QCD.

Science.gov (United States)

Gao, Jun; Li, Chong Sheng; Zhu, Hua Xing

2013-01-25

We present the complete calculation of the top-quark decay width at next-to-next-to-leading order in QCD, including next-to-leading electroweak corrections as well as finite bottom quark mass and W boson width effects. In particular, we also show the first results of the fully differential decay rates for the top-quark semileptonic decay t → W(+)(l(+)ν)b at next-to-next-to-leading order in QCD. Our method is based on the understanding of the invariant mass distribution of the final-state jet in the singular limit from effective field theory. Our result can be used to study arbitrary infrared-safe observables of top-quark decay with the highest perturbative accuracy.

Data-driven automatic parking constrained control for four-wheeled mobile vehicles

Directory of Open Access Journals (Sweden)

Wenxu Yan

2016-11-01

Full Text Available In this article, a novel data-driven constrained control scheme is proposed for automatic parking systems. The design of the proposed scheme only depends on the steering angle and the orientation angle of the car, and it does not involve any model information of the car. Therefore, the proposed scheme-based automatic parking system is applicable to different kinds of cars. In order to further reduce the desired trajectory coordinate tracking errors, a coordinates compensation algorithm is also proposed. In the design procedure of the controller, a novel dynamic anti-windup compensator is used to deal with the change magnitude and rate saturations of automatic parking control input. It is theoretically proven that all the signals in the closed-loop system are uniformly ultimately bounded based on Lyapunov stability analysis method. Finally, a simulation comparison among the proposed scheme with coordinates compensation and Proportion Integration Differentiation (PID control algorithm is given. It is shown that the proposed scheme with coordinates compensation has smaller tracking errors and more rapid responses than PID scheme.

Weak Second Order Explicit Stabilized Methods for Stiff Stochastic Differential Equations

KAUST Repository

Abdulle, Assyr; Vilmart, Gilles; Zygalakis, Konstantinos C.

2013-01-01

We introduce a new family of explicit integrators for stiff Itô stochastic differential equations (SDEs) of weak order two. These numerical methods belong to the class of one-step stabilized methods with extended stability domains and do not suffer

DEVELOPMENT OF AUTOMATIC EXTRACTION METHOD FOR ROAD UPDATE INFORMATION BASED ON PUBLIC WORK ORDER OUTLOOK

Science.gov (United States)

Sekimoto, Yoshihide; Nakajo, Satoru; Minami, Yoshitaka; Yamaguchi, Syohei; Yamada, Harutoshi; Fuse, Takashi

Recently, disclosure of statistic data, representing financial effects or burden for public work, through each web site of national or local government, enables us to discuss macroscopic financial trends. However, it is still difficult to grasp a basic property nationwide how each spot was changed by public work. In this research, our research purpose is to collect road update information reasonably which various road managers provide, in order to realize efficient updating of various maps such as car navigation maps. In particular, we develop the system extracting public work concerned and registering summary including position information to database automatically from public work order outlook, released by each local government, combinating some web mining technologies. Finally, we collect and register several tens of thousands from web site all over Japan, and confirm the feasibility of our method.

A New Factorisation of a General Second Order Differential Equation

Science.gov (United States)

Clegg, Janet

2006-01-01

A factorisation of a general second order ordinary differential equation is introduced from which the full solution to the equation can be obtained by performing two integrations. The method is compared with traditional methods for solving these type of equations. It is shown how the Green's function can be derived directly from the factorisation…

Relative boundedness and compactness theory for second-order differential operators

Directory of Open Access Journals (Sweden)

Don B. Hinton

1997-01-01

Full Text Available The problem considered is to give necessary and sufficient conditions for perturbations of a second-order ordinary differential operator to be either relatively bounded or relatively compact. Such conditions are found for three classes of operators. The conditions are expressed in terms of integral averages of the coefficients of the perturbing operator.

Multivariate Padé Approximation for Solving Nonlinear Partial Differential Equations of Fractional Order

Directory of Open Access Journals (Sweden)

Veyis Turut

2013-01-01

Full Text Available Two tecHniques were implemented, the Adomian decomposition method (ADM and multivariate Padé approximation (MPA, for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in Caputo sense. First, the fractional differential equation has been solved and converted to power series by Adomian decomposition method (ADM, then power series solution of fractional differential equation was put into multivariate Padé series. Finally, numerical results were compared and presented in tables and figures.

Travelling Wave Solutions of Coupled Burger’s Equations of Time-Space Fractional Order by Novel (Gʹ/G-Expansion Method

Directory of Open Access Journals (Sweden)

Rashida Hussain

2017-04-01

Full Text Available In this paper, Novel (Gʹ/G-expansion method is used to find new generalized exact travelling wave solutions of fractional order coupled Burger’s equations in terms of trigonometric functions, rational functions and hyperbolic functions with arbitrary parameters. For the conversion of the partial differential equation to the ordinary differential equation, complex transformation method is used. Novel (Gʹ/G-expansion method is very effective and provides a powerful mathematical tool to solve nonlinear equations. Moreover, for the representation of these exact solutions we have plotted graphs for different values of parameters which were in travelling waveform.

Stability analysis for neutral stochastic differential equation of second order driven by Poisson jumps

Science.gov (United States)

Chadha, Alka; Bora, Swaroop Nandan

2017-11-01

This paper studies the existence, uniqueness, and exponential stability in mean square for the mild solution of neutral second order stochastic partial differential equations with infinite delay and Poisson jumps. By utilizing the Banach fixed point theorem, first the existence and uniqueness of the mild solution of neutral second order stochastic differential equations is established. Then, the mean square exponential stability for the mild solution of the stochastic system with Poisson jumps is obtained with the help of an established integral inequality.

Boundary-value problems for first and second order functional differential inclusions

Directory of Open Access Journals (Sweden)

Shihuang Hong

2003-03-01

Full Text Available This paper presents sufficient conditions for the existence of solutions to boundary-value problems of first and second order multi-valued differential equations in Banach spaces. Our results obtained using fixed point theorems, and lead to new existence principles.

Solving Second-Order Ordinary Differential Equations without Using Complex Numbers

Science.gov (United States)

Kougias, Ioannis E.

2009-01-01

Ordinary differential equations (ODEs) is a subject with a wide range of applications and the need of introducing it to students often arises in the last year of high school, as well as in the early stages of tertiary education. The usual methods of solving second-order ODEs with constant coefficients, among others, rely upon the use of complex…

Cheap arbitrary high order methods for single integrand SDEs

DEFF Research Database (Denmark)

Debrabant, Kristian; Kværnø, Anne

2017-01-01

For a particular class of Stratonovich SDE problems, here denoted as single integrand SDEs, we prove that by applying a deterministic Runge-Kutta method of order $p_d$ we obtain methods converging in the mean-square and weak sense with order $\\lfloor p_d/2\\rfloor$. The reason is that the B-series...

Entanglement entropy and differential entropy for massive flavors

International Nuclear Information System (INIS)

Jones, Peter A.R.; Taylor, Marika

2015-01-01

In this paper we compute the holographic entanglement entropy for massive flavors in the D3-D7 system, for arbitrary mass and various entangling region geometries. We show that the universal terms in the entanglement entropy exactly match those computed in the dual theory using conformal perturbation theory. We derive holographically the universal terms in the entanglement entropy for a CFT perturbed by a relevant operator, up to second order in the coupling; our results are valid for any entangling region geometry. We present a new method for computing the entanglement entropy of any top-down brane probe system using Kaluza-Klein holography and illustrate our results with massive flavors at finite density. Finally we discuss the differential entropy for brane probe systems, emphasising that the differential entropy captures only the effective lower-dimensional Einstein metric rather than the ten-dimensional geometry.

Bessel-like beams modulated by arbitrary radial functions

Science.gov (United States)

Herman; Wiggins

2000-06-01

An approximate method for determining the radial and axial intensity of a Bessel-like beam is presented for the general case in which a radial Bessel distribution of any order is modulated by an arbitrary function. For Bessel-Gauss, generalized Bessel-Gauss, and Bessel-super-Gauss beams, this simple approximation yields results that are very close to the exact values, while they are exact for Bessel beams. A practical beam that can be generated with a combination of simple lenses is also analyzed and illustrated.

Vertical and lateral forces when a permanent magnet above a superconductor traverses in arbitrary directions

Science.gov (United States)

Yang, Yong

2008-12-01

In an actual levitation system composed of high temperature superconductors (HTSs) and permanent magnets (PMs), the levitating bodies may traverse in arbitrary directions. Many previous researchers assumed that the levitating bodies moved in a vertical direction or a lateral direction in order to simplify the problem. In this paper, the vertical and lateral forces acting on the PM are calculated by the modified frozen-image method when a PM above an HTS traverses in arbitrary directions. In order to study the effects of the movement directions on the vertical and lateral forces, comparisons of the forces that act on a PM traversing in a tilted direction with those that act on a PM traversing in a vertical direction or a lateral direction have been presented.

Angular quadrature generator for neutron transport SN calculations in slab geometry with arbitrary arithmetic precision

International Nuclear Information System (INIS)

Dominguez, Dany S.; Oliveira, Francisco B.S.; Barros, Ricardo C.

2003-01-01

We present in this paper a multiplatform computational code to calculate elements of Gauss-Legendre angular quadrature sets of arbitrary order used in slab-geometry discrete ordinates (S N ) formulation of neutron transport equation. In the code, the values can be computed with arbitrary arithmetic precision based on the approach of exact computing floating-point numbers. Calculation routines have been developed in the common language ANSI C using standard compiler gcc and the libraries of the open code GMP (GNU Multi precision Library). The code has a graphical interface in order to facilitate user interaction and numerical results analysis. The code architecture allows it to run on different platforms such as Unix, Linux and Windows. Numerical results and performance measures are also given. (author)

Vertical and lateral forces when a permanent magnet above a superconductor traverses in arbitrary directions

Energy Technology Data Exchange (ETDEWEB)

Yang Yong [Key Laboratory of Applied Superconductivity, Chinese Academy of Sciences, Beijing 100190 (China); Institute of Electrical Engineering, Chinese Academy of Sciences, Beijing 100190 (China)], E-mail: yy@mail.iee.ac.cn

2008-12-15

In an actual levitation system composed of high temperature superconductors (HTSs) and permanent magnets (PMs), the levitating bodies may traverse in arbitrary directions. Many previous researchers assumed that the levitating bodies moved in a vertical direction or a lateral direction in order to simplify the problem. In this paper, the vertical and lateral forces acting on the PM are calculated by the modified frozen-image method when a PM above an HTS traverses in arbitrary directions. In order to study the effects of the movement directions on the vertical and lateral forces, comparisons of the forces that act on a PM traversing in a tilted direction with those that act on a PM traversing in a vertical direction or a lateral direction have been presented.

Vertical and lateral forces when a permanent magnet above a superconductor traverses in arbitrary directions

International Nuclear Information System (INIS)

Yang Yong

2008-01-01

In an actual levitation system composed of high temperature superconductors (HTSs) and permanent magnets (PMs), the levitating bodies may traverse in arbitrary directions. Many previous researchers assumed that the levitating bodies moved in a vertical direction or a lateral direction in order to simplify the problem. In this paper, the vertical and lateral forces acting on the PM are calculated by the modified frozen-image method when a PM above an HTS traverses in arbitrary directions. In order to study the effects of the movement directions on the vertical and lateral forces, comparisons of the forces that act on a PM traversing in a tilted direction with those that act on a PM traversing in a vertical direction or a lateral direction have been presented.

Factorization of a class of almost linear second-order differential equations

International Nuclear Information System (INIS)

Estevez, P G; Kuru, S; Negro, J; Nieto, L M

2007-01-01

A general type of almost linear second-order differential equations, which are directly related to several interesting physical problems, is characterized. The solutions of these equations are obtained using the factorization technique, and their non-autonomous invariants are also found by means of scale transformations

Some oscillation criteria for the second-order linear delay differential equation

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2011-01-01

Roč. 136, č. 2 (2011), s. 195-204 ISSN 0862-7959 Institutional research plan: CEZ:AV0Z10190503 Keywords : second-order linear differential equation with a delay * oscillatory solution Subject RIV: BA - General Mathematics http://www.dml.cz/handle/10338.dmlcz/141582

On detection and automatic tracking of butt weld line in thin wall pipe welding by a mobile robot with visual sensor

International Nuclear Information System (INIS)

Suga, Yasuo; Ishii, Hideaki; Muto, Akifumi

1992-01-01

An automatic pipe welding mobile robot system with visual sensor was constructed. The robot can move along a pipe, and detect the weld line to be welded by visual sensor. Moreover, in order to make an automatic welding, the welding torch can track the butt weld line of the pipes at a constant speed by rotating the robot head. Main results obtained are summarized as follows: 1) Using a proper lighting fixed in front of the CCD camera, the butt weld line of thin wall pipes can be recongnized stably. In this case, the root gap should be approximately 0.5 mm. 2) In order to detect the weld line stably during moving along the pipe, a brightness distribution measured by the CCD camera should be subjected to smoothing and differentiating and then the weld line is judged by the maximum and minimum values of the differentials. 3) By means of the basic robot system with a visual sensor controlled by a personal computer, the detection and in-process automatic tracking of a weld line are possible. The average tracking error was approximately 0.2 mm and maximum error 0.5 mm and the welding speed was held at a constant value with error of about 0.1 cm/min. (author)

Efficient, Differentially Private Point Estimators

OpenAIRE

Smith, Adam

2008-01-01

Differential privacy is a recent notion of privacy for statistical databases that provides rigorous, meaningful confidentiality guarantees, even in the presence of an attacker with access to arbitrary side information. We show that for a large class of parametric probability models, one can construct a differentially private estimator whose distribution converges to that of the maximum likelihood estimator. In particular, it is efficient and asymptotically unbiased. This result provides (furt...

Oscillation of two-dimensional linear second-order differential systems

International Nuclear Information System (INIS)

Kwong, M.K.; Kaper, H.G.

1985-01-01

This article is concerned with the oscillatory behavior at infinity of the solution y: [a, ∞) → R 2 of a system of two second-order differential equations, y''(t) + Q(t) y(t) = 0, t epsilon[a, ∞); Q is a continuous matrix-valued function on [a, ∞) whose values are real symmetric matrices of order 2. It is shown that the solution is oscillatory at infinity if the largest eigenvalue of the matrix integral/sub a//sup t/ Q(s) ds tends to infinity as t → ∞. This proves a conjecture of D. Hinton and R.T. Lewis for the two-dimensional case. Furthermore, it is shown that considerably weaker forms of the condition still suffice for oscillatory behavior at infinity. 7 references

Real-Time Automatic Fetal Brain Extraction in Fetal MRI by Deep Learning

OpenAIRE

Salehi, Seyed Sadegh Mohseni; Hashemi, Seyed Raein; Velasco-Annis, Clemente; Ouaalam, Abdelhakim; Estroff, Judy A.; Erdogmus, Deniz; Warfield, Simon K.; Gholipour, Ali

2017-01-01

Brain segmentation is a fundamental first step in neuroimage analysis. In the case of fetal MRI, it is particularly challenging and important due to the arbitrary orientation of the fetus, organs that surround the fetal head, and intermittent fetal motion. Several promising methods have been proposed but are limited in their performance in challenging cases and in real-time segmentation. We aimed to develop a fully automatic segmentation method that independently segments sections of the feta...

On integration of the first order differential equations in a finite terms

International Nuclear Information System (INIS)

Malykh, M D

2017-01-01

There are several approaches to the description of the concept called briefly as integration of the first order differential equations in a finite terms or symbolical integration. In the report three of them are considered: 1.) finding of a rational integral (Beaune or Poincaré problem), 2.) integration by quadratures and 3.) integration when the general solution of given differential equation is an algebraical function of a constant (Painlevé problem). Their realizations in Sage are presented. (paper)

Trace maps for arbitrary substitution sequences

International Nuclear Information System (INIS)

Avishai, Y.

1993-01-01

The discovery of quasi-crystals and their 1-dimensional modeling have led to a deep mathematical study of Schroedinger operators with an arbitrary deterministic potential sequence. In this work we address this problem and find trace maps for an arbitrary substitution sequence. our trace maps have lower dimensionality than those of Kolar and Nori, which make them quite attractive for actual applications. (authors)

Propagation of arbitrary initial wave packets in a quantum parametric oscillator: Instability zones for higher order moments

Science.gov (United States)

Biswas, Subhadip; Chattopadhyay, Rohitashwa; Bhattacharjee, Jayanta K.

2018-05-01

We consider the dynamics of a particle in a parametric oscillator with a view to exploring any quantum feature of the initial wave packet that shows divergent (in time) behaviour for parameter values where the classical motion dynamics of the mean position is bounded. We use Ehrenfest's theorem to explore the dynamics of nth order moment which reduces exactly to a linear non autonomous differential equation of order n + 1. It is found that while the width and skewness of the packet is unbounded exactly in the zones where the classical motion is unbounded, the kurtosis of an initially non-gaussian wave packet can become infinitely large in certain additional zones. This implies that the shape of the wave packet can change drastically with time in these zones.

Arbitrary Dimension Convection-Diffusion Schemes for Space-Time Discretizations

Energy Technology Data Exchange (ETDEWEB)

Bank, Randolph E. [Univ. of California, San Diego, CA (United States); Vassilevski, Panayot S. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Zikatanov, Ludmil T. [Bulgarian Academy of Sciences, Sofia (Bulgaria)

2016-01-20

This note proposes embedding a time dependent PDE into a convection-diffusion type PDE (in one space dimension higher) with singularity, for which two discretization schemes, the classical streamline-diffusion and the EAFE (edge average finite element) one, are investigated in terms of stability and error analysis. The EAFE scheme, in particular, is extended to be arbitrary order which is of interest on its own. Numerical results, in combined space-time domain demonstrate the feasibility of the proposed approach.

Numerical methods for the solution of ordinary differential equations

International Nuclear Information System (INIS)

Azeem, M.

1999-01-01

The ode 113 code solves non-stiff differential equations and is a fully variable step, variable order, PECE implementation in terms of modified divided differences of Adams-Bashforth-Moulton family of formulas of order 1-12. The main objectives of this project were to modify PECE mode of ode 113 into PEC mode, study the variable step size and variable order strategy of both the modes and finally, develop the switching strategy between both PECE and PEC modes to minimize the cost of solving the ordinary differential equations. Using some test problems (including stiff, mild stiff and non-stiff), it was found that the PEC mode was more efficient for non-stiff problems at crude and intermediate tolerances and the PECE mode for all problems at the stringent tolerance. An automatic switching strategy was developed using the results observed from the step size and order plots of all the test problems for both the modes and gave the optimum results. (author)

Fractional Order Stochastic Differential Equation with Application in European Option Pricing

Directory of Open Access Journals (Sweden)

Qing Li

2014-01-01

Full Text Available Memory effect is an important phenomenon in financial systems, and a number of research works have been carried out to study the long memory in the financial markets. In recent years, fractional order ordinary differential equation is used as an effective instrument for describing the memory effect in complex systems. In this paper, we establish a fractional order stochastic differential equation (FSDE model to describe the effect of trend memory in financial pricing. We, then, derive a European option pricing formula based on the FSDE model and prove the existence of the trend memory (i.e., the mean value function in the option pricing formula when the Hurst index is between 0.5 and 1. In addition, we make a comparison analysis between our proposed model, the classic Black-Scholes model, and the stochastic model with fractional Brownian motion. Numerical results suggest that our model leads to more accurate and lower standard deviation in the empirical study.

Block Hybrid Collocation Method with Application to Fourth Order Differential Equations

Directory of Open Access Journals (Sweden)

Lee Ken Yap

2015-01-01

Full Text Available The block hybrid collocation method with three off-step points is proposed for the direct solution of fourth order ordinary differential equations. The interpolation and collocation techniques are applied on basic polynomial to generate the main and additional methods. These methods are implemented in block form to obtain the approximation at seven points simultaneously. Numerical experiments are conducted to illustrate the efficiency of the method. The method is also applied to solve the fourth order problem from ship dynamics.

The convergence of the order sequence and the solution function sequence on fractional partial differential equation

Science.gov (United States)

Rusyaman, E.; Parmikanti, K.; Chaerani, D.; Asefan; Irianingsih, I.

2018-03-01

One of the application of fractional ordinary differential equation is related to the viscoelasticity, i.e., a correlation between the viscosity of fluids and the elasticity of solids. If the solution function develops into function with two or more variables, then its differential equation must be changed into fractional partial differential equation. As the preliminary study for two variables viscoelasticity problem, this paper discusses about convergence analysis of function sequence which is the solution of the homogenous fractional partial differential equation. The method used to solve the problem is Homotopy Analysis Method. The results show that if given two real number sequences (αn) and (βn) which converge to α and β respectively, then the solution function sequences of fractional partial differential equation with order (αn, βn) will also converge to the solution function of fractional partial differential equation with order (α, β).

A novel ultrawideband FDTD numerical modeling of ground penetrating radar on arbitrary dispersive soils

NARCIS (Netherlands)

Mescia, L.; Bia, P.; Caratelli, D.

2017-01-01

A novel two-dimensional (2-D) finite-difference timedomain algorithm for modeling ultrawideband pulse propagation in arbitrary dispersive soils is presented. The soil dispersion is modeled by general power law series representation, accounting for multiple higher order dispersive relaxation

Teleportation with an Arbitrary Mixed Resource as a Trace-Preserving Quantum Channel

Institute of Scientific and Technical Information of China (English)

Sergio Albeverio; FEI Shao-Ming; YANG Wen-Li

2002-01-01

General conditions are given in order to perform a perfect teleportation process in the case where theHilbert spaces involved have different dimensions. An explicit expression is obtained for the quantum channel associatedwith the standard teleportation protocol To with an arbitrary mixed state resource. The transmission fidelity of thecorresponding quantum channel is given.

Myshkis type oscillation criteria for second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2015-01-01

Roč. 178, č. 1 (2015), s. 143-161 ISSN 0026-9255 Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillation criteria Subject RIV: BA - General Mathematics Impact factor: 0.664, year: 2015 http://link.springer.com/article/10.1007%2Fs00605-014-0719-y

Renormalization-group improved fully differential cross sections for top pair production

International Nuclear Information System (INIS)

Broggio, A.; Papanastasiou, A.S.; Signer, A.; Zuerich Univ.

2014-07-01

We extend approximate next-to-next-to-leading order results for top-pair production to include the semi-leptonic decays of top quarks in the narrow-width approximation. The new hard-scattering kernels are implemented in a fully differential parton-level Monte Carlo that allows for the study of any IR-safe observable constructed from the momenta of the decay products of the top. Our best predictions are given by approximate NNLO corrections in the production matched to a fixed order calculation with NLO corrections in both the production and decay subprocesses. Being fully differential enables us to make comparisons between approximate results derived via different (PIM and 1PI) kinematics for arbitrary distributions. These comparisons reveal that the renormalization-group framework, from which the approximate results are derived, is rather robust in the sense that applying a realistic error estimate allows us to obtain a reliable prediction with a reduced theoretical error for generic observables and analysis cuts.

Fully automatic time-window selection using machine learning for global adjoint tomography

Science.gov (United States)

Chen, Y.; Hill, J.; Lei, W.; Lefebvre, M. P.; Bozdag, E.; Komatitsch, D.; Tromp, J.

2017-12-01

Selecting time windows from seismograms such that the synthetic measurements (from simulations) and measured observations are sufficiently close is indispensable in a global adjoint tomography framework. The increasing amount of seismic data collected everyday around the world demands "intelligent" algorithms for seismic window selection. While the traditional FLEXWIN algorithm can be "automatic" to some extent, it still requires both human input and human knowledge or experience, and thus is not deemed to be fully automatic. The goal of intelligent window selection is to automatically select windows based on a learnt engine that is built upon a huge number of existing windows generated through the adjoint tomography project. We have formulated the automatic window selection problem as a classification problem. All possible misfit calculation windows are classified as either usable or unusable. Given a large number of windows with a known selection mode (select or not select), we train a neural network to predict the selection mode of an arbitrary input window. Currently, the five features we extract from the windows are its cross-correlation value, cross-correlation time lag, amplitude ratio between observed and synthetic data, window length, and minimum STA/LTA value. More features can be included in the future. We use these features to characterize each window for training a multilayer perceptron neural network (MPNN). Training the MPNN is equivalent to solve a non-linear optimization problem. We use backward propagation to derive the gradient of the loss function with respect to the weighting matrices and bias vectors and use the mini-batch stochastic gradient method to iteratively optimize the MPNN. Numerical tests show that with a careful selection of the training data and a sufficient amount of training data, we are able to train a robust neural network that is capable of detecting the waveforms in an arbitrary earthquake data with negligible detection error

High Weak Order Methods for Stochastic Differential Equations Based on Modified Equations

KAUST Repository

Abdulle, Assyr; Cohen, David; Vilmart, Gilles; Zygalakis, Konstantinos C.

2012-01-01

© 2012 Society for Industrial and Applied Mathematics. Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration

A dose assessment method for arbitrary geometries with virtual reality in the nuclear facilities decommissioning

Science.gov (United States)

Chao, Nan; Liu, Yong-kuo; Xia, Hong; Ayodeji, Abiodun; Bai, Lu

2018-03-01

During the decommissioning of nuclear facilities, a large number of cutting and demolition activities are performed, which results in a frequent change in the structure and produce many irregular objects. In order to assess dose rates during the cutting and demolition process, a flexible dose assessment method for arbitrary geometries and radiation sources was proposed based on virtual reality technology and Point-Kernel method. The initial geometry is designed with the three-dimensional computer-aided design tools. An approximate model is built automatically in the process of geometric modeling via three procedures namely: space division, rough modeling of the body and fine modeling of the surface, all in combination with collision detection of virtual reality technology. Then point kernels are generated by sampling within the approximate model, and when the material and radiometric attributes are inputted, dose rates can be calculated with the Point-Kernel method. To account for radiation scattering effects, buildup factors are calculated with the Geometric-Progression formula in the fitting function. The effectiveness and accuracy of the proposed method was verified by means of simulations using different geometries and the dose rate results were compared with that derived from CIDEC code, MCNP code and experimental measurements.

A higher-order conservation element solution element method for solving hyperbolic differential equations on unstructured meshes

Science.gov (United States)

Bilyeu, David

This dissertation presents an extension of the Conservation Element Solution Element (CESE) method from second- to higher-order accuracy. The new method retains the favorable characteristics of the original second-order CESE scheme, including (i) the use of the space-time integral equation for conservation laws, (ii) a compact mesh stencil, (iii) the scheme will remain stable up to a CFL number of unity, (iv) a fully explicit, time-marching integration scheme, (v) true multidimensionality without using directional splitting, and (vi) the ability to handle two- and three-dimensional geometries by using unstructured meshes. This algorithm has been thoroughly tested in one, two and three spatial dimensions and has been shown to obtain the desired order of accuracy for solving both linear and non-linear hyperbolic partial differential equations. The scheme has also shown its ability to accurately resolve discontinuities in the solutions. Higher order unstructured methods such as the Discontinuous Galerkin (DG) method and the Spectral Volume (SV) methods have been developed for one-, two- and three-dimensional application. Although these schemes have seen extensive development and use, certain drawbacks of these methods have been well documented. For example, the explicit versions of these two methods have very stringent stability criteria. This stability criteria requires that the time step be reduced as the order of the solver increases, for a given simulation on a given mesh. The research presented in this dissertation builds upon the work of Chang, who developed a fourth-order CESE scheme to solve a scalar one-dimensional hyperbolic partial differential equation. The completed research has resulted in two key deliverables. The first is a detailed derivation of a high-order CESE methods on unstructured meshes for solving the conservation laws in two- and three-dimensional spaces. The second is the code implementation of these numerical methods in a computer code. For

A note on monotone solutions for a nonconvex second-order functional differential inclusion

Directory of Open Access Journals (Sweden)

Aurelian Cernea

2011-12-01

Full Text Available The existence of monotone solutions for a second-order functional differential inclusion with Carath\\'{e}odory perturbation is obtained in the case when the multifunction that define the inclusion is upper semicontinuous compact valued and contained in the Fr\\'{e}chet subdifferential of a $\\phi $-convex function of order two.

Stability Analysis for Fractional-Order Linear Singular Delay Differential Systems

Directory of Open Access Journals (Sweden)

Hai Zhang

2014-01-01

Full Text Available We investigate the delay-independently asymptotic stability of fractional-order linear singular delay differential systems. Based on the algebraic approach, the sufficient conditions are presented to ensure the asymptotic stability for any delay parameter. By applying the stability criteria, one can avoid solving the roots of transcendental equations. An example is also provided to illustrate the effectiveness and applicability of the theoretical results.

Design of fractional order differentiator using type-III and type-IV discrete cosine transform

Directory of Open Access Journals (Sweden)

Manjeet Kumar

2017-02-01

Full Text Available In this paper, an interpolation method based on discrete cosine transform (DCT is employed for digital finite impulse response-fractional order differentiator (FIR-FOD design. Here, a fractional order digital differentiator is modeled as finite impulse response (FIR system to get an optimized frequency response that approximates the ideal response of a fractional order differentiator. Next, DCT-III and DCT-IV are utilized to determine the filter coefficients of FIR filter that compute the Fractional derivative of a given signal. To improve the frequency response of the proposed FIR-FOD, the filter coefficients are further modified using windows. Several design examples are presented to demonstrate the superiority of the proposed method. The simulation results have also been compared with the existing FIR-FOD design methods such as DFT interpolation, radial basis function (RBF interpolation, DCT-II interpolation and DST interpolation methods. The result reveals that the proposed FIR-FOD design technique using DCT-III and DCT-IV outperforms DFT interpolation, RBF interpolation, DCT-II interpolation and DST interpolation methods in terms of magnitude error.

Implementation of a Three-Dimensional Pedometer Automatic Accumulating Walking or Jogging Motions in Arbitrary Placement

Directory of Open Access Journals (Sweden)

Jia-Shing Sheu

2014-01-01

Full Text Available This study proposes a method for using a three-axis accelerometer and a single-chip microcontrol unit to implement a three-dimensional (3D pedometer that can automatically identify walking and running motions. The proposed design can calculate the number of walking and running steps down to small numbers of steps and can be easily worn, thus remedying defects of generic mechanical and 3D pedometers. The user’s motion state is calculated using a walk/run mode switching algorithm.

Controllability Results For First Order Impulsive Stochastic Functional Differential Systems with State-Dependent Delay

Directory of Open Access Journals (Sweden)

C. Parthasarathy

2013-03-01

Full Text Available In this paper, we study the controllability results of first order impulsive stochastic differential and neutral differential systems with state-dependent delay by using semigroup theory. The controllability results are derived by the means of Leray-SchauderAlternative fixed point theorem. An example is provided to illustrate the theory.

Stability and square integrability of solutions of nonlinear fourth order differential equations

Directory of Open Access Journals (Sweden)

Moussadek Remili

2016-05-01

Full Text Available The aim of the present paper is to establish a new result, which guarantees the asymptotic stability of zero solution and square integrability of solutions and their derivatives to nonlinear differential equations of fourth order.

Higher-order stochastic differential equations and the positive Wigner function

Science.gov (United States)

Drummond, P. D.

2017-12-01

General higher-order stochastic processes that correspond to any diffusion-type tensor of higher than second order are obtained. The relationship of multivariate higher-order stochastic differential equations with tensor decomposition theory and tensor rank is explained. Techniques for generating the requisite complex higher-order noise are proved to exist either using polar coordinates and γ distributions, or from products of Gaussian variates. This method is shown to allow the calculation of the dynamics of the Wigner function, after it is extended to a complex phase space. The results are illustrated physically through dynamical calculations of the positive Wigner distribution for three-mode parametric downconversion, widely used in quantum optics. The approach eliminates paradoxes arising from truncation of the higher derivative terms in Wigner function time evolution. Anomalous results of negative populations and vacuum scattering found in truncated Wigner quantum simulations in quantum optics and Bose-Einstein condensate dynamics are shown not to occur with this type of stochastic theory.

Periodic solutions of first-order functional differential equations in population dynamics

CERN Document Server

Padhi, Seshadev; Srinivasu, P D N

2014-01-01

This book provides cutting-edge results on the existence of multiple positive periodic solutions of first-order functional differential equations. It demonstrates how the Leggett-Williams fixed-point theorem can be applied to study the existence of two or three positive periodic solutions of functional differential equations with real-world applications, particularly with regard to the Lasota-Wazewska model, the Hematopoiesis model, the Nicholsons Blowflies model, and some models with Allee effects. Many interesting sufficient conditions are given for the dynamics that include nonlinear characteristics exhibited by population models. The last chapter provides results related to the global appeal of solutions to the models considered in the earlier chapters. The techniques used in this book can be easily understood by anyone with a basic knowledge of analysis. This book offers a valuable reference guide for students and researchers in the field of differential equations with applications to biology, ecology, a...

The second-order differential phase contrast and its retrieval for imaging with x-ray Talbot interferometry

International Nuclear Information System (INIS)

Yang Yi; Tang Xiangyang

2012-01-01

Purpose: The x-ray differential phase contrast imaging implemented with the Talbot interferometry has recently been reported to be capable of providing tomographic images corresponding to attenuation-contrast, phase-contrast, and dark-field contrast, simultaneously, from a single set of projection data. The authors believe that, along with small-angle x-ray scattering, the second-order phase derivative Φ ″ s (x) plays a role in the generation of dark-field contrast. In this paper, the authors derive the analytic formulae to characterize the contribution made by the second-order phase derivative to the dark-field contrast (namely, second-order differential phase contrast) and validate them via computer simulation study. By proposing a practical retrieval method, the authors investigate the potential of second-order differential phase contrast imaging for extensive applications. Methods: The theoretical derivation starts at assuming that the refractive index decrement of an object can be decomposed into δ=δ s +δ f , where δ f corresponds to the object's fine structures and manifests itself in the dark-field contrast via small-angle scattering. Based on the paraxial Fresnel-Kirchhoff theory, the analytic formulae to characterize the contribution made by δ s , which corresponds to the object's smooth structures, to the dark-field contrast are derived. Through computer simulation with specially designed numerical phantoms, an x-ray differential phase contrast imaging system implemented with the Talbot interferometry is utilized to evaluate and validate the derived formulae. The same imaging system is also utilized to evaluate and verify the capability of the proposed method to retrieve the second-order differential phase contrast for imaging, as well as its robustness over the dimension of detector cell and the number of steps in grating shifting. Results: Both analytic formulae and computer simulations show that, in addition to small-angle scattering, the

Groups of integral transforms generated by Lie algebras of second-and higher-order differential operators

International Nuclear Information System (INIS)

Steinberg, S.; Wolf, K.B.

1979-01-01

The authors study the construction and action of certain Lie algebras of second- and higher-order differential operators on spaces of solutions of well-known parabolic, hyperbolic and elliptic linear differential equations. The latter include the N-dimensional quadratic quantum Hamiltonian Schroedinger equations, the one-dimensional heat and wave equations and the two-dimensional Helmholtz equation. In one approach, the usual similarity first-order differential operator algebra of the equation is embedded in the larger one, which appears as a quantum-mechanical dynamic algebra. In a second approach, the new algebra is built as the time evolution of a finite-transformation algebra on the initial conditions. In a third approach, the algebra to inhomogeneous similarity algebra is deformed to a noncompact classical one. In every case, we can integrate the algebra to a Lie group of integral transforms acting effectively on the solution space of the differential equation. (author)

The Algorithm for Algorithms: An Evolutionary Algorithm Based on Automatic Designing of Genetic Operators

Directory of Open Access Journals (Sweden)

Dazhi Jiang

2015-01-01

Full Text Available At present there is a wide range of evolutionary algorithms available to researchers and practitioners. Despite the great diversity of these algorithms, virtually all of the algorithms share one feature: they have been manually designed. A fundamental question is “are there any algorithms that can design evolutionary algorithms automatically?” A more complete definition of the question is “can computer construct an algorithm which will generate algorithms according to the requirement of a problem?” In this paper, a novel evolutionary algorithm based on automatic designing of genetic operators is presented to address these questions. The resulting algorithm not only explores solutions in the problem space like most traditional evolutionary algorithms do, but also automatically generates genetic operators in the operator space. In order to verify the performance of the proposed algorithm, comprehensive experiments on 23 well-known benchmark optimization problems are conducted. The results show that the proposed algorithm can outperform standard differential evolution algorithm in terms of convergence speed and solution accuracy which shows that the algorithm designed automatically by computers can compete with the algorithms designed by human beings.

Maillet type theorem for singular first order nonlinear partial differential equations of totally characteristic type. Part II

Directory of Open Access Journals (Sweden)

Akira Shirai

2015-01-01

Full Text Available In this paper, we study the following nonlinear first order partial differential equation: \\[f(t,x,u,\\partial_t u,\\partial_x u=0\\quad\\text{with}\\quad u(0,x\\equiv 0.\\] The purpose of this paper is to determine the estimate of Gevrey order under the condition that the equation is singular of a totally characteristic type. The Gevrey order is indicated by the rate of divergence of a formal power series. This paper is a continuation of the previous papers [Convergence of formal solutions of singular first order nonlinear partial differential equations of totally characteristic type, Funkcial. Ekvac. 45 (2002, 187-208] and [Maillet type theorem for singular first order nonlinear partial differential equations of totally characteristic type, Surikaiseki Kenkyujo Kokyuroku, Kyoto University 1431 (2005, 94-106]. Especially the last-mentioned paper is regarded as part I of this paper.

Identification of fractional-order systems via a switching differential evolution subject to noise perturbations

Energy Technology Data Exchange (ETDEWEB)

Zhu, Wu, E-mail: dtzhuwu@gmail.com [College of Information Science and Technology, Donghua University, Shanghai 201620 (China); Fang, Jian-an [College of Information Science and Technology, Donghua University, Shanghai 201620 (China); Tang, Yang, E-mail: yang.tang@pik-potsdam.de [Institute of Physics, Humboldt University, Berlin 12489 (Germany); Potsdam Institute for Climate Impact Research, Potsdam 14415 (Germany); Research Institute for Intelligent Control and System, Harbin Institute of Technology, Harbin 150006 (China); Zhang, Wenbing [Institute of Textiles and Clothing, The Hong Kong Polytechnic University, Hong Kong (China); Xu, Yulong [College of Information Science and Technology, Donghua University, Shanghai 201620 (China)

2012-10-01

In this Letter, a differential evolution variant, called switching DE (SDE), has been employed to estimate the orders and parameters in incommensurate fractional-order chaotic systems. The proposed algorithm includes a switching population utilization strategy, where the population size is adjusted dynamically based on the solution-searching status. Thus, this adaptive control method realizes the identification of fractional-order Lorenz, Lü and Chen systems in both deterministic and stochastic environments, respectively. Numerical simulations are provided, where comparisons are made with five other State-of-the-Art evolutionary algorithms (EAs) to verify the effectiveness of the proposed method. -- Highlights: ► Switching population utilization strategy is applied for differential evolution. ► The parameters are estimated in both deterministic and stochastic environments. ► Comparisons with five other EAs verify the effectiveness of the proposed method.

Identification of fractional-order systems via a switching differential evolution subject to noise perturbations

International Nuclear Information System (INIS)

Zhu, Wu; Fang, Jian-an; Tang, Yang; Zhang, Wenbing; Xu, Yulong

2012-01-01

In this Letter, a differential evolution variant, called switching DE (SDE), has been employed to estimate the orders and parameters in incommensurate fractional-order chaotic systems. The proposed algorithm includes a switching population utilization strategy, where the population size is adjusted dynamically based on the solution-searching status. Thus, this adaptive control method realizes the identification of fractional-order Lorenz, Lü and Chen systems in both deterministic and stochastic environments, respectively. Numerical simulations are provided, where comparisons are made with five other State-of-the-Art evolutionary algorithms (EAs) to verify the effectiveness of the proposed method. -- Highlights: ► Switching population utilization strategy is applied for differential evolution. ► The parameters are estimated in both deterministic and stochastic environments. ► Comparisons with five other EAs verify the effectiveness of the proposed method.

Oscillation Criteria of First Order Neutral Delay Differential Equations with Variable Coefficients

Directory of Open Access Journals (Sweden)

Fatima N. Ahmed

2013-01-01

Full Text Available Some new oscillation criteria are given for first order neutral delay differential equations with variable coefficients. Our results generalize and extend some of the well-known results in the literature. Some examples are considered to illustrate the main results.

Interior and exterior resonances in acoustic scattering. pt. 2 - Targets of arbitrary shape (T-matrix approach)

International Nuclear Information System (INIS)

Uberall, H.; Gaunaurd, G.C.; Tanglis, E.

1983-01-01

The T-matrix approach, which describes the scattering of acoustic waves (or of other waves) from objects of arbitrary shape and geometry, is here 'married' to the resonance scattering theory in order to obtain the (complex) resonance frequencies of an arbitrary shaped target. For the case of nearly impenetrable targets the partial-wave scattering amplitudes are splitted into terms corresponding to 'internal' resonances, plus an apparently nonresonant background amplitude which, however, contains the broad resonances caused by 'external' diffracted (or Franz-type, creeping) waves, in addition to geometrically reflected and refracted (ray) contributions

Initial-value problems for first-order differential recurrence equations with auto-convolution

Directory of Open Access Journals (Sweden)

Mircea Cirnu

2011-01-01

Full Text Available A differential recurrence equation consists of a sequence of differential equations, from which must be determined by recurrence a sequence of unknown functions. In this article, we solve two initial-value problems for some new types of nonlinear (quadratic first order homogeneous differential recurrence equations, namely with discrete auto-convolution and with combinatorial auto-convolution of the unknown functions. In both problems, all initial values form a geometric progression, but in the second problem the first initial value is exempted and has a prescribed form. Some preliminary results showing the importance of the initial conditions are obtained by reducing the differential recurrence equations to algebraic type. Final results about solving the considered initial value problems, are shown by mathematical induction. However, they can also be shown by changing the unknown functions, or by the generating function method. So in a remark, we give a proof of the first theorem by the generating function method.

Integral transformation of the Navier-Stokes equations for laminar flow in channels of arbitrary two-dimensional geometry

International Nuclear Information System (INIS)

Perez Guerrero, Jesus Salvador

1995-01-01

Laminar developing flow in channels of arbitrary geometry was studied by solving the Navier-Stokes equations in the stream function-only formulation through the Generalized Integral Transform Technique (GITT). The stream function is expanded in an infinite system based on eigenfunctions obtained by considering solely the diffusive terms of the original formulation. The Navier-Stokes equations are transformed into an infinite system of ordinary differential equations, by using the transformation and inversion formulae. For computational purposes, the infinite series is truncated, according to an automatic error control procedure. The ordinary differential is solved through well-established scientific subroutines from widely available mathematical libraries. The classical problem of developing flow between parallel-plates is analysed first, as for both uniform and irrotational inlet conditions. The effect of truncating the duct length in the accuracy of the obtained solution is studied. A convergence analysis of the results obtained by the GITT is performed and compared with results obtained by finite difference and finite element methods, for different values of Reynolds number. The problem of flow over a backward-facing step then follows. Comparisons with experimental results in the literature indicate an excellent agreement. The numerical co-validation was established for a test case, and perfect agreement is reached against results considered as benchmarks in the recent literature. The results were shown to be physically more reasonable than others obtained by purely numerical methods, in particular for situations where three-dimensional effects are identified. Finally, a test problem for an irregular by shoped duct was studied and compared against results found in the literature, with good agreement and excellent convergence rates for the stream function field along the whole channel, for different values of Reynolds number. (author)

On oscillations of solutions to second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2013-01-01

Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001.xml?format=INT

FORCED OSCILLATIONS OF SECOND ORDER SUPER-LINEAR DIFFERENTIAL EQUATION WITH IMPULSES

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

At first,by means of Kartsatos technique,we reduce the impulsive differential equation to a second order nonlinear impulsive homogeneous equation.We find some suitable impulse functions such that all the solutions to the equation are oscillatory.Several criteria on the oscillations of solutions are given.At last,we give an example to demonstrate our results.

Triple positive  solutions of nth order impulsive integro-differential equations

Directory of Open Access Journals (Sweden)

Zeyong Qiu

2011-07-01

Full Text Available In this paper, we prove the existence of at least three positive solutions of boundary value problems for nth order nonlinear impulsive integro-differential equations of mixed type on infinite interval with infinite number of impulsive times. Our results are obtained by applying a new fixed point theorem introduced by Avery and Peterson.

On oscillations of solutions to second-order linear delay differential equations

Czech Academy of Sciences Publication Activity Database

Opluštil, Z.; Šremr, Jiří

2013-01-01

Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001. xml ?format=INT

Color ordering in QCD

OpenAIRE

Schuster, Theodor

2013-01-01

We derive color decompositions of arbitrary tree and one-loop QCD amplitudes into color ordered objects called primitive amplitudes. Furthermore, we derive general fermion flip and reversion identities spanning the null space among the primitive amplitudes and use them to prove that all color ordered tree amplitudes of massless QCD can be written as linear combinations of color ordered tree amplitudes of $\\mathcal{N}=4$ super Yang-Mills theory.

Numerical solution of three-dimensional magnetic differential equations

International Nuclear Information System (INIS)

Reiman, A.H.; Greenside, H.S.

1987-02-01

A computer code is described that solves differential equations of the form B . del f = h for a single-valued solution f, given a toroidal three-dimensional divergence-free field B and a single-valued function h. The code uses a new algorithm that Fourier decomposes a given function in a set of flux coordinates in which the field lines are straight. The algorithm automatically adjusts the required integration lengths to compensate for proximity to low order rational surfaces. Applying this algorithm to the Cartesian coordinates defines a transformation to magnetic coordinates, in which the magnetic differential equation can be accurately solved. Our method is illustrated by calculating the Pfirsch-Schlueter currents for a stellarator

Automatic text summarization

CERN Document Server

Torres Moreno, Juan Manuel

2014-01-01

This new textbook examines the motivations and the different algorithms for automatic document summarization (ADS). We performed a recent state of the art. The book shows the main problems of ADS, difficulties and the solutions provided by the community. It presents recent advances in ADS, as well as current applications and trends. The approaches are statistical, linguistic and symbolic. Several exemples are included in order to clarify the theoretical concepts.  The books currently available in the area of Automatic Document Summarization are not recent. Powerful algorithms have been develop

EXISTENCE OF SOLUTION TO NONLINEAR SECOND ORDER NEUTRAL STOCHASTIC DIFFERENTIAL EQUATIONS WITH DELAY

Institute of Scientific and Technical Information of China (English)

无

2010-01-01

This paper is concerned with nonlinear second order neutral stochastic differential equations with delay in a Hilbert space. Sufficient conditions for the existence of solution to the system are obtained by Picard iterations.

Projection operator and propagator for an arbitrary integral spin

CERN Document Server

Huang Shi Zhong; Wu Ning; Zheng Zhi Peng

2002-01-01

Based on the solution of the Bargmann-Wigner equation for an arbitrary integral spin, a direct derivation of the projection operator and propagator for an arbitrary integral spin is presented. The explicit form for the spin projection operators constructed by Behrends and Fronsdal is confirmed. The commutation rules and a general expression for the Feynman propagator for a free particle of arbitrary integral spin are deduced

LSODE, 1. Order Stiff or Non-Stiff Ordinary Differential Equations System Initial Value Problems

International Nuclear Information System (INIS)

Hindmarsh, A.C.; Petzold, L.R.

2005-01-01

1 - Description of program or function: LSODE (Livermore Solver for Ordinary Differential Equations) solves stiff and non-stiff systems of the form dy/dt = f. In the stiff case, it treats the Jacobian matrix df/dy as either a dense (full) or a banded matrix, and as either user-supplied or internally approximated by difference quotients. It uses Adams methods (predictor-corrector) in the non-stiff case, and Backward Differentiation Formula (BDF) methods (the Gear methods) in the stiff case. The linear systems that arise are solved by direct methods (LU factor/solve). The LSODE source is commented extensively to facilitate modification. Both a single-precision version and a double-precision version are available. 2 - Methods: It is assumed that the ODEs are given explicitly, so that the system can be written in the form dy/dt = f(t,y), where y is the vector of dependent variables, and t is the independent variable. LSODE contains two variable-order, variable- step (with interpolatory step-changing) integration methods. The first is the implicit Adams or non-stiff method, of orders one through twelve. The second is the backward differentiation or stiff method (or BDF method, or Gear's method), of orders one through five. 3 - Restrictions on the complexity of the problem: The differential equations must be given in explicit form, i.e., dy/dt = f(y,t). Problems with intermittent high-speed transients may cause inefficient or unstable performance

Composition between mecd and runge-Kutta algorithm method for large system of second order differential equations

International Nuclear Information System (INIS)

Supriyono; Miyoshi, T.

1997-01-01

NECD Method and runge-Kutta method for large system of second order ordinary differential equations in comparing algorithm. The paper introduce a extrapolation method used for solving the large system of second order ordinary differential equation. We call this method the modified extrapolated central difference (MECD) method. for the accuracy and efficiency MECD method. we compare the method with 4-th order runge-Kutta method. The comparison results show that, this method has almost the same accuracy as the 4-th order runge-Kutta method, but the computation time is about half of runge-Kutta. The MECD was declare by the author and Tetsuhiko Miyoshi of the Dept. Applied Science Yamaguchi University Japan

Differential evolution algorithm based automatic generation control for interconnected power systems with

Directory of Open Access Journals (Sweden)

Banaja Mohanty

2014-09-01

Full Text Available This paper presents the design and performance analysis of Differential Evolution (DE algorithm based Proportional–Integral (PI and Proportional–Integral–Derivative (PID controllers for Automatic Generation Control (AGC of an interconnected power system. Initially, a two area thermal system with governor dead-band nonlinearity is considered for the design and analysis purpose. In the proposed approach, the design problem is formulated as an optimization problem control and DE is employed to search for optimal controller parameters. Three different objective functions are used for the design purpose. The superiority of the proposed approach has been shown by comparing the results with a recently published Craziness based Particle Swarm Optimization (CPSO technique for the same interconnected power system. It is noticed that, the dynamic performance of DE optimized PI controller is better than CPSO optimized PI controllers. Additionally, controller parameters are tuned at different loading conditions so that an adaptive gain scheduling control strategy can be employed. The study is further extended to a more realistic network of two-area six unit system with different power generating units such as thermal, hydro, wind and diesel generating units considering boiler dynamics for thermal plants, Generation Rate Constraint (GRC and Governor Dead Band (GDB non-linearity.

Differential effects of exogenous and endogenous attention on second-order texture contrast sensitivity

Science.gov (United States)

Barbot, Antoine; Landy, Michael S.; Carrasco, Marisa

2012-01-01

The visual system can use a rich variety of contours to segment visual scenes into distinct perceptually coherent regions. However, successfully segmenting an image is a computationally expensive process. Previously we have shown that exogenous attention—the more automatic, stimulus-driven component of spatial attention—helps extract contours by enhancing contrast sensitivity for second-order, texture-defined patterns at the attended location, while reducing sensitivity at unattended locations, relative to a neutral condition. Interestingly, the effects of exogenous attention depended on the second-order spatial frequency of the stimulus. At parafoveal locations, attention enhanced second-order contrast sensitivity to relatively high, but not to low second-order spatial frequencies. In the present study we investigated whether endogenous attention—the more voluntary, conceptually-driven component of spatial attention—affects second-order contrast sensitivity, and if so, whether its effects are similar to those of exogenous attention. To that end, we compared the effects of exogenous and endogenous attention on the sensitivity to second-order, orientation-defined, texture patterns of either high or low second-order spatial frequencies. The results show that, like exogenous attention, endogenous attention enhances second-order contrast sensitivity at the attended location and reduces it at unattended locations. However, whereas the effects of exogenous attention are a function of the second-order spatial frequency content, endogenous attention affected second-order contrast sensitivity independent of the second-order spatial frequency content. This finding supports the notion that both exogenous and endogenous attention can affect second-order contrast sensitivity, but that endogenous attention is more flexible, benefitting performance under different conditions. PMID:22895879

Differential effects of exogenous and endogenous attention on second-order texture contrast sensitivity.

Science.gov (United States)

Barbot, Antoine; Landy, Michael S; Carrasco, Marisa

2012-08-15

The visual system can use a rich variety of contours to segment visual scenes into distinct perceptually coherent regions. However, successfully segmenting an image is a computationally expensive process. Previously we have shown that exogenous attention--the more automatic, stimulus-driven component of spatial attention--helps extract contours by enhancing contrast sensitivity for second-order, texture-defined patterns at the attended location, while reducing sensitivity at unattended locations, relative to a neutral condition. Interestingly, the effects of exogenous attention depended on the second-order spatial frequency of the stimulus. At parafoveal locations, attention enhanced second-order contrast sensitivity to relatively high, but not to low second-order spatial frequencies. In the present study we investigated whether endogenous attention-the more voluntary, conceptually-driven component of spatial attention--affects second-order contrast sensitivity, and if so, whether its effects are similar to those of exogenous attention. To that end, we compared the effects of exogenous and endogenous attention on the sensitivity to second-order, orientation-defined, texture patterns of either high or low second-order spatial frequencies. The results show that, like exogenous attention, endogenous attention enhances second-order contrast sensitivity at the attended location and reduces it at unattended locations. However, whereas the effects of exogenous attention are a function of the second-order spatial frequency content, endogenous attention affected second-order contrast sensitivity independent of the second-order spatial frequency content. This finding supports the notion that both exogenous and endogenous attention can affect second-order contrast sensitivity, but that endogenous attention is more flexible, benefitting performance under different conditions.

Third-order ordinary differential equations Y”' = f(x, y, y'', y′”) with ...

African Journals Online (AJOL)

dimensional symmetry algebra. Mathematics Subject Classication (2010): 34A05, 34A25, 53A55, 76M60. Key words: Linearization, third order ODEs, point transformation, contact transformation, Lie symmetries, relative differential invariants.

Perceptually stable regions for arbitrary polygons.

Science.gov (United States)

Rocha, J

2003-01-01

Zou and Yan have recently developed a skeletonization algorithm of digital shapes based on a regularity/singularity analysis; they use the polygon whose vertices are the boundary pixels of the image to compute a constrained Delaunay triangulation (CDT) in order to find local symmetries and stable regions. Their method has produced good results but it is slow since its complexity depends on the number of contour pixels. This paper presents an extension of their technique to handle arbitrary polygons, not only polygons of short edges. Consequently, not only can we achieve results as good as theirs for digital images, but we can also compute skeletons of polygons of any number of edges. Since we can handle polygonal approximations of figures, the skeletons are more resilient to noise and faster to process.

DAF: differential ACE filtering image quality assessment by automatic color equalization

Science.gov (United States)

Ouni, S.; Chambah, M.; Saint-Jean, C.; Rizzi, A.

2008-01-01

Ideally, a quality assessment system would perceive and measure image or video impairments just like a human being. But in reality, objective quality metrics do not necessarily correlate well with perceived quality [1]. Plus, some measures assume that there exists a reference in the form of an "original" to compare to, which prevents their usage in digital restoration field, where often there is no reference to compare to. That is why subjective evaluation is the most used and most efficient approach up to now. But subjective assessment is expensive, time consuming and does not respond, hence, to the economic requirements [2,3]. Thus, reliable automatic methods for visual quality assessment are needed in the field of digital film restoration. The ACE method, for Automatic Color Equalization [4,6], is an algorithm for digital images unsupervised enhancement. It is based on a new computational approach that tries to model the perceptual response of our vision system merging the Gray World and White Patch equalization mechanisms in a global and local way. Like our vision system ACE is able to adapt to widely varying lighting conditions, and to extract visual information from the environment efficaciously. Moreover ACE can be run in an unsupervised manner. Hence it is very useful as a digital film restoration tool since no a priori information is available. In this paper we deepen the investigation of using the ACE algorithm as a basis for a reference free image quality evaluation. This new metric called DAF for Differential ACE Filtering [7] is an objective quality measure that can be used in several image restoration and image quality assessment systems. In this paper, we compare on different image databases, the results obtained with DAF and with some subjective image quality assessments (Mean Opinion Score MOS as measure of perceived image quality). We study also the correlation between objective measure and MOS. In our experiments, we have used for the first image

A differential equation approach to minor loops in the Jiles-Atherton hysteresis model

International Nuclear Information System (INIS)

Carpenter, K.H.

1991-01-01

Jiles and Atherton, in a series of papers, present physically based differential equations for magnetization in ferromagnetic materials. however, if one directly solves their differential equations, the minor loops obtained can be negative slopes, which is a nonphysical behavior. Only one of their papers gives a method for obtaining minor loops, and the method does not use a differential equation, but requires a priori knowledge of the loop turning points in order to obtain a scale factor and offset which allow a portion of a major loop to serve as a portion of a minor one. In this paper, the reason for the failure of the differential equations to yield physical minor loops is explained, and a modified solution for minor loops is presented which retains the features of Jiles and Atherton's original minor loops, but only requires knowledge of the initial point on each portion of the loop to obtain the solution. This yields a general differential equation formulation for the Jiles-Atherton theory that can be used with circuit simulations having arbitrary excitations and initial conditions for ferromagnetic components

Growth and Zeros of Meromorphic Solutions to Second-Order Linear Differential Equations

Directory of Open Access Journals (Sweden)

Maamar Andasmas

2016-04-01

Full Text Available The main purpose of this article is to investigate the growth of meromorphic solutions to homogeneous and non-homogeneous second order linear differential equations f00+Af0+Bf = F, where A(z, B (z and F (z are meromorphic functions with finite order having only finitely many poles. We show that, if there exist a positive constants σ > 0, α > 0 such that |A(z| ≥ eα|z|σ as |z| → +∞, z ∈ H, where dens{|z| : z ∈ H} > 0 and ρ = max{ρ(B, ρ(F} < σ, then every transcendental meromorphic solution f has an infinite order. Further, we give some estimates of their hyper-order, exponent and hyper-exponent of convergence of distinct zeros.

Buckling analysis of laminated plates using the extended Kantorovich method and a system of first-order differential equations

Energy Technology Data Exchange (ETDEWEB)

Singhatanadgid, Pairod; Jommalai, Panupan [Chulalongkorn University, Bangkok (Thailand)

2016-05-15

The extended Kantorovich method using multi-term displacement functions is applied to the buckling problem of laminated plates with various boundary conditions. The out-of-plane displacement of the buckled plate is written as a series of products of functions of parameter x and functions of parameter y. With known functions in parameter x or parameter y, a set of governing equations and a set of boundary conditions are obtained after applying the variational principle to the total potential energy of the system. The higher order differential equations are then transformed into a set of first-order differential equations and solved for the buckling load and mode. Since the governing equations are first-order differential equations, solutions can be obtained analytically with the out-of-plane displacement written in the form of an exponential function. The solutions from the proposed technique are verified with solutions from the literature and FEM solutions. The bucking loads correspond very well to other available solutions in most of the comparisons. The buckling modes also compare very well with the finite element solutions. The proposed solution technique transforms higher-order differential equations to first-order differential equations, and they are analytically solved for out-of-plane displacement in the form of an exponential function. Therefore, the proposed solution technique yields a solution which can be considered as an analytical solution.

The Reduction of Chazy Classes and Other Third-Order Differential Equations Related to Boundary Layer Flow Models

International Nuclear Information System (INIS)

Fakhar, K.; Kara, A. H.

2012-01-01

We study the symmetries, conservation laws and reduction of third-order equations that evolve from a prior reduction of models that arise in fluid phenomena. These could be the ordinary differential equations (ODEs) that are reductions of partial differential equations (PDEs) or, alternatively, PDEs related to given ODEs. In this class, the analysis includes the well-known Blasius, Chazy, and other associated third-order ODEs. (general)

A differential transformation approach for solving functional differential equations with multiple delays

Science.gov (United States)

Rebenda, Josef; Šmarda, Zdeněk

2017-07-01

In the paper an efficient semi-analytical approach based on the method of steps and the differential transformation is proposed for numerical approximation of solutions of functional differential models of delayed and neutral type on a finite interval of arbitrary length, including models with several constant delays. Algorithms for both commensurate and non-commensurate delays are described, applications are shown in examples. Validity and efficiency of the presented algorithms is compared with the variational iteration method, the Adomian decomposition method and the polynomial least squares method numerically. Matlab package DDE23 is used to produce reference numerical values.

Algebraic limit cycles in polynomial systems of differential equations

International Nuclear Information System (INIS)

Llibre, Jaume; Zhao Yulin

2007-01-01

Using elementary tools we construct cubic polynomial systems of differential equations with algebraic limit cycles of degrees 4, 5 and 6. We also construct a cubic polynomial system of differential equations having an algebraic homoclinic loop of degree 3. Moreover, we show that there are polynomial systems of differential equations of arbitrary degree that have algebraic limit cycles of degree 3, as well as give an example of a cubic polynomial system of differential equations with two algebraic limit cycles of degree 4

EXISTENCE OF POSITIVE SOLUTION TO TWO-POINT BOUNDARY VALUE PROBLEM FOR A SYSTEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper, we consider a two-point boundary value problem for a system of second order ordinary differential equations. Under some conditions, we show the existence of positive solution to the system of second order ordinary differential equa-tions.

Limit cycles via higher order perturbations for some piecewise differential systems

Science.gov (United States)

Buzzi, Claudio A.; Lima, Maurício Firmino Silva; Torregrosa, Joan

2018-05-01

A classical perturbation problem is the polynomial perturbation of the harmonic oscillator, (x‧ ,y‧) =(- y + εf(x , y , ε) , x + εg(x , y , ε)) . In this paper we study the limit cycles that bifurcate from the period annulus via piecewise polynomial perturbations in two zones separated by a straight line. We prove that, for polynomial perturbations of degree n , no more than Nn - 1 limit cycles appear up to a study of order N. We also show that this upper bound is reached for orders one and two. Moreover, we study this problem in some classes of piecewise Liénard differential systems providing better upper bounds for higher order perturbation in ε, showing also when they are reached. The Poincaré-Pontryagin-Melnikov theory is the main technique used to prove all the results.

EXISTENCE OF PERIODIC SOLUTION TO HIGHER ORDER DIFFERENTIAL EQUATIONS WITH DEVIATING ARGUMENT

Institute of Scientific and Technical Information of China (English)

无

2009-01-01

In this paper,using the coincidence degree theory of Mawhin,we investigate the existence of periodic solutions to higher order differential equations with deviating argument. Some new results on the existence of periodic solutions to the equations are obtained. In addition,we give an example to illustrate the main results.

Engineering arbitrary pure and mixed quantum states

International Nuclear Information System (INIS)

Pechen, Alexander

2011-01-01

Controlled manipulation by atomic- and molecular-scale quantum systems has attracted a lot of research attention in recent years. A fundamental problem is to provide deterministic methods for controlled engineering of arbitrary quantum states. This work proposes a deterministic method for engineering arbitrary pure and mixed states of a wide class of quantum systems. The method exploits a special combination of incoherent and coherent controls (incoherent and coherent radiation) and has two properties which are specifically important for manipulating by quantum systems: it realizes the strongest possible degree of their state control, complete density matrix controllability, meaning the ability to steer arbitrary pure and mixed initial states into any desired pure or mixed final state, and it is all-to-one, such that each particular control transfers all initial system states into one target state.

Analytic Expression of Arbitrary Matrix Elements for Boson Exponential Quadratic Polynomial Operators

Institute of Scientific and Technical Information of China (English)

XU Xiu-Wei; REN Ting-Qi; LIU Shu-Yan; MA Qiu-Ming; LIU Sheng-Dian

2007-01-01

Making use of the transformation relation among usual, normal, and antinormal ordering for the multimode boson exponential quadratic polynomial operators (BEQPO's), we present the analytic expression of arbitrary matrix elements for BEQPO's. As a preliminary application, we obtain the exact expressions of partition function about the boson quadratic polynomial system, matrix elements in particle-number, coordinate, and momentum representation, and P representation for the BEQPO's.

Oscillation Criteria in First Order Neutral Delay Impulsive Differential Equations with Constant Coefficients

International Nuclear Information System (INIS)

Dimitrova, M. B.; Donev, V. I.

2008-01-01

This paper is dealing with the oscillatory properties of first order neutral delay impulsive differential equations and corresponding to them inequalities with constant coefficients. The established sufficient conditions ensure the oscillation of every solution of this type of equations.

Automatic simplification of systems of reaction-diffusion equations by a posteriori analysis.

Science.gov (United States)

Maybank, Philip J; Whiteley, Jonathan P

2014-02-01

Many mathematical models in biology and physiology are represented by systems of nonlinear differential equations. In recent years these models have become increasingly complex in order to explain the enormous volume of data now available. A key role of modellers is to determine which components of the model have the greatest effect on a given observed behaviour. An approach for automatically fulfilling this role, based on a posteriori analysis, has recently been developed for nonlinear initial value ordinary differential equations [J.P. Whiteley, Model reduction using a posteriori analysis, Math. Biosci. 225 (2010) 44-52]. In this paper we extend this model reduction technique for application to both steady-state and time-dependent nonlinear reaction-diffusion systems. Exemplar problems drawn from biology are used to demonstrate the applicability of the technique. Copyright © 2014 Elsevier Inc. All rights reserved.

Some aspects of the inverse problem for general first order systems

International Nuclear Information System (INIS)

Sarlet, W.; Cantrijn, F.

1978-01-01

In this paper we investigate the most general systems of first-order ordinary differential equations which satisfy the integrability conditions of the inverse problem for canonical formulations, i.e., are derivable from a variational principle. It is shown that they have a structure which is invariant under arbitrary coordinate-transformations. A special class of transformations, called identity-isotopic transformations, is described. It is illustrated how these transformations share all properties with classical canonical transformations except one: the solutions of the system, viewed as transformations from the set of initial values, generally are not identity-isotopic. This failure is shown to be due to the explicit time-dependence of the functions which characterise the system, and the special role of this time-dependence is further clarified using the possibility of a reduction to conventional Hamiltonian systems. All properties are derived here within the framework of an analytic treatment of a system of differential equations

Differential and Difference Boundary Value Problem for Loaded Third-Order Pseudo-Parabolic Differential Equations and Difference Methods for Their Numerical Solution

Science.gov (United States)

Beshtokov, M. Kh.

2017-12-01

Boundary value problems for loaded third-order pseudo-parabolic equations with variable coefficients are considered. A priori estimates for the solutions of the problems in the differential and difference formulations are obtained. These a priori estimates imply the uniqueness and stability of the solution with respect to the initial data and the right-hand side on a layer, as well as the convergence of the solution of each difference problem to the solution of the corresponding differential problem.

The determination of the Dirac density matrix of the d-dimensional harmonic oscillator for an arbitrary number of closed shells

International Nuclear Information System (INIS)

Howard, I.A.; March, N.H.; Nieto, L.M.

2002-01-01

In 1959, March and Young (Nucl. Phys. 12 237) rewrote the equation of motion for the Dirac density matrix γ(x, x 0 ) in terms of sum and difference variables. Here, γ(r-bar, r-bar 0 ) for the d-dimensional isotropic harmonic oscillator for an arbitrary number of closed shells is shown to satisfy, using the variables vertical bar r-bar + r-bar 0 vertical bar/2 and vertical bar r-bar - r-bar 0 vertical bar/2, a generalized partial differential equation embracing the March-Young equation for d=1. As applications, we take in turn the cases d=1, 2, 3 and 4, and obtain both the density matrix γ (r-bar, r-bar 0 ) and the diagonal density ρ(r)=γ(r-bar, r-bar 0 ) vertical bar r-bar 0 =r-bar, this diagonal element already being known to satisfy a third-order linear homogeneous differential equation for d=1 through 3. Some comments are finally made on the d-dimensional kinetic energy density, which is important for first-principles density functional theory in allowing one to bypass one-particle Schroedinger equations (the so-called Slater-Kohn-Sham equations). (author)

Optimized difference schemes for multidimensional hyperbolic partial differential equations

Directory of Open Access Journals (Sweden)

Adrian Sescu

2009-04-01

Full Text Available In numerical solutions to hyperbolic partial differential equations in multidimensions, in addition to dispersion and dissipation errors, there is a grid-related error (referred to as isotropy error or numerical anisotropy that affects the directional dependence of the wave propagation. Difference schemes are mostly analyzed and optimized in one dimension, wherein the anisotropy correction may not be effective enough. In this work, optimized multidimensional difference schemes with arbitrary order of accuracy are designed to have improved isotropy compared to conventional schemes. The derivation is performed based on Taylor series expansion and Fourier analysis. The schemes are restricted to equally-spaced Cartesian grids, so the generalized curvilinear transformation method and Cartesian grid methods are good candidates.

Effective modeling and reverse-time migration for novel pure acoustic wave in arbitrary orthorhombic anisotropic media

Science.gov (United States)

Xu, Shigang; Liu, Yang

2018-03-01

The conventional pseudo-acoustic wave equations (PWEs) in arbitrary orthorhombic anisotropic (OA) media usually have coupled P- and SV-wave modes. These coupled equations may introduce strong SV-wave artifacts and numerical instabilities in P-wave simulation results and reverse-time migration (RTM) profiles. However, pure acoustic wave equations (PAWEs) completely decouple the P-wave component from the full elastic wavefield and naturally solve all the aforementioned problems. In this article, we present a novel PAWE in arbitrary OA media and compare it with the conventional coupled PWEs. Through decomposing the solution of the corresponding eigenvalue equation for the original PWE into an ellipsoidal differential operator (EDO) and an ellipsoidal scalar operator (ESO), the new PAWE in time-space domain is constructed by applying the combination of these two solvable operators and can effectively describe P-wave features in arbitrary OA media. Furthermore, we adopt the optimal finite-difference method (FDM) to solve the newly derived PAWE. In addition, the three-dimensional (3D) hybrid absorbing boundary condition (HABC) with some reasonable modifications is developed for reducing artificial edge reflections in anisotropic media. To improve computational efficiency in 3D case, we adopt graphic processing unit (GPU) with Compute Unified Device Architecture (CUDA) instead of traditional central processing unit (CPU) architecture. Several numerical experiments for arbitrary OA models confirm that the proposed schemes can produce pure, stable and accurate P-wave modeling results and RTM images with higher computational efficiency. Moreover, the 3D numerical simulations can provide us with a comprehensive and real description of wave propagation.

Construction of a Smooth Lyapunov Function for the Robust and Exact Second-Order Differentiator

Directory of Open Access Journals (Sweden)

Tonametl Sanchez

2016-01-01

Full Text Available Differentiators play an important role in (continuous feedback control systems. In particular, the robust and exact second-order differentiator has shown some very interesting properties and it has been used successfully in sliding mode control, in spite of the lack of a Lyapunov based procedure to design its gains. As contribution of this paper, we provide a constructive method to determine a differentiable Lyapunov function for such a differentiator. Moreover, the Lyapunov function is used to provide a procedure to design the differentiator’s parameters. Also, some sets of such parameters are provided. The determination of the positive definiteness of the Lyapunov function and negative definiteness of its derivative is converted to the problem of solving a system of inequalities linear in the parameters of the Lyapunov function candidate and also linear in the gains of the differentiator, but bilinear in both.

Third-order differential ladder operators and supersymmetric quantum mechanics

International Nuclear Information System (INIS)

Mateo, J; Negro, J

2008-01-01

Hierarchies of one-dimensional Hamiltonians in quantum mechanics admitting third-order differential ladder operators are studied. Each Hamiltonian has associated three-step Darboux (pseudo)-cycles and Painleve IV equations as a closure condition. The whole hierarchy is generated applying some operations on the cycles. These operations are investigated in the frame of supersymmetric quantum mechanics and mainly involve algebraic manipulations. A consistent geometric representation for the hierarchy and cycles is built that also helps in understanding the operations. Three kinds of hierarchies are distinguished and a realization based on the harmonic oscillator Hamiltonian is supplied, giving an interpretation for the spectral properties of the Hamiltonians of each hierarchy

Numerical solutions of multi-order fractional differential equations by Boubaker polynomials

Directory of Open Access Journals (Sweden)

Bolandtalat A.

2016-01-01

Full Text Available In this paper, we have applied a numerical method based on Boubaker polynomials to obtain approximate numerical solutions of multi-order fractional differential equations. We obtain an operational matrix of fractional integration based on Boubaker polynomials. Using this operational matrix, the given problem is converted into a set of algebraic equations. Illustrative examples are are given to demonstrate the efficiency and simplicity of this technique.

N-point g-loop vertex for a free fermionic theory with arbitrary spin

International Nuclear Information System (INIS)

Di Vecchia, P.; Pezzella, F.; Frau, M.; Hornfeck, K.

1990-01-01

We use the sewing procedure of the operator formalism to construct explicitly the N-point g-loop vertex V N;g for a free fermionic (b, c)-system with conformal weight (λ, 1-λ). We show that this vertex has the structure we expect from geometrical arguments. We obtain also several geometrical objects, e.g. the holomorphic λ-differentials on an arbitrary Riemann surface, which turn out to be expressed as a Poincare θ-series over all elements of the Schottky group. From V N;g we compute explicitly correlation functions for our system, finding agreement with the geometrical procedure. (orig.)

N-point g-loop vertex for a free fermionic theory with arbitrary spin

International Nuclear Information System (INIS)

Di Vecchia, P.; Pezzella, F.; Frau, M.; Hornfeck, K.

1989-07-01

We use the sewing procedure of the operator fomalism to construct explicitly the N-Point g-Loop Vertex V N;g for a free fermionic (b, c)-system with conformal weight (λ, 1-λ). We show that this Vertex has the structure we expect from geometrical arguments. We obtain also several geometrical objects, e.g. the holomorphic λ differentials on an arbitrary Riemann surface, which turn out to be expressed as a Poincare θ series over all elements of the Schottky group. From V N;g we compute explicitly correlation functions for our system, finding agreement with the geometrical procedure. (orig.)

Automatic Program Development

DEFF Research Database (Denmark)

Automatic Program Development is a tribute to Robert Paige (1947-1999), our accomplished and respected colleague, and moreover our good friend, whose untimely passing was a loss to our academic and research community. We have collected the revised, updated versions of the papers published in his...... honor in the Higher-Order and Symbolic Computation Journal in the years 2003 and 2005. Among them there are two papers by Bob: (i) a retrospective view of his research lines, and (ii) a proposal for future studies in the area of the automatic program derivation. The book also includes some papers...... by members of the IFIP Working Group 2.1 of which Bob was an active member. All papers are related to some of the research interests of Bob and, in particular, to the transformational development of programs and their algorithmic derivation from formal specifications. Automatic Program Development offers...

On a higher order multi-term time-fractional partial differential equation involving Caputo-Fabrizio derivative

OpenAIRE

Pirnapasov, Sardor; Karimov, Erkinjon

2017-01-01

In the present work we discuss higher order multi-term partial differential equation (PDE) with the Caputo-Fabrizio fractional derivative in time. We investigate a boundary value problem for fractional heat equation involving higher order Caputo-Fabrizio derivatives in time-variable. Using method of separation of variables and integration by parts, we reduce fractional order PDE to the integer order. We represent explicit solution of formulated problem in particular case by Fourier series.

Periodic solutions of singular second order differential equations : upper and lower functions

Czech Academy of Sciences Publication Activity Database

Hakl, Robert; Torres, P.J.; Zamora, M.

2011-01-01

Roč. 74, č. 18 (2011), s. 7078-7093 ISSN 0362-546X Institutional research plan: CEZ:AV0Z10190503 Keywords : second-order differential equation * singularity at the phase variable * Rayleigh-Plesset equation Subject RIV: BA - General Mathematics Impact factor: 1.536, year: 2011 http://www.sciencedirect.com/science/article/pii/S0362546X11005049

Interval oscillation criteria for second-order forced impulsive delay differential equations with damping term.

Science.gov (United States)

Thandapani, Ethiraju; Kannan, Manju; Pinelas, Sandra

2016-01-01

In this paper, we present some sufficient conditions for the oscillation of all solutions of a second order forced impulsive delay differential equation with damping term. Three factors-impulse, delay and damping that affect the interval qualitative properties of solutions of equations are taken into account together. The results obtained in this paper extend and generalize some of the the known results for forced impulsive differential equations. An example is provided to illustrate the main result.

Virtual reality based adaptive dose assessment method for arbitrary geometries in nuclear facility decommissioning.

Science.gov (United States)

Liu, Yong-Kuo; Chao, Nan; Xia, Hong; Peng, Min-Jun; Ayodeji, Abiodun

2018-05-17

This paper presents an improved and efficient virtual reality-based adaptive dose assessment method (VRBAM) applicable to the cutting and dismantling tasks in nuclear facility decommissioning. The method combines the modeling strength of virtual reality with the flexibility of adaptive technology. The initial geometry is designed with the three-dimensional computer-aided design tools, and a hybrid model composed of cuboids and a point-cloud is generated automatically according to the virtual model of the object. In order to improve the efficiency of dose calculation while retaining accuracy, the hybrid model is converted to a weighted point-cloud model, and the point kernels are generated by adaptively simplifying the weighted point-cloud model according to the detector position, an approach that is suitable for arbitrary geometries. The dose rates are calculated with the Point-Kernel method. To account for radiation scattering effects, buildup factors are calculated with the Geometric-Progression formula in the fitting function. The geometric modeling capability of VRBAM was verified by simulating basic geometries, which included a convex surface, a concave surface, a flat surface and their combination. The simulation results show that the VRBAM is more flexible and superior to other approaches in modeling complex geometries. In this paper, the computation time and dose rate results obtained from the proposed method were also compared with those obtained using the MCNP code and an earlier virtual reality-based method (VRBM) developed by the same authors. © 2018 IOP Publishing Ltd.

Estimation of periodic solutions number of first-order differential equations

Science.gov (United States)

Ivanov, Gennady; Alferov, Gennady; Gorovenko, Polina; Sharlay, Artem

2018-05-01

The paper deals with first-order differential equations under the assumption that the right-hand side is a periodic function of time and continuous in the set of arguments. Pliss V.A. obtained the first results for a particular class of equations and showed that a number of theorems can not be continued. In this paper, it was possible to reduce the restrictions on the degree of smoothness of the right-hand side of the equation and obtain upper and lower bounds on the number of possible periodic solutions.

Transient potentials in dendritic systems of arbitrary geometry.

Science.gov (United States)

Butz, E G; Cowan, J D

1974-09-01

A simple graphical calculus is developed that generates analytic solutions for membrane potential transforms at any point on the dendritic tree of neurons with arbitrary dendritic geometries, in response to synaptic "current" inputs. Such solutions permit the computation of transients in neurons with arbitrary geometry and may facilitate analysis of the role of dendrites in such cells.

The geometric approach to sets of ordinary differential equations and Hamiltonian dynamics

Science.gov (United States)

Estabrook, F. B.; Wahlquist, H. D.

1975-01-01

The calculus of differential forms is used to discuss the local integration theory of a general set of autonomous first order ordinary differential equations. Geometrically, such a set is a vector field V in the space of dependent variables. Integration consists of seeking associated geometric structures invariant along V: scalar fields, forms, vectors, and integrals over subspaces. It is shown that to any field V can be associated a Hamiltonian structure of forms if, when dealing with an odd number of dependent variables, an arbitrary equation of constraint is also added. Families of integral invariants are an immediate consequence. Poisson brackets are isomorphic to Lie products of associated CT-generating vector fields. Hamilton's variational principle follows from the fact that the maximal regular integral manifolds of a closed set of forms must include the characteristics of the set.

Maximum principles for boundary-degenerate second-order linear elliptic differential operators

OpenAIRE

Feehan, Paul M. N.

2012-01-01

We prove weak and strong maximum principles, including a Hopf lemma, for smooth subsolutions to equations defined by linear, second-order, partial differential operators whose principal symbols vanish along a portion of the domain boundary. The boundary regularity property of the smooth subsolutions along this boundary vanishing locus ensures that these maximum principles hold irrespective of the sign of the Fichera function. Boundary conditions need only be prescribed on the complement in th...

An implicit measure of associations with mental illness versus physical illness: response latency decomposition and stimuli differential functioning in relation to IAT order of associative conditions and accuracy.

Science.gov (United States)

Mannarini, Stefania; Boffo, Marilisa

2014-01-01

The present study aimed at the definition of a latent measurement dimension underlying an implicit measure of automatic associations between the concept of mental illness and the psychosocial and biogenetic causal explanatory attributes. To this end, an Implicit Association Test (IAT) assessing the association between the Mental Illness and Physical Illness target categories to the Psychological and Biologic attribute categories, representative of the causal explanation domains, was developed. The IAT presented 22 stimuli (words and pictures) to be categorized into the four categories. After 360 university students completed the IAT, a Many-Facet Rasch Measurement (MFRM) modelling approach was applied. The model specified a person latency parameter and a stimulus latency parameter. Two additional parameters were introduced to denote the order of presentation of the task associative conditions and the general response accuracy. Beyond the overall definition of the latent measurement dimension, the MFRM was also applied to disentangle the effect of the task block order and the general response accuracy on the stimuli response latency. Further, the MFRM allowed detecting any differential functioning of each stimulus in relation to both block ordering and accuracy. The results evidenced: a) the existence of a latency measurement dimension underlying the Mental Illness versus Physical Illness - Implicit Association Test; b) significant effects of block order and accuracy on the overall latency; c) a differential functioning of specific stimuli. The results of the present study can contribute to a better understanding of the functioning of an implicit measure of semantic associations with mental illness and give a first blueprint for the examination of relevant issues in the development of an IAT.

Do Automatic Self-Associations Relate to Suicidal Ideation?

NARCIS (Netherlands)

Glashouwer, Klaske A.; de Jong, Peter J.; Penninx, Brenda W. J. H.; Kerkhof, Ad J. F. M.; van Dyck, Richard; Ormel, Johan

Dysfunctional self-schemas are assumed to play an important role in suicidal ideation. According to recent information-processing models, it is important to differentiate between 'explicit' beliefs and automatic associations. Explicit beliefs stem from the weighting of propositions and their

Integrable dissipative nonlinear second order differential equations via factorizations and Abel equations

Energy Technology Data Exchange (ETDEWEB)

Mancas, Stefan C. [Department of Mathematics, Embry–Riddle Aeronautical University, Daytona Beach, FL 32114-3900 (United States); Rosu, Haret C., E-mail: hcr@ipicyt.edu.mx [IPICYT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Apdo Postal 3-74 Tangamanga, 78231 San Luis Potosí, SLP (Mexico)

2013-09-02

We emphasize two connections, one well known and another less known, between the dissipative nonlinear second order differential equations and the Abel equations which in their first-kind form have only cubic and quadratic terms. Then, employing an old integrability criterion due to Chiellini, we introduce the corresponding integrable dissipative equations. For illustration, we present the cases of some integrable dissipative Fisher, nonlinear pendulum, and Burgers–Huxley type equations which are obtained in this way and can be of interest in applications. We also show how to obtain Abel solutions directly from the factorization of second order nonlinear equations.

AIDEN: A Density Conscious Artificial Immune System for Automatic Discovery of Arbitrary Shape Clusters in Spatial Patterns

Directory of Open Access Journals (Sweden)

Vishwambhar Pathak

2012-11-01

Full Text Available Recent efforts in modeling of dynamics of the natural immune cells leading to artificial immune systems (AIS have ignited contemporary research interest in finding out its analogies to real world problems. The AIS models have been vastly exploited to develop dependable robust
solutions to clustering. Most of the traditional clustering methods bear limitations in their capability to detect clusters of arbitrary shapes in a fully unsupervised manner. In this paper the recognition and communication dynamics of T Cell Receptors, the recognizing elements in innate immune
system, has been modeled with a kernel density estimation method. The model has been shown to successfully discover non spherical clusters in spatial patterns. Modeling the cohesion of the antibodies and pathogens with ‘local influence’ measure inducts comprehensive extension of the
antibody representation ball (ARB, which in turn corresponds to controlled expansion of clusters and prevents overfitting.

A generic approach for the automatic verification of featured, parameterised systems

OpenAIRE

Miller, A.; Calder, M.

2005-01-01

A general technique is presented that allows property based feature analysis of systems consisting of an arbitrary number of components. Each component may have an arbitrary set of safe features. The components are defined in a guarded command form and the technique combines model checking and abstraction. Features must fulfill certain criteria in order to be safe, the criteria express constraints on the variables which occur in feature guards. The main result is a generalisation theorem whic...

Path integral solution of linear second order partial differential equations I: the general construction

International Nuclear Information System (INIS)

LaChapelle, J.

2004-01-01

A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schroedinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette

Automatic yield-line analysis of slabs using discontinuity layout optimization.

Science.gov (United States)

Gilbert, Matthew; He, Linwei; Smith, Colin C; Le, Canh V

2014-08-08

The yield-line method of analysis is a long established and extremely effective means of estimating the maximum load sustainable by a slab or plate. However, although numerous attempts to automate the process of directly identifying the critical pattern of yield-lines have been made over the past few decades, to date none has proved capable of reliably analysing slabs of arbitrary geometry. Here, it is demonstrated that the discontinuity layout optimization (DLO) procedure can successfully be applied to such problems. The procedure involves discretization of the problem using nodes inter-connected by potential yield-line discontinuities, with the critical layout of these then identified using linear programming. The procedure is applied to various benchmark problems, demonstrating that highly accurate solutions can be obtained, and showing that DLO provides a truly systematic means of directly and reliably automatically identifying yield-line patterns. Finally, since the critical yield-line patterns for many problems are found to be quite complex in form, a means of automatically simplifying these is presented.

Novel Numerical Methods for Optimal Control Problems Involving Fractional-Order Differential Equations

Science.gov (United States)

2018-03-14

UNIVERSITY OF TECHNOLOGY Final Report 03/14/2018 DISTRIBUTION A: Distribution approved for public release. AF Office Of Scientific Research (AFOSR...optimal control problems involving fractional-order differential equations Wang, Song Curtin University of Technology Kent Street, Bentley WA6102...Article history : Received 3 October 2016 Accepted 26 March 2017 Available online 29 April 2017 Keywords: Hamilton–Jacobi–Bellman equation Financial

An Arbitrary Benchmark CAPM: One Additional Frontier Portfolio is Sufficient

OpenAIRE

Ekern, Steinar

2008-01-01

First draft: July 16, 2008 This version: October 7, 2008 The benchmark CAPM linearly relates the expected returns on an arbitrary asset, an arbitrary benchmark portfolio, and an arbitrary MV frontier portfolio. The benchmark is not required to be on the frontier and may be non-perfectly correlated with the frontier portfolio. The benchmark CAPM extends and generalizes previous CAPM formulations, including the zero beta, two correlated frontier portfolios, riskless augmented frontier, an...

m-POINT BOUNDARY VALUE PROBLEM FOR SECOND ORDER IMPULSIVE DIFFERENTIAL EQUATION AT RESONANCE

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In his paper,we obtain a general theorem concerning the existence of solutions to an m-point boundary value problem for the second-order differential equation with impulses.Moreover,the result can also be applied to study the usual m-point boundary value problem at resonance without impulses.

SELFADJUSTING AUTOMATIC CONTROL OF SOWING UNIT

Directory of Open Access Journals (Sweden)

A. Yu. Izmaylov

2015-01-01

Full Text Available The selfadjusting automatic control of sowing unit and differentiated introduction of mineral fertilizers doses according to agrochemical indicators of the soil (precision agriculture are used wider nowadays. It was defined that the main requirement to the differentiated seeding and fertilizing is an accuracy and duration of transition from one norm to another. Established that at a speed of unit of 10 km/h object moves for 0.5 s about on 1.5 m and more. Thus in this device the radio channel originated differentiated correction is updated in 10 s, and in the RTK mode - 0.5-2 s that breaks the accuracy of introduction of seeds and fertilizers. The block schematic diagram of system of automatic control of technological process of seeding and mineral fertilizing with use of navigation means of machine-tractor aggregates orientation in the field and technical means for realization of technology of precision agriculture at sowing and fertilizers application due to electronic maps of soil fertility and navigation satellite systems was worked out. It was noted that for regulation of a fertilizing dose it is necessary to complete the unit with the electric drive, and for error reduction use navigation GLONASS, GPS, Galileo receivers. To tracking of four leading navigation systems GPS/GLONASS/Galileo/Compass receiver with 32 canals developed by domestic-owned firm «KB NAVIS» was suggested. It was established that the automated device created by All-Russia Research Institute of Mechanization for Agriculture information based on NAVSTAR and GLONASS/GPS system successfully operates seeding and make possible the differentiate fertilizing.

The solutions of second-order linear differential systems with constant delays

Science.gov (United States)

Diblík, Josef; Svoboda, Zdeněk

2017-07-01

The representations of solutions to initial problems for non-homogenous n-dimensional second-order differential equations with delays x″(t )-2 A x'(t -τ )+(A2+B2)x (t -2 τ )=f (t ) by means of special matrix delayed functions are derived. Square matrices A and B are commuting and τ > 0. Derived representations use what is called a delayed exponential of a matrix and results generalize some of known results previously derived for homogenous systems.

Conditioned craving cues elicit an automatic approach tendency

NARCIS (Netherlands)

van Gucht, D.; Vansteenwegen, D.; Van den Bergh, O.; Beckers, T.

2008-01-01

In two experiments, we used a Pavlovian differential conditioning procedure to induce craving for chocolate. As a result of repeated pairing with chocolate intake, initially neutral cues came to elicit an automatic approach tendency in a speeded stimulus-response compatibility reaction time task.

Multivariate η-μ fading distribution with arbitrary correlation model

Science.gov (United States)

Ghareeb, Ibrahim; Atiani, Amani

2018-03-01

An extensive analysis for the multivariate ? distribution with arbitrary correlation is presented, where novel analytical expressions for the multivariate probability density function, cumulative distribution function and moment generating function (MGF) of arbitrarily correlated and not necessarily identically distributed ? power random variables are derived. Also, this paper provides exact-form expression for the MGF of the instantaneous signal-to-noise ratio at the combiner output in a diversity reception system with maximal-ratio combining and post-detection equal-gain combining operating in slow frequency nonselective arbitrarily correlated not necessarily identically distributed ?-fading channels. The average bit error probability of differentially detected quadrature phase shift keying signals with post-detection diversity reception system over arbitrarily correlated and not necessarily identical fading parameters ?-fading channels is determined by using the MGF-based approach. The effect of fading correlation between diversity branches, fading severity parameters and diversity level is studied.

Contact symmetries of general linear second-order ordinary differential equations: letter to the editor

NARCIS (Netherlands)

Martini, Ruud; Kersten, P.H.M.

1983-01-01

Using 1-1 mappings, the complete symmetry groups of contact transformations of general linear second-order ordinary differential equations are determined from two independent solutions of those equations, and applied to the harmonic oscillator with and without damping.

The fractional calculus theory and applications of differentiation and integration to arbitrary order

CERN Document Server

Oldham, Keith B

1974-01-01

In this book, we study theoretical and practical aspects of computing methods for mathematical modelling of nonlinear systems. A number of computing techniques are considered, such as methods of operator approximation with any given accuracy; operator interpolation techniques including a non-Lagrange interpolation; methods of system representation subject to constraints associated with concepts of causality, memory and stationarity; methods of system representation with an accuracy that is the best within a given class of models; methods of covariance matrix estimation;methods for low-rank mat

Loss of nonclassical properties of quantum states in linear phase-insensitive processes with arbitrary time-dependent parameters

International Nuclear Information System (INIS)

Dodonov, V.V.

2009-01-01

Conditions of disappearance of different 'nonclassical' properties (usual and high-order squeezing, sub-Poissonian statistics, negativity of s-parametrized quasidistributions) are derived for a quantum oscillator, whose evolution is governed by the standard master equation of quantum optics with arbitrary time-dependent coefficients.

Sturm-Picone type theorems for second-order nonlinear differential equations

Directory of Open Access Journals (Sweden)

Aydin Tiryaki

2014-06-01

Full Text Available The aim of this article is to give Sturm-Picone type theorems for the pair of second-order nonlinear differential equations $$\\displaylines{ (p_1(t|x'|^{\\alpha-1}x''+q_1(tf_1(x=0 \\cr (p_2(t|y'|^{\\alpha-1}y''+q_2(tf_2(y=0,\\quad t_1

ON THE BOUNDEDNESS AND THE STABILITY OF SOLUTION TO THIRD ORDER NON-LINEAR DIFFERENTIAL EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2008-01-01

In this paper we investigate the global asymptotic stability,boundedness as well as the ultimate boundedness of solutions to a general third order nonlinear differential equation,using complete Lyapunov function.

The nonlocal boundary value problems for strongly singular higher-order nonlinear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan

-, č. 35 (2015), s. 23-50 ISSN 1126-8042 Institutional support: RVO:67985840 Keywords : higher order functional differential equations * Dirichlet boundary value problem * strong singularity Subject RIV: BA - General Mathematics http://ijpam.uniud.it/online_issue/201535/03-Mukhigulashvili.pdf

Geometric phases in astigmatic optical modes of arbitrary order

International Nuclear Information System (INIS)

Habraken, Steven J. M.; Nienhuis, Gerard

2010-01-01

The transverse spatial structure of a paraxial beam of light is fully characterized by a set of parameters that vary only slowly under free propagation. They specify bosonic ladder operators that connect modes of different orders, in analogy to the ladder operators connecting harmonic-oscillator wave functions. The parameter spaces underlying sets of higher-order modes are isomorphic to the parameter space of the ladder operators. We study the geometry of this space and the geometric phase that arises from it. This phase constitutes the ultimate generalization of the Gouy phase in paraxial wave optics. It reduces to the ordinary Gouy phase and the geometric phase of nonastigmatic optical modes with orbital angular momentum in limiting cases. We briefly discuss the well-known analogy between geometric phases and the Aharonov-Bohm effect, which provides some complementary insights into the geometric nature and origin of the generalized Gouy phase shift. Our method also applies to the quantum-mechanical description of wave packets. It allows for obtaining complete sets of normalized solutions of the Schroedinger equation. Cyclic transformations of such wave packets give rise to a phase shift, which has a geometric interpretation in terms of the other degrees of freedom involved.

Higher-Order Integral Equation Methods in Computational Electromagnetics

DEFF Research Database (Denmark)

Jørgensen, Erik; Meincke, Peter

Higher-order integral equation methods have been investigated. The study has focused on improving the accuracy and efficiency of the Method of Moments (MoM) applied to electromagnetic problems. A new set of hierarchical Legendre basis functions of arbitrary order is developed. The new basis...

Fully automatic adjoints: a robust and efficient mechanism for generating adjoint ocean models

Science.gov (United States)

Ham, D. A.; Farrell, P. E.; Funke, S. W.; Rognes, M. E.

2012-04-01

The problem of generating and maintaining adjoint models is sufficiently difficult that typically only the most advanced and well-resourced community ocean models achieve it. There are two current technologies which each suffer from their own limitations. Algorithmic differentiation, also called automatic differentiation, is employed by models such as the MITGCM [2] and the Alfred Wegener Institute model FESOM [3]. This technique is very difficult to apply to existing code, and requires a major initial investment to prepare the code for automatic adjoint generation. AD tools may also have difficulty with code employing modern software constructs such as derived data types. An alternative is to formulate the adjoint differential equation and to discretise this separately. This approach, known as the continuous adjoint and employed in ROMS [4], has the disadvantage that two different model code bases must be maintained and manually kept synchronised as the model develops. The discretisation of the continuous adjoint is not automatically consistent with that of the forward model, producing an additional source of error. The alternative presented here is to formulate the flow model in the high level language UFL (Unified Form Language) and to automatically generate the model using the software of the FEniCS project. In this approach it is the high level code specification which is differentiated, a task very similar to the formulation of the continuous adjoint [5]. However since the forward and adjoint models are generated automatically, the difficulty of maintaining them vanishes and the software engineering process is therefore robust. The scheduling and execution of the adjoint model, including the application of an appropriate checkpointing strategy is managed by libadjoint [1]. In contrast to the conventional algorithmic differentiation description of a model as a series of primitive mathematical operations, libadjoint employs a new abstraction of the simulation

Conformal array design on arbitrary polygon surface with transformation optics

Energy Technology Data Exchange (ETDEWEB)

Deng, Li, E-mail: dengl@bupt.edu.cn; Hong, Weijun, E-mail: hongwj@bupt.edu.cn; Zhu, Jianfeng; Peng, Biao; Li, Shufang [Beijing Key Laboratory of Network System Architecture and Convergence, School of Information and Communication Engineering, Beijing University of Posts and Telecommunications, 100876 Beijing (China); Wu, Yongle, E-mail: wuyongle138@gmail.com [Beijing Key Laboratory of Work Safety Intelligent Monitoring, School of Electronic Engineering, Beijing University of Posts and Telecommunications, 100876 Beijing (China)

2016-06-15

A transformation-optics based method to design a conformal antenna array on an arbitrary polygon surface is proposed and demonstrated in this paper. This conformal antenna array can be adjusted to behave equivalently as a uniformly spaced linear array by applying an appropriate transformation medium. An typical example of general arbitrary polygon conformal arrays, not limited to circular array, is presented, verifying the proposed approach. In summary, the novel arbitrary polygon surface conformal array can be utilized in array synthesis and beam-forming, maintaining all benefits of linear array.

Conformal array design on arbitrary polygon surface with transformation optics

International Nuclear Information System (INIS)

Deng, Li; Hong, Weijun; Zhu, Jianfeng; Peng, Biao; Li, Shufang; Wu, Yongle

2016-01-01

A transformation-optics based method to design a conformal antenna array on an arbitrary polygon surface is proposed and demonstrated in this paper. This conformal antenna array can be adjusted to behave equivalently as a uniformly spaced linear array by applying an appropriate transformation medium. An typical example of general arbitrary polygon conformal arrays, not limited to circular array, is presented, verifying the proposed approach. In summary, the novel arbitrary polygon surface conformal array can be utilized in array synthesis and beam-forming, maintaining all benefits of linear array.

A Quartic Conformally Covariant Differential Operator for Arbitrary Pseudo-Riemannian Manifolds (Summary

Directory of Open Access Journals (Sweden)

Stephen M. Paneitz

2008-03-01

Full Text Available This is the original manuscript dated March 9th 1983, typeset by the Editors for the Proceedings of the Midwest Geometry Conference 2007 held in memory of Thomas Branson. Stephen Paneitz passed away on September 1st 1983 while attending a conference in Clausthal and the manuscript was never published. For more than 20 years these few pages were circulated informally. In November 2004, as a service to the mathematical community, Tom Branson added a scan of the manuscript to his website. Here we make it available more formally. It is surely one of the most cited unpublished articles. The differential operator defined in this article plays a key rôle in conformal differential geometry in dimension 4 and is now known as the Paneitz operator.

Predicting top-of-atmosphere radiance for arbitrary viewing geometries from the visible to thermal infrared: generalization to arbitrary average scene temperatures

Science.gov (United States)

Florio, Christopher J.; Cota, Steve A.; Gaffney, Stephanie K.

2010-08-01

In a companion paper presented at this conference we described how The Aerospace Corporation's Parameterized Image Chain Analysis & Simulation SOftware (PICASSO) may be used in conjunction with a limited number of runs of AFRL's MODTRAN4 radiative transfer code, to quickly predict the top-of-atmosphere (TOA) radiance received in the visible through midwave IR (MWIR) by an earth viewing sensor, for any arbitrary combination of solar and sensor elevation angles. The method is particularly useful for large-scale scene simulations where each pixel could have a unique value of reflectance/emissivity and temperature, making the run-time required for direct prediction via MODTRAN4 prohibitive. In order to be self-consistent, the method described requires an atmospheric model (defined, at a minimum, as a set of vertical temperature, pressure and water vapor profiles) that is consistent with the average scene temperature. MODTRAN4 provides only six model atmospheres, ranging from sub-arctic winter to tropical conditions - too few to cover with sufficient temperature resolution the full range of average scene temperatures that might be of interest. Model atmospheres consistent with intermediate temperature values can be difficult to come by, and in any event, their use would be too cumbersome for use in trade studies involving a large number of average scene temperatures. In this paper we describe and assess a method for predicting TOA radiance for any arbitrary average scene temperature, starting from only a limited number of model atmospheres.

On method of solving third-order ordinary differential equations directly using Bernstein polynomials

Science.gov (United States)

Khataybeh, S. N.; Hashim, I.

2018-04-01

In this paper, we propose for the first time a method based on Bernstein polynomials for solving directly a class of third-order ordinary differential equations (ODEs). This method gives a numerical solution by converting the equation into a system of algebraic equations which is solved directly. Some numerical examples are given to show the applicability of the method.

Differential geometry on Hopf algebras and quantum groups

International Nuclear Information System (INIS)

Watts, P.

1994-01-01

The differential geometry on a Hopf algebra is constructed, by using the basic axioms of Hopf algebras and noncommutative differential geometry. The space of generalized derivations on a Hopf algebra of functions is presented via the smash product, and used to define and discuss quantum Lie algebras and their properties. The Cartan calculus of the exterior derivative, Lie derivative, and inner derivation is found for both the universal and general differential calculi of an arbitrary Hopf algebra, and, by restricting to the quasitriangular case and using the numerical R-matrix formalism, the aforementioned structures for quantum groups are determined

Functional contributions and interactions between the human hippocampus and subregions of the striatum during arbitrary associative learning and memory

Science.gov (United States)

Mattfeld, Aaron T.; Stark, Craig E. L.

2015-01-01

The hippocampus and striatum are thought to have different functional roles in learning and memory. It is unknown under what experimental conditions their contributions are dissimilar or converge, and the extent to which they interact over the course of learning. In order to evaluate both the functional contributions of as well as the interactions between the human hippocampus and striatum, the present study used high-resolution functional magnetic resonance imaging (fMRI) and variations of a conditional visuomotor associative learning task that either taxed arbitrary associative learning (Experiment 1) or stimulus-response learning (Experiment 2). In the first experiment we observed changes in activity in the hippocampus and anterior caudate that reflect differences between the two regions consistent with distinct computational principles. In the second experiment we observed activity in the putamen that reflected content specific representations during the learning of arbitrary conditional visuomotor associations. In both experiments the hippocampus and ventral striatum demonstrated dynamic functional coupling during the learning of new arbitrary associations, but not during retrieval of well-learned arbitrary associations using control variants of the tasks that did not preferentially tax one system versus the other. These findings suggest that both the hippocampus and subregions of the dorsal striatum contribute uniquely to the learning of arbitrary associations while the hippocampus and ventral striatum interact over the course of learning. PMID:25560298

Controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with delay and Poisson jumps

Directory of Open Access Journals (Sweden)

Diem Dang Huan

2015-12-01

Full Text Available The current paper is concerned with the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps in Hilbert spaces. Using the theory of a strongly continuous cosine family of bounded linear operators, stochastic analysis theory and with the help of the Banach fixed point theorem, we derive a new set of sufficient conditions for the controllability of nonlocal second-order impulsive neutral stochastic functional integro-differential equations with infinite delay and Poisson jumps. Finally, an application to the stochastic nonlinear wave equation with infinite delay and Poisson jumps is given.

A general comparison theorem for backward stochastic differential equations

OpenAIRE

Cohen, Samuel N.; Elliott, Robert J.; Pearce, Charles E. M.

2010-01-01

A useful result when dealing with backward stochastic differential equations is the comparison theorem of Peng (1992). When the equations are not based on Brownian motion, the comparison theorem no longer holds in general. In this paper we present a condition for a comparison theorem to hold for backward stochastic differential equations based on arbitrary martingales. This theorem applies to both vector and scalar situations. Applications to the theory of nonlinear expectat...

Differential Higgs boson pair production at next-to-next-to-leading order in QCD

International Nuclear Information System (INIS)

Florian, Daniel de; Mazzitelli, Javier; Grazzini, Massimiliano; Hanga, Catalin; Lindert, Jonas M.; Kallweit, Stefan; Maierhoefer, Philipp; Rathlev, Dirk

2016-06-01

We report on the first fully differential calculation for double Higgs boson production through gluon fusion in hadron collisions up to next-to-next-to-leading order (NNLO) in QCD perturbation theory. The calculation is performed in the heavy-top limit of the Standard Model, and in the phenomenological results we focus on pp collisions at √(s)=14 TeV. We present differential distributions through NNLO for various observables including the transverse-momentum and rapidity distributions of the two Higgs bosons. NNLO corrections are at the level of 10%-25% with respect to the next-to-leading order (NLO) prediction with a residual scale uncertainty of 5%-15% and an overall mild phase-space dependence. Only at NNLO the perturbative expansion starts to converge yielding overlapping scale uncertainty bands between NNLO and NLO in most of the phase-space. The calculation includes NLO predictions for pp→HH+jet+X. Corrections to the corresponding distributions exceed 50% with a residual scale dependence of 20%-30%.

Rethinking pedagogy for second-order differential equations: a simplified approach to understanding well-posed problems

Science.gov (United States)

Tisdell, Christopher C.

2017-07-01

Knowing an equation has a unique solution is important from both a modelling and theoretical point of view. For over 70 years, the approach to learning and teaching 'well posedness' of initial value problems (IVPs) for second- and higher-order ordinary differential equations has involved transforming the problem and its analysis to a first-order system of equations. We show that this excursion is unnecessary and present a direct approach regarding second- and higher-order problems that does not require an understanding of systems.

On the mild solutions of higher-order differential equations in Banach spaces

Directory of Open Access Journals (Sweden)

Nguyen Thanh Lan

2003-01-01

Full Text Available For the higher-order abstract differential equation u(n(t=Au(t+f(t, t∈ℝ, we give a new definition of mild solutions. We then characterize the regular admissibility of a translation-invariant subspace ℳ of BUC(ℝ,E with respect to the above-mentioned equation in terms of solvability of the operator equation AX−Xn=C. As applications, periodicity and almost periodicity of mild solutions are also proved.

Symmetric multiparty-controlled teleportation of an arbitrary two-particle entanglement

International Nuclear Information System (INIS)

Deng Fuguo; Zhou Hongyu; Li Chunyan; Wang Yan; Li Yansong

2005-01-01

We present a way for symmetric multiparty-controlled teleportation of an arbitrary two-particle entangled state based on Bell-basis measurements by using two Greenberger-Horne-Zeilinger states, i.e., a sender transmits an arbitrary two-particle entangled state to a distant receiver, an arbitrary one of the n+1 agents, via the control of the others in a network. It will be shown that the outcomes in the cases that n is odd or is even are different in principle as the receiver has to perform a controlled-NOT operation on his particles for reconstructing the original arbitrary entangled state in addition to some local unitary operations in the former. Also we discuss the applications of this controlled teleporation for quantum secret sharing of classical and quantum information. As all the instances can be used to carry useful information, its efficiency for qubit approaches the maximal value

New numerical approximation for solving fractional delay differential equations of variable order using artificial neural networks

Science.gov (United States)

Zúñiga-Aguilar, C. J.; Coronel-Escamilla, A.; Gómez-Aguilar, J. F.; Alvarado-Martínez, V. M.; Romero-Ugalde, H. M.

2018-02-01

In this paper, we approximate the solution of fractional differential equations with delay using a new approach based on artificial neural networks. We consider fractional differential equations of variable order with the Mittag-Leffler kernel in the Liouville-Caputo sense. With this new neural network approach, an approximate solution of the fractional delay differential equation is obtained. Synaptic weights are optimized using the Levenberg-Marquardt algorithm. The neural network effectiveness and applicability were validated by solving different types of fractional delay differential equations, linear systems with delay, nonlinear systems with delay and a system of differential equations, for instance, the Newton-Leipnik oscillator. The solution of the neural network was compared with the analytical solutions and the numerical simulations obtained through the Adams-Bashforth-Moulton method. To show the effectiveness of the proposed neural network, different performance indices were calculated.

High-order fractional partial differential equation transform for molecular surface construction.

Science.gov (United States)

Hu, Langhua; Chen, Duan; Wei, Guo-Wei

2013-01-01

Fractional derivative or fractional calculus plays a significant role in theoretical modeling of scientific and engineering problems. However, only relatively low order fractional derivatives are used at present. In general, it is not obvious what role a high fractional derivative can play and how to make use of arbitrarily high-order fractional derivatives. This work introduces arbitrarily high-order fractional partial differential equations (PDEs) to describe fractional hyperdiffusions. The fractional PDEs are constructed via fractional variational principle. A fast fractional Fourier transform (FFFT) is proposed to numerically integrate the high-order fractional PDEs so as to avoid stringent stability constraints in solving high-order evolution PDEs. The proposed high-order fractional PDEs are applied to the surface generation of proteins. We first validate the proposed method with a variety of test examples in two and three-dimensional settings. The impact of high-order fractional derivatives to surface analysis is examined. We also construct fractional PDE transform based on arbitrarily high-order fractional PDEs. We demonstrate that the use of arbitrarily high-order derivatives gives rise to time-frequency localization, the control of the spectral distribution, and the regulation of the spatial resolution in the fractional PDE transform. Consequently, the fractional PDE transform enables the mode decomposition of images, signals, and surfaces. The effect of the propagation time on the quality of resulting molecular surfaces is also studied. Computational efficiency of the present surface generation method is compared with the MSMS approach in Cartesian representation. We further validate the present method by examining some benchmark indicators of macromolecular surfaces, i.e., surface area, surface enclosed volume, surface electrostatic potential and solvation free energy. Extensive numerical experiments and comparison with an established surface model

Automatic Assessment of Pathological Voice Quality Using Higher-Order Statistics in the LPC Residual Domain

Directory of Open Access Journals (Sweden)

JiYeoun Lee

2009-01-01

Full Text Available A preprocessing scheme based on linear prediction coefficient (LPC residual is applied to higher-order statistics (HOSs for automatic assessment of an overall pathological voice quality. The normalized skewness and kurtosis are estimated from the LPC residual and show statistically meaningful distributions to characterize the pathological voice quality. 83 voice samples of the sustained vowel /a/ phonation are used in this study and are independently assessed by a speech and language therapist (SALT according to the grade of the severity of dysphonia of GRBAS scale. These are used to train and test classification and regression tree (CART. The best result is obtained using an optima l decision tree implemented by a combination of the normalized skewness and kurtosis, with an accuracy of 92.9%. It is concluded that the method can be used as an assessment tool, providing a valuable aid to the SALT during clinical evaluation of an overall pathological voice quality.

Positive nondecreasing solutions for a multi-term fractional-order functional differential equation with integral conditions

OpenAIRE

Ahmed M. A. El-Sayed; Ebtisam O. Bin-Taher

2011-01-01

In this article, we prove the existence of positive nondecreasing solutions for a multi-term fractional-order functional differential equations. We consider Cauchy boundary problems with: nonlocal conditions, two-point boundary conditions, integral conditions, and deviated arguments.

Compound words prompt arbitrary semantic associations in conceptual memory

Directory of Open Access Journals (Sweden)

Bastien eBoutonnet

2014-03-01

Full Text Available Linguistic relativity theory has received empirical support in domains such as colour perception and object categorisation. It is unknown however, whether relations between words idiosyncratic to language impact nonverbal representations and conceptualisations. For instance, would one consider the concepts of horse and sea as related were it not for the existence of the compound seahorse? Here, we investigated such arbitrary conceptual relationships using a non-linguistic picture relatedness task in participants undergoing event-related brain potential recordings. Picture pairs arbitrarily related because of a compound and presented in the compound order elicited N400 amplitudes similar to unrelated pairs. Surprisingly, however, pictures presented in the reverse order (as in the sequence horse – sea reduced N400 amplitudes significantly, demonstrating the existence of a link in memory between these two concepts otherwise unrelated. These results break new ground in the domain of linguistic relativity by revealing predicted semantic associations driven by lexical relations intrinsic to language.

Uniqueness of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations

International Nuclear Information System (INIS)

Tran Duc Van

1994-01-01

The notion of global quasi-classical solutions of the Cauchy problems for first-order nonlinear partial differential equations is presented, some uniqueness theorems and a stability result are established by the method based on the theory of differential inclusions. In particular, the answer to an open problem of S.N. Kruzhkov is given. (author). 10 refs, 1 fig

On one two-point BVP for the fourth order linear ordinary differential equation

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Manjikashvili, M.

2017-01-01

Roč. 24, č. 2 (2017), s. 265-275 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : fourth order linear ordinary differential equations * two-point boundary value problems Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0077/gmj-2016-0077. xml

On one two-point BVP for the fourth order linear ordinary differential equation

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Manjikashvili, M.

2017-01-01

Roč. 24, č. 2 (2017), s. 265-275 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : fourth order linear ordinary differential equations * two-point boundary value problems Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2017.24.issue-2/gmj-2016-0077/gmj-2016-0077.xml

Interacting fields of arbitrary spin and N > 4 supersymmetric self-dual Yang-Mills equations

International Nuclear Information System (INIS)

Devchand, Ch.; Ogievetsky, V.

1996-06-01

We show that the self-dual Yang-Mills equations afford supersymmetrization to systems of equations invariant under global N-extended super-Poincare transformations for arbitrary values of N, without the limitation (N ≤ 4) applicable to standard non-self-dual Yang-Mills theories. These systems of equations provide novel classically consistent interactions for vector supermultiplets containing fields of spin up to N-2/2. The equations of motion of the component fields of spin greater than 1/2 are interacting variants of the first-order Dirac-Fierz equations for zero rest-mass fields of arbitrary spin. The interactions are governed by conserved currents which are constructed by an iterative procedure. In (arbitrarily extended) chiral superspace, the equations of motion for the (arbitrarily large) self-dual supermultiplet are shown to be completely equivalent to the set of algebraic supercurvature defining the self-dual superconnection. (author). 25 refs

A novel measuring method for arbitrary optical vortex by three spiral spectra

Energy Technology Data Exchange (ETDEWEB)

Ni, Bo [School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006 (China); Guo, Lana [School of Electronics and Information, Guangdong Polytechnic Normal University, Guangzhou 510665 (China); Yue, Chengfeng [School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006 (China); Tang, Zhilie, E-mail: tangzhl@scnu.edu.cn [School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006 (China)

2017-02-26

In this letter, the topological charge of non-integer vortices determined by three arbitrary spiral spectra is theoretically demonstrated for the first time. Based on the conclusion, a novel method to measure non-integer vortices is presented. This method is applicable not only to arbitrary non-integer vortex but also to arbitrary integer vortex. - Highlights: • Different non-integer vortices cannot have three spiral spectra is demonstrated. • Relationship between the non-integer topological charge and the spiral spectra is presented. • Topological charge of non-integer vortices can be determined by three arbitrary spiral spectra.

An Efficient Side-Channel Protected AES Implementation with Arbitrary Protection Order

OpenAIRE

Groß, Hannes; Mangard, Stefan; Korak, Thomas

2017-01-01

Passive physical attacks, like power analysis, pose a serious threat to the security of digital circuits. In this work, we introduce an efficient sidechannel protected Advanced Encryption Standard (AES) hardware design that is completely scalable in terms of protection order. Therefore, we revisit the private circuits scheme of Ishai et al. [13] which is known to be vulnerable to glitches. We demonstrate how to achieve resistance against multivariate higher-order attacks in the presence of gl...

Generalized solutions of nonlinear partial differential equations

CERN Document Server

Rosinger, EE

1987-01-01

During the last few years, several fairly systematic nonlinear theories of generalized solutions of rather arbitrary nonlinear partial differential equations have emerged. The aim of this volume is to offer the reader a sufficiently detailed introduction to two of these recent nonlinear theories which have so far contributed most to the study of generalized solutions of nonlinear partial differential equations, bringing the reader to the level of ongoing research.The essence of the two nonlinear theories presented in this volume is the observation that much of the mathematics concernin

Mellin-Barnes representations of Feynman diagrams, linear systems of differential equations, and polynomial solutions

International Nuclear Information System (INIS)

Kalmykov, Mikhail Yu.; Kniehl, Bernd A.

2012-05-01

We argue that the Mellin-Barnes representations of Feynman diagrams can be used for obtaining linear systems of homogeneous differential equations for the original Feynman diagrams with arbitrary powers of propagators without recourse to the integration-by-parts technique. These systems of differential equation can be used (i) for the differential reductions to sets of basic functions and (ii) for counting the numbers of master-integrals.

Lie symmetries of systems of second-order linear ordinary differential equations with constant coefficients.

Science.gov (United States)

Boyko, Vyacheslav M; Popovych, Roman O; Shapoval, Nataliya M

2013-01-01

Lie symmetries of systems of second-order linear ordinary differential equations with constant coefficients are exhaustively described over both the complex and real fields. The exact lower and upper bounds for the dimensions of the maximal Lie invariance algebras possessed by such systems are obtained using an effective algebraic approach.

OSCILLATION OF A SECOND-ORDER HALF-LINEAR NEUTRAL DAMPED DIFFERENTIAL EQUATION WITH TIME-DELAY

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In this paper,the oscillation for a class of second-order half-linear neutral damped differential equation with time-delay is studied.By means of Yang-inequality,the generalized Riccati transformation and a certain function,some new sufficient conditions for the oscillation are given for all solutions to the equation.

An Arbitrary Lagrangian-Eulerian Discretization of MHD on 3D Unstructured Grids

Energy Technology Data Exchange (ETDEWEB)

Rieben, R N; White, D A; Wallin, B K; Solberg, J M

2006-06-12

We present an arbitrary Lagrangian-Eulerian (ALE) discretization of the equations of resistive magnetohydrodynamics (MHD) on unstructured hexahedral grids. The method is formulated using an operator-split approach with three distinct phases: electromagnetic diffusion, Lagrangian motion, and Eulerian advection. The resistive magnetic dynamo equation is discretized using a compatible mixed finite element method with a 2nd order accurate implicit time differencing scheme which preserves the divergence-free nature of the magnetic field. At each discrete time step, electromagnetic force and heat terms are calculated and coupled to the hydrodynamic equations to compute the Lagrangian motion of the conducting materials. By virtue of the compatible discretization method used, the invariants of Lagrangian MHD motion are preserved in a discrete sense. When the Lagrangian motion of the mesh causes significant distortion, that distortion is corrected with a relaxation of the mesh, followed by a 2nd order monotonic remap of the electromagnetic state variables. The remap is equivalent to Eulerian advection of the magnetic flux density with a fictitious mesh relaxation velocity. The magnetic advection is performed using a novel variant of constrained transport (CT) that is valid for unstructured hexahedral grids with arbitrary mesh velocities. The advection method maintains the divergence free nature of the magnetic field and is second order accurate in regions where the solution is sufficiently smooth. For regions in which the magnetic field is discontinuous (e.g. MHD shocks) the method is limited using a novel variant of algebraic flux correction (AFC) which is local extremum diminishing (LED) and divergence preserving. Finally, we verify each stage of the discretization via a set of numerical experiments.

Automatic definition of prescription isodose for stereotactic irradiations of arteriovenous malformations

International Nuclear Information System (INIS)

Dejean, C.; Lefkopoulos, D.; Foulquier, J.N.; Schlienger, M.; Touboul, E.

2001-01-01

To evaluate dosimetric consequences generated by the automatic definition based on lesion coverage of prescription isodose. A clinical series of 124 arteriovenous malformations was analysed. Plan quality was quantified by the standard deviation of the differential dose volume histogram calculated in the lesion. We define two quantitative protocols based on lesion coverage for the automatic definition of prescription isodose using a volumetric definition of coverage (90% of lesion volume), and an isodose-based definition proposed) by RTOG (prescription isodose equals minimum isodose in the lesion divided by 0.9). We have evaluated the plans obtained for these two protocols, calculating several dose-volume indices. These indices are presented as a function of dose-volume histogram standard deviation in order to quantify the consequences of their variations for this representative series of plans. The margin our team tolerates is such that the sum of under-dosed lesion and overdosed healthy tissues factors remains lower than one. Protocol based on volumetric coverage gives results situated within this margin. Protocol based on RTOG definition produces conformation indices that could be greater than 1. The absolute dose would be decided taking into account examined dose-volume indices and clinical data. A protocol for automatic definition of prescription isodose using volumetric lesion coverage seems to be more judiciously adapted to arteriovenous malformation conformal plans in stereotactic conditions because of variations observed in the overdosage of healthy tissues. (authors)

Information guided channel hopping with an arbitrary number of transmit antennas

KAUST Repository

Yang, Yuli; Aï ssa, Sonia

2012-01-01

In order to realize the information guided channel hopping, also known as spatial modulation, with more design flexibility, in this paper we propose a novel scheme that allows operation with an arbitrary number of transmit antennas. Once the number of transmit antennas is not a power of two, the antennas' symbols are mapped by different numbers of bits. Subsequently, constellations with different orders are exploited for the modulation of radiated symbols so as to guarantee that the total number of bits transmitted at each time slot remains the same. Furthermore, we introduce a decoding algorithm with low complexity for this design. Numerical results on bit error rate performance are provided and substantiate that the proposed scheme turns out to be a promising alternative to the design of information guided channel hopping. © 2012 IEEE.

Information guided channel hopping with an arbitrary number of transmit antennas

KAUST Repository

Yang, Yuli

2012-10-01

In order to realize the information guided channel hopping, also known as spatial modulation, with more design flexibility, in this paper we propose a novel scheme that allows operation with an arbitrary number of transmit antennas. Once the number of transmit antennas is not a power of two, the antennas\\' symbols are mapped by different numbers of bits. Subsequently, constellations with different orders are exploited for the modulation of radiated symbols so as to guarantee that the total number of bits transmitted at each time slot remains the same. Furthermore, we introduce a decoding algorithm with low complexity for this design. Numerical results on bit error rate performance are provided and substantiate that the proposed scheme turns out to be a promising alternative to the design of information guided channel hopping. © 2012 IEEE.

Research on wireless communication technology based on automatic logistics system of welder

Directory of Open Access Journals (Sweden)

Sun Xuan

2018-01-01

Full Text Available In order to meet the requirements of high real-time and high stability of data transmission in automatic welding system, RTU data format and real-time communication mechanism are adopted in this system. In the automatic logistics system through the Ethernet and wireless WIFI technology will palletizer, stacker, AGV car organically together to complete the palletizer automatic crawling the goods, AGV car automatic delivery, stacking machine automatically out of the Dimensional warehouse. .

Research on wireless communication technology based on automatic logistics system of welder

OpenAIRE

Sun Xuan; Wang Zhi-yong; Ma Zhe-dong

2018-01-01

In order to meet the requirements of high real-time and high stability of data transmission in automatic welding system, RTU data format and real-time communication mechanism are adopted in this system. In the automatic logistics system through the Ethernet and wireless WIFI technology will palletizer, stacker, AGV car organically together to complete the palletizer automatic crawling the goods, AGV car automatic delivery, stacking machine automatically out of the Dimensional warehouse. .

Almost-Periodic Weak Solutions of Second-Order Neutral Delay-Differential Equations with Piecewise Constant Argument

Directory of Open Access Journals (Sweden)

Wang Li

2008-01-01

Full Text Available We investigate the existence of almost-periodic weak solutions of second-order neutral delay-differential equations with piecewise constant argument of the form , where denotes the greatest integer function, is a real nonzero constant, and is almost periodic.

Positive nondecreasing solutions for a multi-term fractional-order functional differential equation with integral conditions

Directory of Open Access Journals (Sweden)

Ahmed M. A. El-Sayed

2011-12-01

Full Text Available In this article, we prove the existence of positive nondecreasing solutions for a multi-term fractional-order functional differential equations. We consider Cauchy boundary problems with: nonlocal conditions, two-point boundary conditions, integral conditions, and deviated arguments.

Universal sequence map (USM of arbitrary discrete sequences

Directory of Open Access Journals (Sweden)

Almeida Jonas S

2002-02-01

Full Text Available Abstract Background For over a decade the idea of representing biological sequences in a continuous coordinate space has maintained its appeal but not been fully realized. The basic idea is that any sequence of symbols may define trajectories in the continuous space conserving all its statistical properties. Ideally, such a representation would allow scale independent sequence analysis – without the context of fixed memory length. A simple example would consist on being able to infer the homology between two sequences solely by comparing the coordinates of any two homologous units. Results We have successfully identified such an iterative function for bijective mappingψ of discrete sequences into objects of continuous state space that enable scale-independent sequence analysis. The technique, named Universal Sequence Mapping (USM, is applicable to sequences with an arbitrary length and arbitrary number of unique units and generates a representation where map distance estimates sequence similarity. The novel USM procedure is based on earlier work by these and other authors on the properties of Chaos Game Representation (CGR. The latter enables the representation of 4 unit type sequences (like DNA as an order free Markov Chain transition table. The properties of USM are illustrated with test data and can be verified for other data by using the accompanying web-based tool:http://bioinformatics.musc.edu/~jonas/usm/. Conclusions USM is shown to enable a statistical mechanics approach to sequence analysis. The scale independent representation frees sequence analysis from the need to assume a memory length in the investigation of syntactic rules.

Automatic metal parts inspection: Use of thermographic images and anomaly detection algorithms

Science.gov (United States)

Benmoussat, M. S.; Guillaume, M.; Caulier, Y.; Spinnler, K.

2013-11-01

A fully-automatic approach based on the use of induction thermography and detection algorithms is proposed to inspect industrial metallic parts containing different surface and sub-surface anomalies such as open cracks, open and closed notches with different sizes and depths. A practical experimental setup is developed, where lock-in and pulsed thermography (LT and PT, respectively) techniques are used to establish a dataset of thermal images for three different mockups. Data cubes are constructed by stacking up the temporal sequence of thermogram images. After the reduction of the data space dimension by means of denoising and dimensionality reduction methods; anomaly detection algorithms are applied on the reduced data cubes. The dimensions of the reduced data spaces are automatically calculated with arbitrary criterion. The results show that, when reduced data cubes are used, the anomaly detection algorithms originally developed for hyperspectral data, the well-known Reed and Xiaoli Yu detector (RX) and the regularized adaptive RX (RARX), give good detection performances for both surface and sub-surface defects in a non-supervised way.

An optimization method for the distance between exits of buildings considering uncertainties based on arbitrary polynomial chaos expansion

NARCIS (Netherlands)

Xie, Qimiao; Wang, Jinhui; Lu, Shouxiang; Hensen, J.L.M.

2016-01-01

The distance between exits is an important design parameter in fire safety design of buildings. In order to find the optimal distance between exits under uncertainties with a low computational cost, the surrogate model (i.e. approximation model) of evacuation time is constructed by the arbitrary

Local linearization methods for the numerical integration of ordinary differential equations: An overview

International Nuclear Information System (INIS)

Jimenez, J.C.

2009-06-01

Local Linearization (LL) methods conform a class of one-step explicit integrators for ODEs derived from the following primary and common strategy: the vector field of the differential equation is locally (piecewise) approximated through a first-order Taylor expansion at each time step, thus obtaining successive linear equations that are explicitly integrated. Hereafter, the LL approach may include some additional strategies to improve that basic affine approximation. Theoretical and practical results have shown that the LL integrators have a number of convenient properties. These include arbitrary order of convergence, A-stability, linearization preserving, regularity under quite general conditions, preservation of the dynamics of the exact solution around hyperbolic equilibrium points and periodic orbits, integration of stiff and high-dimensional equations, low computational cost, and others. In this paper, a review of the LL methods and their properties is presented. (author)

Third-order operator-differential equations with discontinuous coefficients and operators in the boundary conditions

Directory of Open Access Journals (Sweden)

Araz R. Aliev

2013-10-01

Full Text Available We study a third-order operator-differential equation on the semi-axis with a discontinuous coefficient and boundary conditions which include an abstract linear operator. Sufficient conditions for the well-posed and unique solvability are found by means of properties of the operator coefficients in a Sobolev-type space.

Parasupersymmetry and N-fold supersymmetry in quantum many-body systems. I: General formalism and second order

International Nuclear Information System (INIS)

Tanaka, Toshiaki

2007-01-01

We propose an elegant formulation of parafermionic algebra and parasupersymmetry of arbitrary order in quantum many-body systems without recourse to any specific matrix representation of parafermionic operators and any kind of deformed algebra. Within our formulation, we show generically that every parasupersymmetric quantum system of order p consists of N-fold supersymmetric pairs with N≤p and thus has weak quasi-solvability and isospectral property. We also propose a new type of non-linear supersymmetries, called quasi-parasupersymmetry, which is less restrictive than parasupersymmetry and is different from N-fold supersymmetry even in one-body systems though the conserved charges are represented by higher-order linear differential operators. To illustrate how our formulation works, we construct second-order parafermionic algebra and three simple examples of parasupersymmetric quantum systems of order 2, one is essentially equivalent to the one-body Rubakov-Spiridonov type and the others are two-body systems in which two supersymmetries are folded. In particular, we show that the first model admits a generalized 2-fold superalgebra

A Clustering-Based Automatic Transfer Function Design for Volume Visualization

Directory of Open Access Journals (Sweden)

Tianjin Zhang

2016-01-01

Full Text Available The two-dimensional transfer functions (TFs designed based on intensity-gradient magnitude (IGM histogram are effective tools for the visualization and exploration of 3D volume data. However, traditional design methods usually depend on multiple times of trial-and-error. We propose a novel method for the automatic generation of transfer functions by performing the affinity propagation (AP clustering algorithm on the IGM histogram. Compared with previous clustering algorithms that were employed in volume visualization, the AP clustering algorithm has much faster convergence speed and can achieve more accurate clustering results. In order to obtain meaningful clustering results, we introduce two similarity measurements: IGM similarity and spatial similarity. These two similarity measurements can effectively bring the voxels of the same tissue together and differentiate the voxels of different tissues so that the generated TFs can assign different optical properties to different tissues. Before performing the clustering algorithm on the IGM histogram, we propose to remove noisy voxels based on the spatial information of voxels. Our method does not require users to input the number of clusters, and the classification and visualization process is automatic and efficient. Experiments on various datasets demonstrate the effectiveness of the proposed method.

Stabilization at almost arbitrary points for chaotic systems

International Nuclear Information System (INIS)

Huang, C.-S.; Lian, K.-Y.; Su, C.-H.; Wu, J.-W.

2008-01-01

We consider how to design a feasible control input for chaotic systems via a suitable input channel to achieve the stabilization at arbitrary points. Regarding the nonlinear systems without naturally defined input vectors, we propose a local stabilization controller which works for almost arbitrary points. Subsequently, according to topologically transitive property for chaotic systems, the feedback control force is activated only when the trajectory passes through the neighboring region of the regulated point. Hence the global stabilization is achieved whereas the control effort of the hybrid controller is extremely low

Bisimulation on Markov Processes over Arbitrary Measurable Spaces

DEFF Research Database (Denmark)

Bacci, Giorgio; Bacci, Giovanni; Larsen, Kim Guldstrand

2014-01-01

We introduce a notion of bisimulation on labelled Markov Processes over generic measurable spaces in terms of arbitrary binary relations. Our notion of bisimulation is proven to coincide with the coalgebraic definition of Aczel and Mendler in terms of the Giry functor, which associates with a mea......We introduce a notion of bisimulation on labelled Markov Processes over generic measurable spaces in terms of arbitrary binary relations. Our notion of bisimulation is proven to coincide with the coalgebraic definition of Aczel and Mendler in terms of the Giry functor, which associates...

The focal boundary value problem for strongly singular higher-order nonlinear functional-differential equations

Czech Academy of Sciences Publication Activity Database

Mukhigulashvili, Sulkhan; Půža, B.

2015-01-01

Roč. 2015, January (2015), s. 17 ISSN 1687-2770 Institutional support: RVO:67985840 Keywords : higher order nonlinear functional-differential equations * two-point right-focal boundary value problem * strong singularity Subject RIV: BA - General Mathematics Impact factor: 0.642, year: 2015 http://link.springer.com/article/10.1186%2Fs13661-014-0277-1

Arbitrary protein−protein docking targets biologically relevant interfaces

International Nuclear Information System (INIS)

Martin, Juliette; Lavery, Richard

2012-01-01

Protein-protein recognition is of fundamental importance in the vast majority of biological processes. However, it has already been demonstrated that it is very hard to distinguish true complexes from false complexes in so-called cross-docking experiments, where binary protein complexes are separated and the isolated proteins are all docked against each other and scored. Does this result, at least in part, reflect a physical reality? False complexes could reflect possible nonspecific or weak associations. In this paper, we investigate the twilight zone of protein-protein interactions, building on an interesting outcome of cross-docking experiments: false complexes seem to favor residues from the true interaction site, suggesting that randomly chosen partners dock in a non-random fashion on protein surfaces. Here, we carry out arbitrary docking of a non-redundant data set of 198 proteins, with more than 300 randomly chosen "probe" proteins. We investigate the tendency of arbitrary partners to aggregate at localized regions of the protein surfaces, the shape and compositional bias of the generated interfaces, and the potential of this property to predict biologically relevant binding sites. We show that the non-random localization of arbitrary partners after protein-protein docking is a generic feature of protein structures. The interfaces generated in this way are not systematically planar or curved, but tend to be closer than average to the center of the proteins. These results can be used to predict biological interfaces with an AUC value up to 0.69 alone, and 0.72 when used in combination with evolutionary information. An appropriate choice of random partners and number of docking models make this method computationally practical. It is also noted that nonspecific interfaces can point to alternate interaction sites in the case of proteins with multiple interfaces. We illustrate the usefulness of arbitrary docking using PEBP (Phosphatidylethanolamine binding

Arbitrary protein−protein docking targets biologically relevant interfaces

Directory of Open Access Journals (Sweden)

Martin Juliette

2012-05-01

Full Text Available Abstract Background Protein-protein recognition is of fundamental importance in the vast majority of biological processes. However, it has already been demonstrated that it is very hard to distinguish true complexes from false complexes in so-called cross-docking experiments, where binary protein complexes are separated and the isolated proteins are all docked against each other and scored. Does this result, at least in part, reflect a physical reality? False complexes could reflect possible nonspecific or weak associations. Results In this paper, we investigate the twilight zone of protein-protein interactions, building on an interesting outcome of cross-docking experiments: false complexes seem to favor residues from the true interaction site, suggesting that randomly chosen partners dock in a non-random fashion on protein surfaces. Here, we carry out arbitrary docking of a non-redundant data set of 198 proteins, with more than 300 randomly chosen "probe" proteins. We investigate the tendency of arbitrary partners to aggregate at localized regions of the protein surfaces, the shape and compositional bias of the generated interfaces, and the potential of this property to predict biologically relevant binding sites. We show that the non-random localization of arbitrary partners after protein-protein docking is a generic feature of protein structures. The interfaces generated in this way are not systematically planar or curved, but tend to be closer than average to the center of the proteins. These results can be used to predict biological interfaces with an AUC value up to 0.69 alone, and 0.72 when used in combination with evolutionary information. An appropriate choice of random partners and number of docking models make this method computationally practical. It is also noted that nonspecific interfaces can point to alternate interaction sites in the case of proteins with multiple interfaces. We illustrate the usefulness of arbitrary docking

Usage of aids monitoring in automatic braking systems of modern cars

OpenAIRE

Dembitskyi V.; Mazylyuk P.; Sitovskyi O.

2016-01-01

Increased safety can be carried out at the expense the installation on vehicles of automatic braking systems, that monitor the traffic situation and the actions of the driver. In this paper considered the advantages and disadvantages of automatic braking systems, were analyzed modern tracking tools that are used in automatic braking systems. Based on the statistical data on accidents, are set the main dangers, that the automatic braking system will be reduced. In order to ensure the acc...

Unsplit complex frequency shifted perfectly matched layer for second-order wave equation using auxiliary differential equations.

Science.gov (United States)

Gao, Yingjie; Zhang, Jinhai; Yao, Zhenxing

2015-12-01

The complex frequency shifted perfectly matched layer (CFS-PML) can improve the absorbing performance of PML for nearly grazing incident waves. However, traditional PML and CFS-PML are based on first-order wave equations; thus, they are not suitable for second-order wave equation. In this paper, an implementation of CFS-PML for second-order wave equation is presented using auxiliary differential equations. This method is free of both convolution calculations and third-order temporal derivatives. As an unsplit CFS-PML, it can reduce the nearly grazing incidence. Numerical experiments show that it has better absorption than typical PML implementations based on second-order wave equation.

Generalized Sudan's List Decoding for Order Domain Codes

DEFF Research Database (Denmark)

Geil, Hans Olav; Matsumoto, Ryutaroh

2007-01-01

We generalize Sudan's list decoding algorithm without multiplicity to evaluation codes coming from arbitrary order domains. The number of correctable errors by the proposed method is larger than the original list decoding without multiplicity....

Towards automatic pulmonary nodule management in lung cancer screening with deep learning.

Science.gov (United States)

Ciompi, Francesco; Chung, Kaman; van Riel, Sarah J; Setio, Arnaud Arindra Adiyoso; Gerke, Paul K; Jacobs, Colin; Scholten, Ernst Th; Schaefer-Prokop, Cornelia; Wille, Mathilde M W; Marchianò, Alfonso; Pastorino, Ugo; Prokop, Mathias; van Ginneken, Bram

2017-04-19

The introduction of lung cancer screening programs will produce an unprecedented amount of chest CT scans in the near future, which radiologists will have to read in order to decide on a patient follow-up strategy. According to the current guidelines, the workup of screen-detected nodules strongly relies on nodule size and nodule type. In this paper, we present a deep learning system based on multi-stream multi-scale convolutional networks, which automatically classifies all nodule types relevant for nodule workup. The system processes raw CT data containing a nodule without the need for any additional information such as nodule segmentation or nodule size and learns a representation of 3D data by analyzing an arbitrary number of 2D views of a given nodule. The deep learning system was trained with data from the Italian MILD screening trial and validated on an independent set of data from the Danish DLCST screening trial. We analyze the advantage of processing nodules at multiple scales with a multi-stream convolutional network architecture, and we show that the proposed deep learning system achieves performance at classifying nodule type that surpasses the one of classical machine learning approaches and is within the inter-observer variability among four experienced human observers.

Closed description of arbitrariness in resolving quantum master equation

Energy Technology Data Exchange (ETDEWEB)

Batalin, Igor A., E-mail: batalin@lpi.ru [P.N. Lebedev Physical Institute, Leninsky Prospect 53, 119 991 Moscow (Russian Federation); Tomsk State Pedagogical University, Kievskaya St. 60, 634061 Tomsk (Russian Federation); Lavrov, Peter M., E-mail: lavrov@tspu.edu.ru [Tomsk State Pedagogical University, Kievskaya St. 60, 634061 Tomsk (Russian Federation); National Research Tomsk State University, Lenin Av. 36, 634050 Tomsk (Russian Federation)

2016-07-10

In the most general case of the Delta exact operator valued generators constructed of an arbitrary Fermion operator, we present a closed solution for the transformed master action in terms of the original master action in the closed form of the corresponding path integral. We show in detail how that path integral reduces to the known result in the case of being the Delta exact generators constructed of an arbitrary Fermion function.

Higher-Order and Symbolic Computation. LISP and Symbolic Computationditorial

DEFF Research Database (Denmark)

Danvy, Olivier; Dybvig, R. Kent; Lawall, Julia

2008-01-01

system for these static checks and a corresponding type-inference algorithm. In "An Investigation of Jones Optimality and BTI-Universal Specializers," Robert Glueck establishes a connection between Jones optimal-program specializers and binding-time improvers. This article completes a study started...... at ASIA-PEPM 2002 [1]. In "On the Implementation of Automatic Differentiation Tools," Christian H. Bischof, Paul D. Hovland, and Boyana Norris present a survey of some recent tools for the Automatic Differentiation technology (concentrating mainly on ADIC, ADIFOR and sketching XAIF). They also offer...... for removing tuple constructions and tuple selections. This technique solves the problem of efficiently passing tuples to polymorphic functions by avoiding extra memory operations in selecting components of the tuple....

CPDES2: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in two dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on two-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect indices which is vectorizable on some of the newer scientific computers.

CPDES3: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on three-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect induces which is vectorizable on some of the newer scientific computers.

New Traveling Wave Solutions of the Higher Dimensional Nonlinear Partial Differential Equation by the Exp-Function Method

Directory of Open Access Journals (Sweden)

Hasibun Naher

2012-01-01

Full Text Available We construct new analytical solutions of the (3+1-dimensional modified KdV-Zakharov-Kuznetsev equation by the Exp-function method. Plentiful exact traveling wave solutions with arbitrary parameters are effectively obtained by the method. The obtained results show that the Exp-function method is effective and straightforward mathematical tool for searching analytical solutions with arbitrary parameters of higher-dimensional nonlinear partial differential equation.

Accelerating flight: Edge with arbitrary acceleration

CSIR Research Space (South Africa)

Gledhill, Irvy MA

2011-11-01

Full Text Available This study concludes the possession of a theoretical framework for arbitrary manoeuvre which allows us to keep an eye on transformations. In the theory, relative frame equations are useful in guiding us in what to look for. The code...

Conversion of KEGG metabolic pathways to SBGN maps including automatic layout.

Science.gov (United States)

Czauderna, Tobias; Wybrow, Michael; Marriott, Kim; Schreiber, Falk

2013-08-16

Biologists make frequent use of databases containing large and complex biological networks. One popular database is the Kyoto Encyclopedia of Genes and Genomes (KEGG) which uses its own graphical representation and manual layout for pathways. While some general drawing conventions exist for biological networks, arbitrary graphical representations are very common. Recently, a new standard has been established for displaying biological processes, the Systems Biology Graphical Notation (SBGN), which aims to unify the look of such maps. Ideally, online repositories such as KEGG would automatically provide networks in a variety of notations including SBGN. Unfortunately, this is non-trivial, since converting between notations may add, remove or otherwise alter map elements so that the existing layout cannot be simply reused. Here we describe a methodology for automatic translation of KEGG metabolic pathways into the SBGN format. We infer important properties of the KEGG layout and treat these as layout constraints that are maintained during the conversion to SBGN maps. This allows for the drawing and layout conventions of SBGN to be followed while creating maps that are still recognizably the original KEGG pathways. This article details the steps in this process and provides examples of the final result.

Automatic associations with the sensory aspects of smoking: Positive in habitual smokers but negative in non-smokers

OpenAIRE

Huijding, Jorg; Jong, Peter

2006-01-01

textabstractTo test whether pictorial stimuli that focus on the sensory aspects of smoking elicit different automatic affective associations in smokers than in non-smokers, 31 smoking and 33 non-smoking students completed a single target IAT. Explicit attitudes were assessed using a semantic differential. Automatic affective associations were positive in smokers but negative in non-smokers. Only automatic affective associations but not self-reported attitudes were significantly correlated wit...

Modelling magnetic laminations under arbitrary starting state and flux waveform

International Nuclear Information System (INIS)

Bottauscio, Oriano; Chiampi, Mario; Ragusa, Carlo

2005-01-01

A numerical model able to predict the behaviour of a magnetic sheet under arbitrary supply conditions has been developed. The electromagnetic field problem is formulated in terms of an electric vector potential, which provides the magnetic field strength evolution. The hysteretic behaviour of the material is represented through the dynamic Preisach model where the activation law of the bi-state operators is modified in order to guarantee a smooth response. The problem has been solved through a time step procedure using the fixed Point technique for handling nonlinearity. The model has been validated by comparison with suitable experiments and it is applied to the investigation of the influence of the materials' starting state on the magnetic behaviour

ADI splitting schemes for a fourth-order nonlinear partial differential equation from image processing

KAUST Repository

Calatroni, Luca

2013-08-01

We present directional operator splitting schemes for the numerical solution of a fourth-order, nonlinear partial differential evolution equation which arises in image processing. This equation constitutes the H -1-gradient flow of the total variation and represents a prototype of higher-order equations of similar type which are popular in imaging for denoising, deblurring and inpainting problems. The efficient numerical solution of this equation is very challenging due to the stiffness of most numerical schemes. We show that the combination of directional splitting schemes with implicit time-stepping provides a stable and computationally cheap numerical realisation of the equation.

Linear or linearizable first-order delay ordinary differential equations and their Lie point symmetries

Science.gov (United States)

Dorodnitsyn, Vladimir A.; Kozlov, Roman; Meleshko, Sergey V.; Winternitz, Pavel

2018-05-01

A recent article was devoted to an analysis of the symmetry properties of a class of first-order delay ordinary differential systems (DODSs). Here we concentrate on linear DODSs, which have infinite-dimensional Lie point symmetry groups due to the linear superposition principle. Their symmetry algebra always contains a two-dimensional subalgebra realized by linearly connected vector fields. We identify all classes of linear first-order DODSs that have additional symmetries, not due to linearity alone, and we present representatives of each class. These additional symmetries are then used to construct exact analytical particular solutions using symmetry reduction.