On reduction of Hilbert-Blumenthal varieties

Let $O_F$ be the ring of integers of a totally real field $F$ of degree $g$. We study the reduction of the moduli space of separably polarized abelian $O_F$-varieties of dimension $g$ modulo $p$ for a fixed prime $p$. The invariants and related conditions for the objects in the moduli space are discussed. We construct a scheme-theoretic stratification by $a$-numbers on the Rapoport locus and study the relation with the slope stratification. In particular, we recover the main results of Goren and Oort [GO, J. Alg. Geom. 2000] on the stratifications when $p$ is unramified in $O_F$. We also prove the strong Grothendieck conjecture for the moduli space in some restricted cases, particularly when $p$ is totally ramified in $O_F$.Comment: A shortened revised versio

Algorithms for Positive Semidefinite Factorization

This paper considers the problem of positive semidefinite factorization (PSD factorization), a generalization of exact nonnegative matrix factorization. Given an $m$-by-$n$ nonnegative matrix $X$ and an integer $k$, the PSD factorization problem consists in finding, if possible, symmetric $k$-by-$k$ positive semidefinite matrices $\{A^1,...,A^m\}$ and $\{B^1,...,B^n\}$ such that $X_{i,j}=\text{trace}(A^iB^j)$ for $i=1,...,m$, and $j=1,...n$. PSD factorization is NP-hard. In this work, we introduce several local optimization schemes to tackle this problem: a fast projected gradient method and two algorithms based on the coordinate descent framework. The main application of PSD factorization is the computation of semidefinite extensions, that is, the representations of polyhedrons as projections of spectrahedra, for which the matrix to be factorized is the slack matrix of the polyhedron. We compare the performance of our algorithms on this class of problems. In particular, we compute the PSD extensions of size $k=1+ \lceil \log_2(n) \rceil$ for the regular $n$-gons when $n=5$, $8$ and $10$. We also show how to generalize our algorithms to compute the square root rank (which is the size of the factors in a PSD factorization where all factor matrices $A^i$ and $B^j$ have rank one) and completely PSD factorizations (which is the special case where the input matrix is symmetric and equality $A^i=B^i$ is required for all $i$).Comment: 21 pages, 3 figures, 3 table

On the equivalence between the cell-based smoothed finite element method and the virtual element method

We revisit the cell-based smoothed finite element method (SFEM) for quadrilateral elements and extend it to arbitrary polygons and polyhedrons in 2D and 3D, respectively. We highlight the similarity between the SFEM and the virtual element method (VEM). Based on the VEM, we propose a new stabilization approach to the SFEM when applied to arbitrary polygons and polyhedrons. The accuracy and the convergence properties of the SFEM are studied with a few benchmark problems in 2D and 3D linear elasticity. Later, the SFEM is combined with the scaled boundary finite element method to problems involving singularity within the framework of the linear elastic fracture mechanics in 2D