Gradient based and least squares based iterative algorithms for matrix equations AXB + CXD = F

This paper develops a gradient based and a least squares based iterative algorithms for solving matrix equation AXB + CX^TD = F. The basic idea is to decompose the matrix equation (system) under consideration into two subsystems by applying the hierarchical identification principle and to derive the iterative algorithms by extending the iterative methods for solving Ax = b and AXB = F. The analysis shows that when the matrix equation has a unique solution (under the sense of least squares), the iterative solution converges to the exact solution for any initial values. A numerical example verifies the proposed theorems.     

Extending compact topologies to compact Hausdorff topologies in ZF

Within the framework of Zermelo-Fraenkel set theory ZF, we investigate the set-theoretical strength of the following statements:

(1)

For every family(Ai)_i∈Iof sets there exists a family(Ti)_i∈Isuch that for everyi∈I(Ai,Ti)is a compactT₂space.

(2)

For every family(Ai)_i∈Iof sets there exists a family(Ti)_i∈Isuch that for everyi∈I(Ai,Ti)is a compact, scattered, T₂space.

(3)

For every set X, every compactR₁topology (itsT₀-reflection isT₂) on X can be enlarged to a compactT₂topology.

We show:

(a)

(1) implies every infinite set can be split into two infinite sets.

(b)

(2) iff AC.

(c)

(3) and “there exists a free ultrafilter” iff AC.

We also show that if the topology of certain compact T₁ spaces can be enlarged to a compact T₂ topology then (1) holds true. But in general, compact T₁ topologies do not extend to compact T₂ ones.     

Exponentiation in V-categories

For a Heyting algebra V which, as a category, is monoidal closed, we obtain characterizations of exponentiable objects and morphisms in the category of V-categories and apply them to some well-known examples. In the case these characterizations of exponentiable morphisms and objects in the categories (P)Met of (pre)metric spaces and non-expansive maps show in particular that exponentiable metric spaces are exactly the almost convex metric spaces, while exponentiable complete metric spaces are the complete totally convex ones.     

Squarefree P-modules and the cd-index

In this paper, we introduce a new algebraic concept, which we call squarefree P-modules. This concept is inspired from Karu's proof of the non-negativity of the cd-indices of Gorenstein* posets, and supplies a way to study cd-indices from the viewpoint of commutative algebra. Indeed, by using the theory of squarefree P-modules, we give several new algebraic and combinatorial results on CW-posets. First, we define an analogue of the cd-index for any CW-poset and prove its non-negativity when a CW-poset is Cohen–Macaulay. This result proves that the h-vector of the barycentric subdivision of a Cohen–Macaulay regular CW-complex is unimodal. Second, we prove that the Stanley–Reisner ring of the barycentric subdivision of an odd dimensional Cohen–Macaulay polyhedral complex has the weak Lefschetz property. Third, we obtain sharp upper bounds of the cd-indices of Gorenstein* posets for a fixed rank generating function.     

A Construction of Compactly Supported Biorthogonal Scaling Vectors and Multiwavelets on R

In Kessler (Appl. Comput. Harmonic Anal.9 (2000), 146-165), a construction was given for a class of orthogonal compactly supported scaling vectors on R², called short scaling vectors, and their associated multiwavelets. The span of the translates of the scaling functions along a triangular lattice includes continuous piecewise linear functions on the lattice, although the scaling functions are fractal interpolation functions and possibly nondifferentiable. In this paper, a similar construction will be used to create biorthogonal scaling vectors and their associated multiwavelets. The additional freedom will allow for one of the dual spaces to consist entirely of the continuous piecewise linear functions on a uniform subdivision of the original triangular lattice.     

Large time behavior to the system of incompressible non-Newtonian fluids in R

In this paper, the authors study the large time behavior for the weak solutions to a class system of the incompressible non-Newtonian fluids in R². It is proved that the weak solutions decay in L² norm at (1+t)^−1/2 and the estimate for the decay rate is sharp in the sense that it coincides with the decay rate of a solution to the heat equation.     

