Adaptive Algorithm for Constrained Least-Squares Problems

This paper is concerned with the implementation and testing of an algorithm for solving constrained least-squares problems. The algorithm is an adaptation to the least-squares case of sequential quadratic programming (SQP) trust-region methods for solving general constrained optimization problems. At each iteration, our local quadratic subproblem includes the use of the Gauss–Newton approximation but also encompasses a structured secant approximation along with tests of when to use this approximation. This method has been tested on a selection of standard problems. The results indicate that, for least-squares problems, the approach taken here is a viable alternative to standard general optimization methods such as the Byrd–Omojokun trust-region method and the Powell damped BFGS line search method.     

Global Phase Time and Wave Function for the Kantowski-Sachs Anisotropic Universe

A consistent quantization with a clear notion of time and evolution is given for the anisotropic Kantowski–Sachs cosmological model. It is shown that a suitable coordinate choice allows to obtain a solution of the Wheeler–DeWitt equation in the form of definite energy states, and that the results can be associated to two disjoint equivalent theories, one for each sheet of the constraint surface.     

A primer on perfect simulation for spatial point processes

This primer provides a self-contained exposition of the case where spatial birth-and-death processes are used for perfect simulation of locally stable point processes. Particularly, a simple dominating coupling from the past (CFTP) algorithm and the CFTP algorithms introduced in [13], [14], and [5] are studied. Some empirical results for the algorithms are discussed. Received: 30 June 2002     

Large set of P3‐decompositions

Let G=(V(G),E(G)) be a graph. A (n,G, λ)‐GD is a partition of the edges of λK_n into subgraphs (G‐blocks), each of which is isomorphic to G. The (n,G,λ)‐GD is named as graph design for G or G‐decomposition. The large set of (n,G,λ)‐GD is denoted by (n,G,λ)‐LGD. In this work, we obtain the existence spectrum of (n,P₃,λ)‐LGD. © 2002 Wiley Periodicals, Inc. J Combin Designs 10: 151–159, 2002; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/jcd.10008     

Integral Inequalities and Global Solutions of Semilinear Evolution Equations

A criterion for the nonexplosion of solutions to semilinear evolution equations on Banach spaces is proved. The result is obtained by applying a modification of the Bihari type inequality to the case of a weakly singular nonlinear integral inequality.     

具有模糊数的模糊多目标群体决策优选模型与方法

多目标群体决策问题是运筹学的一个重要研究领域，目前已经提出了一些有效的决策方法。但对目标值和权重均为模糊数的模糊多目标群体决策问题却研究不多，本对此类模糊多目标群体决策问题进行了探讨，利用相对正理想方案与相对负理想方案概念定义了相对差异距离，进而建立了模糊多目标群体决策优选模型与方法，并通过战役决心方案的评价说明了该方法是可行、有效的，可作为军事决策与决策支持系统的备选方法。     

Simple, expedient methods for the determination of water and electrolyte contents of cellulose solvent systems

The concentrations of water, W, and electrolytes present in solutions of LiCl in N,N-dimethylacetamide, LiCl/DMAc, and of tetrabutylammonium fluoride. x-hydrate in DMSO, TBAF.xW/DMSO can be accurately and expediently determined by three independent methods, UV–vis, FTIR and EMF measurement. The first relies on the use of solvatochromic probes whose spectra are sensitive to solution composition. It is applicable to W/LiCl/DMAc solutions but not to TBAF.xW/DMSO, because the charge-transfer complex bands of the probes are suppressed by strong interactions with the latter electrolyte. Integration of ν_OH band of water may be employed in order to determine [W], hence [electrolyte] by weight difference. EMF measurement uses ion-selective electrodes in order to determine [electrolyte], hence [W] by weight difference. Results of the latter method were in excellent agreement with those of FTIR. The reason for the failure of Karl Fischer titration is addressed, and the relevance of the results obtained to functionalization of cellulose under homogenous solution conditions is briefly commented on.     

Patterns and Symmetries in the Visual Cortex and in Natural Images

As borders between different entities, lines are an important element of natural images. Indeed, the neurons of the mammalian visual cortex are tuned to respond best to lines of a given orientation. This preferred orientation varies continuously across most of the cortex, but also has vortex-like singularities known as pinwheels. In attempting to describe such patterns of orientation preference, we are led to consider underlying rotation symmetries: Oriented segments in natural images tend to be collinear; neurons are more likely to be connected if their preferred orientations are aligned to their topographic separation. These are indications of a reduced symmetry requiring joint rotations of both orientation preference and the underlying topography. This is verified by direct statistical tests in both natural images and in cortical maps. Using the statistics of natural scenes we construct filters that are best suited to extracting information from such images, and find qualitative similarities to mammalian vision. PACS84.35+i 89.70.+c 87.57.Nk     

A dynamical test of phase transition order: New things in old places or old wine in new bottles

We first discuss nonlinear aspects of phase transition theory applied to a particular liquid crystal phase transition. A simple derivation is given to show how two coupled Goldstone modes (one appearing as gauge fluctuations of the ordered phase) can force a phase transition, against all expectations, to take place discontinuously (theory of Halperin, Lubensky, and Ma)-but the discontinuity may be immeasurably small. Then, we describe a new dynamical test of phase transition order, developed by Cladiset al., that turns out to be more sensitive than x-ray diffraction and adiabatic calorimetry. Quantitative data found by this new method are in excellent agreement with the measurements of adiabatic calorimetry and x-ray diffraction as well as expectations implicit in the predictions of HLM.This is the text of an after-banquet talk given at the CNLS Workshop on the Dynamics of Concentrated Systems.     

Rehabilitation of the Gauss-Jordan algorithm

Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.