No abstract uploaded!

No abstract uploaded!

We propose a semi-discrete scheme for 2D Keller-Segel equations based on a symmetrization reformation, which is equivalent to the convex splitting method and is free of any nonlinear solver. We show that, this new scheme is unconditionally stable as long as the initial condition does not exceed certain threshold, and it asymptotically preserves the quasi-static limit in the transient regime. Furthermore, we prove that the fully discrete scheme is conservative and positivity preserving, which makes it ideal for simulations. The analogical schemes for the radial symmetric cases and the subcritical degenerate cases are also presented and analyzed. With extensive numerical tests, we verify the claimed properties of the methods and demonstrate their superiority in various challenging applications.

Parameterizations of manifolds are widely applied to the fields of numerical partial differential equations and computer graphics. To this end, in recent years several efficient and reliable numerical algorithms have been developed by different research groups for the computation of triangular and tetrahedral mesh parameterizations. However, it is still challenging when the topology of manifolds is nontrivial, e.g., the 3-manifold of a topological solid torus. In this paper, we propose a novel volumetric stretch energy minimization algorithm for volume-preserving parameterizations of toroidal polyhedra with a single boundary being mapped to a standard torus. In addition, the algorithm can also be used to compute the equiareal mapping between a genus-one closed surface and the standard torus. Numerical experiments indicate that the developed algorithm is effective and performs well on the bijectivity of the mapping. Applications on manifold registrations and partitions are demonstrated to show the robustness of our algorithms.

Denote by Mnthe set of n ×ncomplex matrices. Let f:Mn→[0, ∞)be a continuous map such that f(μUAU∗) =f(A)for any complex unit μ, A ∈Mnand unitary U∈Mn, f(X) =0if and only if X=0and the induced map t →f(tX)is monotonically increasing on [0, ∞)for any rank onenilpotent X∈Mn. Characterization is given for surjective maps φon Mnsatisfying f(AB−BA) =f(φ(A)φ(B) −φ(B)φ(A)). The general theorem isthen used to deduce results on special cases when the function is the pseudo spectrum and the pseudo spectral radius.