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.     

Expected Reflection Distance in G(r, 1, n) after a Fixed Number of Reflections

Extending to r > 1 a formula of the authors, we compute the expected reflection distance of a product of t random reflections in the complex reflection group G(r, 1, n). The result relies on an explicit decomposition of the reflection distance function into irreducible G(r, 1, n)-characters and on the eigenvalues of certain adjacency matrices.Received December 8, 2003     

The Interaction of Shocks with Dispersive Waves II. Incompressible-Integrable Limit

This is the second in a two-part series of articles in which we analyze a system similar in structure to the well-known Zakharov equations from weak plasma turbulence theory, but with a nonlinear conservation equation allowing finite time shock formation. In this article we analyze the incompressible limit in which the shock speed is large compared to the underlying group velocity of the dispersive wave (a situation typically encountered in applications). After presenting some exact solutions of the full system, a multiscale perturbation method is used to resolve several basic wave interactions. The analysis breaks down into two categories: the nonlinear limit and the linear limit, corresponding to the form of the equations when the group velocity to shock speed ratio, denoted by ε, is zero. The former case is an integrable limit in which the model reduces to the cubic nonlinear Schrödinger equation governing the dispersive wave envelope. We focus on the interaction of a “fast” shock wave and a single hump soliton. In the latter case, the ε=0 problem reduces to the linear Schrödinger equation, and the focus is on a fast shock interacting with a dispersive wave whose amplitude is cusped and exponentially decaying. To motivate the time scales and structure of the shock-dispersive wave interactions at lowest orders, we first analyze a simpler system of ordinary differential equations structurally similar to the original system. Then we return to the fully coupled partial differential equations and develop a multiscale asymptotic method to derive the effective leading-order shock equations and the leading-order modulation equations governing the phase and amplitude of the dispersive wave envelope. The leading-order interaction equations admit a fairly complete analysis based on characteristic methods. Conditions are derived in which: (a) the shock passes through the soliton, (b) the shock is completely blocked by the soliton, or (c) the shock reverses direction. In the linear limit, a phenomenon is described in which the dispersive wave induces the formation of a second, transient shock front in the rapidly moving hyperbolic wave. In all cases, we can characterize the long-time dynamics of the shock. The influence of the shock on the dispersive wave is manifested, to leading order, in the generalized frequency of the dispersive wave: the fast-time part of the frequency is the shock wave itself. Hence, the frequency undergoes a sudden jump across the shock layer.In the last section, a sequence of numerical experiments depicting some of the interesting interactions predicted by the analysis is performed on the leading-order shock equations.     

From monomeric to polymeric ferroelectric liquid crystals A comparative study of ferroelectric properties

Two new ferroelectric oligosiloxanes, a cyclic tetramer and a twin, have been synthesized. By a comparative study with their corresponding monomer and side chain polysiloxanes, the influence of oligo- and polymerization on the liquid crystalline and ferroelectric properties have been investigated. Polymerization leads to a stabilization of LC phases through increase of the clearing temperatures and suppression of crystallization. Oligomerization also leads to mesophase broadening, but, due to the low degree of polymerization, the effect is inferior to the linear polysiloxanes. The low viscosity of the oligosiloxanes ensures response times in the microsecond region, thus being comparable with their monomer and conventional LMWFLCs. It is found that polymerization increases the spontaneous polarization P_s. This is attributed to the density increase after polymerization, enhancing the inter-mesogenic interactions. The collective and local dynamics of the OFLCs are influenced differently with respect to their molecular structures. Each oligomer is already a good model for its corresponding polymer concerning the soft mode dynamics. For the local β-relaxation a similar temperature dependence of the relaxation times τ for the cyclic tetramer and for the side chain polysiloxanes is observed. The long axial rotation of the twin, having a very efficient decoupling, is significantly faster, thus resembling the monomer.     

Non-destructive analysis of archeological samples by energy-dispersive X-ray fluorescence

Summary A non-destructive method is described for the determination of major and minor constituents in archeological specimens by energy-dispersive X-ray fluorescence. Homogeneity tests are made by measuring at various sites of the specimen. In the same way, mean values are obtained for inhomogeneous specimen without taking samples. For calibration, powder standards are used. In case of the determination of elements with numbers up to 14 (Si) a vacuum chamber is used and the dimensions of the specimens are limited by the dimensions of that vacuum chamber, whereas for the determination of elements from K up to U specimens of any size, form or weight are suitable.

---

Zerstörungsfreie Analyse von archäologischen Proben mit Hilfe der Energie-dispersiven Röntgenfluorescenzanalyse

Zusammenfassung Eine zerstörungsfreie Methode für die Bestimmung von Haupt- und Nebenbestandteilen in archäologischen Proben mit Hilfe der Energie-dispersiven Röntgenfluorescenzanalyse wird beschrieben. Für Homogenitätstests wird an mehreren Stellen der Probe gemessen. In der gleichen Weise werden für inhomogene Proben Mittelwerte erhalten ohne Probenahme. Für die Eichung werden Pulverstandards verwendet. Im Falle der Bestimmung von Elementen mit Ordnungszahlen bis 14 (Si) wird eine Vakuumkammer eingesetzt, und die Dimensionen der Proben sind durch die Dimensionen dieser Vakuumkammer begrenzt, während für die Bestimmung der Elemente K bis U Proben jeder Größe, jeder Form oder jeden Gewichts verwendbar sind.