Analysis of trigonometric implicit Runge-Kutta methods

Using generalized collocation techniques based on fitting functions that are trigonometric (rather than algebraic as in classical integrators), we develop a new class of multistage, one-step, variable stepsize, and variable coefficients implicit Runge-Kutta methods to solve oscillatory ODE problems. The coefficients of the methods are functions of the frequency and the stepsize. We refer to this class as trigonometric implicit Runge-Kutta (TIRK) methods. They integrate an equation exactly if its solution is a trigonometric polynomial with a known frequency. We characterize the order and A-stability of the methods and establish results similar to that of classical algebraic collocation RK methods. (c) 2006 Elsevier B.V. All rights reserved.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

The cost of life expectancy and the implicit social valuation of life

A new method of estimating the economic value of life is proposed. Using cross-country data, an equation is estimated to explain life expectancy as a function of real consumption of goods and services. The associated cost function for life expectancy in terms of the prices of specific goods and services is used to estimate the cost of a reduction in age-specific mortality rates sufficient to save the life of one person. The cost of saving a life in OECD countries is as much as 1000 times that in the poorest countries. Ethical implications are discussed.

Chemical synthesis, characterization and activity of RK-1, a novel alpha-defensin-related peptide

The 32-residue peptide, RK-1, a novel kidney-derived three disulfide-bonded member of the antimicrobial alpha-defensin family, was synthesized by the continuous now Fmoc-solid phase method. The crude, cleaved and S-reduced Linear peptide was both efficiently folded and oxidized in an acidic solution of aqueous dimethyl sulfoxide. Following purification of the resulting product, it was shown by a variety of analytical techniques, including matrix assisted laser desorption time of flight mass spectrometry, to possess a very high degree of purity. The disulfide bond pairing of the synthetic peptide was determined by H-1-NMR spectroscopy and confirmed to be a Cys(1)-Cys(6), Cys(2)-Cys(4), Cys(3)-Cys(5) arrangement similar to other mammalian alpha-defensin peptides. The synthetic RK-1 was also shown to inhibit the growth of Escherichia coli type strain NCTC 10418, Copyright (C) 2000 European Peptide Society and John Wiley & Sons, Ltd.

Three-dimensional structure of RK-1: A novel alpha-defensin peptide

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of RK-1, an antimicrobial peptide from rabbit kidney recently discovered from homology screening based on the distinctive physicochemical properties of the corticostatins/defensins. RK-1 consists of 32 residues, including six cysteines arranged into three disulfide bonds. It exhibits antimicrobial activity against Escherichia coli and activates Ca2+ channels in vitro. Through its physicochemical similarity, identical cysteine spacing, and linkage to the corticostatins/defensins, it was presumed to be a member of this family. However, RK-1 lacks both a large number of arginines in the primary sequence and a high overall positive charge, which are characteristic of this family of peptides. The three-dimensional solution structure, determined by NMR, consists of a triple-stranded antiparallel beta -sheet and a series of turns and is similar to the known structures of other alpha -defensins. This has enabled the definitive classification of RK-1 as a member of this family of antimicrobial peptides. Ultracentrifuge measurements confirmed that like rabbit neutrophil defensins, RK-1 is monomeric in solution, in contrast to human neutrophil defensins, which are dimeric.

Stereotypic explanatory bias: Implicit stereotyping as a predictor of discrimination

Two experiments examined whether a measure of implicit stereotyping based on the tendency to explain Black stereotype-incongruent events more often than Black stereotype-congruent events (Stereotypic Explanatory Bias or SEB) is predictive of behavior toward a partner in an interracial interaction. In Experiment I SEB predicted White males' choice to ask stereotypic questions of a Black female (but not a White male or White female) in an interview. In Experiment 2 the type of explanation (internal or external attribution) made for stereotype-inconsistency was examined. Results showed that White participants who made internal attributions for Black stereotype-incongruent behavior were rated more positively and those who made external attributions were rated more negatively by a Black male confederate. These results point to the potential of implicit stereotyping as an important predictor of behavior in an interracial interaction. (C) 2002 Elsevier Science (USA). All rights reserved.

Implicit Taylor methods for stiff stochastic differential equations

In this paper we discuss implicit Taylor methods for stiff Ito stochastic differential equations. Based on the relationship between Ito stochastic integrals and backward stochastic integrals, we introduce three implicit Taylor methods: the implicit Euler-Taylor method with strong order 0.5, the implicit Milstein-Taylor method with strong order 1.0 and the implicit Taylor method with strong order 1.5. The mean-square stability properties of the implicit Euler-Taylor and Milstein-Taylor methods are much better than those of the corresponding semi-implicit Euler and Milstein methods and these two implicit methods can be used to solve stochastic differential equations which are stiff in both the deterministic and the stochastic components. Numerical results are reported to show the convergence properties and the stability properties of these three implicit Taylor methods. The stability analysis and numerical results show that the implicit Euler-Taylor and Milstein-Taylor methods are very promising methods for stiff stochastic differential equations.

