Search NASA⌕ Search

SEARCH · Search NASA

Results for “irregular”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

Investigating the Impact of Power-Take-Off System Parameters and Control Law on a Rotational Wave Energy Converter’s Peak-to-Average Power Ratio Reduction

Due to the irregular nature of real waves, the power captured in a wave energy converter (WEC) system is highly variable. This is an important barrier to the effective use of WECs. To address this challenge, this study focuses on a rotational WEC power-take-off system in which high-speed and high-efficiency generators along with a torque/power smoothing inertia element can be effectively utilized. In the first phase of this study, the U.S. Department of Energy’s reference model 3 (WEC-Sim RM3; two-body point absorber), along with a slider-crank WEC, were integrated for linear to rotational conversion. Relative motion between the float and spar in RM3 was the driving force for this slider-crank WEC, which is connected to a motor/generator set through a gearbox. RM3 geometry was scaled down by 25 times to work within the limits of the physical motor/generator set used in the experimentation. Once the integration in a hardware-in-the-loop simulation environment was successfully completed, data on the peak-to-average power ratio was collected for various wave conditions including regular and irregular waves. The control algorithm designed to keep the system in resonance with waves was able to maintain relatively high speed depending on the specific gear ratio and wave period. Initial results with hardware-in-the-loop simulations reveal that gear ratio and crank radius have a strong impact on the peak-to-average power ratio. In addition, it was found that output power from the generator was maximized at a larger gear ratio, as the crank radius was increased.

