Related papers: Robust and efficient validation of the linear hexa…
At present, mass matrix of solid fifteen node wedge element is computed by means of eighteen-point (Gauss points) numerical integration scheme. Herein, this widely accepted scheme is being challenged. We derive a novel, easy-to-implement,…
In this paper we describe a procedure for calculating the number of regular octahedrons that have vertices with coordinates in the set {0,1,...,n}. As a result, we introduce a new sequence in ``The Online Encyclopedia of Integer Sequences"…
By a direct computation, we show that the $P_2$ interpolation of a $P_3$ function is also a local $H^1$-projection on uniform tetrahedral meshes, i.e., the difference is $H^1$-orthogonal to the $P_2$ Lagrange basis function on the support…
A simple cut-and-patch method is presented for the construction and classification for fullerenes belonging to the octahedral point groups, $O$ or $O_h$. In order to satisfy the symmetry requirement of the octahedral group, suitable numbers…
This paper presents a new mixed finite element method for the Cahn-Hilliard equation. The well-posedness of the mixed formulation is established and the error estimates for its linearized fully discrete scheme are provided. The new mixed…
We propose one finite element method for both second order linear uniformly elliptic PDE in non-divergence form and the uniformly elliptic Hamilton-Jacobi-Bellman (HJB) equation. For both linear elliptic PDE in non-divergence form and the…
In this paper, we propose a numerical method for verifying the positiveness of solutions to semilinear elliptic equations. We provide a sufficient condition for a solution to an elliptic equation to be positive in the domain of the…
Matching codes are stabilizer codes based on Kitaev's honeycomb lattice model. The hexagonal form of these codes are particularly well-suited to the heavy-hexagon device layouts currently pursued in the hardware of IBM Quantum. Here we show…
This paper describes in detail the implementation of a finite element technique for solving the compressible Navier-Stokes equations that is provably robust and demonstrates excellent performance on modern computer hardware. The method is…
A rational triangle has rational edge-lengths and area; a rational tetrahedron has rational faces and volume; either is Heronian when its edge-lengths are integer, and proper when its content is nonzero. A variant proof is given, via…
Knowing the symmetries of a polyhedron can be very useful for the analysis of its structure as well as for practical polyhedral computations. In this note, we study symmetry groups preserving the linear, projective and combinatorial…
Polyhedral projection is a main operation of the polyhedron abstract domain.It can be computed via parametric linear programming (PLP), which is more efficient than the classic Fourier-Motzkin elimination method.In prior work, PLP was done…
The classical Cauchy rigidity theorem for convex polytopes reads that if two convex polytopes have isometric developments then they are congruent. In other words, we can decide whether two polyhedra are isometric or not by using their…
The finite-element method is a preferred numerical method when electromagnetic fields at high accuracy are to be computed in nano-optics design. Here, we demonstrate a finite-element method using hp-adaptivity on tetrahedral meshes for…
This paper proposes an efficient algorithm for testing copositivity of homogeneous polynomials over the positive semidefinite cone. The algorithm is based on a novel matrix optimization reformulation and requires solving a hierarchy of…
Finite elements of higher continuity, say conforming in $H^2$ instead of $H^1$, require a mapping from reference cells to mesh cells which is continuously differentiable across cell interfaces. In this article, we propose an algorithm to…
In this paper we bring together tropical linear algebra and convex 3-dimensional bodies. We show how certain convex 3-dimensional bodies having 20 vertices and 12 facets can be encoded in a $4\times 4$ integer zero-diagonal matrix $A$. A…
Number Decision Diagrams (NDD) provide a natural finite symbolic representation for regular set of integer vectors encoded as strings of digit vectors (least or most significant digit first). The convex hull of the set of vectors…
We consider the discretization of a boundary value problem for a general linear second-order elliptic operator with smooth coefficients using the Virtual Element approach. As in [59] the problem is supposed to have a unique solution, but…
This paper introduces a code generator designed for node-level optimized, extreme-scalable, matrix-free finite element operators on hybrid tetrahedral grids. It optimizes the local evaluation of bilinear forms through various techniques…