Local modification of Pythagorean-hodograph quintic spline curves using the B-spline form

The problems of determining the B–spline form of a C ² Pythagorean–hodograph (PH) quintic spline curve interpolating given points, and of using this form to make local modifications, are addressed. To achieve the correct order of continuity, a quintic B–spline basis constructed on a knot sequence in which each (interior) knot is of multiplicity 3 is required. C ² quintic bases on uniform triple knots are constructed for both open and closed C ² curves, and are used to derive simple explicit formulae for the B–spline control points of C ² PH quintic spline curves. These B-spline control points are verified, and generalized to the case of non–uniform knots, by applying a knot removal scheme to the Bézier control points of the individual PH quintic spline segments, associated with a set of six–fold knots. Based on the B–spline form, a scheme for the local modification of planar PH quintic splines, in response to a control point displacement, is proposed. Only two contiguous spline segments are modified, but to preserve the PH nature of the modified segments, the continuity between modified and unmodified segments must be relaxed from C ² to C ¹. A number of computed examples are presented, to compare the shape quality of PH quintic and “ordinary” cubic splines subject to control point modifications.     

A geometric approach to Euler’s addition formula

Plane quartic curves given by equations of the form y ²=P(x) with polynomials P of degree 4 represent singular models of elliptic curves which are directly related to elliptic integrals in the form studied by Euler and for which he developed his famous addition formulas. For cubic curves, the well-known secant and tangent construction establishes an immediate connection of addition formulas for the corresponding elliptic integrals with the structure of an algebraic group. The situation for quartic curves is considerably more complicated due to the presence of the singularity. We present a geometric construction, similar in spirit to the secant method for cubic curves, which defines an addition law on a quartic elliptic curve given by rational functions. Furthermore, we show how this addition on the curve itself corresponds to the addition in the (generalized) Jacobian variety of the curve, and we show how any addition formula for elliptic integrals of the form \(\int (1/\sqrt{P(x)})\,\mathrm{d}x\) with a quartic polynomial P can be derived directly from this addition law.     

Congruent elliptic curves with non-trivial Shafarevich-Tate groups: Distribution part

Given a large positive number x and a positive integer k, we denote by Qk(x) the set of congruent elliptic curves E(n): y2= z3- n2 z with positive square-free integers n x congruent to one modulo eight,having k prime factors and each prime factor congruent to one modulo four. We obtain the asymptotic formula for the number of congruent elliptic curves E(n)∈ Qk(x) with Mordell-Weil ranks zero and 2-primary part of Shafarevich-Tate groups isomorphic to(Z/2Z)2. We also get a lower bound for the number of E(n)∈ Qk(x)with Mordell-Weil ranks zero and 2-primary part of Shafarevich-Tate groups isomorphic to(Z/2Z)4. The key ingredient of the proof of these results is an independence property of residue symbols. This property roughly says that the number of positive square-free integers n x with k prime factors and residue symbols(quadratic and quartic) among its prime factors being given compatible values does not depend on the actual values.     

On the algorithm of smoothing by a spline with bilateral constraints

In this paper, the problem of constructing a spline σ in the Hilbert space, which satisfies the bilateral constraints z ^? ≤ Aσ ≤ z ⁺ with a linear operator A and minimizes the squared Hilbert seminorm is studied. A solution to this problem can be obtained with convex programming iterative methods, in particular, with the gradient projection method. A modification of the gradient projection method is proposed, which allows one to find a set of active constraints with a smaller number of iterations. The efficiency of the modification proposed is demonstrated in the problem of approximation with a pseudolinear bivariate spline.     

The Nöther and Riemann-Roch type theorems for piecewise algebraic curve

