This paper presents a novel scheme for constructing bivariate spline surfaces over triangular meshes which are topologically equivalent to a disk. The core part of the scheme is a set of knot selection rules that define local configurations of a triangulation called the directed-one-ring-cycle (D1RC) configurations and bivariate splines defined over a D1RC configuration that are new non-tensor-product splines and possess many nice properties of a univariate B-spline. Using D1RC splines, we take an input triangular mesh as a control mesh and define a bivariate spline surface from the control mesh, which mimics the standard NURBS modeling. Moreover, we can introduce sharp features into the overall smooth spline surface by simply setting special D1RC configurations. As a result, the proposed scheme can define spline surfaces in a way similar to that of NURBS, but has less restriction on the connectivity of the input control mesh.

It is known that the maximum diameter for the rupture-risk assessment of the abdominal aortic aneurysm is a generally good method, but not sufficient. Alternative features obtained with computational modeling may provide additional useful criteria. Though computational approaches are noninvasive, they are often time-consuming because of the high computational complexity. In this paper, we present a highly parallel algorithm for the numerical simulation of unsteady blood flows in the patient-specific abdominal aorta. We model the blood flow with the unsteady incompressible Navier-Stokes equations, and solve the discretized system with a highly scalable domain decomposition method. With this approach, the complete flow field can be obtained in less than an hour, instead of days with older methods. We show experimentally that the proposed method offers accurate solutions of the pressure, the velocity and the wall shear stress, and the parallel efficiency is higher than 70% on a parallel computer with more than 1,000 processor cores.

Given an open book decomposition (Σ,τ) of a three manifold Y , Thurston and Winkelnkemper [TW] construct a specific contact form a on Y . Given a spin-c Dirac operator D on Y , the contact form naturally associates a one parameter family of Dirac operators Dr = D− ir 2 cl(a) for r ≥ 0. When r >> 1, we prove that the spectrum of Dr = D0 − ir 2 cl(a) within [−1 2r 1 2 , 1 2r 1 2 ] are almost uniformly distributed. With the result in Part I [Ts1], it implies that the subleading order term of the spectral flow from D0 to Dr is of order r(logr) 9 2 . Besides the interests of the spectral flow, the method of this paper provide a tool to analyze the Dirac operator on an open book decomposition.

No abstract uploaded!

Genome comparison is a vital research area of bioinformatics. For large-scale genome comparisons, the Multiple Sequence Alignment (MSA) methods have been impractical to use due to its algorithmic complexity. In this study, we propose a novel alignment-free method based on the one-to-one correspondence between a DNA sequence and its complete central moment vector of the cumulative Fourier power and phase spectra. In addition, the covariance between the four nucleotides in the power and phase spectra is included. We use the cumulative Fourier power and phase spectra to define a 28-dimensional vector for each DNA sequence. Euclidean distances between the vectors can measure the dissimilarity between DNA sequences. We perform testing with datasets of different sizes and types including simulated DNA sequences, exon-intron and complete genomes. The results show that our method is more accurate and efficient for performing hierarchical clustering than other alignment-free methods and MSA methods.