50 EE - Wind and Water Power Program - Water (EE-4↗

Accelerating GNNs on GPU Sparse Tensor Cores through N:M Sparsity-Oriented Graph Reordering

Recent advancements in GPU hardware support have introduced the capability to leverage N:M sparse patterns for substantial performance gains. Graphs in Graph Neural Networks (GNNs) are typically sparse, but the sparsity is often irregular, not conforming to such sparse patterns. In this paper, we propose a novel graph reordering algorithm, the first of its kind, to reshape irregular graph data into the N:M structured sparse pattern at the tile level, allowing linear-algebra-based graph operations in GNNs to benefit from the N:M sparse hardware. The optimization is lossless, maintaining the accuracy of GNN. It can remove 98-100\% violations of the N:M sparse patterns at the vector level, and increase the proportion of conforming graphs in SuiteSparse collection from 5-9\% to 88.7-93.5\%. On A100 GPUs, the optimization accelerates Sparse Matrix Matrix (SpMM) by up to 43X (2.3X -- 7.5X on average) and speeds up the key graph operations in GNNs on real graphs by as much as 8.6X (3.5X on average).

artificial intelligence, graph neural networks↗

Development of NDE/NDT Tools for High-Volume & High-Speed Inspection of CFRP Structures in Automotive Manufacturing

Main advantages of the air-coupled ultrasound testing (ACUT) and electromagnetic testing (EMT) techniques for NDE of CFRP composites were non-contact sensing, scalability for high-speed inspection, cost-effectiveness, and non-hazardous operation. Despite these advantages, no systems that would satisfy the project requirements were commercially available. Hence, one of the major efforts of the Michigan State University (MSU) team at the initial stage of the project was to close this technological gap by developing, optimizing, and validating array sensors that would provide sufficient sensitivity, spatial coverage, and resolution for robust defect detection. Optimization of the ACUT and EMT sensor designs was performed using experimentally validated finite element models. Initial experiments using array probes were conducted on relatively flat CFRP samples. In parallel, the MSU team designed and assembled a portable platform with two robotic arms. The robots were equipped with newly designed sensors that enabled high-speed NDE of curved CFRP parts. Presently, the developed robotic platform can be used as a demo/template NDE system, which is easily adaptable to manufacturing environments and in-line NDE. The ACUT NDE system developed by the MSU team used a high-power 4-channel pulser receiver for parallel data acquisition. The array probes were designed by stacking commercially available ACUT transducers, which operated in the frequency range between 100 kHz and 500 kHz. MSU optimized the excitation procedure and developed wave focusing cones so as to reduce the crosstalk between the transducers and to provide higher pulse repletion frequency (PRF). The through-transmission (TT) and single-side access (SSA) inspection modes were successfully implemented. In the TT-ACUT, structural defects in CFRP were detected by passing ultrasonic waves through the test part. Hence, the ACUT transmitters and receivers needed to be placed on the opposite sides of the test part. In the SSA-ACUT, guided waves (GW) were excited in the test part using the transmitters and were sensed by the receivers from the same side. Multi-channel TT-ACUT and SSA-ACUT provided high-speed NDE, and were successfully validated on CFRP test samples with interlaminar delaminations and other embedded defects The EM techniques developed by the MSU team included: 1) eddy current testing (ECT), 2) capacitive imaging (CI) and hybrid dual-mode imaging. In ECT, structural damage was detected in CFRP using coils sensor arrays. In ECT, the excitation magnetic field is generated by passing an alternating current through a coil, which is placed above the test sample. The excitation field penetrates the conductive sample and induces the eddy currents in its transect. In turn, the eddy currents generate the reaction field, which affects the total field sensed by a coil. Hence, the presence of structural flaws will alter the eddy current flow and the picked-up signal. ECT is mostly sensitive to local changes of the electric conductivity of the test sample, and CFRPs are mostly conductive in the direction of carbon fibers. Hence, ECT was well suited for the detection of fiber damage/fiber irregularities. The MSU team developed printed circuit boards (PCB) with coil sensor arrays optimized for NDE of CFRP. Unlike most commercial probes designed for ECT of metallic structures, the MSU array probes were designed for operation in [1-10] MHz frequency range, which was optimal for low-conductive CFRP. Multiple sensing topologies (coil groups excitation/sensing arrangements) were implemented and successfully validated. Capacitive Imaging (CI) technique developed by MSU was complementary to ECT. In contrast to ECT, which was sensitive to local changes of the electrical conductivity, the CI was sensitive to local changes of the dielectric constant. Therefore, CI could provide information about matrix damage/matrix irregularities in CFRP. The MSU CI sensor arrays were made of multiple circular or rectangular open-plate capacitors printed on PCB. Sensors of this type are not commercially available. In addition to ECT and CI, the MSU team developed a hybrid (dual-mode) inductive/capacitive measurement technique that synergistically combined the benefits of inductive and capacitive sensing for rapid NDE of fiber reinforced polymer (FRP) composite structures. Fiber damage and fiber irregularities in FRPs were detected by configuring hybrid sensors as coil sensors. Similarly, matrix damage, matrix irregularities and interlaminar delaminations were detected by configuring hybrid sensors as capacitive sensors. ECT and CI were performed sequentially by means of electronic switching. Hence, eliminating the need for mounting two separate sensor arrays on the probe. Portable robotic platform was developed by MSU for multi-technique high-speed NDE of CFRP test parts. The platform had two 6-axis robots, which enabled inspection of curved parts in approximately a 6×6×6 ft 3 active scan area. On the software side, the MSU team integrated scripts for NDE hardware control with scripts for robot motion control. MSU also implemented automated path planning for the robots, reconstruction of part’s surfaces via stereovision, 3D rendering of inspection data, and image processing algorithms for enhanced defect detection. Automotive composite parts manufactured by Plasan Composites from Phase I were used to validate the ACUT and EMT techniques on representative testbeds. Among those parts were three X-braces for a Dodge Viper, one composite calibration plaque with known defects at known locations, and four other test sections, including sections from a front splitter, a corner section from a composite hood, and a high-pressure RTM panel made using non crimp fabric. Other test samples included CFRP and GFRP calibration plates with fiber/matrix defects fabricated at MSU/CVRC.

36 MATERIALS SCIENCE↗

Radar studies of the equatorial electrojet at three frequencies.

Radar measurements at 16.25, 49.92, and 146.25 MHz of plasma irregularities in the equatorial electrojet were carried out with the objective to investigate the dependence upon wavelength of the properties of the type I and type II irregularities and to see to what extent this dependence is consistent with present theoretical ideas. The equipment and experimental procedure employed is discussed in detail, and the data obtained are analyzed. The following conclusions have been drawn: the Doppler velocity of the type I echoes corresponds to the acoustic velocity rather than the electron streaming velocity; the acoustic velocity varies with wavelength in at least approximate agreement with theory; and the type II irregularities are probably not excited directly, but rather result from nonlinear interactions among larger irregularities.

