skip to main content

Title: A Simple Discretization of the Vector Dirichlet Energy

We present a simple and concise discretization of the covariant derivative vector Dirichlet energy for triangle meshes in 3D using Crouzeix‐Raviart finite elements. The discretization is based on linear discontinuous Galerkin elements, and is simple to implement, without compromising on quality: there are two degrees of freedom for each mesh edge, and the sparse Dirichlet energy matrix can be constructed in a single pass over all triangles using a short formula that only depends on the edge lengths, reminiscent of the scalar cotangent Laplacian. Our vector Dirichlet energy discretization can be used in a variety of applications, such as the calculation of Killing fields, parallel transport of vectors, and smooth vector field design. Experiments suggest convergence and suitability for applications similar to other discretizations of the vector Dirichlet energy.

more » « less
Award ID(s):
1717178 1717268
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Date Published:
Journal Name:
Computer Graphics Forum
Page Range / eLocation ID:
p. 81-92
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We develop a new finite difference method for the wave equation in second order form. The finite difference operators satisfy a summation-by-parts (SBP) property. With boundary conditions and material interface conditions imposed weakly by the simultaneous-approximation-term (SAT) method, we derive energy estimates for the semi-discretization. In addition, error estimates are derived by the normal mode analysis. The proposed method is termed as energy-based because of its similarity with the energy-based discontinuous Galerkin method. When imposing the Dirichlet boundary condition and material interface conditions, the traditional SBP-SAT discretization uses a penalty term with a mesh-dependent parameter, which is not needed in our method. Furthermore, numerical dissipation can be added to the discretization through the boundary and interface conditions. We present numerical experiments that verify convergence and robustness of the proposed method.

    more » « less
  2. Binary neural networks (BNNs) substitute complex arithmetic operations with simple bit-wise operations. The binarized weights and activations in BNNs can drastically reduce memory requirement and energy consumption, making it attractive for edge ML applications with limited resources. However, the severe memory capacity and energy constraints of low-power edge devices call for further reduction of BNN models beyond binarization. Weight pruning is a proven solution for reducing the size of many neural network (NN) models, but the binary nature of BNN weights make it difficult to identify insignificant weights to remove. In this paper, we present a pruning method based on latent weight with layer-level pruning sensitivity analysis which reduces the over-parameterization of BNNs, allowing for accuracy gains while drastically reducing the model size. Our method advocates for a heuristics that distinguishes weights by their latent weights, a real-valued vector used to compute the pseudogradient during backpropagation. It is tested using three different convolutional NNs on the MNIST, CIFAR-10, and Imagenette datasets with results indicating a 33%--46% reduction in operation count, with no accuracy loss, improving upon previous works in accuracy, model size, and total operation count. 
    more » « less
  3. Abstract

    A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrumf0(,p)with the Green’s functionG(r,p;p), which describes the monoenergetic spectrum solution in whichf0δ(pp)asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution forG(r,p;p). In this paper, we explore for the first time, solutions for more general and realistic forms forf0(,p). The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering timeτ(r,p)=τ0(p/p0)αin the shear flow region 0 <r<r2, andτ(r,p)=τ0(p/p0)α(r/r2)s, wheres> 0 in the regionr>r2is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distributionψp(r,p;p)that particles observed at (r,p) originated fromr→ ∞ with momentump. The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described.

    more » « less
  4. Abstract

    In this paper, we study the superconvergence of the semi-discrete discontinuous Galerkin (DG) method for linear hyperbolic equations in one spatial dimension. The asymptotic errors in cell averages, downwind point values, and the postprocessed solution are derived for the initial discretization by Gaussian projection (for periodic boundary condition) or Cao projection Cao et al. (SIAM J. Numer. Anal.5, 2555–2573 (2014)) (for Dirichlet boundary condition). We proved that the error constant in the superconvergence of order$$2k+1$$2k+1for DG methods based on upwind-biased fluxes depends on the parity of the orderk. The asymptotic errors are demonstrated by various numerical experiments for scalar and vector hyperbolic equations.

    more » « less
  5. Abstract

    Analytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones.

    more » « less