Error analysis of a Monte Carlo algorithm for computing bilinear forms of matrix powers

In this paper we present error analysis for a Monte Carlo algorithm for evaluating bilinear forms of matrix powers. An almost Optimal Monte Carlo (MAO) algorithm for solving this problem is formulated. Results for the structure of the probability error are presented and the construction of robust and interpolation Monte Carlo algorithms are discussed. Results are presented comparing the performance of the Monte Carlo algorithm with that of a corresponding deterministic algorithm. The two algorithms are tested on a well balanced matrix and then the effects of perturbing this matrix, by small and large amounts, is studied.

Comparison of the computational cost of a Monte Carlo and deterministic algorithm for computing bilinear forms of matrix powers

In this paper we consider bilinear forms of matrix polynomials and show that these polynomials can be used to construct solutions for the problems of solving systems of linear algebraic equations, matrix inversion and finding extremal eigenvalues. An almost Optimal Monte Carlo (MAO) algorithm for computing bilinear forms of matrix polynomials is presented. Results for the computational costs of a balanced algorithm for computing the bilinear form of a matrix power is presented, i.e., an algorithm for which probability and systematic errors are of the same order, and this is compared with the computational cost for a corresponding deterministic method.

Miniversal deformations of matrices of bilinear forms

Arnold [V.I. Arnold, On matrices depending on parameters, Russian Math. Surveys 26 (2) (1971) 29-43] constructed miniversal deformations of square complex matrices under similarity; that is, a simple normal form to which not only a given square matrix A but all matrices B close to it can be reduced by similarity transformations that smoothly depend on the entries of B. We construct miniversal deformations of matrices under congruence. (C) 2011 Elsevier Inc. All rights reserved.

Quadratic and Symmetric Bilinear Forms on Modules with Unique Base Over a Semiring

27 pages

Quadratic and symmetric bilinear forms on modules with unique base over a semiring

Peer reviewed

Tridiagonal canonical matrices of bilinear or sesquilinear forms and of pairs of symmetric, skew-symmetric, or Hermitian forms

Tridiagonal canonical forms of square matrices under congruence or *congruence, pairs of symmetric or skew-symmetric matrices under congruence, and pairs of Hermitian matrices under *congruence are given over an algebraically closed field of characteristic different from 2. (C) 2008 Elsevier Inc. All rights reserved.

Monte Carlo numerical treatment of large linear algebra problems

In this paper we deal with performance analysis of Monte Carlo algorithm for large linear algebra problems. We consider applicability and efficiency of the Markov chain Monte Carlo for large problems, i.e., problems involving matrices with a number of non-zero elements ranging between one million and one billion. We are concentrating on analysis of the almost Optimal Monte Carlo (MAO) algorithm for evaluating bilinear forms of matrix powers since they form the so-called Krylov subspaces. Results are presented comparing the performance of the Robust and Non-robust Monte Carlo algorithms. The algorithms are tested on large dense matrices as well as on large unstructured sparse matrices.

Is the Helmholtz equation really sign-indefinite?

The usual variational (or weak) formulations of the Helmholtz equation are sign-indefinite in the sense that the bilinear forms cannot be bounded below by a positive multiple of the appropriate norm squared. This is often for a good reason, since in bounded domains under certain boundary conditions the solution of the Helmholtz equation is not unique at wavenumbers that correspond to eigenvalues of the Laplacian, and thus the variational problem cannot be sign-definite. However, even in cases where the solution is unique for all wavenumbers, the standard variational formulations of the Helmholtz equation are still indefinite when the wavenumber is large. This indefiniteness has implications for both the analysis and the practical implementation of finite element methods. In this paper we introduce new sign-definite (also called coercive or elliptic) formulations of the Helmholtz equation posed in either the interior of a star-shaped domain with impedance boundary conditions, or the exterior of a star-shaped domain with Dirichlet boundary conditions. Like the standard variational formulations, these new formulations arise just by multiplying the Helmholtz equation by particular test functions and integrating by parts.

On a formula for the spectral flow and its applications

We consider a continuous path of bounded symmetric Fredholm bilinear forms with arbitrary endpoints on a real Hilbert space, and we prove a formula that gives the spectral flow of the path in terms of the spectral flow of the restriction to a finite codimensional closed subspace. We also discuss the case of restrictions to a continuous path of finite codimensional closed subspaces. As an application of the formula, we introduce the notion of spectral flow for a periodic semi-Riemannian geodesic, and we compute its value in terms of the Maslov index. (C) 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

An algebraic theory for generalized Jordan chains and partial signatures in the Lagrangian Grassmannian

We discuss an algebraic theory for generalized Jordan chains and partial signatures, that are invariants associated to sequences of symmetric bilinear forms on a vector space. We introduce an intrinsic notion of partial signatures in the Lagrangian Grassmannian of a symplectic space that does not use local coordinates, and we give a formula for the Maslov index of arbitrary real analytic paths in terms of partial signatures.

Classificação de cônicas e quádricas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

On Finite Element Methods for 2nd order (semi–) periodic Eigenvalue Problems