Two, more readily computable equivariant Nielsen numbers I. Nielsen theory for M-ads

In this, the first of two papers outlining a Nielsen theory for “two, more readily computable equivariant numbers”, we define and study two Nielsen type numbers N(f,k;X−{Xν}_ν∈M) and N(f,k;X,{Xν}_ν∈M), where f and k are M-ad maps. While a Nielsen theory of M-ads is of interest in its own right, our main motivation lies in the fact that maps of M-ads accurately mirror one of two fundamental structures of equivariant maps. Being simpler however, M-ad Nielsen numbers are easier to study and to compute than equivariant Nielsen numbers. In the sequel, we show our M-ad numbers can be used to form both upper and lower bounds on their equivariant counterparts.The numbers N(f,k;X−{Xν}_ν∈M) and N(f,k;X,{Xν}_ν∈M), generalize the generalizations to coincidences, of Zhao's Nielsen number on the complement N(f;X−A), respectively Schirmer's relative Nielsen number N(f;X,A). Our generalizations are from the category of pairs, to the category of M-ads. The new numbers are lower bounds for the number of coincidence points of all maps f^′ and k^′ which are homotopic as maps ofM-ads to f, respectively k firstly on the complement of the union of the subspaces Xν in the domain M-ad X, and secondly on all of X. The second number is shown to be greater than or equal to a sum of the first of our numbers. Conditions are given which allow for both equality, and Möbius inversion. Finally we show that the fixed point case of our second number generalizes Schirmer's triad Nielsen number N(f;X₁∪X₂).Our work is very different from what at first sight appears to be similar partial results due to P. Wong. The differences, while in some sense subtle in terms of definition, are profound in terms of commutability. In order to work in a variety of both fixed point and coincidence points contexts, we introduce in this first paper and extend in the second, the concept of an essentiality on a topological category. This allows us to give computational theorems within this diversity. Finally we include an introduction to both papers here.     

A geometric realization of quantum groups of type D

We geometrize quantum groups of type D in the spirit of Beilinson et al. (1990) [1].     

Topological sums and products in ZF-set theory

In ZF, i.e., Zermelo-Fraenkel set theory without the axiom of choice, the category Top of topological spaces and continuous maps is well-behaved. In particular, Top has sums (=coproducts) and products. However, it may happen that for families (Xi)_i∈I and (Yi)_i∈I with the property that each Xi is homeomorphic to the corresponding Yi neither their sums _⊕i∈IXi and _⊕i∈IYi nor their products _∏i∈IXi and _∏i∈IYi are homeomorphic. It will be shown that the axiom of choice is not only sufficient but also necessary to rectify this defect.     

Polynomials satisfied by two linked matrices

Polynomials in two variables, evaluated at A and with A being a square complex matrix and being its transform belonging to the set {A⁼, A^†, A^∗}, in which A⁼, A^†, and A^∗ denote, respectively, any reflexive generalized inverse, the Moore-Penrose inverse, and the conjugate transpose of A, are considered. An essential role, in characterizing when such polynomials are satisfied by two matrices linked as above, is played by the condition that the column space of A is the column space of . The results given unify a number of prior, isolated results.     

A Test Set for stiff Initial Value Problem Solvers in the open source software R: Package deTestSet

In this paper we present the R package deTestSet that includes challenging test problems written as ordinary differential equations (ODEs), differential algebraic equations (DAEs) of index up to 3 and implicit differential equations (IDEs). In addition it includes 6 new codes to solve initial value problems (IVPs). The R package is derived from the Test Set for Initial Value Problem Solvers available at http://www.dm.uniba.it/~testset which includes documentation of the test problems, experimental results from a number of proven solvers, and Fortran subroutines providing a common interface to the defining problem functions. Many of these facilities are now available in the R package deTestSet, which comprises an R interface to the test problems and to most of the Fortran solvers. The package deTestSet is free software which is distributed under the GNU General Public License, as part of the R open source software project.     

