Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Static Performance of a Wing-Mounted Thrust Reverser Concept

An experimental investigation was conducted in the Jet-Exit Test Facility at NASA Langley Research Center to study the static aerodynamic performance of a wing-mounted thrust reverser concept applicable to subsonic transport aircraft. This innovative engine powered thrust reverser system is designed to utilize wing-mounted flow deflectors to produce aircraft deceleration forces. Testing was conducted using a 7.9%-scale exhaust system model with a fan-to-core bypass ratio of approximately 9.0, a supercritical left-hand wing section attached via a pylon, and wing-mounted flow deflectors attached to the wing section. Geometric variations of key design parameters investigated for the wing-mounted thrust reverser concept included flow deflector angle and chord length, deflector edge fences, and the yaw mount angle of the deflector system (normal to the engine centerline or parallel to the wing trailing edge). All tests were conducted with no external flow and high pressure air was used to simulate core and fan engine exhaust flows. Test results indicate that the wing-mounted thrust reverser concept can achieve overall thrust reverser effectiveness levels competitive with (parallel mount), or better than (normal mount) a conventional cascade thrust reverser system. By removing the thrust reverser system from the nacelle, the wing-mounted concept offers the nacelle designer more options for improving nacelle aero dynamics and propulsion-airframe integration, simplifying nacelle structural designs, reducing nacelle weight, and improving engine maintenance access.

Asbury, Scott C.↗
NASA NTRSDate not supplied

Sparse Linear Algebra Toolkit for Computational Aerodynamics

Finding solutions to sparse linear systems of equations is an essential step in Computational Engineering applications of interest to NASA. Linear systems of equations are composed and solved in almost every computational engineering application. The characteristics of linear systems vary greatly from one application to another. Accordingly, there are a wide variety of methods for the solution of linear systems of equations. The operations and methods prepared by the authors are focused on linear systems of interest to NASA, primarily those associated with Computational Fluid Dynamics (CFD), Aeroelasticity, and Aeroacoustics. The Sparse Linear Algebra Toolkit (SLAT) is a coordinated collection of software featuring operations, methods, and data structures that are useful when solving sparse linear systems of equations on modern computer architectures. The implemented operations and methods are designed and tuned for parallelism in shared memory, in distributed memory, and across the hybrid combination of distributed-shared memory. The toolkit includes novel methods and implementations for modern architectures and facilitates development of new approaches for meeting NASA’s evolving computational engineering challenges using evolving computer architectures that are not available in vendor libraries. In this paper, significant features and interfaces within SLAT are presented and verified for simulations performed with NASA’s CFD solver, FUN3D. The runtime and scaling performance of the Generalized Minimum Residual (GMRES) method implemented in SLAT is analyzed for the linear subproblems within the solution of turbulent Navier-Stokes equations employed in the simulation of high-lift configurations. Prior to this work, the SPARSKIT GMRES implementation was the only Krylov subspace method available within FUN3D. A strong scaling study shows the SLAT GMRES implementation facilitates accurate Reynolds-averaged Navier-Stokes CFD solutions between 15% and 56% faster than the SPARSKIT GMRES implementation.

Stephen L Wood↗
NASA NTRS2019

Thermal design and validation of Mars 202 gas Dust Removal Tool (gDRT)

As part of the science goals for the planned Mars 2020 mission, two instruments, PIXL and SHERLOC, intend to study the Mars surface at a close scale. These two instruments are planned to be used on smooth rock surfaces that are free of dust and other particles. Such surfaces are prepared by using a drill with an abrading drill bit; residual dust and particles are then blown away using a compressed gas system—the Gas Dust Removal Tool (gDRT). Early in the design process a risk of valve leakage below the vendor’s -20C rating was identified. To mitigate this risk a parallel path was implemented: 1) qualify the valves to operate at -135C, and 2) develop thermal control capability to elevated temperatures in case qualification efforts are unsuccessful. While qualification efforts are ongoing, it is expected that the valves will be qualified to an operating temperature of -135C. In the event that qualification is not successful, thermal control has been realized via a thermally isolated valve configuration and thermostatically controlled heaters. Thermal testing of gDRT has validated this design and enabled thermal model correlation in order to provide more accurate survival energy predictions. The resulting energy consumption for heating the valves at the Jezero landing site represents a minimal impact to the Rover energy budget.

Novak, Keith↗
NASA NTRS2021

Collaborative Pose Estimation of An Unknown Target Using Multiple Spacecraft

