A Mixed Real and Floating-Point Solver

Reasoning about mixed real and floating-point constraints is essential for developing accurate analysis tools for floating-point pro- grams. This paper presents FPRoCK, a prototype tool for solving mixed real and floating-point formulas. FPRoCK transforms a mixed formula into an equisatisfiable one over the reals. This formula is then solved using an off-the-shelf SMT solver. FPRoCK is also integrated with the PRECiSA static analyzer, which computes a sound estimation of the round-off error of a floating-point program. It is used to detect infeasible computational paths, thereby improving the accuracy of PRECiSA

Provably Correct Floating-Point Implementation of a Point-In-Polygon Algorithm

The problem of determining whether or not a point lies inside a given polygon occurs in many applications. In air traffic management concepts, a correct solution to the point-in-polygon problem is critical to geofencing systems for Unmanned Aerial Vehicles and in weather avoidance applications. Many mathematical methods can be used to solve the point-in-polygon problem. Unfortunately, a straightforward floating- point implementation of these methods can lead to incorrect results due to round-off errors. In particular, these errors may cause the control flow of the program to diverge with respect to the ideal real-number algorithm. This divergence potentially results in an incorrect point-in- polygon determination even when the point is far from the edges of the polygon. This paper presents a provably correct implementation of a point-in-polygon method that is based on the computation of the winding number. This implementation is mechanically generated from a source- to-source transformation of the ideal real-number specification of the algorithm. The correctness of this implementation is formally verified within the Frama-C analyzer, where the proof obligations are discharged using the Prototype Verification System (PVS)

Accelerating Dust Temperature Calculations with Graphics Processing Units

When calculating the infrared spectral energy distributions (SEDs) of galaxies in radiation-transfer models, the calculation of dust grain temperatures is generally the most time-consuming part of the calculation. Because of its highly parallel nature, this calculation is perfectly suited for massively parallel general-purpose Graphics Processing Units (GPUs). This paper presents an implementation of the calculation of dust grain equilibrium temperatures on GPUs in the Monte-Carlo radiation transfer code Sunrise, using the CUDA API. The GPU can perform this calculation 69 times faster than the 8 CPU cores, showing great potential for accelerating calculations of galaxy SEDs.Comment: 7 pages, 2 figures, accepted to New Astronomy. Minor updates to text and performance based on feedback from refere

Engineering model 8-cm thruster subsystem

An Engineering Model (EM) 8 cm Ion Thruster Propulsion Subsystem was developed for operation at a thrust level 5 mN (1.1 mlb) at a specific impulse 1 sub sp = 2667 sec with a total system input power P sub in = 165 W. The system dry mass is 15 kg with a mercury-propellant-reservoir capacity of 8.75 kg permitting uninterrupted operation for about 12,500 hr. The subsystem can be started from a dormant condition in a time less than or equal to 15 min. The thruster has a design lifetime of 20,000 hr with 10,000 startup cycles. A gimbal unit is included to provide a thrust vector deflection capability of + or - 10 degrees in any direction from the zero position. The EM subsystem development program included thruster optimization, power-supply circuit optimization and flight packaging, subsystem integration, and subsystem acceptance testing including a cyclic test of the total propulsion package

The Political Economy of Exchange Rate Policies in Latin America and the Caribbean: The Case of Peru

The main objective of this document is to identify structural long-term factors that can help to understand the exchange rate policy in Peru during the last fifty years. Rather than identifying particular causes that could help to understand specific decisions at some point in time, the emphasis is on those general political economy arguments that explain the main trends and changes in exchange rate policy over long periods of time. Thus, two main things need to be understood: the choice of the exchange rate regime and its level.

Implementation of the Hierarchical Reference Theory for simple one-component fluids

Combining renormalization group theoretical ideas with the integral equation approach to fluid structure and thermodynamics, the Hierarchical Reference Theory is known to be successful even in the vicinity of the critical point and for sub-critical temperatures. We here present a software package independent of earlier programs for the application of this theory to simple fluids composed of particles interacting via spherically symmetrical pair potentials, restricting ourselves to hard sphere reference systems. Using the hard-core Yukawa potential with z=1.8/sigma for illustration, we discuss our implementation and the results it yields, paying special attention to the core condition and emphasizing the decoupling assumption's role.Comment: RevTeX, 16 pages, 2 figures. Minor changes, published versio