Balsley, B. B.↗

Observations of artificially produced scintillations using satellite transmissions

The ionospheric modification experiment, utilizing a high-powered transmitter, provides an opportunity to study ionospheric irregularities under relatively known conditions. The irregularities were studied by means of transionospheric signals from the polar-orbiting satellite ESSA 8 transmitting at 137 MHz. These observations show that scintillations occur when the satellite to ground station geometry is such that the ray from the satellite passes through the region in the ionosphere under modification. In general, a cut across the illuminated volume is obtained; thus an active region of about 100 km in diameter is indicated. For the most part, the artificial scintillations appear to be similar to those obtained from naturally occurring irregularities, but a systematic change from 3 to 1 sec in the fluctuation period is usually observed as the satellite traverses from the northern to the southern portions of the active region. The change in period appears to be produced by a systematic change in irregularity scale size from about 4 to about 1 km.

Pope, J. H.↗

Elastohydrodynamic film thickness measurements of artificially produced nonsmooth surfaces

Optical interferometry is used to measure the elastohydrodynamic (EHD) film thickness associated with artificially produced nonsmooth surfaces. The nonsmooth surfaces are produced by modifying the surfaces of highly-polished balls with irregularities in the form of multiple grooves and dents. By closely spacing these irregularities it is possible not only to produce depressions on the surface of the balls but also to generate pseudo asperities. The average roughness wavelength of this artificially-produced, nonsmooth, surface approximates the average fundamental roughness wavelength found on surfaces of some mechanical elements operating under concentrated contact. By comparing the measured film thickness profiles to the stylus traces of the irregularities, it was possible to observe the local deformations associated with micro-EHD pressure generation. In both pure rolling and pure sliding conditions the artificially-produced asperities are deformed and complete separation exists between them and the mating surface. Such findings demonstrate the importance of local surface topography and resulting micro-EHD effects on the film thickness between rough surfaces in concentrated contact. Sliding data are presented which demonstrate a severe constriction, caused by the irregularities, at the exit of the Hertzian region.

Cusano, C.↗

Elastohydrodynamic film thickness measurements of artificially-produced nonsmooth surfaces

Optical interferometry is used to measure the elastohydrodynamic (EHD) film thickness associated with artificially produced, nonsmooth surfaces. The nonsmooth surfaces are produced by modifying the surfaces of highly-polished balls with irregularities in the form of multiple grooves and dents. By closely spacing these irregularities, it is possible not only to produce depressions on the surface of the balls but also to generate pseudo asperities. The average roughness wavelength of this artificially-produced, nonsmooth surface approximates the average fundamental roughness wavelength found on surfaces of some mechanical elements operating under concentrated contact. By comparing the measured film thickness profiles to the stylus traces of the irregularities, it is possible to observe the local deformations associated with micro-EHD pressure generation. In both pure rolling and pure sliding conditions, the artificially-produced 'asperities' are deformed and complete separation exists between them and the mating surface. Such findings demonstrate the importance of local surface topography and resulting micro-EHD effects on the film thickness between rough surfaces in concentrated contact. In addition, sliding data are presented which demonstrate a severe constriction, caused by the irregularities, at the exit of the Hertzian region.

Cusano, C.↗

Phase perturbation measurements through a heated ionosphere

High frequency radiowaves incident on an overdense (i.e., HF-frequency penetration frequency) ionosphere produce electron density irregularities. The effect of such ionospheric irregularities on the phase of UHF-radiowaves was determined. For that purpose the phase of radiowaves originating from celestial radio sources was observed with two antennas. The radiosources were chosen such that the line of sight to at least one of the antennas (usually both) passed through the modified volume of the ionosphere. Observations at 430 MHz and at 2380 MHz indicate that natural irregularities have a much stronger effect on the UHF phase fluctuations than the HF-induced irregularities for presently achieved HF-power densities of 20-80 uW/sq m. It is not clear whether some of the effects observed are the result of HF-modification of the ionosphere. Upper limits on the phase perturbations produced by HF-modification are 10 deg at 2380 MHz and 80 deg at 430 MHz.