A piecewise algebraic curve is a curve determined by the zero set of a bivariate spline function. In this paper, the Nöther type theorems for C ^µ piecewise algebraic curves are obtained. The theory of the linear series of sets of places on the piecewise algebraic curve is also established. In this theory, singular cycles are put into the linear series, and a complete series of the piecewise algebraic curves consists of all effective ordinary cycles in an equivalence class and all effective singular cycles which are equivalent specifically to any effective ordinary cycle in the equivalence class. This theory is a generalization of that of linear series of the algebraic curve. With this theory and the fundamental theory of multivariate splines on smoothing cofactors and global conformality conditions, and the results on the general expression of multivariate splines, we get a formula on the index, the order and the dimension of a complete series of the irreducible C ^µ piecewise algebraic curves and the degree, the genus and the smoothness of the curves, hence the Riemann-Roch type theorem of the C ^µ piecewise algebraic curve is established.     

Non-polynomial spline method for the solution of two-dimensional linear wave equations with a nonlinear source term

In this paper, two classes of methods are developed for the solution of two space dimensional wave equations with a nonlinear source term. We have used non-polynomial cubic spline function approximations in both space directions. The methods involve some parameters, by suitable choices of the parameters, a new high accuracy three time level scheme of order O(h ⁴ + k ⁴ + τ ² + τ ² h ² + τ ² k ²) has been obtained. Stability analysis of the methods have been carried out. The results of some test problems are included to demonstrate the practical usefulness of the proposed methods. The numerical results for the solution of two dimensional sine-Gordon equation are compared with those already available in literature.     

On using cubic spline for the solution of problems in calculus of variations

Two different approaches based on cubic B-spline are developed to approximate the solution of problems in calculus of variations. Both direct and indirect methods will be described. It is known that, when using cubic spline for interpolating a function g∈C⁴[a,b] on a uniform partition with the step size h, the optimal order of convergence derived is O(h⁴). In Zarebnia and Birjandi (J. Appl. Math. 1–10, 2012) a non-optimal O(h²) method based on cubic B-spline has been used to solve the problems in calculus of variations. In this paper at first we will obtain an optimal O(h⁴) indirect method using cubic B-spline to approximate the solution. The convergence analysis will be discussed in details. Also a locally superconvergent O(h⁶) indirect approximation will be describe. Finally the direct method based on cubic spline will be developed. Some test problems are given to demonstrate the efficiency and applicability of the numerical methods.     

On the problem of determining the parameters of a layered elastic medium and an impulse source

Considering the linear system of elasticity equations describing the wave propagation in the half-space ? ₊ ³ = {x ∈ ?³ | x ₃ > 0} we address the problem of determining the density and elastic parameters which are piecewise constant functions of x ₃. The shape is unknown of a point-like impulse source that excites elastic oscillations in the half-space. We show that under certain assumptions on the source shape and the parameters of the elastic medium the displacements of the boundary points of the half-space for some finite time interval (0, T) uniquely determine the normalized density (with respect to the first layer) and the elastic Lamé parameters for x ₃ ∈ [0, H], where H = H(T). We give an algorithmic procedure for constructing the required parameters.     

The Elekes–Szabó Theorem in four dimensions

Let F ∈ C[x, y, s, t] be an irreducible constant-degree polynomial, and let A,B,C,D ? C be finite sets of size n. We show that F vanishes on at most O(n^8/3) points of the Cartesian product A × B × C × D, unless F has a special group-related form. A similar statement holds for A,B,C,D of unequal sizes, with a suitably modified bound on the number of zeros. This is a four-dimensional extension of our recent improved analysis of the original Elekes–Szabó theorem in three dimensions. We give three applications: an expansion bound for three-variable real polynomials that do not have a special form, a bound on the number of coplanar quadruples on a space curve that is neither planar nor quartic, and a bound on the number of four-point circles on a plane curve that has degree at least five.     

Holomorphic differentials and Laguerre deformation of surfaces

A Laguerre geometric local characterization is given of L-minimal surfaces and Laguerre deformations (T-transforms) of L-minimal isothermic surfaces in terms of the holomorphicity of a quartic and a quadratic differential. This is used to prove that, via their L-Gauss maps, the T-transforms of L-minimal isothermic surfaces have constant mean curvature \(H=r\) in some translate of hyperbolic 3-space \({\mathbb {H}}^3(-r^2)\subset \mathbb {R}^4_1\), de Sitter 3-space \({\mathbb {S}}^3_1(r^2)\subset \mathbb {R}^4_1\), or have mean curvature \(H=0\) in some translate of a time-oriented lightcone in \(\mathbb {R}^4_1\). As an application, we show that various instances of the Lawson isometric correspondence can be viewed as special cases of the T-transformation of L-isothermic surfaces with holomorphic quartic differential.     

