Some results on holomorphic generalized functions

In this note, we present three independent results within generalized complex analysis (in the Colombeau sense). The first of them deals with non-removable singularities; we construct a generalized function u on an open subset Omega of C(n), which is not a holomorphic generalized function on Omega but it is a holomorphic generalized function on Omega\S, where S is a hypersurface contained in Omega. The second result shows the existence of a holomorphic generalized function with prescribed values in the zero-set of a classical holomorphic function. The last result states the existence of a compactly supported solution to the (partial derivative) over bar operator.

A generalized notion of compliance

It is a known fact in structural optimization that for structures subject to prescribed non-zero displacements the work done by the loads is not agood measure of compliance, neither is the stored elastic energy. We briefly discuss a possible alternative measure of compliance, valid for general boundary conditions. We also present the adjoint states (necessary for the computation of the structural derivative) for the three functionals under consideration. (C) 2011 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

On derivatives and norms of generalized matrix functions and respective symmetric powers

In recent papers, the authors obtained formulas for directional derivatives of all orders, of the immanant and of the m-th xi-symmetric tensor power of an operator and a matrix, when xi is a character of the full symmetric group. The operator norm of these derivatives was also calculated. In this paper, similar results are established for generalized matrix functions and for every symmetric tensor power.

Mathematical models of generalized diffusion

We discuss in this paper equations describing processes involving non-linear and higher-order diffusion. We focus on a particular case (u(t) = 2 lambda (2)(uu(x))(x) + lambda (2)u(xxxx)), which is put into analogy with the KdV equation. A balance of nonlinearity and higher-order diffusion enables the existence of self-similar solutions, describing diffusive shocks. These shocks are continuous solutions with a discontinuous higher-order derivative at the shock front. We argue that they play a role analogous to the soliton solutions in the dispersive case. We also discuss several physical instances where such equations are relevant.

HIGHER-DERIVATIVE SCHWINGER MODEL

Using the operator formalism, we obtain the bosonic representation for the free fermion field satisfying an equation of motion with higher-order derivatives. Then, we consider the operator solution of a generalized Schwinger model with higher-derivative coupling. Since the increasing of the derivative order implies the introduction of an equivalent number of extra fermionic degrees of freedom, the mass acquired by the gauge field is bigger than the one for the standard two-dimensional QED. An analysis of the problem from the functional integration point of view corroborates the findings of canonical quantization, and corrects certain results previously announced in the literature on the basis of Fujikawa's technique.

Fermion helicity flip in higher-derivative electromagnetism

It is shown that massive fermions have their helicity flipped on account of their interaction with an electromagnetic field described by Podolsky's generalized electrodynamics. Massless fermions, in turn, seem to be unaffected by the electromagnetic field as far as their helicity is concerned. © Springer-Verlag 1997.

Composition for a class of generalized functions in Colombeau's theory

In Colombeau's theory, given an open subset Ω of ℝn, there is a differential algebra G(Ω) of generalized functions which contains in a natural way the space D′(Ω) of distributions as a vector subspace. There is also a simpler version of the algebra G,(Ω). Although this subalgebra does not contain, in canonical way, the space D′(Ω) is enough for most applications. This work is developed in the simplified generalized functions framework. In several applications it is necessary to compute higher intrinsic derivatives of generalized functions, and since these derivatives are multilinear maps, it is necessary to define the space of generalized functions in Banach spaces. In this article we introduce the composite function for a special class of generalized mappings (defined in open subsets of Banach spaces with values in Banach spaces) and we compute the higher intrinsic derivative of this composite function.

Higher-derivative non-Abelian gauge fields via the Faddeev-Jackiw formalism

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

On abstract differential equations with nonlocal conditions involving the temporal derivative of the solution

In this paper we discuss the existence of solutions for a class of abstract differential equations with nonlocal conditions for which the nonlocal term involves the temporal derivative of the solution. Some concrete applications to parabolic differential equations with nonlocal conditions are considered. (C) 2012 Royal Dutch Mathematical Society (KWG). Published by Elsevier B.V. All rights reserved.

R-estimator of location of the generalized secant hyperbolic distribution

The generalized secant hyperbolic distribution (GSHD) proposed in Vaughan (2002) includes a wide range of unimodal symmetric distributions, with the Cauchy and uniform distributions being the limiting cases, and the logistic and hyperbolic secant distributions being special cases. The current article derives an asymptotically efficient rank estimator of the location parameter of the GSHD and suggests the corresponding one- and two-sample optimal rank tests. The rank estimator derived is compared to the modified MLE of location proposed in Vaughan (2002). By combining these two estimators, a computationally attractive method for constructing an exact confidence interval of the location parameter is developed. The statistical procedures introduced in the current article are illustrated by examples.