Frey, A.↗

The role of in-situ measurements in scintillation modelling

It is shown that the conflicting equatorial scintillation observations obtained from greatly separated ground stations can be organized in the framework of a longitudinal variation of irregularity occurrence, given satellite-borne, in situ measurements of irregularity amplitude of the global F-region irregularity morphology's general features. High-inclination satellite data are used to delineate the morphological features of the polar cap by means of such a method. The lack of diurnal and magnetic control of the irregularity morphology within the low solar flux, northern winter polar cap distinguishes this region from the auroral oval regime. A polar-orbiting communication system sensitive to phase perturbations may observe large differences in the phase-to-amplitude scintillation ratio, as it traverses the auroral oval and proceeds into the polar cap, with its sun-aligned arc system.

Basu, S.↗

A theoretical model of artificial spread F echoes

A theoretical model is developed for artificial spread F echoes elicited by irradiating the ionospheric F region with signals from ground-based HF transmitters. Account is taken of the irregularity polarizations, scale length and the magnetic dip angle of the echo. Ray tracing equations are defined for wave propagation in a horizontally stratified ionosphere which has been bathed with HF signals and therefore contains wavelike structures. Irregularities polarized within the meridian plane are found to cause the spread F echoes, while perpendicularly polarized irregularities do not. A magnetic dip angle of 5 deg must be exceeded for the spread F to become strong. The irregularities need scale lengths exceeding 100 m. Modes of operation are identified for inducing the required spread F echoes using ground-based radiotelescopes.

Kuo, S. P.↗

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a small band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary layer flow. The location, and width, of this most receptive band of wavenumbers corresponds to a relative detuning of O(R sub l.b.(exp -3/8)) with respect to the lower-neutral instability wavenumber at the frequency under consideration, R sub l.b. being the Reynolds number based on a typical boundary-layer thickness at the lower branch of the neutral stability curve. Surface imperfections in the form of discrete mode waviness in this range of wavenumbers lead to initial instability amplitudes which are O(R sub l.b.(exp 3/8)) larger than those caused by a single, isolated roughness element. In contrast, irregularities with a continuous spatial spectrum produce much smaller instability amplitudes, even compared to the isolated case, since the increase due to the resonant nature of the response is more than that compensated for by the asymptotically small band-width of the receptivity process. Analytical expressions for the maximum possible instability amplitudes, as well as their expectation for an ensemble of statistically irregular surfaces with random phase distributions, are also presented.

Choudhari, Meelan↗

Vestibular afferent responses to linear accelerations in the alert squirrel monkey

The spontaneous activity of 40 otolith afferents and 44 canal afferents was recorded in 4 alert, intact squirrel monkeys. Polarization vectors and response properties of otolith afferents were determined during static re-orientations relative to gravity and during Earth-horizontal, sinusoidal, linear oscillations. Canal afferents were tested for sensitivity to linear accelerations. For regular otolith afferents, a significant correlation between upright discharge rate and sensitivity to dynamic acceleration in the horizontal plane was observed. This correlation was not present in irregular units. The sensitivity of otolith afferents to both static tilts and dynamic linear acceleration was much greater in irregularly discharging units than in regularly discharging units. The spontaneous activity and static and dynamic response properties of regularly discharging otolith afferents were similar to those reported in barbiturate-anesthetized squirrel monkeys. Irregular afferents also had similar dynamic response properties when compared to anesthetized monkeys. However, this sample of irregular afferents in alert animals had higher resting discharge rates and greater sensitivity to static tilts. The majority of otolith polarization vectors were oriented near the horizontal in the plane of the utricular maculae; however, directions of maximum sensitivity were different during dynamic and static testing. Canal afferents were not sensitive to static tilts or linear oscillations of the head.

Somps, Christopher J.↗

NGC 5291: Implications for the Formation of Dwarf Galaxies

