On a nonlinear multigrid algorithm with primal relaxation for the image total variation minimisation

Digital image restoration has drawn much attention in the recent years and a lot of research has been done on effective variational partial differential equation models and their theoretical studies. However there remains an urgent need to develop fast and robust iterative solvers, as the underlying problem sizes are large. This paper proposes a fast multigrid method using primal relaxations. The basic primal relaxation is known to get stuck at a ‘local’ non-stationary minimum of the solution, which is usually believed to be ‘non-smooth’. Our idea is to utilize coarse level corrections, overcoming the deadlock of a basic primal relaxation scheme. A further refinement is to allow non-regular coarse levels to correct the solution, which helps to improve the multilevel method. Numerical experiments on both 1D and 2D images are presented.     

Monte Carlo approximation through Gibbs output in generalized linear mixed models

Geyer (J. Roy. Statist. Soc. 56 (1994) 291) proposed Monte Carlo method to approximate the whole likelihood function. His method is limited to choosing a proper reference point. We attempt to improve the method by assigning some prior information to the parameters and using the Gibbs output to evaluate the marginal likelihood and its derivatives through a Monte Carlo approximation. Vague priors are assigned to the parameters as well as the random effects within the Bayesian framework to represent a non-informative setting. Then the maximum likelihood estimates are obtained through the Newton Raphson method. Thus, out method serves as a bridge between Bayesian and classical approaches. The method is illustrated by analyzing the famous salamander mating data by generalized linear mixed models.     

Generalization of Strang's Preconditioner with Applications to Toeplitz Least Squares Problems

In this paper, we propose a method to generalize Strang's circulant preconditioner for arbitrary n-by-n matrices A_n. The th column of our circulant preconditioner Sn is equal to the th column of the given matrix A_n. Thus if A_n is a square Toeplitz matrix, then Sn is just the Strang circulant preconditioner. When Sn is not Hermitian, our circulant preconditioner can be defined as . This construction is similar to the forward-backward projection method used in constructing preconditioners for tomographic inversion problems in medical imaging. We show that if the matrix A_n has decaying coefficients away from the main diagonal, then is a good preconditioner for An. Comparisons of our preconditioner with other circulant-based preconditioners are carried out for some 1-D Toeplitz least squares problems: min ∥ b - Ax∥₂. Preliminary numerical results show that our preconditioner performs quite well, in comparison to other circulant preconditioners. Promising test results are also reported for a 2-D deconvolution problem arising in ground-based atmospheric imaging.     

Behavior of a PZT ring under non-uniform mechanical stress

Ring-shaped lead zirconate titanate (PZT) piezoelectric vibrators were subjected to non-uniform mechanical stress applied by bolt clamping. The effect of mechanical stress on the effective electromechanical coupling factor (k_eff) and mechanical quality factor (Q_m) of the thickness and wall thickness modes was studied by an equivalent electric circuit analysis. The initiation and propagation of cracks under mechanical stress were also discussed based on the resonance method and the indentation technique. k_eff for both the thickness and wall thickness modes decreased with increase in mechanical stress due to de-poling of the PZT. Q_m of the thickness mode dropped sharply with increase in mechanical stress while Q_m of the wall thickness mode remained almost unchange. Existence of microcracks in a PZT vibrator can be detected by the occurrence of spurious vibrations at the wall thickness mode in the electrical impedance vs. frequency spectra.     

Definite regular quadratic forms over

Let be a power of an odd prime, and be the ring of polynomials over a finite field of elements. A quadratic form over is said to be regular if globally represents all polynomials that are represented by the genus of . In this paper, we study definite regular quadratic forms over . It is shown that for a fixed , there are only finitely many equivalence classes of regular definite primitive quadratic forms over , regardless of the number of variables. Characterizations of those which are universal are also given.

     

Optimal Partition Trees

We revisit one of the most fundamental classes of data structure problems in computational geometry: range searching. Matoušek (Discrete Comput. Geom. 10:157–182, 1993) gave a partition tree method for d-dimensional simplex range searching achieving O(n) space and O(n ^1−1/d ) query time. Although this method is generally believed to be optimal, it is complicated and requires O(n ^1+ε ) preprocessing time for any fixed ε>0. An earlier method by Matoušek (Discrete Comput. Geom. 8:315–334, 1992) requires O(nlogn) preprocessing time but O(n ^1−1/d log ^O(1) n) query time. We give a new method that achieves simultaneously O(nlogn) preprocessing time, O(n) space, and O(n ^1−1/d ) query time with high probability. Our method has several advantages:

On Levels in Arrangements of Surfaces in Three Dimensions

A favorite open problem in combinatorial geometry is to determine the worst-case complexity of a level in an arrangement. Up to now, nontrivial upper bounds in three dimensions are known only for the linear cases of planes and triangles. We propose the first technique that can deal with more general surfaces in three dimensions. For example, in an arrangement of n ??pseudo-planes?? or ??pseudo-spherical patches?? (where the main criterion is that each triple of surfaces has at most two common intersections), we prove that there are at most O(n ^2.997) vertices at any given level.     

Air stability of PTCDI‐C13‐based n‐OFETs on polymer interfacial layers

This Letter reports the novel use of poly(9‐vinylcarbazole) (PVK) as a dielectric interfacial layer for n‐type organic field‐effect transistors (n‐OFETs). With PVK, both the air stability and electron mobility of N,N′‐ditridecylperylene‐3,4,9,10‐tetracarboxylic diimide (PTCDI‐C13)‐based OFETs were improved. Among the PVKs with different weight‐average molecular weight (M_w), PVK with high M_w showed good performance. The high glass transition temperature of PVK enabled thermal post annealing of the active layer, which resulted in a high electron mobility of 0.61 cm²/Vs. This mobility was maintained at 90% and 59% after 4 days and 105 days in air, respectively. The PVK interfacial layer reduced the trapped charges in the PTCDI‐C13‐based n‐OFET for air‐exposure and caused a decrease in the threshold voltage shift.