skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Numerical Fuzz: A Type System for Rounding Error Analysis
Algorithms operating on real numbers are implemented as floating-point computations in practice, but floatingpoint operations introduceroundoff errorsthat can degrade the accuracy of the result. We propose Λ num , a functional programming language with a type system that can express quantitative bounds on roundoff error. Our type system combines a sensitivity analysis, enforced through a linear typing discipline, with a novel graded monad to track the accumulation of roundoff errors. We prove that our type system is sound by relating the denotational semantics of our language to the exact and floating-point operational semantics. To demonstrate our system, we instantiate Λ num with error metrics proposed in the numerical analysis literature and we show how to incorporate rounding operations that faithfully model aspects of the IEEE 754 floating-point standard. To show that Λ num can be a useful tool for automated error analysis, we develop a prototype implementation for Λ num that infers error bounds that are competitive with existing tools, while often running significantly faster. Finally, we consider semantic extensions of our graded monad to bound error under more complex rounding behaviors, such as non-deterministic and randomized rounding.  more » « less
Award ID(s):
2153916 2219757
PAR ID:
10612719
Author(s) / Creator(s):
;
Publisher / Repository:
Association for Computing Machinery (ACM)
Date Published:
Journal Name:
Proceedings of the ACM on Programming Languages
Volume:
8
Issue:
PLDI
ISSN:
2475-1421
Format(s):
Medium: X Size: p. 1954-1978
Size(s):
p. 1954-1978
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We generalize a magnetogram-matching Biot–Savart law (BSl) from planar to spherical geometry. For a given coronal current densityJ, this law determines the magnetic field B ˜ whose radial component vanishes at the surface. The superposition of B ˜ with a potential field defined by a given surface radial field,Br, provides the entire configuration whereBrremains unchanged by the currents. Using this approach, we (1) upgrade our regularized BSls for constructing coronal magnetic flux ropes (MFRs) and (2) propose a new method for decomposing a measured photospheric magnetic field as B = B pot + B T + B S ˜ , where the potential,Bpot, toroidal,BT, and poloidal, B S ˜ , fields are determined byBr,Jr, and the surface divergence ofB–Bpot, respectively, all derived from magnetic data. OurBTis identical to the one in the alternative Gaussian decomposition by P. W. Schuck et al., whileBpotand B S ˜ are different from their poloidal fields B P < and B P > , which arepotentialin the infinitesimal proximity to the upper and lower side of the surface, respectively. In contrast, our B S ˜ has no such constraints and, asBpotandBT, refers to thesameupper side of the surface. In spite of these differences, for a continuousJdistribution across the surface,Bpotand B S ˜ are linear combinations of B P < and B P > . We demonstrate that, similar to the Gaussian method, our decomposition allows one to identify the footprints and projected surface-location of MFRs in the solar corona, as well as the direction and connectivity of their currents. 
    more » « less
  2. Abstract The sensitivity of urban canopy air temperature ( T a ) to anthropogenic heat flux ( Q A H ) is known to vary with space and time, but the key factors controlling such spatiotemporal variabilities remain elusive. To quantify the contributions of different physical processes to the magnitude and variability of Δ T a / Δ Q A H (where Δ represents a change), we develop a forcing-feedback framework based on the energy budget of air within the urban canopy layer and apply it to diagnosing Δ T a / Δ Q A H simulated by the Community Land Model Urban over the contiguous United States (CONUS). In summer, the median Δ T a / Δ Q A H is around 0.01 K  W  m 2 1 over the CONUS. Besides the direct effect of Q A H on T a , there are important feedbacks through changes in the surface temperature, the atmosphere–canopy air heat conductance ( c a ), and the surface–canopy air heat conductance. The positive and negative feedbacks nearly cancel each other out and Δ T a / Δ Q A H is mostly controlled by the direct effect in summer. In winter, Δ T a / Δ Q A H becomes stronger, with the median value increased by about 20% due to weakened negative feedback associated with c a . The spatial and temporal (both seasonal and diurnal) variability of Δ T a / Δ Q A H as well as the nonlinear response of Δ T a to Δ Q A H are strongly related to the variability of c a , highlighting the importance of correctly parameterizing convective heat transfer in urban canopy models. 
    more » « less
  3. We determine for which exotic tori T \mathcal {T} of dimension d ≠<#comment/> 4 d\neq 4 the homomorphism from the group of isotopy classes of orientation-preserving diffeomorphisms of T \mathcal {T} to S L d ( Z ) \mathrm {SL}_d(\mathbf {Z}) given by the action on the first homology group is split surjective. As part of the proof we compute the mapping class group of all exotic tori T \mathcal {T} that are obtained from the standard torus by a connected sum with an exotic sphere. Moreover, we show that any nontrivial S L d ( Z ) \mathrm {SL}_d(\mathbf {Z}) -action on T \mathcal {T} agrees on homology with the standard action, up to an automorphism of S L d ( Z ) \mathrm {SL}_d(\mathbf {Z}) . When combined, these results in particular show that many exotic tori do not admit any nontrivial differentiable action by S L d ( Z ) \mathrm {SL}_d(\mathbf {Z})
    more » « less
  4. Abstract We study the radial evolution of the inertial-range solar wind plasma turbulence and its anisotropy in the outer heliosphere. We use magnetic field (B) measurements from the Voyager 2 spacecraft for heliocentric distancesRfrom 1 to 33 au. We find that the perpendicular and trace power spectral densities (PSDs) of the magnetic field ( E B and E B Tr ) still follow a Kolmogorov-like spectrum until 33 au. The parallel magnetic field PSD, E B , transits from a power-law index of −2 to −5/3 as the distance crossesR∼ 10 au. The PSD at frequencies 0.01 Hz <f< 0.2 Hz flattens atR> 20 au, gradually approaching anf−1spectrum, probably due to instrument noise. At 0.002 Hz <f< 0.1 Hz, quasi-parallel propagation dominates at 1 au <R< 7 au, with quasi-perpendicular propagation gradually emerging atR> 5 au. ForR> 7 au, oblique propagation becomes the primary mode of propagation. At smaller frequencies off< 0.01 Hz, E B increases with propagation angle at 1 au <R< 5 au, and in contrast decreases with propagation angle atR> 5 au due to the enhanced power level at propagation angles smaller than 20°. Such enhancement may derive from the injection of wave energy from the pickup ion source into the background turbulent cascade, and the injected wave energy is transferred across scales without leaving local enhancements in E B or E B Tr
    more » « less
  5. Abstract The family of transition-metal dipnictides has been of theoretical and experimental interest because this family hosts topological states and extremely large magnetoresistance (MR). Recently, T a A s 2 , a member of this family, has been predicted to support a topological crystalline insulating state. Here, by using high-resolution angle-resolved photoemission spectroscopy (ARPES), we reveal both closed and open pockets in the metallic Fermi surface (FS) and linearly dispersive bands on the ( 2 01 ) surface, along with the presence of extreme MR observed from magneto-transport measurements. A comparison of the ARPES results with first-principles computations shows that the linearly dispersive bands on the measured surface of T a A s 2 are trivial bulk bands. The absence of symmetry-protected surface state on the ( 2 01 ) surface indicates its topologically dark nature. The presence of open FS features suggests that the open-orbit fermiology could contribute to the extremely large MR of T a A s 2
    more » « less