A reliable method for pose estimation of an unknown and uncooperative space target using monocular vision remains an open problem. Vision-based pose determination can be challenging in case of unfavorable illumination, time-varying conditions due to rotational motion and relative orbit, and scale ambiguity resolution. To address these challenges, we propose a novel collaborative pose determination algorithm called Multi- Spacecraft Simultaneous Estimation of Pose and Shape algorithm or M-SEPS.Within M-SEPS, a team of chaser spacecraft, each equipped with a monocular camera, exchange information over a local network to jointly estimate the relative kinematic state of the target and its sparse shape landmarks. In this approach, each spacecraft processes its own images and observes particular target landmarks in parallel and in a distributed fashion. Then, the local network is exploited by the spacecraft to share their consensus proposals and aggregate them to achieve the joint estimate. We validate our algorithm using simulations of relative orbits and observations, captured by each chaser spacecraft. To the best of the authors’ knowledge, this is the first cooperative, vision-based algorithm for estimating the pose and shape of a space object for an arbitrary number of spacecraft.

Chung, Soon-Jo↗
NASA NTRS1994

The determination of interplanetary magnetic field polarities around sector boundaries using E greater than 2 keV electrons

The determination of the polarities of interplanetary magnetic fields (whether the field direction is outward from or inward toward the sun) has been based on a comparison of observed field directions with the nominal Parker spiral angle. These polarities can be mapped back to the solar source field polarities. This technique fails when field directions deviate substantially from the Parker angle or when fields are substantially kinked. We introduce a simple new technique to determine the polarities of interplanetary fields using E greater than 2 keV interplanetary electrons which stream along field lines away from the sun. Those electrons usually show distinct unidirectional pitch-angle anisotropies either parallel or anti-parallel to the field. Since the electron flow direction is known to be outward from the sun, the anisotropies parallel to the field indicate outward-pointing, positive-polarity fields, and those anti-parallel indicate inward-pointing, negative-polarity fields. We use data from the UC Berkeley electron experiment on the International Sun Earth Explorer 3 (ISSE-3) spacecraft to compare the field polarities deduced from the electron data, Pe (outward or inward), with the polarities inferred from field directions, Pd, around two sector boundaries in 1979. We show examples of large (greater than 100 deg) changes in azimuthal field direction Phi over short (less than 1 hr) time scales, some with and some without reversals in Pe. The latter cases indicate that such large directional changes can occur in unipolar structures. On the other hand, we found an example of a change in Pe during which the rotation in Phi was less than 30 deg, indicating polarity changes in nearly unidirectional structures. The field directions are poor guides to the polarities in these cases.

Kahler, S.↗
NASA NTRS1994

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗
NASA NTRS2004

Particle Acceleration and Magnetic Field Generation in Electron-Positron Relativistic Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at the comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. Additionally, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportional to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. While some Fermi acceleration may occur at the jet front, the majority of electron and positron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform: small-scale magnetic fields which contribute to the electron's (positron's) transverse deflection behind the jet head. This small scale magnetic field structure is appropriate to the generation of jitter radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron radiation calculated assuming a a uniform magnetic field. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

Nishikawa, K.-L.↗
NASA NTRS1976

Solid-propellant rocket motor internal ballistic performance variation analysis, phase 2

The Monte Carlo method was used to investigate thrust imbalance and its first time derivative throughtout the burning time of pairs of solid rocket motors firing in parallel. Results obtained compare favorably with Titan 3 C flight performance data. Statistical correlations of the thrust imbalance at various times with corresponding nominal trace slopes suggest several alternative methods of predicting thrust imbalance. The effect of circular-perforated grain deformation on internal ballistics is discussed, and a modified design analysis computer program which permits such an evaluation is presented. Comparisons with SRM firings indicate that grain deformation may account for a portion of the so-called scale factor on burning rate between large motors and strand burners or small ballistic test motors. Thermoelastic effects on burning rate are also investigated. Burning surface temperature is calculated by coupling the solid phase energy equation containing a strain rate term with a model of gas phase combustion zone using the Zeldovich-Novozhilov technique. Comparisons of solutions with and without the strain rate term indicate a small but possibly significant effect of the thermoelastic coupling.

Sforzini, R. H.↗
NASA NTRS2005

Particle Acceleration and Magnetic Field Generation in Electron-Positron Relativistic Shocks

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel, and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the nonlinear fluctuation amplitudes of densities, currents, and electric and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at a comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. In addition, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by tine Weibel instability scale proportionally to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. While some Fermi acceleration may occur at the jet front, the majority of electron and positron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields, which contribute to the electron s (positron s) transverse deflection behind the jet head. This small- scale magnetic field structure is appropriate to the generation of "jitter" radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron radiation calculated assuming a uniform magnetic field. The jitter radiation resulting from small-scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

Nishikawa, K.-I.↗
NASA NTRS2005

Effect of Free Jet on Refraction and Noise

