Related papers: Mass Matrix Integration Scheme for Fifteen-node We…
The focus in this work is on interior-point methods for inequality-constrained quadratic programs, and particularly on the system of nonlinear equations to be solved for each value of the barrier parameter. Newton iterations give high…
Based on existing data, we argue for a peculiar structure of the neutrino mass matrix, that has a block of relatively large elements--a dominant block. We analyze this ansatz and extract its predictions, assuming that the O(1) coefficients…
The hierarchical structure in the quark masses and mixings allows its ten physical parameters to be most conveniently encoded in mass matrices of the upper triangular form. We classify these matrices in the hierarchical, minimal parameter…
This paper introduces an efficient algorithm for computing the best approximation of a given matrix onto the intersection of linear equalities, inequalities and the doubly nonnegative cone (the cone of all positive semidefinite matrices…
We study the problem of estimating precision matrices in Gaussian distributions that are multivariate totally positive of order two ($\mathrm{MTP}_2$). The precision matrix in such a distribution is an M-matrix. This problem can be…
In this paper, a direct finite element method is proposed for solving interface problems on unfitted meshes. This new method treats the two interface conditions as an $H^{\frac12}(\Gamma)\times H^{-\frac12}(\Gamma)$ pair for the mutual…
Few Shot Segmentation aims to segment novel object classes given only a handful of labeled examples, enabling rapid adaptation with minimal supervision. Current literature crucially lacks a selection method that goes beyond visual…
We study numerical integration of smooth functions defined over the $s$-dimensional unit cube. A recent work by Dick et al. (2019) has introduced so-called extrapolated polynomial lattice rules, which achieve the almost optimal rate of…
A novel finite element scheme is studied for solving the time-dependent Maxwell's equations on unstructured grids efficiently. Similar to the traditional Yee scheme, the method has one degree of freedom for most edges and a sparse inverse…
A new weak Galerkin (WG) finite element method is introduced and analyzed in this paper for the biharmonic equation in its primary form. This method is highly robust and flexible in the element construction by using discontinuous piecewise…
Casting nonlocal problems in variational form and discretizing them with the finite element (FE) method facilitates the use of nonlocal vector calculus to prove well-posedeness, convergence, and stability of such schemes. Employing an FE…
Highly oscillatory integrals of composite type arise in electronic engineering and their calculations is a challenging problem. In this paper, we propose two Gaussian quadrature rules for computing such integrals. The first one is…
As inelastic structures are ubiquitous in many engineering fields, a central task in computational mechanics is to develop accurate, robust and efficient tools for their analysis. Motivated by the poor performances exhibited by standard…
This paper introduces a numerical scheme for time harmonic Maxwell's equations by using weak Galerkin (WG) finite element methods. The WG finite element method is based on two operators: discrete weak curl and discrete weak gradient, with…
The mass matrix for Gauss-Lobatto grid points is usually approximated by Gauss-Lobatto quadrature because this leads to a diagonal matrix that is easy to invert. The exact mass matrix and its inverse are full. We show that the exact mass…
The scalar wave equation is solved using higher order immersed finite elements. We demonstrate that higher order convergence can be obtained. Small cuts with the background mesh are stabilized by adding penalty terms to the weak…
A new approach to building explicit time-marching stencil computation schemes for the transient 2D acoustic wave equation is implemented. It is based on using Poisson's formula and its three time level modification combined with polynomial…
Quasi-Monte Carlo (QMC) methods are equal weight quadrature rules to approximate integrals over the unit cube with respect to the uniform measure. In this paper we discuss QMC integration with respect to general product measures defined on…
We present a multi-scale approach to efficiently embed an ab initio correlated chemical fragment described by its energy-weighted density matrices, and entangled with a wider mean-field many-electron system. This approach, first presented…
The renormalization group equations for the general 2 by 2, complex, neutrino mass matrix are shown to have exact, analytic solutions. Simple formulas are given for the physical mixing angle, complex phase and mass ratio in terms of their…