Approximation of conic section by quartic Bzier curve with endpoints continuity condition

A new method for approximation of conic section by quartic B′ezier curve is presented, based on the quartic B′ezier approximation of circular arcs. Here we give an upper bound of the Hausdorff distance between the conic section and the approximation curve, and show that the error bounds have the approximation order of eight. Furthermore, our method yields quartic G2 continuous spline approximation of conic section when using the subdivision scheme,and the effectiveness of this method is demonstrated by some numerical examples.     

Stability of Cubic and Quartic Functional Equations in Non-Archimedean Spaces

We prove generalized Hyers-Ulam–Rassias stability of the cubic functional equation f(kx+y)+f(kx?y)=k[f(x+y)+f(x?y)]+2(k ³?k)f(x) for all \(k\in \Bbb{N}\) and the quartic functional equation f(kx+y)+f(kx?y)=k ²[f(x+y)+f(x?y)]+2k ²(k ²?1)f(x)?2(k ²?1)f(y) for all \(k\in \Bbb{N}\) in non-Archimedean normed spaces.     

The graph Kre(4) does not exist

Suppose that a strongly regular graph Γ with parameters (v, k, λ, μ) has eigenvalues k, r, and s. If the graphs Γ and \(\bar \Gamma \) are connected, then the following inequalities, known as Krein’s conditions, hold: (i) (r + 1)(k + r + 2rs) ≤ (k + r)(s + 1)² and (ii) (s + 1)(k + s + 2rs) ≤ (k + s)(r + 1)². We say that Γ is a Krein graph if one of Krein’s conditions (i) and (ii) is an equality for this graph. A triangle-free Krein graph has parameters ((r ² + 3r)², r ³ + 3r ² + r, 0, r ² + r). We denote such a graph by Kre(r). It is known that, in the cases r = 1 and r = 2, the graphs Kre(r) exist and are unique; these are the Clebsch and Higman–Sims graphs, respectively. The latter was constructed in 1968 together with the Higman–Sims sporadic simple group. A.L. Gavrilyuk and A.A. Makhnev have proved that the graph Kre(3) does not exist. In this paper, it is proved that the graph Kre(4) (a strongly regular graph with parameters (784, 116, 0, 20)) does not exist either.     

The cubic Fermat equation and complex multiplication on the Deuring normal form

The arithmetic on elliptic curves in Deuring normal form is shown to be related to solutions of the Fermat equation 27X ³+27Y ³=X ³ Y ³. This arithmetic is used to give conditions for the existence of multipliers μ on supersingular elliptic curves in characteristic p for which μ ²=?3p. Together with an explicit factorization of a certain class equation, these conditions imply that the number of irreducible binomial quadratic factors (mod p) of the Legendre polynomial P _(p?e)/3(x) of degree (p?e)/3 is a simple linear function of the class number of the quadratic field \(\mathbb{Q}(\sqrt{-3p})\).     

Steepest descent curves of convex functions on surfaces of constant curvature

Let S be a complete surface of constant curvature K = ±1, i.e., S ² or л ², and Ω ? S a bounded convex subset. If S = S ², assume also diameter(Ω) < π/2. It is proved that the length of any steepest descent curve of a quasi-convex function in Ω is less than or equal to the perimeter of Ω. This upper bound is actually proved for the class of G-curves, a family of curves that naturally includes all steepest descent curves. In case S = S ², the existence of G-curves, whose length is equal to the perimeter of their convex hull, is also proved, showing that the above estimate is indeed optimal. The results generalize theorems by Manselli and Pucci on steepest descent curves in the Euclidean plane.     

Frobenius nonclassicality of Fermat curves with respect to cubics