This article investigates the role of a free jet on the sound radiated from a jet. In particular, the role of an infinite wind tunnel, which simulates the forward flight condition, is compared to that of a finite wind tunnel. The second configuration is usually used in experiments, where the microphones are located in a static ambient medium far outside the free jet. To study the effect of the free jet on noise, both propagation and source strength need to be addressed. In this work, the exact Green's function in a locally parallel flow is derived for a simulated flight case. Numerical examples are presented that show a reduction in the magnitude of the Green's function in the aft arc and an increase in the forward arc for the simulated flight condition. The effect of finite wind tunnel on refraction is sensitive to the source location and is most pronounced in the aft arc. A Reynolds-averaged Navier-Stokes solution (RANS) yields the required mean flow and turbulence scales that are used in the jet mixing noise spectrum calculations. In addition to the sound/flow interaction, the separate effect of source strength and elongation of the noise-generating region of the jet in a forward flight is studied. Comparisons are made with experiments for the static and finite tunnel cases. Finally, the standard free-jet shear corrections that convert the finite wind tunnel measurements to an ideal wind tunnel arrangement are evaluated.

Khavaran, Abbas↗
NASA NTRS2012

Surface Textures and Features Indicative of Endogenous Growth at the McCartys Flow Field, NM, as an Analog to Martian Volcanic Plains

Basaltic lavas typically form channels or tubes, which are recognized on the Earth and Mars. Although largely unrecognized in the planetary community, terrestrial inflated sheet flows also display morphologies that share many commonalities with lava plains on Mars. The McCartys lava flow field is among the youngest (approx.3000 yrs) basaltic flows in the continental United States. The southwest sections of the flow displays smooth, flat-topped plateaus with irregularly shaped pits and hummocky inter-plateau units that form a polygonal surface. Plateaus are typically elongate in map view, up to 20 m high and display lineations within the glassy crust. Lineated surfaces occasionally display small < 1m diameter lava coils. Lineations are generally straight and parallel each other, sometimes for over 100 meters. The boundaries between plateaus and depressions are also lineated and tilted to angles sometimes approaching vertical. Plateau-parallel cracks, sometimes containing squeeze-ups, mark the boundary between tilted crust and plateau. Some plateau depressions display level floors with hummocky surfaces, while some are bowl shaped with floors covered in broken lava slabs. The lower walls of pits sometimes display lateral, sagged lava wedges. Infrequently, pit floors display the upper portion of a tumulus from an older flow. In some places the surface crust has been disrupted forming a slabby texture. Slabs are typically on the scale of a meter or less across and no less than 7-10 cm thick. The slabs preserve the lineated textures of the undisturbed plateau crust. It appears that this style of terrain represents the emplacement of an extensive sheet that experiences inflation episodes within preferred regions where lateral spreading of the sheet is inhibited, thereby forming plateaus. Rough surfaces represent inflation-related disruption of pahoehoe lava and not a a lava. Depressions are often the result of non-inflation and can be clearly identified by lateral squeeze-outs along the pit walls that form when the rising crust exposes the still liquid core of the sheet. The plains of Tharsis and Elysium, Mars, display many analogous features

Bleacher, Jacob E.↗
NASA NTRS2016

Van Allen Probes Observations of Magnetic Field Dipolarization and Its Associated O+ Flux Variations in the Inner Magnetosphere at L 6.6

We investigate the magnetic field dipolarization in the inner magnetosphere and its associated ion flux variations, using the magnetic field and energetic ion flux data acquired by the Van Allen Probes. From a study of 74 events that appeared at L= 4.5-6.6 between 1 October 2012 and 31 October 2013, we reveal the following characteristics of the dipolarization in the inner magnetosphere: (1) its time scale is approximately 5 min; (2) it is accompanied by strong magnetic fluctuations that have a dominant frequency close to the O+ gyrofrequency; (3) ion fluxes at 20-50 keV are simultaneously enhanced with larger magnitudes for O+ than for H+; (4) after a few minutes of the dipolarization, the flux enhancement at 0.1-5 keV appears with a clear energy-dispersion signature only for O+; and (5) the energy-dispersed O+ flux enhancement appears in directions parallel or antiparallel to the magnetic field. From these characteristics, we discuss possible mechanisms that can provide selective acceleration to O+ ions at > 20 keV. We conclude that O+ ions at L= 5.4-6.6 undergo nonadiabatic local acceleration caused by oscillating electric field associated with the magnetic fluctuations and/or adiabatic convective transport from the plasma sheet to the inner magnetosphere by the impulsive electric field. At L= 4.5-5.4, however, only the former acceleration is plausible. We also conclude that the field-aligned energy-dispersed O+ ions at 0.1-5 keV originate from the ionosphere and are extracted nearly simultaneously to the onset of the dipolarization.

