Approximate factorization constraint preconditioners for saddle-point matrices

We consider the application of the conjugate gradient method to the solution of large, symmetric indefinite linear systems. Special emphasis is put on the use of constraint preconditioners and a new factorization that can reduce the number of flops required by the preconditioning step. Results concerning the eigenvalues of the preconditioned matrix and its minimum polynomial are given. Numerical experiments validate these conclusions.

Implicit-factorization preconditioning and iterative solvers for regularized saddle-point systems

We consider conjugate-gradient like methods for solving block symmetric indefinite linear systems that arise from saddle-point problems or, in particular, regularizations thereof. Such methods require preconditioners that preserve certain sub-blocks from the original systems but allow considerable flexibility for the remaining blocks. We construct a number of families of implicit factorizations that are capable of reproducing the required sub-blocks and (some) of the remainder. These generalize known implicit factorizations for the unregularized case. Improved eigenvalue clustering is possible if additionally some of the noncrucial blocks are reproduced. Numerical experiments confirm that these implicit-factorization preconditioners can be very effective in practice.

A benchmark method for global optimization problems in structure determination from powder diffraction data

Quasi-Newton-Raphson minimization and conjugate gradient minimization have been used to solve the crystal structures of famotidine form B and capsaicin from X-ray powder diffraction data and characterize the chi(2) agreement surfaces. One million quasi-Newton-Raphson minimizations found the famotidine global minimum with a frequency of ca 1 in 5000 and the capsaicin global minimum with a frequency of ca 1 in 10 000. These results, which are corroborated by conjugate gradient minimization, demonstrate the existence of numerous pathways from some of the highest points on these chi(2) agreement surfaces to the respective global minima, which are passable using only downhill moves. This important observation has significant ramifications for the development of improved structure determination algorithms.

Simple adaptive momentum: new algorithm for training multilayer perceptrons

The speed of convergence while training is an important consideration in the use of neural nets. The authors outline a new training algorithm which reduces both the number of iterations and training time required for convergence of multilayer perceptrons, compared to standard back-propagation and conjugate gradient descent algorithms.

Iterative Methods for Equality-Constrained Least Squares Problems

We consider the linear equality-constrained least squares problem (LSE) of minimizing ${\|c - Gx\|}_2 $, subject to the constraint $Ex = p$. A preconditioned conjugate gradient method is applied to the Kuhn–Tucker equations associated with the LSE problem. We show that our method is well suited for structural optimization problems in reliability analysis and optimal design. Numerical tests are performed on an Alliant FX/8 multiprocessor and a Cray-X-MP using some practical structural analysis data.

The counter-propagating Rossby-wave perspective on baroclinic instability. II: Application to the Charney model