For Fermat curves F: aX ⁿ + bY ⁿ = Z ⁿ defined over F _q , we establish necessary and sufficient conditions for F to be F _q -Frobenius nonclassical with respect to the linear system of plane cubics. In the new F _q -Frobenius nonclassical cases, we determine explicit formulas for the number N _q (F) of F _q -rational points on F. For the remaining Fermat curves, nice upper bounds for N _q (F) are immediately given by the Stöhr–Voloch Theory.     

8-ranks of class groups of quadratic number fields and their densities

For \(F = \mathbb{Q}\left( {\sqrt {\varepsilon pq} } \right)\), ? ∈ {±1, ±2}, primes ?p ≡ q ≡ 1 mod 4, we give the necessary and sufficient conditions for 8-ranks of narrow class groups of F equal to 1 or 2 such that we can calculate their densities. All results are stated in terms of congruence relations of p, q modulo 2 ⁿ , the quartic residue symbol \(\left( {\frac{p}{q}} \right)_4\) and binary quadratic forms such as q ^h(?2p)/4 = x ² + 2py ², where h(?2p) is the class number of \(\mathbb{Q}\left( {\sqrt { - 2p} } \right)\). The results are very useful for numerical computations.     

Bernstein-Bézier techniques for divergence of polynomial spline vector fields in ℝ<Superscript><Emphasis Type="Italic">n</Emphasis></Superscript>

Bernstein-Bézier techniques for analyzing polynomial spline fields in n variables and their divergence are developed. Dimension and a minimal determining set for continuous piecewise divergence-free spline fields on the Alfeld split of a simplex in ? ⁿ are obtained using the new techniques, as well as the dimension formula for continuous piecewise divergence-free splines on the Alfeld refinement of an arbitrary simplicial partition in ? ⁿ .     

Quantization causes waves: Smooth finitely computable functions are affine

Every automaton (a letter-to-letter transducer) A whose both input and output alphabets are F _p = {0, 1,..., p - 1} produces a 1-Lipschitz map f _A from the space Z _p of p-adic integers to Z _p . The map fA can naturally be plotted in a unit real square I² ? R²: To an m-letter non-empty word v = γ _m-1γ _m-2... γ0 there corresponds a number 0.v ∈ R with base-p expansion 0.γ _m-1γ _m-2... γ0; so to every m-letter input word w = α _m-1α _m-2 ··· α0 of A and to the respective m-letter output word a(w) = β _m-1β _m-2 ··· β0 of A there corresponds a point (0.w; 0.a(w)) ∈ R². Denote P(A) a closure of the point set (0.w; 0.a(w)) where w ranges over all non-empty words.We prove that once some points of P(A) constitute a C ²-smooth curve in R², the curve is a segment of a straight line with a rational slope. Moreover, when identifying P(A) with a subset of a 2-dimensional torus T² ∈ R³, the smooth curves from P(A) constitute a collection of torus windings which can be ascribed to complex-valued functions ψ(x, t) = e ^i(Ax-2πBt) (x, t ∈ R), i.e., to matter waves. As automata are causal discrete systems, the main result may serve a mathematical reasoning why wave phenomena are inherent in quantum systems: This is just because of causality principle and discreteness of matter.     

Novel polynomial Bernstein bases and Bézier curves based on a general notion of polynomial blossoming

We introduce the G-blossom of a polynomial by altering the diagonal property of the classical blossom, replacing the identity function by arbitrary linear functions G=G(t). By invoking the G-blossom, we construct G-Bernstein bases and G-Bézier curves and study their algebraic and geometric properties. We show that the G-blossom provides the dual functionals for the G-Bernstein basis functions and we use this dual functional property to prove that G-Bernstein basis functions form a partition of unity and satisfy a Marsden identity. We also show that G-Bézier curves share several other properties with classical Bézier curves, including affine invariance, interpolation of end points, and recursive algorithms for evaluation and subdivision. We investigate the effect of the linear functions G on the shape of the corresponding G-Bézier curves, and we derive some necessary and sufficient conditions on the linear functions G which guarantee that the corresponding G-Bézier curves are of Pólya type and variation diminishing. Finally we prove that the control polygons generated by recursive subdivision converge to the original G-Bézier curve, and we derive the geometric rate of convergence of this algorithm.