Global and local optimization using radial basis function response surface models

The focus of this paper is the optimization of complex multi-parameter systems. We consider systems in which the objective function is not known explicitly, and can only be evaluated through computationally intensive numerical simulation or through costly physical experiments. The objective function may also contain many local extrema which may be of interest. Given objective function values at a scattered set of parameter values, we develop a response surface model that can dramatically reduce the required computation time for parameter optimization runs. The response surface model is developed using radial basis functions, producing a model whose objective function values match those of the original system at all sampled data points. Interpolation to any other point is easily accomplished and generates a model which represents the system over the entire parameter space. This paper presents the details of the use of radial basis functions to transform scattered data points, obtained from a complex continuum mechanics simulation of explosive materials, into a response surface model of a function over the given parameter space. Response surface methodology and radial basis functions are discussed in general and are applied to a global optimization problem for an explosive oil well penetrator.     

Interaural intensity discrimination: insensitivity at 1000 Hz

Recent data from three laboratories have replicated Mills' [J. Acoust. Soc. Am. 32, 132-134 (1960)] finding that interaural intensity discrimination is relatively poorer for tones of 1000 Hz than for tones of either higher or lower frequencies. To get a finer look at this frequency effect, interaural intensity difference thresholds were obtained from four subjects for tones of several frequencies around 1000 Hz. An adaptive two-interval forced-choice procedure was employed, in which the overall intensity of the signals was varied randomly in order to prevent subjects from listening to monaural loudness changes. Despite large intersubject differences in overall sensitivity to interaural intensity differences, all four subjects showed a local peak in their threshold functions at or near 1000 Hz. This curious "1000-Hz effect" might be explained by imagining that an interaural intensity comparator operates more efficiently as frequency increases, but that a peripheral interaural intensity difference to interaural-time difference conversion contributes to laterality judgments for low-frequency tones, thus acting to lower thresholds again for frequencies below 1000 Hz.     

Binaural modulation masking

Modulation thresholds were measured in three subjects for a sinusoidally amplitude-modulated (SAM) wideband noise (the signal) in the presence of a second amplitude-modulated wideband noise (the masker). In monaural conditions (Mm-Sm) masker and signal were presented to only one ear; in binaural conditions (M0-S pi) the masker was presented diotically while the phase of modulation of the SAM noise signal was inverted in one ear relative to the other. In experiment 1 masker modulation frequency (fm) was fixed at 16 Hz, and signal modulation frequency (fs) was varied from 2-512 Hz. For monaural presentation, masking generally decreased as fs diverged from fm, although there was a secondary increase in masking for very low signal modulation frequencies, as reported previously [Bacon and Grantham, J. Acoust. Soc. Am. 85, 2575-2580 (1989)]. The binaural masking patterns did not show this low-frequency upturn: binaural thresholds continued to improve as fs decreased from 16 to 2 Hz. Thus, comparing masked monaural and masked binaural thresholds, there was an average binaural advantage, or masking-level difference (MLD) of 9.4 dB at fs = 2 Hz and 5.3 dB at fs = 4 Hz. In addition, there were positive MLDs for the on-frequency condition (fm = fs = 16 Hz: average MLD = 4.4 dB) and for the highest signal frequency tested (fs = 512 Hz: average MLD = 7.3 dB). In experiment 2 the signal was a SAM noise (fs = 16 Hz), and the masker was a wideband noise, amplitude-modulated by a narrow band of noise centered at fs. There was no effect on monaural or binaural thresholds as masker modulator bandwidth was varied from 4 to 20 Hz (the average MLD remained constant at 8.0 dB), which suggests that the observed "tuning" for modulation may be based on temporal pattern discrimination and not on a critical-band-like filtering mechanism. In a final condition the masker modulator was a 10-Hz-wide band of noise centered at the 64-Hz signal modulation frequency. The average MLD in this case was 7.4 dB. The results are discussed in terms of various binaural capacities that probably play a role in binaural release from modulation masking, including detection of varying interaural intensity differences (IIDs) and discrimination of interaural correlation.     

Frobenius pseudoprimes

The proliferation of probable prime tests in recent years has produced a plethora of definitions with the word ``pseudoprime' in them. Examples include pseudoprimes, Euler pseudoprimes, strong pseudoprimes, Lucas pseudoprimes, strong Lucas pseudoprimes, extra strong Lucas pseudoprimes and Perrin pseudoprimes. Though these tests represent a wealth of ideas, they exist as a hodge-podge of definitions rather than as examples of a more general theory. It is the goal of this paper to present a way of viewing many of these tests as special cases of a general principle, as well as to re-formulate them in the context of finite fields.

One aim of the reformulation is to enable the creation of stronger tests; another is to aid in proving results about large classes of pseudoprimes.

     

An improved error bound for a finite element approximation of a model for phase separation of a multi-component alloy

Using the approach of Rulla (1996 SIAM J. Numer. Anal. 33, 68-87)for analysing the time discretization error and assuming moreregularity on the initial data, we improve on the error boundderived by Barrett and Blowey (1996 IMA J. Numer. Anal. 16,257-287) for a fully practical piecewise linear finite elementapproximation with a backward Euler time discretization of amodel for phase separation of a multi-component alloy.     

Synthesis and properties of unusual π-networks : Cycloocta[1,2,3,4]-def]benzo[3,4]cyclobuta[6,7[biphenylene

The title compound, an example of a hydrocarbon containing three linearly fused 4n-membered rings, has been synthesized and characterized as a thermally stable species possessing no unusual chemical reactivity. Despite the presence of these conventionally antiaromatic structural units this hydrocarbon is predicted to have a net diamagnetic ring current. The NMR results, which are discussed in some detail, are shown to be consistent with this expectation. The impact of this result on the NMR criterion of aromaticity is discussed briefly.     

Forty years of metrology with synchrotron radiation at SURF

The Canadian Light Source, the University of Western Ontario, IBM Canada and BigBangwidth during 2006 developed the Remote Beamline Access (RBA) system. This system was demonstrated using a standalone soft X-ray machine as well on the VESPERS beamline. The RBA system serves as a proof of concept and basis for ScienceStudio, which is now under development.