Can anticipatory skills be learned through implicit video-based perceptual training?

The aim of this experiment was to determine the effectiveness of two video-based perceptual training approaches designed to improve the anticipatory skills of junior tennis players. Players were assigned equally to an explicit learning group, an implicit learning group, a placebo group or a control group. A progressive temporal occlusion paradigm was used to examine, before and after training, the ability of the players to predict the direction of an opponent's service in an in-vivo on-court setting. The players responded either through hitting a return stroke or making a verbal prediction of stroke direction. Results revealed that the implicit learning group, whose training required them to predict serve speed direction while viewing temporally occluded video footage of the return-of-serve scenario, significantly improved their prediction accuracy after the training intervention. However, this training effect dissipated after a 32 day unfilled retention interval. The explicit learning group, who received instructions about the specific aspects of the pre-contact service kinematics that are informative with respect to service direction, did not demonstrate any significant performance improvements after the intervention. This, together with the absence of any significant improvements for the placebo and control groups, demonstrated that the improvement observed for the implicit learning group was not a consequence of either expectancy or familiarity effects.

A four-stage index 2 Diagonally Implicit Runge-Kutta method

High index Differential Algebraic Equations (DAEs) force standard numerical methods to lower order. Implicit Runge-Kutta methods such as RADAU5 handle high index problems but their fully implicit structure creates significant overhead costs for large problems. Singly Diagonally Implicit Runge-Kutta (SDIRK) methods offer lower costs for integration. This paper derives a four-stage, index 2 Explicit Singly Diagonally Implicit Runge-Kutta (ESDIRK) method. By introducing an explicit first stage, the method achieves second order stage calculations. After deriving and solving appropriate order conditions., numerical examples are used to test the proposed method using fixed and variable step size implementations. (C) 2001 IMACS. Published by Elsevier Science B.V. All rights reserved.

Implicit vector equilibrium problems with applications

Let X and Y be Hausdorff topological vector spaces, K a nonempty, closed, and convex subset of X, C: K--> 2(Y) a point-to-set mapping such that for any x is an element of K, C(x) is a pointed, closed, and convex cone in Y and int C(x) not equal 0. Given a mapping g : K --> K and a vector valued bifunction f : K x K - Y, we consider the implicit vector equilibrium problem (IVEP) of finding x* is an element of K such that f (g(x*), y) is not an element of - int C(x) for all y is an element of K. This problem generalizes the (scalar) implicit equilibrium problem and implicit variational inequality problem. We propose the dual of the implicit vector equilibrium problem (DIVEP) and establish the equivalence between (IVEP) and (DIVEP) under certain assumptions. Also, we give characterizations of the set of solutions for (IVP) in case of nonmonotonicity, weak C-pseudomonotonicity, C-pseudomonotonicity, and strict C-pseudomonotonicity, respectively. Under these assumptions, we conclude that the sets of solutions are nonempty, closed, and convex. Finally, we give some applications of (IVEP) to vector variational inequality problems and vector optimization problems. (C) 2003 Elsevier Science Ltd. All rights reserved.

Explicit/implicit partitioning and a new explicit form of the generalized alpha method

Most finite element packages use the Newmark algorithm for time integration of structural dynamics. Various algorithms have been proposed to better optimize the high frequency dissipation of this algorithm. Hulbert and Chung proposed both implicit and explicit forms of the generalized alpha method. The algorithms optimize high frequency dissipation effectively, and despite recent work on algorithms that possess momentum conserving/energy dissipative properties in a non-linear context, the generalized alpha method remains an efficient way to solve many problems, especially with adaptive timestep control. However, the implicit and explicit algorithms use incompatible parameter sets and cannot be used together in a spatial partition, whereas this can be done for the Newmark algorithm, as Hughes and Liu demonstrated, and for the HHT-alpha algorithm developed from it. The present paper shows that the explicit generalized alpha method can be rewritten so that it becomes compatible with the implicit form. All four algorithmic parameters can be matched between the explicit and implicit forms. An element interface between implicit and explicit partitions can then be used, analogous to that devised by Hughes and Liu to extend the Newmark method. The stability of the explicit/implicit algorithm is examined in a linear context and found to exceed that of the explicit partition. The element partition is significantly less dissipative of intermediate frequencies than one using the HHT-alpha method. The explicit algorithm can also be rewritten so that the discrete equation of motion evaluates forces from displacements and velocities found at the predicted mid-point of a cycle. Copyright (C) 2003 John Wiley Sons, Ltd.