We deal with a class of elliptic eigenvalue problems (EVPs) on a rectangle Ω ⊂ R^2 , with periodic or semi–periodic boundary conditions (BCs) on ∂Ω. First, for both types of EVPs, we pass to a proper variational formulation which is shown to fit into the general framework of abstract EVPs for symmetric, bounded, strongly coercive bilinear forms in Hilbert spaces, see, e.g., [13, §6.2]. Next, we consider finite element methods (FEMs) without and with numerical quadrature. The aim of the paper is to show that well–known error estimates, established for the finite element approximation of elliptic EVPs with classical BCs, hold for the present types of EVPs too. Some attention is also paid to the computational aspects of the resulting algebraic EVP. Finally, the analysis is illustrated by two non-trivial numerical examples, the exact eigenpairs of which can be determined.

Chronic Telogen Effluvium And Female Pattern Hair Loss Are Separate And Distinct Forms Of Alopecia: A Histomorphometric And Immunohistochemical Analysis.

Chronic telogen effluvium (CTE), a poorly understood condition, can be confused with or may be a prodrome to female pattern hair loss (FPHL). The pathogenesis of both is related to follicle cycle shortening and possibly to blood supply changes. To analyze a number of histomorphometric and immunohistochemical findings through vascular endothelial growth factor (VEGF), Ki-67, and CD31 immunostaining in scalp biopsies of 20 patients with CTE, 17 patients with mild FPHL and 9 controls. Ki-67 index and VEGF optical density were analyzed at the follicular outer sheath using ImageJ software. CD31 microvessel density was assessed by a Chalkley grid. Significant follicle miniaturization and higher density of nonanagen follicles were found in FPHL, compared with patients with CTE and controls. Ki-67+ index correlated positively with FPHL histological features. The FPHL group showed the highest VEGF optical density, followed by the CTE and control groups. No differences were found in CD31 microvessel density between the three groups. Histomorphometric results establish CTE as a distinct disorder, separate from FPHL from its outset. Its pathogenic mechanisms are also distinct. These findings support the proposed mechanism of 'immediate telogen release' for CTE, leading to cycle synchronization. For FPHL, accelerated anagen follicular mitotic rates and, thus, higher Ki-67 and VEGF values, would leave less time for differentiation, resulting in hair miniaturization.

Use Of Experimental Design To Optimize A Triple-potential Waveform To Develop A Method For The Determination Of Streptomycin And Dihydrostreptomycin In Pharmaceutical Veterinary Dosage Forms By Hplc-pad.

An HPLC-PAD method using a gold working electrode and a triple-potential waveform was developed for the simultaneous determination of streptomycin and dihydrostreptomycin in veterinary drugs. Glucose was used as the internal standard, and the triple-potential waveform was optimized using a factorial and a central composite design. The optimum potentials were as follows: amperometric detection, E1=-0.15V; cleaning potential, E2=+0.85V; and reactivation of the electrode surface, E3=-0.65V. For the separation of the aminoglycosides and the internal standard of glucose, a CarboPac™ PA1 anion exchange column was used together with a mobile phase consisting of a 0.070 mol L(-1) sodium hydroxide solution in the isocratic elution mode with a flow rate of 0.8 mL min(-1). The method was validated and applied to the determination of streptomycin and dihydrostreptomycin in veterinary formulations (injection, suspension and ointment) without any previous sample pretreatment, except for the ointments, for which a liquid-liquid extraction was required before HPLC-PAD analysis. The method showed adequate selectivity, with an accuracy of 98-107% and a precision of less than 3.9%.

Interleukin-22 forms dimers that are recognized by two interleukin-22R1 receptor chains

Interleukin-22 (IL-22) is a class 2 cytokine whose primary structure is similar to that of interleukin 10 (IL-10) and interferon-gamma (IFN-gamma). IL-22 induction during acute phase immune response indicates its involvement in mechanisms of inflammation. Structurally different from IL-10 and a number of other members of IL-10 family, which form intertwined inseparable V-shaped dimers of two identical polypeptide chains, a single polypeptide chain of IL-22 folds on itself in a relatively globular structure. Here we present evidence, based on native gel electrophoresis, glutaraldehyde cross-linking, dynamic light scattering, and small angle x-ray scattering experiments, that human IL-22 forms dimers and tetramers in solution under protein concentrations assessable by these experiments. Unexpectedly, low-resolution molecular shape of IL-22 dimers is strikingly similar to that of IL-10 and other intertwined cytokine dimeric forms. Furthermore, we determine an ab initio molecular shape of the IL-22/IL-22R1 complex which reveals the V-shaped IL-22 dimer interacting with two cognate IL-22R1 molecules. Based on this collective evidence, we argue that dimerization might be a common mechanism of all class 2 cytokines for the molecular recognition with their respective membrane receptor. We also speculate that the IL-22 tetramer formation could represent a way to store the cytokine in nonactive form at high concentrations that could be readily converted into functionally active monomers and dimers upon interaction with the cognate cellular receptors.