The possible formation and evolution of dwarf irregular galaxies from material derived from perturbed evolved galaxies is addressed via an H I study of a likely example, the peculiar system NGC 5291. This system, located in the western outskirts of the cluster Abell 3574, contains the lenticular galaxy NGC 5291 which is in close proximity to a disturbed companion and is flanked by an extensive complex of numerous knots extending roughly 4 min north and 4 min south of the galaxy. In an initial optical and radio study, Longmore et al. (1979, MNRAS, 188, 285) showed that these knots have the spectra of vigorous star-forming regions, and suggested that some may in fact be young dwarf irregular galaxies. High resolution 21-cm line observations taken with the VLA are presented here and reveal that the H I distribution associated with this system encompasses not only the entire N-S complex of optical knots, but also forms an incomplete ring or tail that extends approximately 3 min to the west. The H I associated with NGC 5291 itself shows a high velocity range; the Seashell is not detected. The formation mechanism for this unusual system is unclear and two models - a large, low-luminosity ram-swept disk, and a ram-swept interaction-are discussed. The H I in the system contains numerous concentrations, mostly along the N-S arc of the star-forming complexes, which generally coincide with one or more optical knots; the larger H I features contain several x 10(exp 9) solar mass of gas. Each of the knots is compared to a set of criteria designed to determine if these objects are bound against their own internal kinetic energy and are tidally stable relative to the host galaxy. An analysis of the properties of the H I concentrations surrounding the optical star-forming complexes indicates that at least the largest of these is a bound system; it also possesses a stellar component. It is suggested that this object is a genuinely young dwarf irregular galaxy that has evolved from the material associated with the system and that this entire complex contains several proto- or young dwarf irregular galaxies in various stages of development. We are therefore witnessing the early evolution of a number of genuinely young galaxies. Given the evident importance of the NGC 5291 system as a 'nursery' for young galaxies, careful modeling is required if we are to understand this remarkable galaxy.

Malphrus, Benjamin K.↗

Investigations of Remote Plasma Irregularites by Radio Sounding: Applications of the Radio Plasma Imager on IMAGE

The Radio Plasma Imager (RPI) on the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission operates like a radar by transmitting and receiving coherent electromagnetic pulses. Long-range echoes of electromagnetic sounder waves are reflected at remote plasma cutoffs. Thus, analyses of RPI observations will yield the plasma parameters and distances to the remote reflection points. These analyses assume that the reflecting plasma surfaces are cold and are sufficiently smooth that they effectively behave as plane mirrors to the incoming sounder waves, i.e., that geometric optics can be used. The RPI will employ pulse compression and spectral integration techniques, perfected in ground-based ionospheric digital sounders, in order to enhance the signal-to-noise ratio in long-range magnetospheric sounding. When plasma irregularities exist in the remote magnetospheric plasmas that are being probed by the sounder waves, echo signatures may become complicated. Ionospheric sounding experience indicates that while topside sounding echo strengths can actually be enhanced by the presence of irregularities, ground-based sounding indicates that coherent detection techniques can still be employed. In this paper we investigate the plasma conditions that will allow coherent signals to be detected by the RPI and the signatures to be expected, such as scattering and plasma resonances, in the presence of multi-scale irregularities, may possibly have on RPI signals. Sounding of irregular plasma structures in the plasmasphere, plasmapause and magnetopause are also discussed.

Fung, Shing F.↗

Initial Results from the Vector Electric Field Investigation on the C/NOFS Satellite

Initial results are presented from the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite, a mission designed to understand, model, and forecast the presence of equatorial ionospheric irregularities. The VEFI instrument includes a vector DC electric field detector, a fixed-bias Langmuir probe operating in the ion saturation regime, a flux gate magnetometer, an optical lightning detector, and associated electronics including a burst memory. The DC electric field detector has revealed zonal and meridional electric fields that undergo a diurnal variation, typically displaying eastward and outward-directed fields during the day and westward and downward-directed fields at night. In general, the measured DC electric field amplitudes are in the 0.5-2 mV/m range, corresponding to I3 x B drifts of the order of 30-150 m/s. What is surprising is the high degree of large-scale (10's to 100's of km) structure in the DC electric field, particularly at night, regardless of whether well-defined spread-F plasma density depletions are present. The spread-F density depletions and corresponding electric fields that have been detected thus far have displayed a preponderance to appear between midnight and dawn. Associated with the narrow plasma depletions that are detected are broad spectra of electric field and plasma density irregularities for which a full vector set of measurements is available for detailed study. On some occasions, localized regions of low frequency (< 8 Hz) magnetic field broadband irregularities have been detected, suggestive of filamentary currents, although there is no one-to-one correspondence of these waves with the observed plasma density depletions, at least within the data examined thus far. Finally, the data set includes a wide range of ELF/VLF/HF waves corresponding to a variety of plasma waves, in particular banded ELF hiss, whistlers, and lower hybrid wave turbulence triggered by lightning-induced sferics. The VEFI data set represents a treasure trove of measurements that are germane to numerous fundamental aspects of the electrodynamics and irregularities inherent to the Earth's low latitude ionosphere.