Magnetosphere↗
NASA NTRS1990

Structural properties of laminated Douglas fir/epoxy composite material

This publication contains a compilation of static and fatigue strength data for laminated-wood material made from Douglas fir and epoxy. Results of tests conducted by several organizations are correlated to provide insight into the effects of variables such as moisture, size, lamina-to-lamina joint design, wood veneer grade, and the ratio of cyclic stress to steady stress during fatigue testing. These test data were originally obtained during development of wood rotor blades for large-scale wind turbines of the horizontal-axis (propeller) configuration. Most of the strength property data in this compilation are not found in the published literature. Test sections ranged from round cylinders 2.25 in. in diameter to rectangular slabs 6 by 24 in. in cross section and approximately 30 ft. long. All specimens were made from Douglas fir veneers 0.10 in. thick, bonded together with the WEST epoxy system developed for fabrication and repair of wood boats. Loading was usually parallel to the grain. Size effects (reduction in strength with increase in test volume) are observed in some of the test data, and a simple mathematical model is presented that includes the probability of failure. General characteristics of the wood/epoxy laminate are discussed, including features that make it useful for a wide variety of applications.

Spera, David A.↗
NASA NTRS1998

A Theoretical and Observational Study of Large and Small-Scale Dynamics in the Mesosphere and Lower Thermosphere in Support of TIMED

The long-term objectives of this IDS component of the TIMED mission are to define, through a combination of modeling and data analysis efforts, the influences of large- and small-scale waves (tidal, planetary, and gravity waves) on the circulation, structure, composition, and variability of the mesosphere and lower thermosphere. Particular emphasis is aimed at the TIMED core science objectives of defining the structure and energetics of the TIMED core region, approx. 60 to 180 km. To do so, however, it is also essential to understand the dynamical influences and contributions because of their considerable direct and indirect influences on structure and energetics. To accomplish these goals, we anticipate two parallel modeling efforts and a significant data analysis effort following TIMED launch. The modeling efforts will address both: the residual circulation response to mean and variable wave forcing of the MLT through development of a spherical model addressing balanced flow dynamics and transience and a more comprehensive effort addressing gravity wave propagation, interactions, nonlinearity, and spectral evolution and transports in canonical and realistic environments. The data analysis effort is envisioned to involve application of ground-based measurements and TIMED data analyses to a quantification of gravity wave forcing and effects at one or more sites with comprehensive measurement suites and extension of the inferred forcing globally, based on distributed ground-based measurements. Research results and publications resulting from the research are included.

Fritts, David C.↗
NASA NTRS2015

20-mN Variable Specific Impulse (Isp) Colloid Thruster

Busek Company, Inc., has designed and manufactured an electrospray emitter capable of generating 20 mN in a compact package (7x7x1.7 in). The thruster consists of nine porous-surface emitters operating in parallel from a common propellant supply. Each emitter is capable of supporting over 70,000 electrospray emission sites with the plume from each emitter being accelerated through a single aperture, eliminating the need for individual emission site alignment to an extraction grid. The total number of emission sites during operation is expected to approach 700,000. This Phase II project optimized and characterized the thruster fabricated during the Phase I effort. Additional porous emitters also were fabricated for full-scale testing. Propellant is supplied to the thruster via existing feed-system and microvalve technology previously developed by Busek, under the NASA Space Technology 7's Disturbance Reduction System (ST7-DRS) mission and via follow-on electric propulsion programs. This project investigated methods for extending thruster life beyond the previously demonstrated 450 hours. The life-extending capabilities will be demonstrated on a subscale version of the thruster.

Demmons, Nathaniel↗
NASA NTRS1989

Computational chaos in massively parallel neural networks

A fundamental issue which directly impacts the scalability of current theoretical neural network models to massively parallel embodiments, in both software as well as hardware, is the inherent and unavoidable concurrent asynchronicity of emerging fine-grained computational ensembles and the possible emergence of chaotic manifestations. Previous analyses attributed dynamical instability to the topology of the interconnection matrix, to parasitic components or to propagation delays. However, researchers have observed the existence of emergent computational chaos in a concurrently asynchronous framework, independent of the network topology. Researcher present a methodology enabling the effective asynchronous operation of large-scale neural networks. Necessary and sufficient conditions guaranteeing concurrent asynchronous convergence are established in terms of contracting operators. Lyapunov exponents are computed formally to characterize the underlying nonlinear dynamics. Simulation results are presented to illustrate network convergence to the correct results, even in the presence of large delays.

Barhen, Jacob↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the loop unrolling technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