Mach bands and multiscale models of spatial vision:the role of first, second, and third derivative operators in encoding bars and edges

Ernst Mach observed that light or dark bands could be seen at abrupt changes of luminance gradient in the absence of peaks or troughs in luminance. Many models of feature detection share the idea that bars, lines, and Mach bands are found at peaks and troughs in the output of even-symmetric spatial filters. Our experiments assessed the appearance of Mach bands (position and width) and the probability of seeing them on a novel set of generalized Gaussian edges. Mach band probability was mainly determined by the shape of the luminance profile and increased with the sharpness of its corners, controlled by a single parameter (n). Doubling or halving the size of the images had no significant effect. Variations in contrast (20%-80%) and duration (50-300 ms) had relatively minor effects. These results rule out the idea that Mach bands depend simply on the amplitude of the second derivative, but a multiscale model, based on Gaussian-smoothed first- and second-derivative filtering, can account accurately for the probability and perceived spatial layout of the bands. A key idea is that Mach band visibility depends on the ratio of second- to first-derivative responses at peaks in the second-derivative scale-space map. This ratio is approximately scale-invariant and increases with the sharpness of the corners of the luminance ramp, as observed. The edges of Mach bands pose a surprisingly difficult challenge for models of edge detection, but a nonlinear third-derivative operation is shown to predict the locations of Mach band edges strikingly well. Mach bands thus shed new light on the role of multiscale filtering systems in feature coding. © 2012 ARVO.

Relaxation procedures for an iterative MFS algorithm for the stable reconstruction of elastic fields from Cauchy data in two-dimensional isotropic linear elasticity

We investigate two numerical procedures for the Cauchy problem in linear elasticity, involving the relaxation of either the given boundary displacements (Dirichlet data) or the prescribed boundary tractions (Neumann data) on the over-specified boundary, in the alternating iterative algorithm of Kozlov et al. (1991). The two mixed direct (well-posed) problems associated with each iteration are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen via the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise becomes dominant and the errors in predicting the exact solutions increase, is also presented. The MFS-based iterative algorithms with relaxation are tested for Cauchy problems for isotropic linear elastic materials in various geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the proposed method.

Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral equations

We propose and investigate a method for the stable determination of a harmonic function from knowledge of its value and its normal derivative on a part of the boundary of the (bounded) solution domain (Cauchy problem). We reformulate the Cauchy problem as an operator equation on the boundary using the Dirichlet-to-Neumann map. To discretize the obtained operator, we modify and employ a method denoted as Classic II given in [J. Helsing, Faster convergence and higher accuracy for the Dirichlet–Neumann map, J. Comput. Phys. 228 (2009), pp. 2578–2576, Section 3], which is based on Fredholm integral equations and Nyström discretization schemes. Then, for stability reasons, to solve the discretized integral equation we use the method of smoothing projection introduced in [J. Helsing and B.T. Johansson, Fast reconstruction of harmonic functions from Cauchy data using integral equation techniques, Inverse Probl. Sci. Eng. 18 (2010), pp. 381–399, Section 7], which makes it possible to solve the discretized operator equation in a stable way with minor computational cost and high accuracy. With this approach, for sufficiently smooth Cauchy data, the normal derivative can also be accurately computed on the part of the boundary where no data is initially given.

Fast reconstruction of harmonic functions from Cauchy data using integral equation techniques

We consider the problem of stable determination of a harmonic function from knowledge of the solution and its normal derivative on a part of the boundary of the (bounded) solution domain. The alternating method is a procedure to generate an approximation to the harmonic function from such Cauchy data and we investigate a numerical implementation of this procedure based on Fredholm integral equations and Nyström discretization schemes, which makes it possible to perform a large number of iterations (millions) with minor computational cost (seconds) and high accuracy. Moreover, the original problem is rewritten as a fixed point equation on the boundary, and various other direct regularization techniques are discussed to solve that equation. We also discuss how knowledge of the smoothness of the data can be used to further improve the accuracy. Numerical examples are presented showing that accurate approximations of both the solution and its normal derivative can be obtained with much less computational time than in previous works.

An alternating boundary integral based method for a Cauchy problem for the Laplace equation in a quadrant

We study the Cauchy problem for the Laplace equation in a quadrant (quarter-plane) containing a bounded inclusion. Given the values of the solution and its derivative on the edges of the quadrant the solution is reconstructed on the boundary of the inclusion. This is achieved using an alternating iterative method where at each iteration step mixed boundary value problems are being solved. A numerical method is also proposed and investigated for the direct mixed problems reducing these to integral equations over the inclusion. Numerical examples verify the efficiency of the proposed scheme.