A nested Krylov subspace method to compute the sign function of large complex matrices

We present an acceleration of the well-established Krylov-Ritz methods to compute the sign function of large complex matrices, as needed in lattice QCD simulations involving the overlap Dirac operator at both zero and nonzero baryon density. Krylov-Ritz methods approximate the sign function using a projection on a Krylov subspace. To achieve a high accuracy this subspace must be taken quite large, which makes the method too costly. The new idea is to make a further projection on an even smaller, nested Krylov subspace. If additionally an intermediate preconditioning step is applied, this projection can be performed without affecting the accuracy of the approximation, and a substantial gain in efficiency is achieved for both Hermitian and non-Hermitian matrices. The numerical efficiency of the method is demonstrated on lattice configurations of sizes ranging from 4^4 to 10^4, and the new results are compared with those obtained with rational approximation methods.Comment: 17 pages, 12 figures, minor corrections, extended analysis of the preconditioning ste

Smooth parametric surfaces and n-sided patches

The theory of 'geometric continuity' within the subject of CAGD is reviewed. In particular, we are concerned with how parametric surface patches for CAGD can be pieced together to form a smooth Ck surface. The theory is applied to the problem of filling an n-sided hole occurring within a smooth rectangular patch complex. A number of solutions to this problem are surveyed

Strong μ\mu-Bases for Rational Tensor Product Surfaces and Extraneous Factors Associated to Bad Base Points and Anomalies at Infinity

We investigate conditions under which the resultant of a $\mu$-basis for a rational tensor product surface is the implicit equation of the surface without any extraneous factors. In this case, we also derive a formula for the implicit degree of the rational surface based only on the bidegree of the rational parametrization and the bidegrees of the elements of the $\mu$-basis without any knowledge of the number or multiplicities of the base points, assuming only that all the base points are local complete intersections. We conclude that in this case the implicit degree of a rational surface of bidegree $(m,n)$ is at most $mn$, so the rational surface must have at least $mn$ base points counting multiplicity. When the resultant of a $\mu$-basis generates extraneous factors, we show how to predict and compute these extraneous factors from either the existence of bad base points or anomalies occurring in the parametrization at infinity. Examples are provided to flesh out the theory

A survey of the representations of rational ruled surfaces

The rational ruled surface is a typical modeling surface in computer aided geometric design. A rational ruled surface may have different representations with respective advantages and disadvantages. In this paper, the authors revisit the representations of ruled surfaces including the parametric form, algebraic form, homogenous form and Pl¨ucker form. Moreover, the transformations between these representations are proposed such as parametrization for an algebraic form, implicitization for a parametric form, proper reparametrization of an improper one and standardized reparametrization for a general parametrization. Based on these transformation algorithms, one can give a complete interchange graph for the different representations of a rational ruled surface. For rational surfaces given in algebraic form or parametric form not in the standard form of ruled surfaces, the characterization methods are recalled to identify the ruled surfaces from them.Agencia Estatal de Investigació

Bounds for degrees of syzygies of polynomials defining a grade two ideal

We make explicit the exponential bound on the degrees of the polynomials appearing in the Effective Quillen-Suslin Theorem, and apply it jointly with the Hilbert-Burch Theorem to show that the syzygy module of a sequence of $m$ polynomials in $n$ variables defining a complete intersection ideal of grade two is free, and that a basis of it can be computed with bounded degrees. In the known cases, these bounds improve previous results

Modelling the pension system in an overlapping-generations general equilibrium modelling framework

This article presents a theoretical contribution to the field of overlapping-generations general equilibrium modelling, i.e. an upgrade of this branch of models with a pension system. Within the pension block we model both the first pension pillar, financed on a pay-as-you-go basis, and the fully-funded second pillar of the Slovenian pension system. The modelling of the first pension pillar is based on cash flows of the mandatory pension insurance institution, the relationship between the pension base and the pension, and the process of harmonising pension growth to wage growth. The modelling of the second pillar centres on implementation of the liquidity constraint. Use was made of supplementary pension profiles, and the ratio between premia paid and pensions paid out from supplementary pension insurance. The category of total pension was also introduced, and the model ensured that at every point households adjusted the scope of labour supply and their current consumption towards the target total pension.first pension pillar; general equilibrium models; liquidity constraint; MCP; overlapping generations; PAYG; pension system; second pension pillar

Representing rational curve segments and surface patches using semi-algebraic sets

We provide a framework for representing segments of rational planar curves or patches of rational tensor product surfaces with no singularities using semi-algebraic sets. Given a rational planar curve segment or a rational tensor product surface patch with no singularities, we find the implicit equation of the corresponding unbounded curve or surface and then construct an algebraic box defined by some additional equations and inequalities associated to the implicit equation. This algebraic box is proved to include only the given curve segment or surface patch without any extraneous parts of the unbounded curve or surface. We also explain why it is difficult to construct such an algebraic box if the curve segment or surface patch includes some singular points such as self-intersections. In this case, we show how to isolate a neighborhood of these special points from the corresponding curve segment or surface patch and to represent these special points with small curve segments or surface patches. This framework allows us to dispense with expensive approximation methods such as voxels for representing surface patches.National Natural Science Foundation of ChinaMinisterio de Ciencia, Innovación y Universidade

On the computation of singularities of parametrized ruled surfaces

Given a ruled surface V defined in the standard parametric form P(t1, t2), we present an algorithm that determines the singularities (and their multiplicities) of V from the parametrization P. More precisely, from P we construct an auxiliary parametric curve and we show how the problem can be simplified to determine the singularities of this auxiliary curve. Only one univariate resultant has to be computed and no elimination theory techniques are necessary. These results improve some previous algorithms for detecting singularities for the special case of parametric ruled surfaces.Ministerio de Ciencia, Innovación y Universidade