Pfaff, R.↗

Flight Test of the Aerojet 7KS-6000 T-27 Jato Rocket Motor

A flight test of the Aero jet Engineering Corporation's 7KS-6000 T-27 Jato rocket motor was conducted at the Langley Pilotless Aircraft Research Station at Wallops Island, Va, to determine the flight performance characteristics of the motor. The flight test imposed an absolute longitudinal acceleration of 9.8 g upon the rocket motor at 2.8 seconds after launching. The total impulse developed by the motor was 43,400 pound-seconds, and the thrusting time was 7.58 seconds. The maximum thrust was 7200 pounds and occurred at 4.8 seconds after launching. No thrust irregularities attributable to effects of the flight longitudinal acceleration were observed. Certain small thrust irregularities occurred in the flight test which appear to correspond to irregularities observed in static tests conducted elsewhere. A hypothesis regarding the origin of these small irregularities is presented.

Bond, Aleck C.↗

Characterization of the as Manufactured Variability in a CVI SiC/SiC Woven Composite

The microstructure of a 2D woven ceramic matrix composite displays significant variability and irregularity. For example, a chemical vapor infiltrated (CVI) SiC/SiC composite exhibits significant amount of porosity arranged in irregular patterns. Furthermore, the fiber tows within a ply frequently have irregular shape and spacing, and the stacked plies are often misaligned and irregularly nested within each other. The goal of an ongoing project at NASA Glenn is to investigate the effects of the complex microstructure and its variability on the properties and the durability of the material. One key requirement for this effort is the development of methods to characterize the distribution in as-fabricated ceramic matrix composite (CMC) microstructures with the objective of correlating microstructural distribution parameters with mechanical performance. An initial task in this effort was to perform quantitative image analysis of polished cross sections of CVI SiC/SiC composite specimens. This analysis provided sample distributions of various microstructural composite features, including: inter-tow pore sizes and shapes, transverse sectioned tow sizes and shapes, and within ply tow spacing. This information can then be used to quantify the effect of extreme values of these features on the local stress state with the goal of determining the likelihood of matrix cracking at a given external load.

Bonacuse, Peter J.↗

A Fine-grained Asynchronous Bulk Synchronous parallelism model for PGAS applications

The Partitioned Global Address Space (PGAS) model is well suited for executing irregular applications on cluster-based systems, due to its efficient support for short, one-sided messages. Separately, the actor model has been gaining popularity as a productive asynchronous message-passing approach for distributed objects in enterprise and cloud computing platforms, typically implemented in languages such as Erlang, Scala or Rust. To the best of our knowledge, there has been no past work on using the actor model to deliver both productivity and scalability to irregular PGAS applications with large number of small messages. In this paper, we introduce a new programming system for PGAS applications, in which point-to-point remote operations can be expressed as fine-grained asynchronous actor messages. In our approach, the programmer does not need to worry about programming complexities related to message aggregation and termination detection. Our approach can be viewed as extending the classical Bulk Synchronous Parallelism model with fine-grained asynchronous communications within a phase or superstep. Here, we believe that our approach offers a desirable point in the productivity-performance space for PGAS applications, with more scalable performance and higher productivity relative to past approaches. Specifically, for seven irregular mini-applications from the Bale Kernels and three graph kernels executed using 2048 cores in the NERSC Cori system, our approach shows geometric mean performance improvements of ≥ 20X relative to standard PGAS versions (UPC and OpenSHMEM) while maintaining comparable productivity to those versions.

97 MATHEMATICS AND COMPUTING↗