The characterization of a class of multivariate MRA and semi-orthogonal parseval frame wavelets

Let A be a d × d expansive matrix with ∣detA∣ = 2. This paper addresses Parseval frame wavelets (PFWs) in the setting of reducing subspaces of L²(R^d). We prove that all semi-orthogonal PFWs (semi-orthogonal MRA PFWs) are precisely the ones with their dimension functions being non-negative integer-valued (0 or 1). We also characterize all MRA PFWs. Some examples are provided.     

Two local operators and the BLUE

For a given matrix A, a matrix P such that PA = A is said to be a local identity, and such that P²A = PA is said to be a local idempotent. In the paper, some simple properties of such operators are presented. Their relation to the best linear unbiased estimation in the general Gauss-Markov model is also demonstrated.     

Bounds for the error of linear systems of equations using the theory of moments

Consider the system, of linear equations Ax = b where A is an n × n real symmetric, positive definite matrix and b is a known vector. Suppose we are given an approximation to x, ξ, and we wish to determine upper and lower bounds for ∥ x − ξ ∥ where ∥ ··· ∥ indicates the euclidean norm. Given the sequence of vectors {r_i}_{i^k = 0}, where r_i = Ar_{i − 1} and r₀ = b − Aξ, it is shown how to construct a sequence of upper and lower bounds for ∥ x − ξ ∥ using the theory of moments.     

On commutativity of projectors

The purpose of this paper is to revisit two problems discussed previously in the literature, both related to the commutativity property P₁P₂ = P₂P₁, where P₁ and P₂ denote projectors (i.e., idempotent matrices). The first problem was considered by Baksalary et al. [J.K. Baksalary, O.M. Baksalary, T. Szulc, A property of orthogonal projectors, Linear Algebra Appl. 354 (2002) 35-39], who have shown that if P₁ and P₂ are orthogonal projectors (i.e., Hermitian idempotent matrices), then in all nontrivial cases a product of any length having P₁ and P₂ as its factors occurring alternately is equal to another such product if and only if P₁ and P₂ commute. In the present paper a generalization of this result is proposed and validity of the equivalence between commutativity property and any equality involving two linear combinations of two any length products having orthogonal projectors P₁ and P₂ as their factors occurring alternately is investigated. The second problem discussed in this paper concerns specific generalized inverses of the sum P₁ + P₂ and the difference P₁ − P₂ of (not necessary orthogonal) commuting projectors P₁ and P₂. The results obtained supplement those provided in Section 4 of Baksalary and Baksalary [J.K. Baksalary, O.M. Baksalary, Commutativity of projectors, Linear Algebra Appl. 341 (2002) 129-142].     

Strong linear preservers of rank reverse permutability on triangular matrices

Let F be a field with ∣F∣ > 2 and T_n(F) be the set of all n × n upper triangular matrices, where n ? 2. Let k ? 2 be a given integer. A k-tuple of matrices A₁, …, A_k ∈ T_n(F) is called rank reverse permutable if rank(A₁ A₂ ? A_k) = rank(A_k A_k−1 ? A₁). We characterize the linear maps on T_n(F) that strongly preserve the set of rank reverse permutable matrix k-tuples.     

Orbifolds and groupoids

We define a 2-category structure (Pre-Orb) on the category of reduced complex orbifold atlases. We construct a 2-functor F from (Pre-Orb) to the 2-category (Grp) of proper étale effective groupoid objects over the complex manifolds. Both on (Pre-Orb) and (Grp) there are natural equivalence relations on objects: (a natural extension of) equivalence of orbifold atlases on (Pre-Orb) and Morita equivalence in (Grp). We prove that F induces a bijection between the equivalence classes of its source and target.