The constant-density Charney model describes the simplest unstable basic state with a planetary-vorticity gradient, which is uniform and positive, and baroclinicity that is manifest as a negative contribution to the potential-vorticity (PV) gradient at the ground and positive vertical wind shear. Together, these ingredients satisfy the necessary conditions for baroclinic instability. In Part I it was shown how baroclinic growth on a general zonal basic state can be viewed as the interaction of pairs of ‘counter-propagating Rossby waves’ (CRWs) that can be constructed from a growing normal mode and its decaying complex conjugate. In this paper the normal-mode solutions for the Charney model are studied from the CRW perspective. Clear parallels can be drawn between the most unstable modes of the Charney model and the Eady model, in which the CRWs can be derived independently of the normal modes. However, the dispersion curves for the two models are very different; the Eady model has a short-wave cut-off, while the Charney model is unstable at short wavelengths. Beyond its maximum growth rate the Charney model has a neutral point at finite wavelength (r=1). Thereafter follows a succession of unstable branches, each with weaker growth than the last, separated by neutral points at integer r—the so-called ‘Green branches’. A separate branch of westward-propagating neutral modes also originates from each neutral point. By approximating the lower CRW as a Rossby edge wave and the upper CRW structure as a single PV peak with a spread proportional to the Rossby scale height, the main features of the ‘Charney branch’ (0gradients, are arguably more applicable to the atmosphere than modes of the Eady model where the positive PV gradient exists only at the tropopause

The response of a uniform horizontal temperature gradient to heating

The response of a uniform horizontal temperature gradient to prescribed fixed heating is calculated in the context of an extended version of surface quasigeostrophic dynamics. It is found that for zero mean surface flow and weak cross-gradient structure the prescribed heating induces a mean temperature anomaly proportional to the spatial Hilbert transform of the heating. The interior potential vorticity generated by the heating enhances this surface response. The time-varying part is independent of the heating and satisfies the usual linearized surface quasigeostrophic dynamics. It is shown that the surface temperature tendency is a spatial Hilbert transform of the temperature anomaly itself. It then follows that the temperature anomaly is periodically modulated with a frequency proportional to the vertical wind shear. A strong local bound on wave energy is also found. Reanalysis diagnostics are presented that indicate consistency with key findings from this theory.

Projected gradient approach to the numerical solution of the SCoTLASS

The SCoTLASS problem-principal component analysis modified so that the components satisfy the Least Absolute Shrinkage and Selection Operator (LASSO) constraint-is reformulated as a dynamical system on the unit sphere. The LASSO inequality constraint is tackled by exterior penalty function. A globally convergent algorithm is developed based on the projected gradient approach. The algorithm is illustrated numerically and discussed on a well-known data set. (c) 2004 Elsevier B.V. All rights reserved.

Structural study of BSA/poly(ethylene glycol) lipid conjugate complexes

In this work we report the structural characteristics of bovine serum albumin/poly(ethylene glycol) lipid conjugate (BSA/PEG(2000)-PE) complexes under physiological conditions (37 degrees C and pH 7.4) for particular fractions of BSA to PEG-lipid concentration, CBSA/C-PEG2000-PE. Ultraviolet fluorescence spectroscopy (UV) results shown that PEG(2000)-PE is associated to BSA, leading to;protein unfolding for fixed C-BSA = 0.01 wt % and variable C-PEG2000-PE = 0.0015-0.6 wt %. Tryptophan groups on the BSA surface are in contact with the PEG-lipid at C-PEG2000-PE = 0.0015 wt %, while they are exposed to water at C-PEG2000-PE (>)0.0015 wt %. Dynamic and static light scattering (DLS and SLS) and small-angle neutron scattering (SANS) point out the existence of individual BSAIPEG-lipid complexes in the system for fixed C-BSA = 1 wt % and variable C-PEG2000-PE = 0.15-2 wt %. DLS shows that there is only one BSA molecule per protein/PEG-lipid complex, while SLS shows that the PEG-lipid associates to the BSA without promoting aggregation between adjacent protein/ polymer-lipid conjugate complexes. SANS was used to show that BSA/PEG(2000)-PE complexes adopt an oblate ellipsoidal shape. Partially unfolded BSA is contained in the core of the oblate ellipsoid, which is surrounded by an external shell containing the PEG(2000)-PE.

Interactions of KLVFF-PEG peptide conjugate with fibrinogen in neutral aqueous solutions

In this work we report on the interaction of KLVFF-PEG with fibrinogen (Fbg) in neutral aqueous solutions at 20 degrees C, for particular ratios of KLVFF-PEG to Fbg concentration, Delta = CKLVFF-PEG/C-Fbg- Our results show the formation of Fbg/KLVFF-PEG complexes for Delta > 0, such that there is not an extended network of complexes throughout the solution. In addition, cleaved protein and Fbg dimers are identified in the solution for Delta >= 0. There is a dramatic change in the tertiary structure of the Fbg upon KLVFF-PEG binding, although the KLVFF-PEG binds to the Fbg without affecting the secondary structure elements of the glycoprotein.

Nematic and columnar ordering of a PEG-peptide conjugate in aqueous solution

The self-assembly in aqueous solution of a PEG-peptide conjugate is studied by spectroscopy, electron microscopy, rheology and small-angle Xray and neutron scattering (SAXS and SANS). The peptide fragment, FFKLVFF is based on fragment KLVFF of the amyloid beta-peptide, A beta(16-20), extended by two hydrophobic phenylalanine units. This is conjugated to PEG which confers water solubility and leads to distinct self-assembled structures. Small-angle scattering reveals the formation of cylindrical fibrils comprising a peptide core and PEG corona. This constrained structure leads to a model parallel beta-sheet self-assembled structure with a radial arrangement of beta sheets. Oil increasing concentration, successively nematic and hexagonal columnar phases are formed. The flow-induced alignment of both structures was studied in situ by SANS using a Couette cell. Shear-induced alignment is responsible for the shear thinning behaviour observed by dynamic shear rheometry. Incomplete recovery of moduli after cessation of shear is consistent with the observation from SANS of retained orientation in the sample.

Orientational ordering in the nematic phase of a polyethylene glycol-peptide conjugate in aqueous solution

The orientational ordering of the nematic phase of a polyethylene glycol (PEG)-peptide block copolymer in aqueous solution is probed by small-angle neutron scattering (SANS), with the sample subjected to steady shear in a Couette cell. The PEG-peptide conjugate forms fibrils that behave as semiflexible rodlike chains. The orientational order parameters (P) over bar (2) and (P) over bar (4) are obtained by modeling the data using a series expansion approach to the form factor of uniform cylinders. The method used is independent of assumptions on the form of the singlet orientational distribution function. Good agreement with the anisotropic two-dimensional SANS patterns is obtained. The results show shear alignment starting at very low shear rates, and the orientational order parameters reach a plateau at higher shear rates with a pseudologarithmic dependence on shear rate. The most probable distribution functions correspond to fibrils parallel to the flow direction under shear, but a sample at rest shows a bimodal distribution with some of the rodlike peptide fibrils oriented perpendicular to the flow direction.