New methods for interpretation of magnetic vector and gradient tensor data I: eigenvector analysis and the normalised source strength

NASA Astrophysics Data System (ADS)

Clark, David A.

2012-09-01

Acquisition of magnetic gradient tensor data is likely to become routine in the near future. New methods for inverting gradient tensor surveys to obtain source parameters have been developed for several elementary, but useful, models. These include point dipole (sphere), vertical line of dipoles (narrow vertical pipe), line of dipoles (horizontal cylinder), thin dipping sheet, and contact models. A key simplification is the use of eigenvalues and associated eigenvectors of the tensor. The normalised source strength (NSS), calculated from the eigenvalues, is a particularly useful rotational invariant that peaks directly over 3D compact sources, 2D compact sources, thin sheets and contacts, and is independent of magnetisation direction. In combination the NSS and its vector gradient determine source locations uniquely. NSS analysis can be extended to other useful models, such as vertical pipes, by calculating eigenvalues of the vertical derivative of the gradient tensor. Inversion based on the vector gradient of the NSS over the Tallawang magnetite deposit obtained good agreement between the inferred geometry of the tabular magnetite skarn body and drill hole intersections. Besides the geological applications, the algorithms for the dipole model are readily applicable to the detection, location and characterisation (DLC) of magnetic objects, such as naval mines, unexploded ordnance, shipwrecks, archaeological artefacts, and buried drums.

Pixel-By Estimation of Scene Motion in Video

NASA Astrophysics Data System (ADS)

Tashlinskii, A. G.; Smirnov, P. V.; Tsaryov, M. G.

2017-05-01

The paper considers the effectiveness of motion estimation in video using pixel-by-pixel recurrent algorithms. The algorithms use stochastic gradient decent to find inter-frame shifts of all pixels of a frame. These vectors form shift vectors' field. As estimated parameters of the vectors the paper studies their projections and polar parameters. It considers two methods for estimating shift vectors' field. The first method uses stochastic gradient descent algorithm to sequentially process all nodes of the image row-by-row. It processes each row bidirectionally i.e. from the left to the right and from the right to the left. Subsequent joint processing of the results allows compensating inertia of the recursive estimation. The second method uses correlation between rows to increase processing efficiency. It processes rows one after the other with the change in direction after each row and uses obtained values to form resulting estimate. The paper studies two criteria of its formation: gradient estimation minimum and correlation coefficient maximum. The paper gives examples of experimental results of pixel-by-pixel estimation for a video with a moving object and estimation of a moving object trajectory using shift vectors' field.

Non-singular spherical harmonic expressions of geomagnetic vector and gradient tensor fields in the local north-oriented reference frame

NASA Astrophysics Data System (ADS)

Du, J.; Chen, C.; Lesur, V.; Wang, L.

2014-12-01

General expressions of magnetic vector (MV) and magnetic gradient tensor (MGT) in terms of the first- and second-order derivatives of spherical harmonics at different degrees and orders, are relatively complicated and singular at the poles. In this paper, we derived alternative non-singular expressions for the MV, the MGT and also the higher-order partial derivatives of the magnetic field in local north-oriented reference frame. Using our newly derived formulae, the magnetic potential, vector and gradient tensor fields at an altitude of 300 km are calculated based on a global lithospheric magnetic field model GRIMM_L120 (version 0.0) and the main magnetic field model of IGRF11. The corresponding results at the poles are discussed and the validity of the derived formulas is verified using the Laplace equation of the potential field.

Multi-color incomplete Cholesky conjugate gradient methods for vector computers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Poole, E.L.

1986-01-01

This research is concerned with the solution on vector computers of linear systems of equations. Ax = b, where A is a large, sparse symmetric positive definite matrix with non-zero elements lying only along a few diagonals of the matrix. The system is solved using the incomplete Cholesky conjugate gradient method (ICCG). Multi-color orderings are used of the unknowns in the linear system to obtain p-color matrices for which a no-fill block ICCG method is implemented on the CYBER 205 with O(N/p) length vector operations in both the decomposition of A and, more importantly, in the forward and back solvesmore » necessary at each iteration of the method. (N is the number of unknowns and p is a small constant). A p-colored matrix is a matrix that can be partitioned into a p x p block matrix where the diagonal blocks are diagonal matrices. The matrix is stored by diagonals and matrix multiplication by diagonals is used to carry out the decomposition of A and the forward and back solves. Additionally, if the vectors across adjacent blocks line up, then some of the overhead associated with vector startups can be eliminated in the matrix vector multiplication necessary at each conjugate gradient iteration. Necessary and sufficient conditions are given to determine which multi-color orderings of the unknowns correspond to p-color matrices, and a process is indicated for choosing multi-color orderings.« less

Navier-Stokes dynamics on a differential one-form

NASA Astrophysics Data System (ADS)

Story, Troy L.

2006-11-01

After transforming the Navier-Stokes dynamic equation into a characteristic differential one-form on an odd-dimensional differentiable manifold, exterior calculus is used to construct a pair of differential equations and tangent vector(vortex vector) characteristic of Hamiltonian geometry. A solution to the Navier-Stokes dynamic equation is then obtained by solving this pair of equations for the position x^k and the conjugate to the position bk as functions of time. The solution bk is shown to be divergence-free by contracting the differential 3-form corresponding to the divergence of the gradient of the velocity with a triple of tangent vectors, implying constraints on two of the tangent vectors for the system. Analysis of the solution bk shows it is bounded since it remains finite as | x^k | ->,, and is physically reasonable since the square of the gradient of the principal function is bounded. By contracting the characteristic differential one-form with the vortex vector, the Lagrangian is obtained.

A Genealogy of Convex Solids Via Local and Global Bifurcations of Gradient Vector Fields

NASA Astrophysics Data System (ADS)

Domokos, Gábor; Holmes, Philip; Lángi, Zsolt

2016-12-01

Three-dimensional convex bodies can be classified in terms of the number and stability types of critical points on which they can balance at rest on a horizontal plane. For typical bodies, these are non-degenerate maxima, minima, and saddle points, the numbers of which provide a primary classification. Secondary and tertiary classifications use graphs to describe orbits connecting these critical points in the gradient vector field associated with each body. In previous work, it was shown that these classifications are complete in that no class is empty. Here, we construct 1- and 2-parameter families of convex bodies connecting members of adjacent primary and secondary classes and show that transitions between them can be realized by codimension 1 saddle-node and saddle-saddle (heteroclinic) bifurcations in the gradient vector fields. Our results indicate that all combinatorially possible transitions can be realized in physical shape evolution processes, e.g., by abrasion of sedimentary particles.

Non-singular spherical harmonic expressions of geomagnetic vector and gradient tensor fields in the local north-oriented reference frame

NASA Astrophysics Data System (ADS)

Du, J.; Chen, C.; Lesur, V.; Wang, L.

2015-07-01

General expressions of magnetic vector (MV) and magnetic gradient tensor (MGT) in terms of the first- and second-order derivatives of spherical harmonics at different degrees/orders are relatively complicated and singular at the poles. In this paper, we derived alternative non-singular expressions for the MV, the MGT and also the third-order partial derivatives of the magnetic potential field in the local north-oriented reference frame. Using our newly derived formulae, the magnetic potential, vector and gradient tensor fields and also the third-order partial derivatives of the magnetic potential field at an altitude of 300 km are calculated based on a global lithospheric magnetic field model GRIMM_L120 (GFZ Reference Internal Magnetic Model, version 0.0) with spherical harmonic degrees 16-90. The corresponding results at the poles are discussed and the validity of the derived formulas is verified using the Laplace equation of the magnetic potential field.

Contaminant Gradients in Trees: Directional Tree Coring Reveals Boundaries of Soil and Soil-Gas Contamination with Potential Applications in Vapor Intrusion Assessment.

PubMed

Wilson, Jordan L; Samaranayake, V A; Limmer, Matthew A; Schumacher, John G; Burken, Joel G

2017-12-19

Contaminated sites pose ecological and human-health risks through exposure to contaminated soil and groundwater. Whereas we can readily locate, monitor, and track contaminants in groundwater, it is harder to perform these tasks in the vadose zone. In this study, tree-core samples were collected at a Superfund site to determine if the sample-collection location around a particular tree could reveal the subsurface location, or direction, of soil and soil-gas contaminant plumes. Contaminant-centroid vectors were calculated from tree-core data to reveal contaminant distributions in directional tree samples at a higher resolution, and vectors were correlated with soil-gas characterization collected using conventional methods. Results clearly demonstrated that directional tree coring around tree trunks can indicate gradients in soil and soil-gas contaminant plumes, and the strength of the correlations were directly proportionate to the magnitude of tree-core concentration gradients (spearman's coefficient of -0.61 and -0.55 in soil and tree-core gradients, respectively). Linear regression indicates agreement between the concentration-centroid vectors is significantly affected by in planta and soil concentration gradients and when concentration centroids in soil are closer to trees. Given the existing link between soil-gas and vapor intrusion, this study also indicates that directional tree coring might be applicable in vapor intrusion assessment.

Contaminant gradients in trees: Directional tree coring reveals boundaries of soil and soil-gas contamination with potential applications in vapor intrusion assessment

USGS Publications Warehouse

Wilson, Jordan L.; Samaranayake, V.A.; Limmer, Matthew A.; Schumacher, John G.; Burken, Joel G.

2017-01-01

Contaminated sites pose ecological and human-health risks through exposure to contaminated soil and groundwater. Whereas we can readily locate, monitor, and track contaminants in groundwater, it is harder to perform these tasks in the vadose zone. In this study, tree-core samples were collected at a Superfund site to determine if the sample-collection location around a particular tree could reveal the subsurface location, or direction, of soil and soil-gas contaminant plumes. Contaminant-centroid vectors were calculated from tree-core data to reveal contaminant distributions in directional tree samples at a higher resolution, and vectors were correlated with soil-gas characterization collected using conventional methods. Results clearly demonstrated that directional tree coring around tree trunks can indicate gradients in soil and soil-gas contaminant plumes, and the strength of the correlations were directly proportionate to the magnitude of tree-core concentration gradients (spearman’s coefficient of -0.61 and -0.55 in soil and tree-core gradients, respectively). Linear regression indicates agreement between the concentration-centroid vectors is significantly affected by in-planta and soil concentration gradients and when concentration centroids in soil are closer to trees. Given the existing link between soil-gas and vapor intrusion, this study also indicates that directional tree coring might be applicable in vapor intrusion assessment.

A Low Frequency Electromagnetic Sensor for Underwater Geo-Location

DTIC Science & Technology

2011-05-01

used a set of commercially available fluxgate magnetometers to measure the magnetic field gradients associated with a magnetic dipole transmitter...insight into the operational capabilities of commercial fluxgate sensors. Figure 42. Applied Physics Systems 1540 magnetometer ...a magnetic field gradient receiver array. Highest quality gradient estimates were achieved with three vector magnetometers equally spaced and

Vector 33: A reduce program for vector algebra and calculus in orthogonal curvilinear coordinates

NASA Astrophysics Data System (ADS)

Harper, David

1989-06-01

This paper describes a package with enables REDUCE 3.3 to perform algebra and calculus operations upon vectors. Basic algebraic operations between vectors and between scalars and vectors are provided, including scalar (dot) product and vector (cross) product. The vector differential operators curl, divergence, gradient and Laplacian are also defined, and are valid in any orthogonal curvilinear coordinate system. The package is written in RLISP to allow algebra and calculus to be performed using notation identical to that for operations. Scalars and vectors can be mixed quite freely in the same expression. The package will be of interest to mathematicians, engineers and scientists who need to perform vector calculations in orthogonal curvilinear coordinates.

Numerical simulations of short-mixing-time double-wave-vector diffusion-weighting experiments with multiple concatenations on whole-body MR systems

NASA Astrophysics Data System (ADS)

Finsterbusch, Jürgen

2010-12-01

Double- or two-wave-vector diffusion-weighting experiments with short mixing times in which two diffusion-weighting periods are applied in direct succession, are a promising tool to estimate cell sizes in the living tissue. However, the underlying effect, a signal difference between parallel and antiparallel wave vector orientations, is considerably reduced for the long gradient pulses required on whole-body MR systems. Recently, it has been shown that multiple concatenations of the two wave vectors in a single acquisition can double the modulation amplitude if short gradient pulses are used. In this study, numerical simulations of such experiments were performed with parameters achievable with whole-body MR systems. It is shown that the theoretical model yields a good approximation of the signal behavior if an additional term describing free diffusion is included. More importantly, it is demonstrated that the shorter gradient pulses sufficient to achieve the desired diffusion weighting for multiple concatenations, increase the signal modulation considerably, e.g. by a factor of about five for five concatenations. Even at identical echo times, achieved by a shortened diffusion time, a moderate number of concatenations significantly improves the signal modulation. Thus, experiments on whole-body MR systems may benefit from multiple concatenations.

Detection of blob objects in microscopic zebrafish images based on gradient vector diffusion.

PubMed

Li, Gang; Liu, Tianming; Nie, Jingxin; Guo, Lei; Malicki, Jarema; Mara, Andrew; Holley, Scott A; Xia, Weiming; Wong, Stephen T C

2007-10-01

The zebrafish has become an important vertebrate animal model for the study of developmental biology, functional genomics, and disease mechanisms. It is also being used for drug discovery. Computerized detection of blob objects has been one of the important tasks in quantitative phenotyping of zebrafish. We present a new automated method that is able to detect blob objects, such as nuclei or cells in microscopic zebrafish images. This method is composed of three key steps. The first step is to produce a diffused gradient vector field by a physical elastic deformable model. In the second step, the flux image is computed on the diffused gradient vector field. The third step performs thresholding and nonmaximum suppression based on the flux image. We report the validation and experimental results of this method using zebrafish image datasets from three independent research labs. Both sensitivity and specificity of this method are over 90%. This method is able to differentiate closely juxtaposed or connected blob objects, with high sensitivity and specificity in different situations. It is characterized by a good, consistent performance in blob object detection.

Gradients estimation from random points with volumetric tensor in turbulence

NASA Astrophysics Data System (ADS)

Watanabe, Tomoaki; Nagata, Koji

2017-12-01

We present an estimation method of fully-resolved/coarse-grained gradients from randomly distributed points in turbulence. The method is based on a linear approximation of spatial gradients expressed with the volumetric tensor, which is a 3 × 3 matrix determined by a geometric distribution of the points. The coarse grained gradient can be considered as a low pass filtered gradient, whose cutoff is estimated with the eigenvalues of the volumetric tensor. The present method, the volumetric tensor approximation, is tested for velocity and passive scalar gradients in incompressible planar jet and mixing layer. Comparison with a finite difference approximation on a Cartesian grid shows that the volumetric tensor approximation computes the coarse grained gradients fairly well at a moderate computational cost under various conditions of spatial distributions of points. We also show that imposing the solenoidal condition improves the accuracy of the present method for solenoidal vectors, such as a velocity vector in incompressible flows, especially when the number of the points is not large. The volumetric tensor approximation with 4 points poorly estimates the gradient because of anisotropic distribution of the points. Increasing the number of points from 4 significantly improves the accuracy. Although the coarse grained gradient changes with the cutoff length, the volumetric tensor approximation yields the coarse grained gradient whose magnitude is close to the one obtained by the finite difference. We also show that the velocity gradient estimated with the present method well captures the turbulence characteristics such as local flow topology, amplification of enstrophy and strain, and energy transfer across scales.

IBMISPS (International Brain Mapping & Intraoperative Surgical Planning Symposium)

DTIC Science & Technology

2005-12-01

they received the 2005 Excellence in R, D & E award for their contribution in the feild of prosthetics and brain imaging. Excellence in Educational...specific bipolar magnetic gradient pulses which measure the velocity vector components of motion. Presented here are the development of dynamic MR...movies of quantitative velocity vector components, 30 frames per second. The 3 velocity vector maps with tensor analysis produced maps of the

Electro-gravity via geometric chrononfield

NASA Astrophysics Data System (ADS)

Suchard, Eytan H.

2017-05-01

In De Sitter / Anti De Sitter space-time and in other geometries, reference sub-manifolds from which proper time is measured along integral curves, are described as events. We introduce here a foliation with the help of a scalar field. The scalar field need not be unique but from the gradient of the scalar field, an intrinsic Reeb vector of the foliations perpendicular to the gradient vector is calculated. The Reeb vector describes the acceleration of a physical particle that moves along the integral curves that are formed by the gradient of the scalar field. The Reeb vector appears as a component of an anti-symmetric matrix which is a part of a rank-2, 2-Form. The 2-form is extended into a non-degenerate 4-form and into rank-4 matrix of a 2-form, which when multiplied by a velocity of a particle, becomes the acceleration of the particle. The matrix has one U(1) degree of freedom and an additional SU(2) degrees of freedom in two vectors that span the plane perpendicular to the gradient of the scalar field and to the Reeb vector. In total, there are U(1) x SU(2) degrees of freedom. SU(3) degrees of freedom arise from three dimensional foliations but require an additional symmetry to exist in order to have a valid covariant meaning. Matter in the Einstein Grossmann equation is replaced by the action of the acceleration field, i.e. by a geometric action which is not anticipated by the metric alone. This idea leads to a new formalism that replaces the conventional stress-energy-momentum-tensor. The formalism will be mainly developed for classical physics but will also be discussed for quantized physics based on events instead of particles. The result is that a positive charge manifests small attracting gravity and a stronger but small repelling acceleration field that repels even uncharged particles that have a rest mass. Negative charge manifests a repelling anti-gravity but also a stronger acceleration field that attracts even uncharged particles that have rest mass. Preliminary version: http://sciencedomain.org/abstract/9858

Anisotropic mimetic cosmology

NASA Astrophysics Data System (ADS)

Abbassi, M. H.; Jozani, A.; Sepangi, H. R.

2018-06-01

We consider a mimetic set up in which the mimetic scalar is coupled to a vector field. It is shown that such a field with a timelike component does not contribute to the background equations and yet produces healthy isocurvature perturbations with respect to ghost and gradient instabilities in spite of the absence of any propagating curvature perturbations at the level of the quadratic action. We then consider a vector field with spacelike components, which leads to an anisotropic Bianchi universe, and show that the ghost and gradient instabilities are absent in the limit of high momenta and that the propagating curvature perturbations have healthy UV behavior.

Domain Derivatives in Dielectric Rough Surface Scattering

DTIC Science & Technology

2015-01-01

and require the gradient of the objective function in the unknown model parameter vector at each stage of iteration. For large N, finite...differencing becomes numerically intensive, and an efficient alternative is domain differentiation in which the full gradient is obtained by solving a single...derivative calculation of the gradient for a locally perturbed dielectric interface. The method is non-variational, and algebraic in nature in that it

Minimizing inner product data dependencies in conjugate gradient iteration

NASA Technical Reports Server (NTRS)

Vanrosendale, J.

1983-01-01

The amount of concurrency available in conjugate gradient iteration is limited by the summations required in the inner product computations. The inner product of two vectors of length N requires time c log(N), if N or more processors are available. This paper describes an algebraic restructuring of the conjugate gradient algorithm which minimizes data dependencies due to inner product calculations. After an initial start up, the new algorithm can perform a conjugate gradient iteration in time c*log(log(N)).

Parallel processors and nonlinear structural dynamics algorithms and software

NASA Technical Reports Server (NTRS)

Belytschko, Ted

1990-01-01

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Symbolic computer vector analysis

NASA Technical Reports Server (NTRS)

Stoutemyer, D. R.

1977-01-01

A MACSYMA program is described which performs symbolic vector algebra and vector calculus. The program can combine and simplify symbolic expressions including dot products and cross products, together with the gradient, divergence, curl, and Laplacian operators. The distribution of these operators over sums or products is under user control, as are various other expansions, including expansion into components in any specific orthogonal coordinate system. There is also a capability for deriving the scalar or vector potential of a vector field. Examples include derivation of the partial differential equations describing fluid flow and magnetohydrodynamics, for 12 different classic orthogonal curvilinear coordinate systems.

SPECIAL ISSUE ON OPTICAL PROCESSING OF INFORMATION: Reconstruction of vector physical fields by optical tomography

NASA Astrophysics Data System (ADS)

Kulchin, Yurii N.; Vitrik, O. B.; Kamenev, O. T.; Kirichenko, O. V.; Petrov, Yu S.

1995-10-01

Reconstruction of vector physical fields by optical tomography, with the aid of a system of fibre-optic measuring lines, is considered. The reported experimental results are used to reconstruct the distribution of the square of the gradient of transverse displacements of a flat membrane.

Fractional vector calculus for fractional advection dispersion

NASA Astrophysics Data System (ADS)

Meerschaert, Mark M.; Mortensen, Jeff; Wheatcraft, Stephen W.

2006-07-01

We develop the basic tools of fractional vector calculus including a fractional derivative version of the gradient, divergence, and curl, and a fractional divergence theorem and Stokes theorem. These basic tools are then applied to provide a physical explanation for the fractional advection-dispersion equation for flow in heterogeneous porous media.

High angular resolution diffusion imaging with stimulated echoes: compensation and correction in experiment design and analysis.

PubMed

Lundell, Henrik; Alexander, Daniel C; Dyrby, Tim B

2014-08-01

Stimulated echo acquisition mode (STEAM) diffusion MRI can be advantageous over pulsed-gradient spin-echo (PGSE) for diffusion times that are long compared with T2 . It therefore has potential for biomedical diffusion imaging applications at 7T and above where T2 is short. However, gradient pulses other than the diffusion gradients in the STEAM sequence contribute much greater diffusion weighting than in PGSE and lead to a disrupted experimental design. Here, we introduce a simple compensation to the STEAM acquisition that avoids the orientational bias and disrupted experiment design that these gradient pulses can otherwise produce. The compensation is simple to implement by adjusting the gradient vectors in the diffusion pulses of the STEAM sequence, so that the net effective gradient vector including contributions from diffusion and other gradient pulses is as the experiment intends. High angular resolution diffusion imaging (HARDI) data were acquired with and without the proposed compensation. The data were processed to derive standard diffusion tensor imaging (DTI) maps, which highlight the need for the compensation. Ignoring the other gradient pulses, a bias in DTI parameters from STEAM acquisition is found, due both to confounds in the analysis and the experiment design. Retrospectively correcting the analysis with a calculation of the full B matrix can partly correct for these confounds, but an acquisition that is compensated as proposed is needed to remove the effect entirely. © 2014 The Authors. NMR in Biomedicine published by John Wiley & Sons, Ltd.

Manipulating the Gradient

ERIC Educational Resources Information Center

Gaze, Eric C.

2005-01-01

We introduce a cooperative learning, group lab for a Calculus III course to facilitate comprehension of the gradient vector and directional derivative concepts. The lab is a hands-on experience allowing students to manipulate a tangent plane and empirically measure the effect of partial derivatives on the direction of optimal ascent. (Contains 7…

Gradient-based controllers for timed continuous Petri nets

NASA Astrophysics Data System (ADS)

Lefebvre, Dimitri; Leclercq, Edouard; Druaux, Fabrice; Thomas, Philippe

2015-07-01

This paper is about control design for timed continuous Petri nets that are described as piecewise affine systems. In this context, the marking vector is considered as the state space vector, weighted marking of place subsets are defined as the model outputs and the model inputs correspond to multiplicative control actions that slow down the firing rate of some controllable transitions. Structural and functional sensitivity of the outputs with respect to the inputs are discussed in terms of Petri nets. Then, gradient-based controllers (GBC) are developed in order to adapt the control actions of the controllable transitions according to desired trajectories of the outputs.

Wideband radar cross section reduction using two-dimensional phase gradient metasurfaces

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Yongfeng; Qu, Shaobo; Wang, Jiafu

2014-06-02

Phase gradient metasurface (PGMs) are artificial surfaces that can provide pre-defined in-plane wave-vectors to manipulate the directions of refracted/reflected waves. In this Letter, we propose to achieve wideband radar cross section (RCS) reduction using two-dimensional (2D) PGMs. A 2D PGM was designed using a square combination of 49 split-ring sub-unit cells. The PGM can provide additional wave-vectors along the two in-plane directions simultaneously, leading to either surface wave conversion, deflected reflection, or diffuse reflection. Both the simulation and experiment results verified the wide-band, polarization-independent, high-efficiency RCS reduction induced by the 2D PGM.

A Note on the Application of the Extended Bernoulli Equation

DTIC Science & Technology

1999-02-01

as OV s ... - Vp „ _ = -±L L + VO , (2) Dt p where DIDt denotes the material derivative (discussed in following section); V is the vector...force potential; V is the vector gradient operator; s (J is the deviatoric-stress tensor arising from any type of elasto-viscoplastic constitutive...behavior; and s ^j is index notation for dsy/dxp denoting the following vector condensation of the deviatoric-stress tensor: ds ds ds

Manufacturing in space: Fluid dynamics numerical analysis

NASA Technical Reports Server (NTRS)

Robertson, S. J.; Nicholson, L. A.; Spradley, L. W.

1982-01-01

Numerical computations were performed for natural convection in circular enclosures under various conditions of acceleration. It was found that subcritical acceleration vectors applied in the direction of the temperature gradient will lead to an eventual state of rest regardless of the initial state of motion. Supercritical acceleration vectors will lead to the same steady state condition of motion regardless of the initial state of motion. Convection velocities were computed for acceleration vectors at various angles of the initial temperature gradient. The results for Rayleigh numbers of 1000 or less were found to closely follow Weinbaum's first order theory. Higher Rayleigh number results were shown to depart significantly from the first order theory. Supercritical behavior was confirmed for Rayleigh numbers greater than the known supercritical value of 9216. Response times were determined to provide an indication of the time required to change states of motion for the various cases considered.

ILUBCG2-11: Solution of 11-banded nonsymmetric linear equation systems by a preconditioned biconjugate gradient routine

NASA Astrophysics Data System (ADS)

Chen, Y.-M.; Koniges, A. E.; Anderson, D. V.

1989-10-01

The biconjugate gradient method (BCG) provides an attractive alternative to the usual conjugate gradient algorithms for the solution of sparse systems of linear equations with nonsymmetric and indefinite matrix operators. A preconditioned algorithm is given, whose form resembles the incomplete L-U conjugate gradient scheme (ILUCG2) previously presented. Although the BCG scheme requires the storage of two additional vectors, it converges in a significantly lesser number of iterations (often half), while the number of calculations per iteration remains essentially the same.

On dealing with multiple correlation peaks in PIV

NASA Astrophysics Data System (ADS)

Masullo, A.; Theunissen, R.

2018-05-01

A novel algorithm to analyse PIV images in the presence of strong in-plane displacement gradients and reduce sub-grid filtering is proposed in this paper. Interrogation windows subjected to strong in-plane displacement gradients often produce correlation maps presenting multiple peaks. Standard multi-grid procedures discard such ambiguous correlation windows using a signal to noise (SNR) filter. The proposed algorithm improves the standard multi-grid algorithm allowing the detection of splintered peaks in a correlation map through an automatic threshold, producing multiple displacement vectors for each correlation area. Vector locations are chosen by translating images according to the peak displacements and by selecting the areas with the strongest match. The method is assessed on synthetic images of a boundary layer of varying intensity and a sinusoidal displacement field of changing wavelength. An experimental case of a flow exhibiting strong velocity gradients is also provided to show the improvements brought by this technique.

GPU-accelerated adjoint algorithmic differentiation

NASA Astrophysics Data System (ADS)

Gremse, Felix; Höfter, Andreas; Razik, Lukas; Kiessling, Fabian; Naumann, Uwe

2016-03-01

Many scientific problems such as classifier training or medical image reconstruction can be expressed as minimization of differentiable real-valued cost functions and solved with iterative gradient-based methods. Adjoint algorithmic differentiation (AAD) enables automated computation of gradients of such cost functions implemented as computer programs. To backpropagate adjoint derivatives, excessive memory is potentially required to store the intermediate partial derivatives on a dedicated data structure, referred to as the ;tape;. Parallelization is difficult because threads need to synchronize their accesses during taping and backpropagation. This situation is aggravated for many-core architectures, such as Graphics Processing Units (GPUs), because of the large number of light-weight threads and the limited memory size in general as well as per thread. We show how these limitations can be mediated if the cost function is expressed using GPU-accelerated vector and matrix operations which are recognized as intrinsic functions by our AAD software. We compare this approach with naive and vectorized implementations for CPUs. We use four increasingly complex cost functions to evaluate the performance with respect to memory consumption and gradient computation times. Using vectorization, CPU and GPU memory consumption could be substantially reduced compared to the naive reference implementation, in some cases even by an order of complexity. The vectorization allowed usage of optimized parallel libraries during forward and reverse passes which resulted in high speedups for the vectorized CPU version compared to the naive reference implementation. The GPU version achieved an additional speedup of 7.5 ± 4.4, showing that the processing power of GPUs can be utilized for AAD using this concept. Furthermore, we show how this software can be systematically extended for more complex problems such as nonlinear absorption reconstruction for fluorescence-mediated tomography.

GPU-Accelerated Adjoint Algorithmic Differentiation.

PubMed

Gremse, Felix; Höfter, Andreas; Razik, Lukas; Kiessling, Fabian; Naumann, Uwe

2016-03-01

Many scientific problems such as classifier training or medical image reconstruction can be expressed as minimization of differentiable real-valued cost functions and solved with iterative gradient-based methods. Adjoint algorithmic differentiation (AAD) enables automated computation of gradients of such cost functions implemented as computer programs. To backpropagate adjoint derivatives, excessive memory is potentially required to store the intermediate partial derivatives on a dedicated data structure, referred to as the "tape". Parallelization is difficult because threads need to synchronize their accesses during taping and backpropagation. This situation is aggravated for many-core architectures, such as Graphics Processing Units (GPUs), because of the large number of light-weight threads and the limited memory size in general as well as per thread. We show how these limitations can be mediated if the cost function is expressed using GPU-accelerated vector and matrix operations which are recognized as intrinsic functions by our AAD software. We compare this approach with naive and vectorized implementations for CPUs. We use four increasingly complex cost functions to evaluate the performance with respect to memory consumption and gradient computation times. Using vectorization, CPU and GPU memory consumption could be substantially reduced compared to the naive reference implementation, in some cases even by an order of complexity. The vectorization allowed usage of optimized parallel libraries during forward and reverse passes which resulted in high speedups for the vectorized CPU version compared to the naive reference implementation. The GPU version achieved an additional speedup of 7.5 ± 4.4, showing that the processing power of GPUs can be utilized for AAD using this concept. Furthermore, we show how this software can be systematically extended for more complex problems such as nonlinear absorption reconstruction for fluorescence-mediated tomography.

GPU-Accelerated Adjoint Algorithmic Differentiation

PubMed Central

Gremse, Felix; Höfter, Andreas; Razik, Lukas; Kiessling, Fabian; Naumann, Uwe

2015-01-01

Many scientific problems such as classifier training or medical image reconstruction can be expressed as minimization of differentiable real-valued cost functions and solved with iterative gradient-based methods. Adjoint algorithmic differentiation (AAD) enables automated computation of gradients of such cost functions implemented as computer programs. To backpropagate adjoint derivatives, excessive memory is potentially required to store the intermediate partial derivatives on a dedicated data structure, referred to as the “tape”. Parallelization is difficult because threads need to synchronize their accesses during taping and backpropagation. This situation is aggravated for many-core architectures, such as Graphics Processing Units (GPUs), because of the large number of light-weight threads and the limited memory size in general as well as per thread. We show how these limitations can be mediated if the cost function is expressed using GPU-accelerated vector and matrix operations which are recognized as intrinsic functions by our AAD software. We compare this approach with naive and vectorized implementations for CPUs. We use four increasingly complex cost functions to evaluate the performance with respect to memory consumption and gradient computation times. Using vectorization, CPU and GPU memory consumption could be substantially reduced compared to the naive reference implementation, in some cases even by an order of complexity. The vectorization allowed usage of optimized parallel libraries during forward and reverse passes which resulted in high speedups for the vectorized CPU version compared to the naive reference implementation. The GPU version achieved an additional speedup of 7.5 ± 4.4, showing that the processing power of GPUs can be utilized for AAD using this concept. Furthermore, we show how this software can be systematically extended for more complex problems such as nonlinear absorption reconstruction for fluorescence-mediated tomography. PMID:26941443

Algorithms for parallel and vector computations

NASA Technical Reports Server (NTRS)

Ortega, James M.

1995-01-01

This is a final report on work performed under NASA grant NAG-1-1112-FOP during the period March, 1990 through February 1995. Four major topics are covered: (1) solution of nonlinear poisson-type equations; (2) parallel reduced system conjugate gradient method; (3) orderings for conjugate gradient preconditioners, and (4) SOR as a preconditioner.

Ghost circles in lattice Aubry-Mather theory

NASA Astrophysics Data System (ADS)

Mramor, Blaz; Rink, Bob

Monotone lattice recurrence relations such as the Frenkel-Kontorova lattice, arise in Hamiltonian lattice mechanics, as models for ferromagnetism and as discretization of elliptic PDEs. Mathematically, they are a multi-dimensional counterpart of monotone twist maps. Such recurrence relations often admit a variational structure, so that the solutions x:Z→R are the stationary points of a formal action function W(x). Given any rotation vector ω∈R, classical Aubry-Mather theory establishes the existence of a large collection of solutions of ∇W(x)=0 of rotation vector ω. For irrational ω, this is the well-known Aubry-Mather set. It consists of global minimizers and it may have gaps. In this paper, we study the parabolic gradient flow {dx}/{dt}=-∇W(x) and we will prove that every Aubry-Mather set can be interpolated by a continuous gradient-flow invariant family, the so-called 'ghost circle'. The existence of these ghost circles is known in dimension d=1, for rational rotation vectors and Morse action functions. The main technical result of this paper is therefore a compactness theorem for lattice ghost circles, based on a parabolic Harnack inequality for the gradient flow. This implies the existence of lattice ghost circles of arbitrary rotation vectors and for arbitrary actions. As a consequence, we can give a simple proof of the fact that when an Aubry-Mather set has a gap, then this gap must be filled with minimizers, or contain a non-minimizing solution.

Altered orientation and flight paths of pigeons reared on gravity anomalies: a GPS tracking study.

PubMed

Blaser, Nicole; Guskov, Sergei I; Meskenaite, Virginia; Kanevskyi, Valerii A; Lipp, Hans-Peter

2013-01-01

The mechanisms of pigeon homing are still not understood, in particular how they determine their position at unfamiliar locations. The "gravity vector" theory holds that pigeons memorize the gravity vector at their home loft and deduct home direction and distance from the angular difference between memorized and actual gravity vector. However, the gravity vector is tilted by different densities in the earth crust leading to gravity anomalies. We predicted that pigeons reared on different gravity anomalies would show different initial orientation and also show changes in their flight path when crossing a gravity anomaly. We reared one group of pigeons in a strong gravity anomaly with a north-to-south gravity gradient, and the other group of pigeons in a normal area but on a spot with a strong local anomaly with a west-to-east gravity gradient. After training over shorter distances, pigeons were released from a gravitationally and geomagnetically normal site 50 km north in the same direction for both home lofts. As expected by the theory, the two groups of pigeons showed divergent initial orientation. In addition, some of the GPS-tracked pigeons also showed changes in their flight paths when crossing gravity anomalies. We conclude that even small local gravity anomalies at the birth place of pigeons may have the potential to bias the map sense of pigeons, while reactivity to gravity gradients during flight was variable and appeared to depend on individual navigational strategies and frequency of position updates.

New methods for interpretation of magnetic vector and gradient tensor data II: application to the Mount Leyshon anomaly, Queensland, Australia

NASA Astrophysics Data System (ADS)

Clark, David A.

2013-04-01

Acquisition of magnetic gradient tensor data is anticipated to become routine in the near future. In the meantime, modern ultrahigh resolution conventional magnetic data can be used, with certain important caveats, to calculate magnetic vector components and gradient tensor elements from total magnetic intensity (TMI) or TMI gradient surveys. An accompanying paper presented new methods for inverting gradient tensor data to obtain source parameters for several elementary, but useful, models. These include point dipole (sphere), vertical line of dipoles (narrow vertical pipe), line of dipoles (horizontal cylinder), thin dipping sheet, and contact models. A key simplification is the use of eigenvalues and associated eigenvectors of the tensor. The normalised source strength (NSS), calculated from the eigenvalues, is a particularly useful rotational invariant that peaks directly over 3D compact sources, 2D compact sources, thin sheets, and contacts, independent of magnetisation direction. Source locations can be inverted directly from the NSS and its vector gradient. Some of these new methods have been applied to analysis of the magnetic signature of the Early Permian Mount Leyshon gold-mineralised system, Queensland. The Mount Leyshon magnetic anomaly is a prominent TMI low that is produced by rock units with strong reversed remanence acquired during the Late Palaeozoic Reverse Superchron. The inferred magnetic moment for the source zone of the Mount Leyshon magnetic anomaly is ~1010Am2. Its direction is consistent with petrophysical measurements. Given estimated magnetisation from samples and geological information, this suggests a volume of ~1.5km×1.5km×2km (vertical). The inferred depth of the centre of magnetisation is ~900m below surface, suggesting that the depth extent of the magnetic zone is ~1800m. Some of the deeper, undrilled portion of the magnetic zone could be a mafic intrusion similar to the nearby coeval Fenian Diorite, representing part of the parent magma chamber beneath the Mount Leyshon Intrusive Complex.

Radar observations of density gradients, electric fields, and plasma irregularities near polar cap patches in the context of the gradient-drift instability

NASA Astrophysics Data System (ADS)

Lamarche, Leslie J.; Makarevich, Roman A.

2017-03-01

We present observations of plasma density gradients, electric fields, and small-scale plasma irregularities near a polar cap patch made by the Super Dual Auroral Radar Network radar at Rankin Inlet (RKN) and the northern face of Resolute Bay Incoherent Scatter Radar (RISR-N). RKN echo power and occurrence are analyzed in the context of gradient-drift instability (GDI) theory, with a particular focus on the previously uninvestigated 2-D dependencies on wave propagation, electric field, and gradient vectors, with the latter two quantities evaluated directly from RISR-N measurements. It is shown that higher gradient and electric field components along the wave vector generally lead to the higher observed echo occurrence, which is consistent with the expected higher GDI growth rate, but the relationship with echo power is far less straightforward. The RKN echo power increases monotonically as the predicted linear growth rate approaches zero from negative values but does not continue this trend into positive growth rate values, in contrast with GDI predictions. The observed greater consistency of echo occurrence with GDI predictions suggests that GDI operating in the linear regime can control basic plasma structuring, but measured echo strength may be affected by other processes and factors, such as multistep or nonlinear processes or a shear-driven instability.

Guiding automated left ventricular chamber segmentation in cardiac imaging using the concept of conserved myocardial volume.

PubMed

Garson, Christopher D; Li, Bing; Acton, Scott T; Hossack, John A

2008-06-01

The active surface technique using gradient vector flow allows semi-automated segmentation of ventricular borders. The accuracy of the algorithm depends on the optimal selection of several key parameters. We investigated the use of conservation of myocardial volume for quantitative assessment of each of these parameters using synthetic and in vivo data. We predicted that for a given set of model parameters, strong conservation of volume would correlate with accurate segmentation. The metric was most useful when applied to the gradient vector field weighting and temporal step-size parameters, but less effective in guiding an optimal choice of the active surface tension and rigidity parameters.

A feasibility study of a 3-D finite element solution scheme for aeroengine duct acoustics

NASA Technical Reports Server (NTRS)

Abrahamson, A. L.

1980-01-01

The advantage from development of a 3-D model of aeroengine duct acoustics is the ability to analyze axial and circumferential liner segmentation simultaneously. The feasibility of a 3-D duct acoustics model was investigated using Galerkin or least squares element formulations combined with Gaussian elimination, successive over-relaxation, or conjugate gradient solution algorithms on conventional scalar computers and on a vector machine. A least squares element formulation combined with a conjugate gradient solver on a CDC Star vector computer initially appeared to have great promise, but severe difficulties were encountered with matrix ill-conditioning. These difficulties in conditioning rendered this technique impractical for realistic problems.

Segmentation of touching mycobacterium tuberculosis from Ziehl-Neelsen stained sputum smear images

NASA Astrophysics Data System (ADS)

Xu, Chao; Zhou, Dongxiang; Liu, Yunhui

2015-12-01

Touching Mycobacterium tuberculosis objects in the Ziehl-Neelsen stained sputum smear images present different shapes and invisible boundaries in the adhesion areas, which increases the difficulty in objects recognition and counting. In this paper, we present a segmentation method of combining the hierarchy tree analysis with gradient vector flow snake to address this problem. The skeletons of the objects are used for structure analysis based on the hierarchy tree. The gradient vector flow snake is used to estimate the object edge. Experimental results show that the single objects composing the touching objects are successfully segmented by the proposed method. This work will improve the accuracy and practicability of the computer-aided diagnosis of tuberculosis.

Mosquito communities with trap height and urban-rural gradient in Adelaide, South Australia: implications for disease vector surveillance.

PubMed

Johnston, Emily; Weinstein, Phillip; Slaney, David; Flies, Andrew S; Fricker, Stephen; Williams, Craig

2014-06-01

Understanding the factors influencing mosquito distribution is important for effective surveillance and control of nuisance and disease vector mosquitoes. The goal of this study was to determine how trap height and distance to the city center influenced the abundance and species of mosquitoes collected in Adelaide, South Australia. Mosquito communities were sampled at two heights (<2 m and ~10 m) along an urban-rural gradient. A total of 5,133 mosquitoes was identified over 176 trap nights. Aedes notoscriptus, Ae. vigilax, and Culex molestus were all more abundant in lower traps while Cx. quinquefasciatus (an ornithophilic species) was found to be more abundant in high traps. Distance to city center correlated strongly with the abundance of Ae. vigilax, Ae. camptorhynchus, Cx. globocoxitus, and Cx. molestus, all of which were most common at the sites farthest from the city and closest to the saltmarsh. Overall, the important disease vectors in South Australia (Ae. vigilax, Ae. camptorhynchus, Ae. notoscriptus, and Cx. annulirostris) were more abundant in low traps farthest from the city and closest to the saltmarsh. The current mosquito surveillance practice of setting traps within two meters of the ground is effective for sampling populations of the important disease vector species in South Australia. © 2014 The Society for Vector Ecology.

Vector method for strain estimation in phase-sensitive optical coherence elastography

NASA Astrophysics Data System (ADS)

Matveyev, A. L.; Matveev, L. A.; Sovetsky, A. A.; Gelikonov, G. V.; Moiseev, A. A.; Zaitsev, V. Y.

2018-06-01

A noise-tolerant approach to strain estimation in phase-sensitive optical coherence elastography, robust to decorrelation distortions, is discussed. The method is based on evaluation of interframe phase-variation gradient, but its main feature is that the phase is singled out at the very last step of the gradient estimation. All intermediate steps operate with complex-valued optical coherence tomography (OCT) signals represented as vectors in the complex plane (hence, we call this approach the ‘vector’ method). In comparison with such a popular method as least-square fitting of the phase-difference slope over a selected region (even in the improved variant with amplitude weighting for suppressing small-amplitude noisy pixels), the vector approach demonstrates superior tolerance to both additive noise in the receiving system and speckle-decorrelation caused by tissue straining. Another advantage of the vector approach is that it obviates the usual necessity of error-prone phase unwrapping. Here, special attention is paid to modifications of the vector method that make it especially suitable for processing deformations with significant lateral inhomogeneity, which often occur in real situations. The method’s advantages are demonstrated using both simulated and real OCT scans obtained during reshaping of a collagenous tissue sample irradiated by an IR laser beam producing complex spatially inhomogeneous deformations.

An M-step preconditioned conjugate gradient method for parallel computation

NASA Technical Reports Server (NTRS)

Adams, L.

1983-01-01

This paper describes a preconditioned conjugate gradient method that can be effectively implemented on both vector machines and parallel arrays to solve sparse symmetric and positive definite systems of linear equations. The implementation on the CYBER 203/205 and on the Finite Element Machine is discussed and results obtained using the method on these machines are given.

Toward an improved determination of Earth's lithospheric magnetic field from satellite observations

NASA Astrophysics Data System (ADS)

Kotsiaros, S.

2016-12-01

An analytical and numerical analysis of the spectral properties of the gradient tensor, initially performed by Rummel and van Gelderen (1992) for the gravity potential, shows that when the tensor elements are grouped into sets of semi-tangential and pure-tangential parts, they produce almost identical signal content as the normal element. Moreover, simple eigenvalue relations can be derived between these sets and the spherical harmonic expansion of the potential. This theoretical development generally applies to any potential field. First, the analysis of Rummel and van Gelderen (1992) is adapted to the magnetic field case and then the elements of the magnetic gradient tensor are estimated by 2 years of Swarm data and grouped into Γ(1) = {[∇B]rθ,[∇B]rφ} resp. Γ(2) = {[∇B]θθ-[∇B]φφ, 2[∇B]θφ}. It is shown that the estimated combinations Γ(1) and Γ(2) produce similar signal content as the theoretical radial gradient [∇B]rr. These results demonstrate the ability of multi-satellite missions such as Swarm, which cannot directly measure the radial gradient, to retrieve similar signal content by means of the horizontal gradients. Finally, lithospheric field models are derived using the gradient combinations Γ(1) and Γ(2) and compared with models derived from traditional vector and gradient data. The model resulting from Γ(1) leads to a very similar, and in particular cases improved, model compared to models retrieved by using approximately three times more data, i.e. a full set of vector, North-South and East-West gradients. ReferencesRummel, R., and M. van Gelderen (1992), Spectral analysis of the full gravity tensor, Geophysical Journal International, 111 (1), 159-169.

Retinal Microaneurysms Detection Using Gradient Vector Analysis and Class Imbalance Classification.

PubMed

Dai, Baisheng; Wu, Xiangqian; Bu, Wei

2016-01-01

Retinal microaneurysms (MAs) are the earliest clinically observable lesions of diabetic retinopathy. Reliable automated MAs detection is thus critical for early diagnosis of diabetic retinopathy. This paper proposes a novel method for the automated MAs detection in color fundus images based on gradient vector analysis and class imbalance classification, which is composed of two stages, i.e. candidate MAs extraction and classification. In the first stage, a candidate MAs extraction algorithm is devised by analyzing the gradient field of the image, in which a multi-scale log condition number map is computed based on the gradient vectors for vessel removal, and then the candidate MAs are localized according to the second order directional derivatives computed in different directions. Due to the complexity of fundus image, besides a small number of true MAs, there are also a large amount of non-MAs in the extracted candidates. Classifying the true MAs and the non-MAs is an extremely class imbalanced classification problem. Therefore, in the second stage, several types of features including geometry, contrast, intensity, edge, texture, region descriptors and other features are extracted from the candidate MAs and a class imbalance classifier, i.e., RUSBoost, is trained for the MAs classification. With the Retinopathy Online Challenge (ROC) criterion, the proposed method achieves an average sensitivity of 0.433 at 1/8, 1/4, 1/2, 1, 2, 4 and 8 false positives per image on the ROC database, which is comparable with the state-of-the-art approaches, and 0.321 on the DiaRetDB1 V2.1 database, which outperforms the state-of-the-art approaches.

Efficient boundary hunting via vector quantization

NASA Astrophysics Data System (ADS)

Diamantini, Claudia; Panti, Maurizio

2001-03-01

A great amount of information about a classification problem is contained in those instances falling near the decision boundary. This intuition dates back to the earliest studies in pattern recognition, and in the more recent adaptive approaches to the so called boundary hunting, such as the work of Aha et alii on Instance Based Learning and the work of Vapnik et alii on Support Vector Machines. The last work is of particular interest, since theoretical and experimental results ensure the accuracy of boundary reconstruction. However, its optimization approach has heavy computational and memory requirements, which limits its application on huge amounts of data. In the paper we describe an alternative approach to boundary hunting based on adaptive labeled quantization architectures. The adaptation is performed by a stochastic gradient algorithm for the minimization of the error probability. Error probability minimization guarantees the accurate approximation of the optimal decision boundary, while the use of a stochastic gradient algorithm defines an efficient method to reach such approximation. In the paper comparisons to Support Vector Machines are considered.

Frequency-domain beamformers using conjugate gradient techniques for speech enhancement.

PubMed

Zhao, Shengkui; Jones, Douglas L; Khoo, Suiyang; Man, Zhihong

2014-09-01

A multiple-iteration constrained conjugate gradient (MICCG) algorithm and a single-iteration constrained conjugate gradient (SICCG) algorithm are proposed to realize the widely used frequency-domain minimum-variance-distortionless-response (MVDR) beamformers and the resulting algorithms are applied to speech enhancement. The algorithms are derived based on the Lagrange method and the conjugate gradient techniques. The implementations of the algorithms avoid any form of explicit or implicit autocorrelation matrix inversion. Theoretical analysis establishes formal convergence of the algorithms. Specifically, the MICCG algorithm is developed based on a block adaptation approach and it generates a finite sequence of estimates that converge to the MVDR solution. For limited data records, the estimates of the MICCG algorithm are better than the conventional estimators and equivalent to the auxiliary vector algorithms. The SICCG algorithm is developed based on a continuous adaptation approach with a sample-by-sample updating procedure and the estimates asymptotically converge to the MVDR solution. An illustrative example using synthetic data from a uniform linear array is studied and an evaluation on real data recorded by an acoustic vector sensor array is demonstrated. Performance of the MICCG algorithm and the SICCG algorithm are compared with the state-of-the-art approaches.

Aeromagnetic gradient compensation method for helicopter based on ɛ-support vector regression algorithm

NASA Astrophysics Data System (ADS)

Wu, Peilin; Zhang, Qunying; Fei, Chunjiao; Fang, Guangyou

2017-04-01

Aeromagnetic gradients are typically measured by optically pumped magnetometers mounted on an aircraft. Any aircraft, particularly helicopters, produces significant levels of magnetic interference. Therefore, aeromagnetic compensation is essential, and least square (LS) is the conventional method used for reducing interference levels. However, the LSs approach to solving the aeromagnetic interference model has a few difficulties, one of which is in handling multicollinearity. Therefore, we propose an aeromagnetic gradient compensation method, specifically targeted for helicopter use but applicable on any airborne platform, which is based on the ɛ-support vector regression algorithm. The structural risk minimization criterion intrinsic to the method avoids multicollinearity altogether. Local aeromagnetic anomalies can be retained, and platform-generated fields are suppressed simultaneously by constructing an appropriate loss function and kernel function. The method was tested using an unmanned helicopter and obtained improvement ratios of 12.7 and 3.5 in the vertical and horizontal gradient data, respectively. Both of these values are probably better than those that would have been obtained from the conventional method applied to the same data, had it been possible to do so in a suitable comparative context. The validity of the proposed method is demonstrated by the experimental result.

Vector competence of northern and southern European Culex pipiens pipiens mosquitoes for West Nile virus across a gradient of temperatures.

PubMed

Vogels, C B F; Göertz, G P; Pijlman, G P; Koenraadt, C J M

2017-12-01

In Europe, West Nile virus (WNV) outbreaks have been limited to southern and central European countries. However, competent mosquito vectors and susceptible bird hosts are present in northern Europe. Differences in temperature and vector competence of mosquito populations may explain the absence of WNV outbreaks in northern Europe. The aim of the present study was to directly compare vector competence of northern and southern European Culex pipiens (Cx. p.) pipiens mosquitoes for WNV across a gradient of temperatures. WNV infection and transmission rates were determined for two Cx. p. pipiens populations originating from The Netherlands and Italy, respectively. Mosquitoes were orally exposed by providing an infectious bloodmeal, or by injecting WNV (lineage 2) in the thorax, followed by 14-day incubation at 18, 23, or 28 °C. No differences in infection or transmission rates were found between the Cx. p. pipiens populations with both infection methods, but WNV transmission rates were significantly higher at temperatures above 18 °C. The absence of WNV outbreaks in northern Europe cannot be explained by differences in vector competence between Cx. p. pipiens populations originating from northern and southern Europe. This study suggests that low temperature is a key limiting factor for WNV transmission. © 2017 The Authors. Medical and Veterinary Entomology published by John Wiley & Sons Ltd on behalf of Royal Entomological Society.

Stochastic estimates of gradient from laser measurements for an autonomous Martian roving vehicle

NASA Technical Reports Server (NTRS)

Burger, P. A.

1973-01-01

The general problem of estimating the state vector x from the state equation h = Ax where h, A, and x are all stochastic, is presented. Specifically, the problem is for an autonomous Martian roving vehicle to utilize laser measurements in estimating the gradient of the terrain. Error exists due to two factors - surface roughness and instrumental measurements. The errors in slope depend on the standard deviations of these noise factors. Numerically, the error in gradient is expressed as a function of instrumental inaccuracies. Certain guidelines for the accuracy of permissable gradient must be set. It is found that present technology can meet these guidelines.

Vector calculus in non-integer dimensional space and its applications to fractal media

NASA Astrophysics Data System (ADS)

Tarasov, Vasily E.

2015-02-01

We suggest a generalization of vector calculus for the case of non-integer dimensional space. The first and second orders operations such as gradient, divergence, the scalar and vector Laplace operators for non-integer dimensional space are defined. For simplification we consider scalar and vector fields that are independent of angles. We formulate a generalization of vector calculus for rotationally covariant scalar and vector functions. This generalization allows us to describe fractal media and materials in the framework of continuum models with non-integer dimensional space. As examples of application of the suggested calculus, we consider elasticity of fractal materials (fractal hollow ball and fractal cylindrical pipe with pressure inside and outside), steady distribution of heat in fractal media, electric field of fractal charged cylinder. We solve the correspondent equations for non-integer dimensional space models.

Gradient-based adaptation of general gaussian kernels.

PubMed

Glasmachers, Tobias; Igel, Christian

2005-10-01

Gradient-based optimizing of gaussian kernel functions is considered. The gradient for the adaptation of scaling and rotation of the input space is computed to achieve invariance against linear transformations. This is done by using the exponential map as a parameterization of the kernel parameter manifold. By restricting the optimization to a constant trace subspace, the kernel size can be controlled. This is, for example, useful to prevent overfitting when minimizing radius-margin generalization performance measures. The concepts are demonstrated by training hard margin support vector machines on toy data.

Multi-color incomplete Cholesky conjugate gradient methods for vector computers. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Poole, E. L.

1986-01-01

In this research, we are concerned with the solution on vector computers of linear systems of equations, Ax = b, where A is a larger, sparse symmetric positive definite matrix. We solve the system using an iterative method, the incomplete Cholesky conjugate gradient method (ICCG). We apply a multi-color strategy to obtain p-color matrices for which a block-oriented ICCG method is implemented on the CYBER 205. (A p-colored matrix is a matrix which can be partitioned into a pXp block matrix where the diagonal blocks are diagonal matrices). This algorithm, which is based on a no-fill strategy, achieves O(N/p) length vector operations in both the decomposition of A and in the forward and back solves necessary at each iteration of the method. We discuss the natural ordering of the unknowns as an ordering that minimizes the number of diagonals in the matrix and define multi-color orderings in terms of disjoint sets of the unknowns. We give necessary and sufficient conditions to determine which multi-color orderings of the unknowns correpond to p-color matrices. A performance model is given which is used both to predict execution time for ICCG methods and also to compare an ICCG method to conjugate gradient without preconditioning or another ICCG method. Results are given from runs on the CYBER 205 at NASA's Langley Research Center for four model problems.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Guo, Yiting; Dong, Bin; Wang, Bing

Purpose: Effective and accurate segmentation of the aortic valve (AV) from sequenced ultrasound (US) images remains a technical challenge because of intrinsic factors of ultrasound images that impact the quality and the continuous changes of shape and position of segmented objects. In this paper, a novel shape-constraint gradient Chan-Vese (GCV) model is proposed for segmenting the AV from time serial echocardiography. Methods: The GCV model is derived by incorporating the energy of the gradient vector flow into a CV model framework, where the gradient vector energy term is introduced by calculating the deviation angle between the inward normal force ofmore » the evolution contour and the gradient vector force. The flow force enlarges the capture range and enhances the blurred boundaries of objects. This is achieved by adding a circle-like contour (constructed using the AV structure region as a constraint shape) as an energy item to the GCV model through the shape comparison function. This shape-constrained energy can enhance the image constraint force by effectively connecting separate gaps of the object edge as well as driving the evolution contour to quickly approach the ideal object. Because of the slight movement of the AV in adjacent frames, the initial constraint shape is defined by users, with the other constraint shapes being derived from the segmentation results of adjacent sequence frames after morphological filtering. The AV is segmented from the US images by minimizing the proposed energy function. Results: To evaluate the performance of the proposed method, five assessment parameters were used to compare it with manual delineations performed by radiologists (gold standards). Three hundred and fifteen images acquired from nine groups were analyzed in the experiment. The area-metric overlap error rate was 6.89% ± 2.88%, the relative area difference rate 3.94% ± 2.63%, the average symmetric contour distance 1.08 ± 0.43 mm, the root mean square symmetric contour distance 1.37 ± 0.52 mm, and the maximum symmetric contour distance was 3.57 ± 1.72 mm. Conclusions: Compared with the CV model, as a result of the combination of the gradient vector and neighborhood shape information, this semiautomatic segmentation method significantly improves the accuracy and robustness of AV segmentation, making it feasible for improved segmentation of aortic valves from US images that have fuzzy boundaries.« less

Truncated Conjugate Gradient: An Optimal Strategy for the Analytical Evaluation of the Many-Body Polarization Energy and Forces in Molecular Simulations.

PubMed

Aviat, Félix; Levitt, Antoine; Stamm, Benjamin; Maday, Yvon; Ren, Pengyu; Ponder, Jay W; Lagardère, Louis; Piquemal, Jean-Philip

2017-01-10

We introduce a new class of methods, denoted as Truncated Conjugate Gradient(TCG), to solve the many-body polarization energy and its associated forces in molecular simulations (i.e. molecular dynamics (MD) and Monte Carlo). The method consists in a fixed number of Conjugate Gradient (CG) iterations. TCG approaches provide a scalable solution to the polarization problem at a user-chosen cost and a corresponding optimal accuracy. The optimality of the CG-method guarantees that the number of the required matrix-vector products are reduced to a minimum compared to other iterative methods. This family of methods is non-empirical, fully adaptive, and provides analytical gradients, avoiding therefore any energy drift in MD as compared to popular iterative solvers. Besides speed, one great advantage of this class of approximate methods is that their accuracy is systematically improvable. Indeed, as the CG-method is a Krylov subspace method, the associated error is monotonically reduced at each iteration. On top of that, two improvements can be proposed at virtually no cost: (i) the use of preconditioners can be employed, which leads to the Truncated Preconditioned Conjugate Gradient (TPCG); (ii) since the residual of the final step of the CG-method is available, one additional Picard fixed point iteration ("peek"), equivalent to one step of Jacobi Over Relaxation (JOR) with relaxation parameter ω, can be made at almost no cost. This method is denoted by TCG-n(ω). Black-box adaptive methods to find good choices of ω are provided and discussed. Results show that TPCG-3(ω) is converged to high accuracy (a few kcal/mol) for various types of systems including proteins and highly charged systems at the fixed cost of four matrix-vector products: three CG iterations plus the initial CG descent direction. Alternatively, T(P)CG-2(ω) provides robust results at a reduced cost (three matrix-vector products) and offers new perspectives for long polarizable MD as a production algorithm. The T(P)CG-1(ω) level provides less accurate solutions for inhomogeneous systems, but its applicability to well-conditioned problems such as water is remarkable, with only two matrix-vector product evaluations.

Truncated Conjugate Gradient: An Optimal Strategy for the Analytical Evaluation of the Many-Body Polarization Energy and Forces in Molecular Simulations

PubMed Central

2016-01-01

We introduce a new class of methods, denoted as Truncated Conjugate Gradient(TCG), to solve the many-body polarization energy and its associated forces in molecular simulations (i.e. molecular dynamics (MD) and Monte Carlo). The method consists in a fixed number of Conjugate Gradient (CG) iterations. TCG approaches provide a scalable solution to the polarization problem at a user-chosen cost and a corresponding optimal accuracy. The optimality of the CG-method guarantees that the number of the required matrix-vector products are reduced to a minimum compared to other iterative methods. This family of methods is non-empirical, fully adaptive, and provides analytical gradients, avoiding therefore any energy drift in MD as compared to popular iterative solvers. Besides speed, one great advantage of this class of approximate methods is that their accuracy is systematically improvable. Indeed, as the CG-method is a Krylov subspace method, the associated error is monotonically reduced at each iteration. On top of that, two improvements can be proposed at virtually no cost: (i) the use of preconditioners can be employed, which leads to the Truncated Preconditioned Conjugate Gradient (TPCG); (ii) since the residual of the final step of the CG-method is available, one additional Picard fixed point iteration (“peek”), equivalent to one step of Jacobi Over Relaxation (JOR) with relaxation parameter ω, can be made at almost no cost. This method is denoted by TCG-n(ω). Black-box adaptive methods to find good choices of ω are provided and discussed. Results show that TPCG-3(ω) is converged to high accuracy (a few kcal/mol) for various types of systems including proteins and highly charged systems at the fixed cost of four matrix-vector products: three CG iterations plus the initial CG descent direction. Alternatively, T(P)CG-2(ω) provides robust results at a reduced cost (three matrix-vector products) and offers new perspectives for long polarizable MD as a production algorithm. The T(P)CG-1(ω) level provides less accurate solutions for inhomogeneous systems, but its applicability to well-conditioned problems such as water is remarkable, with only two matrix-vector product evaluations. PMID:28068773

Gradient Evolution-based Support Vector Machine Algorithm for Classification

NASA Astrophysics Data System (ADS)

Zulvia, Ferani E.; Kuo, R. J.

2018-03-01

This paper proposes a classification algorithm based on a support vector machine (SVM) and gradient evolution (GE) algorithms. SVM algorithm has been widely used in classification. However, its result is significantly influenced by the parameters. Therefore, this paper aims to propose an improvement of SVM algorithm which can find the best SVMs’ parameters automatically. The proposed algorithm employs a GE algorithm to automatically determine the SVMs’ parameters. The GE algorithm takes a role as a global optimizer in finding the best parameter which will be used by SVM algorithm. The proposed GE-SVM algorithm is verified using some benchmark datasets and compared with other metaheuristic-based SVM algorithms. The experimental results show that the proposed GE-SVM algorithm obtains better results than other algorithms tested in this paper.

Detection of License Plate using Sliding Window, Histogram of Oriented Gradient, and Support Vector Machines Method

NASA Astrophysics Data System (ADS)

Astawa, INGA; Gusti Ngurah Bagus Caturbawa, I.; Made Sajayasa, I.; Dwi Suta Atmaja, I. Made Ari

2018-01-01

The license plate recognition usually used as part of system such as parking system. License plate detection considered as the most important step in the license plate recognition system. We propose methods that can be used to detect the vehicle plate on mobile phone. In this paper, we used Sliding Window, Histogram of Oriented Gradient (HOG), and Support Vector Machines (SVM) method to license plate detection so it will increase the detection level even though the image is not in a good quality. The image proceed by Sliding Window method in order to find plate position. Feature extraction in every window movement had been done by HOG and SVM method. Good result had shown in this research, which is 96% of accuracy.

Electrokinetic Control of Viscous Fingering

NASA Astrophysics Data System (ADS)

Mirzadeh, Mohammad; Bazant, Martin Z.

2017-10-01

We present a theory of the interfacial stability of two immiscible electrolytes under the coupled action of pressure gradients and electric fields in a Hele-Shaw cell or porous medium. Mathematically, our theory describes a phenomenon of "vector Laplacian growth," in which the interface moves in response to the gradient of a vector-valued potential function through a generalized mobility tensor. Physically, we extend the classical Saffman-Taylor problem to electrolytes by incorporating electrokinetic (EK) phenomena. A surprising prediction is that viscous fingering can be controlled by varying the injection ratio of electric current to flow rate. Beyond a critical injection ratio, stability depends only upon the relative direction of flow and current, regardless of the viscosity ratio. Possible applications include porous materials processing, electrically enhanced oil recovery, and EK remediation of contaminated soils.

Stable solutions of inflation driven by vector fields

NASA Astrophysics Data System (ADS)

Emami, Razieh; Mukohyama, Shinji; Namba, Ryo; Zhang, Ying-li

2017-03-01

Many models of inflation driven by vector fields alone have been known to be plagued by pathological behaviors, namely ghost and/or gradient instabilities. In this work, we seek a new class of vector-driven inflationary models that evade all of the mentioned instabilities. We build our analysis on the Generalized Proca Theory with an extension to three vector fields to realize isotropic expansion. We obtain the conditions required for quasi de-Sitter solutions to be an attractor analogous to the standard slow-roll one and those for their stability at the level of linearized perturbations. Identifying the remedy to the existing unstable models, we provide a simple example and explicitly show its stability. This significantly broadens our knowledge on vector inflationary scenarios, reviving potential phenomenological interests for this class of models.

Fast higher-order MR image reconstruction using singular-vector separation.

PubMed

Wilm, Bertram J; Barmet, Christoph; Pruessmann, Klaas P

2012-07-01

Medical resonance imaging (MRI) conventionally relies on spatially linear gradient fields for image encoding. However, in practice various sources of nonlinear fields can perturb the encoding process and give rise to artifacts unless they are suitably addressed at the reconstruction level. Accounting for field perturbations that are neither linear in space nor constant over time, i.e., dynamic higher-order fields, is particularly challenging. It was previously shown to be feasible with conjugate-gradient iteration. However, so far this approach has been relatively slow due to the need to carry out explicit matrix-vector multiplications in each cycle. In this work, it is proposed to accelerate higher-order reconstruction by expanding the encoding matrix such that fast Fourier transform can be employed for more efficient matrix-vector computation. The underlying principle is to represent the perturbing terms as sums of separable functions of space and time. Compact representations with this property are found by singular-vector analysis of the perturbing matrix. Guidelines for balancing the accuracy and speed of the resulting algorithm are derived by error propagation analysis. The proposed technique is demonstrated for the case of higher-order field perturbations due to eddy currents caused by diffusion weighting. In this example, image reconstruction was accelerated by two orders of magnitude.

Stochastic estimates of gradient from laser measurements for an autonomous Martian Roving Vehicle

NASA Technical Reports Server (NTRS)

Shen, C. N.; Burger, P.

1973-01-01

The general problem presented in this paper is one of estimating the state vector x from the state equation h = Ax, where h, A, and x are all stochastic. Specifically, the problem is for an autonomous Martian Roving Vehicle to utilize laser measurements in estimating the gradient of the terrain. Error exists due to two factors - surface roughness and instrumental measurements. The errors in slope depend on the standard deviations of these noise factors. Numerically, the error in gradient is expressed as a function of instrumental inaccuracies. Certain guidelines for the accuracy of permissable gradient must be set. It is found that present technology can meet these guidelines.-

Thermotropism by primary roots of maize

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fortin, M.-C.; Poff, K.L.

1990-05-01

Sensing in the roots of higher plants has long been recognized to be restricted mainly to gravitropism and thigmotropism. However, root responses to temperature gradients have not been extensively studied. We have designed experiments under controlled conditions to test if and how root direction of maize can be altered by thermal gradients perpendicular to the gravity vector. Primary roots of maize grown on agar plates exhibit positive thermotropism (curvature toward the warmer temperature) when exposed to gradients of 0.5 to 4.2{degree}C cm{sup {minus}1}. The extent of thermotropism depends on the temperature gradient and the temperature at which the root ismore » placed within the gradient. The curvature cannot be accounted for by differential growth as a direct effect of temperature on each side of the root.« less

A kinetic flux vector splitting scheme for shallow water equations incorporating variable bottom topography and horizontal temperature gradients.

PubMed

Saleem, M Rehan; Ashraf, Waqas; Zia, Saqib; Ali, Ishtiaq; Qamar, Shamsul

2018-01-01

This paper is concerned with the derivation of a well-balanced kinetic scheme to approximate a shallow flow model incorporating non-flat bottom topography and horizontal temperature gradients. The considered model equations, also called as Ripa system, are the non-homogeneous shallow water equations considering temperature gradients and non-uniform bottom topography. Due to the presence of temperature gradient terms, the steady state at rest is of primary interest from the physical point of view. However, capturing of this steady state is a challenging task for the applied numerical methods. The proposed well-balanced kinetic flux vector splitting (KFVS) scheme is non-oscillatory and second order accurate. The second order accuracy of the scheme is obtained by considering a MUSCL-type initial reconstruction and Runge-Kutta time stepping method. The scheme is applied to solve the model equations in one and two space dimensions. Several numerical case studies are carried out to validate the proposed numerical algorithm. The numerical results obtained are compared with those of staggered central NT scheme. The results obtained are also in good agreement with the recently published results in the literature, verifying the potential, efficiency, accuracy and robustness of the suggested numerical scheme.

A kinetic flux vector splitting scheme for shallow water equations incorporating variable bottom topography and horizontal temperature gradients

PubMed Central

2018-01-01

This paper is concerned with the derivation of a well-balanced kinetic scheme to approximate a shallow flow model incorporating non-flat bottom topography and horizontal temperature gradients. The considered model equations, also called as Ripa system, are the non-homogeneous shallow water equations considering temperature gradients and non-uniform bottom topography. Due to the presence of temperature gradient terms, the steady state at rest is of primary interest from the physical point of view. However, capturing of this steady state is a challenging task for the applied numerical methods. The proposed well-balanced kinetic flux vector splitting (KFVS) scheme is non-oscillatory and second order accurate. The second order accuracy of the scheme is obtained by considering a MUSCL-type initial reconstruction and Runge-Kutta time stepping method. The scheme is applied to solve the model equations in one and two space dimensions. Several numerical case studies are carried out to validate the proposed numerical algorithm. The numerical results obtained are compared with those of staggered central NT scheme. The results obtained are also in good agreement with the recently published results in the literature, verifying the potential, efficiency, accuracy and robustness of the suggested numerical scheme. PMID:29851978

Skeletonization and Partitioning of Digital Images Using Discrete Morse Theory.

PubMed

Delgado-Friedrichs, Olaf; Robins, Vanessa; Sheppard, Adrian

2015-03-01

We show how discrete Morse theory provides a rigorous and unifying foundation for defining skeletons and partitions of grayscale digital images. We model a grayscale image as a cubical complex with a real-valued function defined on its vertices (the voxel values). This function is extended to a discrete gradient vector field using the algorithm presented in Robins, Wood, Sheppard TPAMI 33:1646 (2011). In the current paper we define basins (the building blocks of a partition) and segments of the skeleton using the stable and unstable sets associated with critical cells. The natural connection between Morse theory and homology allows us to prove the topological validity of these constructions; for example, that the skeleton is homotopic to the initial object. We simplify the basins and skeletons via Morse-theoretic cancellation of critical cells in the discrete gradient vector field using a strategy informed by persistent homology. Simple working Python code for our algorithms for efficient vector field traversal is included. Example data are taken from micro-CT images of porous materials, an application area where accurate topological models of pore connectivity are vital for fluid-flow modelling.

Vector disparity sensor with vergence control for active vision systems.

PubMed

Barranco, Francisco; Diaz, Javier; Gibaldi, Agostino; Sabatini, Silvio P; Ros, Eduardo

2012-01-01

This paper presents an architecture for computing vector disparity for active vision systems as used on robotics applications. The control of the vergence angle of a binocular system allows us to efficiently explore dynamic environments, but requires a generalization of the disparity computation with respect to a static camera setup, where the disparity is strictly 1-D after the image rectification. The interaction between vision and motor control allows us to develop an active sensor that achieves high accuracy of the disparity computation around the fixation point, and fast reaction time for the vergence control. In this contribution, we address the development of a real-time architecture for vector disparity computation using an FPGA device. We implement the disparity unit and the control module for vergence, version, and tilt to determine the fixation point. In addition, two on-chip different alternatives for the vector disparity engines are discussed based on the luminance (gradient-based) and phase information of the binocular images. The multiscale versions of these engines are able to estimate the vector disparity up to 32 fps on VGA resolution images with very good accuracy as shown using benchmark sequences with known ground-truth. The performances in terms of frame-rate, resource utilization, and accuracy of the presented approaches are discussed. On the basis of these results, our study indicates that the gradient-based approach leads to the best trade-off choice for the integration with the active vision system.

Artificial Vector Calibration Method for Differencing Magnetic Gradient Tensor Systems

PubMed Central

Li, Zhining; Zhang, Yingtang; Yin, Gang

2018-01-01

The measurement error of the differencing (i.e., using two homogenous field sensors at a known baseline distance) magnetic gradient tensor system includes the biases, scale factors, nonorthogonality of the single magnetic sensor, and the misalignment error between the sensor arrays, all of which can severely affect the measurement accuracy. In this paper, we propose a low-cost artificial vector calibration method for the tensor system. Firstly, the error parameter linear equations are constructed based on the single-sensor’s system error model to obtain the artificial ideal vector output of the platform, with the total magnetic intensity (TMI) scalar as a reference by two nonlinear conversions, without any mathematical simplification. Secondly, the Levenberg–Marquardt algorithm is used to compute the integrated model of the 12 error parameters by nonlinear least-squares fitting method with the artificial vector output as a reference, and a total of 48 parameters of the system is estimated simultaneously. The calibrated system outputs along the reference platform-orthogonal coordinate system. The analysis results show that the artificial vector calibrated output can track the orientation fluctuations of TMI accurately, effectively avoiding the “overcalibration” problem. The accuracy of the error parameters’ estimation in the simulation is close to 100%. The experimental root-mean-square error (RMSE) of the TMI and tensor components is less than 3 nT and 20 nT/m, respectively, and the estimation of the parameters is highly robust. PMID:29373544

Vector Disparity Sensor with Vergence Control for Active Vision Systems

PubMed Central

Barranco, Francisco; Diaz, Javier; Gibaldi, Agostino; Sabatini, Silvio P.; Ros, Eduardo

2012-01-01

This paper presents an architecture for computing vector disparity for active vision systems as used on robotics applications. The control of the vergence angle of a binocular system allows us to efficiently explore dynamic environments, but requires a generalization of the disparity computation with respect to a static camera setup, where the disparity is strictly 1-D after the image rectification. The interaction between vision and motor control allows us to develop an active sensor that achieves high accuracy of the disparity computation around the fixation point, and fast reaction time for the vergence control. In this contribution, we address the development of a real-time architecture for vector disparity computation using an FPGA device. We implement the disparity unit and the control module for vergence, version, and tilt to determine the fixation point. In addition, two on-chip different alternatives for the vector disparity engines are discussed based on the luminance (gradient-based) and phase information of the binocular images. The multiscale versions of these engines are able to estimate the vector disparity up to 32 fps on VGA resolution images with very good accuracy as shown using benchmark sequences with known ground-truth. The performances in terms of frame-rate, resource utilization, and accuracy of the presented approaches are discussed. On the basis of these results, our study indicates that the gradient-based approach leads to the best trade-off choice for the integration with the active vision system. PMID:22438737

Toward more complete magnetic gradiometry with the Swarm mission

NASA Astrophysics Data System (ADS)

Kotsiaros, Stavros

2016-07-01

An analytical and numerical analysis of the spectral properties of the gradient tensor, initially performed by Rummel and van Gelderen (Geophys J Int 111(1):159-169, 1992) for the gravity potential, shows that when the tensor elements are grouped into sets of semi-tangential and pure-tangential parts, they produce almost identical signal content as the normal element. Moreover, simple eigenvalue relations can be derived between these sets and the spherical harmonic expansion of the potential. This theoretical development generally applies to any potential field. First, the analysis of Rummel and van Gelderen (1992) is adapted to the magnetic field case and then the elements of the magnetic gradient tensor are estimated by 2 years of Swarm data and grouped into \\varvec{Γ }^{(1)} = {[\\varvec{nabla } {{B}}]_{rθ },[\\varvec{nabla } {{B}}]_{r\\varphi }} resp. \\varvec{Γ }^{(2)} = {[\\varvec{nabla } {{B}}]_{θ θ }-[\\varvec{nabla } {{B}}]_{\\varphi \\varphi }, 2[\\varvec{nabla } {{B}}]_{θ \\varphi }}. It is shown that the estimated combinations \\varvec{Γ }^{(1)} and \\varvec{Γ }^{(2)} produce similar signal content as the theoretical radial gradient \\varvec{Γ }^{(0)} = {[\\varvec{nabla } {{B}}]_{rr}}. These results demonstrate the ability of multi-satellite missions such as Swarm, which cannot directly measure the radial gradient, to retrieve similar signal content by means of the horizontal gradients. Finally, lithospheric field models are derived using the gradient combinations \\varvec{Γ }^{(1)} and \\varvec{Γ }^{(2)} and compared with models derived from traditional vector and gradient data. The model resulting from \\varvec{Γ }^{(1)} leads to a very similar, and in particular cases improved, model compared to models retrieved by using approximately three times more data, i.e., a full set of vector, North-South and East-West gradients. This demonstrates the high information content of \\varvec{Γ }^{(1)}.

Multiclass Reduced-Set Support Vector Machines

NASA Technical Reports Server (NTRS)

Tang, Benyang; Mazzoni, Dominic

2006-01-01

There are well-established methods for reducing the number of support vectors in a trained binary support vector machine, often with minimal impact on accuracy. We show how reduced-set methods can be applied to multiclass SVMs made up of several binary SVMs, with significantly better results than reducing each binary SVM independently. Our approach is based on Burges' approach that constructs each reduced-set vector as the pre-image of a vector in kernel space, but we extend this by recomputing the SVM weights and bias optimally using the original SVM objective function. This leads to greater accuracy for a binary reduced-set SVM, and also allows vectors to be 'shared' between multiple binary SVMs for greater multiclass accuracy with fewer reduced-set vectors. We also propose computing pre-images using differential evolution, which we have found to be more robust than gradient descent alone. We show experimental results on a variety of problems and find that this new approach is consistently better than previous multiclass reduced-set methods, sometimes with a dramatic difference.

Consideration of Gravity Gradient Stabilization for Orion

DTIC Science & Technology

1989-03-01

AND ERIC ANDionl STABILIZATION TION. MAY NEED SECOND CONTROL SYSTEM TO CONTROL OVERALL ANGULAR MOMENTUM I MOMENTUM DUMPING I IN RESPONSE TO...FURTHER EXPERIENCE IS GAINED RPEFERS TO ANY DEVICE THAT MAY BEl USED Ift A PRIOCESS TOE ECHANGE ANGULAR MOMENTUM WITH THME SPACIECRAFTI BODY Figure 5...rotating with angular velocity w relative to XYZ. If unit vectors along the X, Y, and Z axes are ij, and k, respectively, the vector r can be written

Altered Orientation and Flight Paths of Pigeons Reared on Gravity Anomalies: A GPS Tracking Study

PubMed Central

Blaser, Nicole; Guskov, Sergei I.; Meskenaite, Virginia; Kanevskyi, Valerii A.; Lipp, Hans-Peter

2013-01-01

The mechanisms of pigeon homing are still not understood, in particular how they determine their position at unfamiliar locations. The “gravity vector” theory holds that pigeons memorize the gravity vector at their home loft and deduct home direction and distance from the angular difference between memorized and actual gravity vector. However, the gravity vector is tilted by different densities in the earth crust leading to gravity anomalies. We predicted that pigeons reared on different gravity anomalies would show different initial orientation and also show changes in their flight path when crossing a gravity anomaly. We reared one group of pigeons in a strong gravity anomaly with a north-to-south gravity gradient, and the other group of pigeons in a normal area but on a spot with a strong local anomaly with a west-to-east gravity gradient. After training over shorter distances, pigeons were released from a gravitationally and geomagnetically normal site 50 km north in the same direction for both home lofts. As expected by the theory, the two groups of pigeons showed divergent initial orientation. In addition, some of the GPS-tracked pigeons also showed changes in their flight paths when crossing gravity anomalies. We conclude that even small local gravity anomalies at the birth place of pigeons may have the potential to bias the map sense of pigeons, while reactivity to gravity gradients during flight was variable and appeared to depend on individual navigational strategies and frequency of position updates. PMID:24194860

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous, Anisotropic Media

NASA Astrophysics Data System (ADS)

Zhou, Bing; Greenhalgh, S. A.

2011-10-01

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called `the perturbation method' and `the matrix method', to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green's function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green's function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green's function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.

Implementation of a Balance Operator in NCOM

DTIC Science & Technology

2016-04-07

the background temperature Tb and salinity Sb fields do), f is the Coriolis parameter, k is the vertical unit vector, ∇ is the horizontal gradient, p... effectively used as a natural metric in the space of cost function gradients. The associated geometry inhibits descent in the unbalanced directions...28) where f is the local Coriolis parameter, ∆yv is the local grid spacing in the y direction at a v point, and the overbars indicates horizontal

Temperature sensing by primary roots of maize

NASA Technical Reports Server (NTRS)

Poff, K. L.

1990-01-01

Zea mays L. seedlings, grown on agar plates at 26 degrees C, reoriented the original vertical direction of their primary root when exposed to a thermal gradient applied perpendicular to the gravity vector. The magnitude and direction of curvature can not be explained simply by either a temperature or a humidity effect on root elongation. It is concluded that primary roots of maize sense temperature gradients in addition to sensing the gravitational force.

Full magnetic gradient tensor from triaxial aeromagnetic gradient measurements: Calculation and application

NASA Astrophysics Data System (ADS)

Luo, Yao; Wu, Mei-Ping; Wang, Ping; Duan, Shu-Ling; Liu, Hao-Jun; Wang, Jin-Long; An, Zhan-Feng

2015-09-01

The full magnetic gradient tensor (MGT) refers to the spatial change rate of the three field components of the geomagnetic field vector along three mutually orthogonal axes. The tensor is of use to geological mapping, resources exploration, magnetic navigation, and others. However, it is very difficult to measure the full magnetic tensor gradient using existing engineering technology. We present a method to use triaxial aeromagnetic gradient measurements for deriving the full MGT. The method uses the triaxial gradient data and makes full use of the variation of the magnetic anomaly modulus in three dimensions to obtain a self-consistent magnetic tensor gradient. Numerical simulations show that the full MGT data obtained with the proposed method are of high precision and satisfy the requirements of data processing. We selected triaxial aeromagnetic gradient data from the Hebei Province for calculating the full MGT. Data processing shows that using triaxial tensor gradient data allows to take advantage of the spatial rate of change of the total field in three dimensions and suppresses part of the independent noise in the aeromagnetic gradient. The calculated tensor components have improved resolution, and the transformed full tensor gradient satisfies the requirement of geological mapping and interpretation.

Implementation of a Balance Operator in NCOM

DTIC Science & Technology

2016-04-07

the background temperature Tb and salinity Sb fields do), f is the Coriolis parameter, k is the vertical unit vector, ∇ is the horizontal gradient, p... effectively used as a natural metric in the space of cost function gradients. The associated geometry inhibits descent in the unbalanced directions and...28) where f is the local Coriolis parameter, ∆yv is the local grid spacing in the y direction at a v point, and the overbars indicates horizontal

Multiple-Point Temperature Gradient Algorithm for Ring Laser Gyroscope Bias Compensation

PubMed Central

Li, Geng; Zhang, Pengfei; Wei, Guo; Xie, Yuanping; Yu, Xudong; Long, Xingwu

2015-01-01

To further improve ring laser gyroscope (RLG) bias stability, a multiple-point temperature gradient algorithm is proposed for RLG bias compensation in this paper. Based on the multiple-point temperature measurement system, a complete thermo-image of the RLG block is developed. Combined with the multiple-point temperature gradients between different points of the RLG block, the particle swarm optimization algorithm is used to tune the support vector machine (SVM) parameters, and an optimized design for selecting the thermometer locations is also discussed. The experimental results validate the superiority of the introduced method and enhance the precision and generalizability in the RLG bias compensation model. PMID:26633401

Evidence that implicit assumptions of ‘no evolution’ of disease vectors in changing environments can be violated on a rapid timescale

PubMed Central

Egizi, Andrea; Fefferman, Nina H.; Fonseca, Dina M.

2015-01-01

Projected impacts of climate change on vector-borne disease dynamics must consider many variables relevant to hosts, vectors and pathogens, including how altered environmental characteristics might affect the spatial distributions of vector species. However, many predictive models for vector distributions consider their habitat requirements to be fixed over relevant time-scales, when they may actually be capable of rapid evolutionary change and even adaptation. We examine the genetic signature of a spatial expansion by an invasive vector into locations with novel temperature conditions compared to its native range as a proxy for how existing vector populations may respond to temporally changing habitat. Specifically, we compare invasions into different climate ranges and characterize the importance of selection from the invaded habitat. We demonstrate that vector species can exhibit evolutionary responses (altered allelic frequencies) to a temperature gradient in as little as 7–10 years even in the presence of high gene flow, and further, that this response varies depending on the strength of selection. We interpret these findings in the context of climate change predictions for vector populations and emphasize the importance of incorporating vector evolution into models of future vector-borne disease dynamics. PMID:25688024

DOE Office of Scientific and Technical Information (OSTI.GOV)

Emami, Razieh; Mukohyama, Shinji; Namba, Ryo

Many models of inflation driven by vector fields alone have been known to be plagued by pathological behaviors, namely ghost and/or gradient instabilities. In this work, we seek a new class of vector-driven inflationary models that evade all of the mentioned instabilities. We build our analysis on the Generalized Proca Theory with an extension to three vector fields to realize isotropic expansion. We obtain the conditions required for quasi de-Sitter solutions to be an attractor analogous to the standard slow-roll one and those for their stability at the level of linearized perturbations. Identifying the remedy to the existing unstable models,more » we provide a simple example and explicitly show its stability. This significantly broadens our knowledge on vector inflationary scenarios, reviving potential phenomenological interests for this class of models.« less

Analysis of Particle Content of Recombinant Adeno-Associated Virus Serotype 8 Vectors by Ion-Exchange Chromatography

PubMed Central

Lock, Martin; Alvira, Mauricio R.

2012-01-01

Abstract Advances in adeno-associated virus (AAV)-mediated gene therapy have brought the possibility of commercial manufacturing of AAV vectors one step closer. To realize this prospect, a parallel effort with the goal of ever-increasing sophistication for AAV vector production technology and supporting assays will be required. Among the important release assays for a clinical gene therapy product, those monitoring potentially hazardous contaminants are most critical for patient safety. A prominent contaminant in many AAV vector preparations is vector particles lacking a genome, which can substantially increase the dose of AAV capsid proteins and lead to possible unwanted immunological consequences. Current methods to determine empty particle content suffer from inconsistency, are adversely affected by contaminants, or are not applicable to all serotypes. Here we describe the development of an ion-exchange chromatography-based assay that permits the rapid separation and relative quantification of AAV8 empty and full vector particles through the application of shallow gradients and a strong anion-exchange monolith chromatography medium. PMID:22428980

A generalized nonlocal vector calculus

NASA Astrophysics Data System (ADS)

Alali, Bacim; Liu, Kuo; Gunzburger, Max

2015-10-01

A nonlocal vector calculus was introduced in Du et al. (Math Model Meth Appl Sci 23:493-540, 2013) that has proved useful for the analysis of the peridynamics model of nonlocal mechanics and nonlocal diffusion models. A formulation is developed that provides a more general setting for the nonlocal vector calculus that is independent of particular nonlocal models. It is shown that general nonlocal calculus operators are integral operators with specific integral kernels. General nonlocal calculus properties are developed, including nonlocal integration by parts formula and Green's identities. The nonlocal vector calculus introduced in Du et al. (Math Model Meth Appl Sci 23:493-540, 2013) is shown to be recoverable from the general formulation as a special example. This special nonlocal vector calculus is used to reformulate the peridynamics equation of motion in terms of the nonlocal gradient operator and its adjoint. A new example of nonlocal vector calculus operators is introduced, which shows the potential use of the general formulation for general nonlocal models.

SU-F-J-95: Impact of Shape Complexity On the Accuracy of Gradient-Based PET Volume Delineation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dance, M; Wu, G; Gao, Y

2016-06-15

Purpose: Explore correlation of tumor complexity shape with PET target volume accuracy when delineated with gradient-based segmentation tool. Methods: A total of 24 clinically realistic digital PET Monte Carlo (MC) phantoms of NSCLC were used in the study. The phantom simulated 29 thoracic lesions (lung primary and mediastinal lymph nodes) of varying size, shape, location, and {sup 18}F-FDG activity. A program was developed to calculate a curvature vector along the outline and the standard deviation of this vector was used as a metric to quantify a shape’s “complexity score”. This complexity score was calculated for standard geometric shapes and MC-generatedmore » target volumes in PET phantom images. All lesions were contoured using a commercially available gradient-based segmentation tool and the differences in volume from the MC-generated volumes were calculated as the measure of the accuracy of segmentation. Results: The average absolute percent difference in volumes between the MC-volumes and gradient-based volumes was 11% (0.4%–48.4%). The complexity score showed strong correlation with standard geometric shapes. However, no relationship was found between the complexity score and the accuracy of segmentation by gradient-based tool on MC simulated tumors (R{sup 2} = 0.156). When the lesions were grouped into primary lung lesions and mediastinal/mediastinal adjacent lesions, the average absolute percent difference in volumes were 6% and 29%, respectively. The former group is more isolated and the latter is more surround by tissues with relatively high SUV background. Conclusion: The complexity shape of NSCLC lesions has little effect on the accuracy of the gradient-based segmentation method and thus is not a good predictor of uncertainty in target volume delineation. Location of lesion within a relatively high SUV background may play a more significant role in the accuracy of gradient-based segmentation.« less

Towards human behavior recognition based on spatio temporal features and support vector machines

NASA Astrophysics Data System (ADS)

Ghabri, Sawsen; Ouarda, Wael; Alimi, Adel M.

2017-03-01

Security and surveillance are vital issues in today's world. The recent acts of terrorism have highlighted the urgent need for efficient surveillance. There is indeed a need for an automated system for video surveillance which can detect identity and activity of person. In this article, we propose a new paradigm to recognize an aggressive human behavior such as boxing action. Our proposed system for human activity detection includes the use of a fusion between Spatio Temporal Interest Point (STIP) and Histogram of Oriented Gradient (HoG) features. The novel feature called Spatio Temporal Histogram Oriented Gradient (STHOG). To evaluate the robustness of our proposed paradigm with a local application of HoG technique on STIP points, we made experiments on KTH human action dataset based on Multi Class Support Vector Machines classification. The proposed scheme outperforms basic descriptors like HoG and STIP to achieve 82.26% us an accuracy value of classification rate.

A Concept of Cross-Ferroic Plasma Turbulence

PubMed Central

Inagaki, S.; Kobayashi, T.; Kosuga, Y.; Itoh, S.-I.; Mitsuzono, T.; Nagashima, Y.; Arakawa, H.; Yamada, T.; Miwa, Y.; Kasuya, N.; Sasaki, M.; Lesur, M.; Fujisawa, A.; Itoh, K.

2016-01-01

The variety of scalar and vector fields in laboratory and nature plasmas is formed by plasma turbulence. Drift-wave fluctuations, driven by density gradients in magnetized plasmas, are known to relax the density gradient while they can generate flows. On the other hand, the sheared flow in the direction of magnetic fields causes Kelvin-Helmholtz type instabilities, which mix particle and momentum. These different types of fluctuations coexist in laboratory and nature, so that the multiple mechanisms for structural formation exist in extremely non-equilibrium plasmas. Here we report the discovery of a new order in plasma turbulence, in which chained structure formation is realized by cross-interaction between inhomogeneities of scalar and vector fields. The concept of cross-ferroic turbulence is developed, and the causal relation in the multiple mechanisms behind structural formation is identified, by measuring the relaxation rate and dissipation power caused by the complex turbulence-driven flux. PMID:26917218

Renormalizability of the gradient flow in the 2D O(N) non-linear sigma model

NASA Astrophysics Data System (ADS)

Makino, Hiroki; Suzuki, Hiroshi

2015-03-01

It is known that the gauge field and its composite operators evolved by the Yang-Mills gradient flow are ultraviolet (UV) finite without any multiplicative wave function renormalization. In this paper, we prove that the gradient flow in the 2D O(N) non-linear sigma model possesses a similar property: The flowed N-vector field and its composite operators are UV finite without multiplicative wave function renormalization. Our proof in all orders of perturbation theory uses a (2+1)-dimensional field theoretical representation of the gradient flow, which possesses local gauge invariance without gauge field. As an application of the UV finiteness of the gradient flow, we construct the energy-momentum tensor in the lattice formulation of the O(N) non-linear sigma model that automatically restores the correct normalization and the conservation law in the continuum limit.

Wideband helicity dependent spoof surface plasmon polaritons coupling metasurface based on dispersion design.

PubMed

Dong, Guoxiang; Shi, Hongyu; He, Yuchen; Zhang, Anxue; Wei, Xiaoyong; Zhuang, Yongyong; Du, Bai; Xia, Song; Xu, Zhuo

2016-12-06

The surface plasmon polaritons (SPPs) have many potential application due to their local field enhancement and sub-wavelength characteristics. Recently, the gradient metasurface is introduced to couple the spoof SPPs in microwave frequency band. One of the most important issue which should be solved is the narrowband of spoof SPPs coupling on the gradient metasurface. Here, the metasurface is proposed to achieve the wideband helicity dependent directional spoof SPPs coupling for circular polarized light. Our research show that the coupling frequency of spoof SPPs on the gradient metasurface is determined by the dispersion of the metasurface, so the coupling frequency can be controlled by dispersion design. The careful design of each cell geometric parameters has provided many appropriate dispersion relations possessed by just one metasurface. The wave vector matching between the propagating wave and the spoof SPPs has been achieved at several frequencies for certain wave vector provided by the metasurface, which leads to wideband spoof SPPs coupling. This work has shown that wideband helicity dependent directional spoof SPPs coupling has been achieved with a high efficiency. Hence, the proposed wideband spoof SPPs coupling presents the improvement in practice applications.

Quantitative structure-retention relationship models for the prediction of the reversed-phase HPLC gradient retention based on the heuristic method and support vector machine.

PubMed

Du, Hongying; Wang, Jie; Yao, Xiaojun; Hu, Zhide

2009-01-01

The heuristic method (HM) and support vector machine (SVM) were used to construct quantitative structure-retention relationship models by a series of compounds to predict the gradient retention times of reversed-phase high-performance liquid chromatography (HPLC) in three different columns. The aims of this investigation were to predict the retention times of multifarious compounds, to find the main properties of the three columns, and to indicate the theory of separation procedures. In our method, we correlated the retention times of many diverse structural analytes in three columns (Symmetry C18, Chromolith, and SG-MIX) with their representative molecular descriptors, calculated from the molecular structures alone. HM was used to select the most important molecular descriptors and build linear regression models. Furthermore, non-linear regression models were built using the SVM method; the performance of the SVM models were better than that of the HM models, and the prediction results were in good agreement with the experimental values. This paper could give some insights into the factors that were likely to govern the gradient retention process of the three investigated HPLC columns, which could theoretically supervise the practical experiment.

Wideband helicity dependent spoof surface plasmon polaritons coupling metasurface based on dispersion design

PubMed Central

Dong, Guoxiang; Shi, Hongyu; He, Yuchen; Zhang, Anxue; Wei, Xiaoyong; Zhuang, Yongyong; Du, Bai; Xia, Song; Xu, Zhuo

2016-01-01

The surface plasmon polaritons (SPPs) have many potential application due to their local field enhancement and sub-wavelength characteristics. Recently, the gradient metasurface is introduced to couple the spoof SPPs in microwave frequency band. One of the most important issue which should be solved is the narrowband of spoof SPPs coupling on the gradient metasurface. Here, the metasurface is proposed to achieve the wideband helicity dependent directional spoof SPPs coupling for circular polarized light. Our research show that the coupling frequency of spoof SPPs on the gradient metasurface is determined by the dispersion of the metasurface, so the coupling frequency can be controlled by dispersion design. The careful design of each cell geometric parameters has provided many appropriate dispersion relations possessed by just one metasurface. The wave vector matching between the propagating wave and the spoof SPPs has been achieved at several frequencies for certain wave vector provided by the metasurface, which leads to wideband spoof SPPs coupling. This work has shown that wideband helicity dependent directional spoof SPPs coupling has been achieved with a high efficiency. Hence, the proposed wideband spoof SPPs coupling presents the improvement in practice applications. PMID:27922132

3PE: A Tool for Estimating Groundwater Flow Vectors

EPA Science Inventory

Evaluation of hydraulic gradients and the associated groundwater flow rates and directions is a fundamental aspect of hydrogeologic characterization. Many methods, ranging in complexity from simple three-point solution techniques to complex numerical models of groundwater flow, ...

RES: Regularized Stochastic BFGS Algorithm

NASA Astrophysics Data System (ADS)

Mokhtari, Aryan; Ribeiro, Alejandro

2014-12-01

RES, a regularized stochastic version of the Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton method is proposed to solve convex optimization problems with stochastic objectives. The use of stochastic gradient descent algorithms is widespread, but the number of iterations required to approximate optimal arguments can be prohibitive in high dimensional problems. Application of second order methods, on the other hand, is impracticable because computation of objective function Hessian inverses incurs excessive computational cost. BFGS modifies gradient descent by introducing a Hessian approximation matrix computed from finite gradient differences. RES utilizes stochastic gradients in lieu of deterministic gradients for both, the determination of descent directions and the approximation of the objective function's curvature. Since stochastic gradients can be computed at manageable computational cost RES is realizable and retains the convergence rate advantages of its deterministic counterparts. Convergence results show that lower and upper bounds on the Hessian egeinvalues of the sample functions are sufficient to guarantee convergence to optimal arguments. Numerical experiments showcase reductions in convergence time relative to stochastic gradient descent algorithms and non-regularized stochastic versions of BFGS. An application of RES to the implementation of support vector machines is developed.

Characterization of thermotropism in primary roots of maize: dependence on temperature and temperature gradient, and interaction with gravitropism

NASA Technical Reports Server (NTRS)

Poff, K. L.

1991-01-01

Thermotropism in primary roots of Zea mays L. was studied with respect to gradient strength (degrees C cm-1), temperature of exposure within a gradient, pre-treatment temperature, and gravitropic stimulation. The magnitude of the response decreased with gradient strength. Maximum thermotropism was independent of gradient strength and pre-treatment temperature. The range of temperature for positive and negative thermotropism did not change with pre-treatment temperature. However, the exact range of temperatures for positive and negative thermotropism varied with gradient strengths. In general, temperatures of exposure lower than 25 degrees C resulted in positive tropic responses while temperatures of exposure of 39 degrees C or more resulted in negative tropic responses. Thermotropism was shown to modify and reverse the normal gravitropic curvature of a horizontal root when thermal gradients were applied opposite the 1 g vector. It is concluded that root thermotropism is a consequence of thermal sensing and that the curvature of the primary root results from the interaction of the thermal and gravitational sensing systems.

Fast temporal neural learning using teacher forcing

NASA Technical Reports Server (NTRS)

Toomarian, Nikzad (Inventor); Bahren, Jacob (Inventor)

1992-01-01

A neural network is trained to output a time dependent target vector defined over a predetermined time interval in response to a time dependent input vector defined over the same time interval by applying corresponding elements of the error vector, or difference between the target vector and the actual neuron output vector, to the inputs of corresponding output neurons of the network as corrective feedback. This feedback decreases the error and quickens the learning process, so that a much smaller number of training cycles are required to complete the learning process. A conventional gradient descent algorithm is employed to update the neural network parameters at the end of the predetermined time interval. The foregoing process is repeated in repetitive cycles until the actual output vector corresponds to the target vector. In the preferred embodiment, as the overall error of the neural network output decreasing during successive training cycles, the portion of the error fed back to the output neurons is decreased accordingly, allowing the network to learn with greater freedom from teacher forcing as the network parameters converge to their optimum values. The invention may also be used to train a neural network with stationary training and target vectors.

Fast temporal neural learning using teacher forcing

NASA Technical Reports Server (NTRS)

Toomarian, Nikzad (Inventor); Bahren, Jacob (Inventor)

1995-01-01

A neural network is trained to output a time dependent target vector defined over a predetermined time interval in response to a time dependent input vector defined over the same time interval by applying corresponding elements of the error vector, or difference between the target vector and the actual neuron output vector, to the inputs of corresponding output neurons of the network as corrective feedback. This feedback decreases the error and quickens the learning process, so that a much smaller number of training cycles are required to complete the learning process. A conventional gradient descent algorithm is employed to update the neural network parameters at the end of the predetermined time interval. The foregoing process is repeated in repetitive cycles until the actual output vector corresponds to the target vector. In the preferred embodiment, as the overall error of the neural network output decreasing during successive training cycles, the portion of the error fed back to the output neurons is decreased accordingly, allowing the network to learn with greater freedom from teacher forcing as the network parameters converge to their optimum values. The invention may also be used to train a neural network with stationary training and target vectors.

Optimization Based Efficiencies in First Order Reliability Analysis

NASA Technical Reports Server (NTRS)

Peck, Jeffrey A.; Mahadevan, Sankaran

2003-01-01

This paper develops a method for updating the gradient vector of the limit state function in reliability analysis using Broyden's rank one updating technique. In problems that use commercial code as a black box, the gradient calculations are usually done using a finite difference approach, which becomes very expensive for large system models. The proposed method replaces the finite difference gradient calculations in a standard first order reliability method (FORM) with Broyden's Quasi-Newton technique. The resulting algorithm of Broyden updates within a FORM framework (BFORM) is used to run several example problems, and the results compared to standard FORM results. It is found that BFORM typically requires fewer functional evaluations that FORM to converge to the same answer.

Mosquito communities and disease risk influenced by land use change and seasonality in the Australian tropics.

PubMed

Meyer Steiger, Dagmar B; Ritchie, Scott A; Laurance, Susan G W

2016-07-07

Anthropogenic land use changes have contributed considerably to the rise of emerging and re-emerging mosquito-borne diseases. These diseases appear to be increasing as a result of the novel juxtapositions of habitats and species that can result in new interchanges of vectors, diseases and hosts. We studied whether the mosquito community structure varied between habitats and seasons and whether known disease vectors displayed habitat preferences in tropical Australia. Using CDC model 512 traps, adult mosquitoes were sampled across an anthropogenic disturbance gradient of grassland, rainforest edge and rainforest interior habitats, in both the wet and dry seasons. Nonmetric multidimensional scaling (NMS) ordinations were applied to examine major gradients in the composition of mosquito and vector communities. We captured ~13,000 mosquitoes from 288 trap nights across four study sites. A community analysis identified 29 species from 7 genera. Even though mosquito abundance and richness were similar between the three habitats, the community composition varied significantly in response to habitat type. The mosquito community in rainforest interiors was distinctly different to the community in grasslands, whereas forest edges acted as an ecotone with shared communities from both forest interiors and grasslands. We found two community patterns that will influence disease risk at out study sites, first, that disease vectoring mosquito species occurred all year round. Secondly, that anthropogenic grasslands adjacent to rainforests may increase the probability of novel disease transmission through changes to the vector community on rainforest edges, as most disease transmitting species predominantly occurred in grasslands. Our results indicate that the strong influence of anthropogenic land use change on mosquito communities could have potential implications for pathogen transmission to humans and wildlife.

Generalized Analysis Tools for Multi-Spacecraft Missions

NASA Astrophysics Data System (ADS)

Chanteur, G. M.

2011-12-01

Analysis tools for multi-spacecraft missions like CLUSTER or MMS have been designed since the end of the 90's to estimate gradients of fields or to characterize discontinuities crossed by a cluster of spacecraft. Different approaches have been presented and discussed in the book "Analysis Methods for Multi-Spacecraft Data" published as Scientific Report 001 of the International Space Science Institute in Bern, Switzerland (G. Paschmann and P. Daly Eds., 1998). On one hand the approach using methods of least squares has the advantage to apply to any number of spacecraft [1] but is not convenient to perform analytical computation especially when considering the error analysis. On the other hand the barycentric approach is powerful as it provides simple analytical formulas involving the reciprocal vectors of the tetrahedron [2] but appears limited to clusters of four spacecraft. Moreover the barycentric approach allows to derive theoretical formulas for errors affecting the estimators built from the reciprocal vectors [2,3,4]. Following a first generalization of reciprocal vectors proposed by Vogt et al [4] and despite the present lack of projects with more than four spacecraft we present generalized reciprocal vectors for a cluster made of any number of spacecraft : each spacecraft is given a positive or nul weight. The non-coplanarity of at least four spacecraft with strictly positive weights is a necessary and sufficient condition for this analysis to be enabled. Weights given to spacecraft allow to minimize the influence of some spacecraft if its location or the quality of its data are not appropriate, or simply to extract subsets of spacecraft from the cluster. Estimators presented in [2] are generalized within this new frame except for the error analysis which is still under investigation. References [1] Harvey, C. C.: Spatial Gradients and the Volumetric Tensor, in: Analysis Methods for Multi-Spacecraft Data, G. Paschmann and P. Daly (eds.), pp. 307-322, ISSI SR-001, 1998. [2] Chanteur, G.: Spatial Interpolation for Four Spacecraft: Theory, in: Analysis Methods for Multi-Spacecraft Data, G. Paschmann and P. Daly (eds.), pp. 371-393, ISSI SR-001, 1998. [3] Chanteur, G.: Accuracy of field gradient estimations by Cluster: Explanation of its dependency upon elongation and planarity of the tetrahedron, pp. 265-268, ESA SP-449, 2000. [4] Vogt, J., Paschmann, G., and Chanteur, G.: Reciprocal Vectors, pp. 33-46, ISSI SR-008, 2008.

Studies on deflection area vectors of QRS and T and ventricular gradient in right ventricular hypertrophy.

PubMed

Kawaguchi, Y

1985-04-01

QRS deflection area vector (Aqrs), T deflection area vector (At) and ventricular gradient (G) in right ventricular hypertrophy were studied in 53 subjects divided on the basis of cardiac catheterization data into four subgroups; normal controls, mild MS group, right ventricular pressure overload group and right ventricular volume overload group. Aqrs, At and G of the four subgroups were calculated using a microcomputer and compared. Aqrs in right ventricular pressure overload group and volume overload group was shifted to the right and slightly anteriorly from that in normal control group. At in right ventricular pressure overload group and volume overload group was shifted slightly upwards and significantly posteriorly from that in the normal control and mild MS groups. G in right ventricular pressure overload group and volume overload group was shifted to the right and significantly posteriorly from that in normal control and mild MS groups. Using multivariative analysis, we developed criteria for diagnosing right ventricular hypertrophy with At: 0.059At(Z) - 0.0145 [At] - 0.2608 less than or equal to 0. Application of this criteria achieved 82.4% (28 of 34) sensitivity in the patients with right ventricular hypertrophy and 90.9% (10 of 11) specificity in the normal control subjects.

Minimum impulse three-body trajectories.

NASA Technical Reports Server (NTRS)

D'Amario, L.; Edelbaum, T. N.

1973-01-01

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

Lattice QCD calculation using VPP500

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kim, Seyong; Ohta, Shigemi

1995-02-01

A new vector parallel supercomputer, Fujitsu VPP500, was installed at RIKEN earlier this year. It consists of 30 vector computers, each with 1.6 GFLOPS peak speed and 256 MB memory, connected by a crossbar switch with 400 MB/s peak data transfer rate each way between any pair of nodes. The authors developed a Fortran lattice QCD simulation code for it. It runs at about 1.1 GFLOPS sustained per node for Metropolis pure-gauge update, and about 0.8 GFLOPS sustained per node for conjugate gradient inversion of staggered fermion matrix.

Cardiac Gene Therapy: Optimization of Gene Delivery Techniques In Vivo

PubMed Central

Katz, Michael G.; Swain, JaBaris D.; White, Jennifer D.; Low, David; Stedman, Hansell

2010-01-01

Abstract Vector-mediated cardiac gene therapy holds tremendous promise as a translatable platform technology for treating many cardiovascular diseases. The ideal technique is one that is efficient and practical, allowing for global cardiac gene expression, while minimizing collateral expression in other organs. Here we survey the available in vivo vector-mediated cardiac gene delivery methods—including transcutaneous, intravascular, intramuscular, and cardiopulmonary bypass techniques—with consideration of the relative merits and deficiencies of each. Review of available techniques suggests that an optimal method for vector-mediated gene delivery to the large animal myocardium would ideally employ retrograde and/or anterograde transcoronary gene delivery,extended vector residence time in the coronary circulation, an increased myocardial transcapillary gradient using physical methods, increased endothelial permeability with pharmacological agents, minimal collateral gene expression by isolation of the cardiac circulation from the systemic, and have low immunogenicity. PMID:19947886

Silencing Genes in the Heart.

PubMed

Fechner, Henry; Vetter, Roland; Kurreck, Jens; Poller, Wolfgang

2017-01-01

Silencing of cardiac genes by RNA interference (RNAi) has developed into a powerful new method to treat cardiac diseases. Small interfering (si)RNAs are the inducers of RNAi, but cultured primary cardiomyocytes and heart are highly resistant to siRNA transfection. This can be overcome by delivery of small hairpin (sh)RNAs or artificial microRNA (amiRNAs) by cardiotropic adeno-associated virus (AAV) vectors. Here we describe as example of the silencing of a cardiac gene, the generation and cloning of shRNA, and amiRNAs directed against the cardiac protein phospholamban. We further describe the generation of AAV shuttle plasmids with self complementary vector genomes, the production of AAV vectors in roller bottles, and their purification via iodixanol gradient centrifugation and concentration with filter systems. Finally we describe the preparation of primary neonatal rat cardiomyocytes (PNRC), the transduction of PNRC with AAV vectors, and the maintenance of the transduced cell culture.

Spectral functions with the density matrix renormalization group: Krylov-space approach for correction vectors

DOE Office of Scientific and Technical Information (OSTI.GOV)

None, None

Frequency-dependent correlations, such as the spectral function and the dynamical structure factor, help illustrate condensed matter experiments. Within the density matrix renormalization group (DMRG) framework, an accurate method for calculating spectral functions directly in frequency is the correction-vector method. The correction vector can be computed by solving a linear equation or by minimizing a functional. Our paper proposes an alternative to calculate the correction vector: to use the Krylov-space approach. This paper also studies the accuracy and performance of the Krylov-space approach, when applied to the Heisenberg, the t-J, and the Hubbard models. The cases we studied indicate that themore » Krylov-space approach can be more accurate and efficient than the conjugate gradient, and that the error of the former integrates best when a Krylov-space decomposition is also used for ground state DMRG.« less

Spectral functions with the density matrix renormalization group: Krylov-space approach for correction vectors

DOE PAGES

None, None

2016-11-21

Frequency-dependent correlations, such as the spectral function and the dynamical structure factor, help illustrate condensed matter experiments. Within the density matrix renormalization group (DMRG) framework, an accurate method for calculating spectral functions directly in frequency is the correction-vector method. The correction vector can be computed by solving a linear equation or by minimizing a functional. Our paper proposes an alternative to calculate the correction vector: to use the Krylov-space approach. This paper also studies the accuracy and performance of the Krylov-space approach, when applied to the Heisenberg, the t-J, and the Hubbard models. The cases we studied indicate that themore » Krylov-space approach can be more accurate and efficient than the conjugate gradient, and that the error of the former integrates best when a Krylov-space decomposition is also used for ground state DMRG.« less

Toward lattice fractional vector calculus

NASA Astrophysics Data System (ADS)

Tarasov, Vasily E.

2014-09-01

An analog of fractional vector calculus for physical lattice models is suggested. We use an approach based on the models of three-dimensional lattices with long-range inter-particle interactions. The lattice analogs of fractional partial derivatives are represented by kernels of lattice long-range interactions, where the Fourier series transformations of these kernels have a power-law form with respect to wave vector components. In the continuum limit, these lattice partial derivatives give derivatives of non-integer order with respect to coordinates. In the three-dimensional description of the non-local continuum, the fractional differential operators have the form of fractional partial derivatives of the Riesz type. As examples of the applications of the suggested lattice fractional vector calculus, we give lattice models with long-range interactions for the fractional Maxwell equations of non-local continuous media and for the fractional generalization of the Mindlin and Aifantis continuum models of gradient elasticity.

Rotating-frame gradient fields for magnetic resonance imaging and nuclear magnetic resonance in low fields

DOEpatents

Bouchard, Louis-Serge; Pines, Alexander; Demas, Vasiliki

2014-01-21

A system and method for Fourier encoding a nuclear magnetic resonance (NMR) signal is disclosed. A static magnetic field B.sub.0 is provided along a first direction. An NMR signal from the sample is Fourier encoded by applying a rotating-frame gradient field B.sub.G superimposed on the B.sub.0, where the B.sub.G comprises a vector component rotating in a plane perpendicular to the first direction at an angular frequency .omega.in a laboratory frame. The Fourier-encoded NMR signal is detected.

SAGRAD: A Program for Neural Network Training with Simulated Annealing and the Conjugate Gradient Method.

PubMed

Bernal, Javier; Torres-Jimenez, Jose

2015-01-01

SAGRAD (Simulated Annealing GRADient), a Fortran 77 program for computing neural networks for classification using batch learning, is discussed. Neural network training in SAGRAD is based on a combination of simulated annealing and Møller's scaled conjugate gradient algorithm, the latter a variation of the traditional conjugate gradient method, better suited for the nonquadratic nature of neural networks. Different aspects of the implementation of the training process in SAGRAD are discussed, such as the efficient computation of gradients and multiplication of vectors by Hessian matrices that are required by Møller's algorithm; the (re)initialization of weights with simulated annealing required to (re)start Møller's algorithm the first time and each time thereafter that it shows insufficient progress in reaching a possibly local minimum; and the use of simulated annealing when Møller's algorithm, after possibly making considerable progress, becomes stuck at a local minimum or flat area of weight space. Outlines of the scaled conjugate gradient algorithm, the simulated annealing procedure and the training process used in SAGRAD are presented together with results from running SAGRAD on two examples of training data.

Integrating the Gradient of the Thin Wire Kernel

NASA Technical Reports Server (NTRS)

Champagne, Nathan J.; Wilton, Donald R.

2008-01-01

A formulation for integrating the gradient of the thin wire kernel is presented. This approach employs a new expression for the gradient of the thin wire kernel derived from a recent technique for numerically evaluating the exact thin wire kernel. This approach should provide essentially arbitrary accuracy and may be used with higher-order elements and basis functions using the procedure described in [4].When the source and observation points are close, the potential integrals over wire segments involving the wire kernel are split into parts to handle the singular behavior of the integrand [1]. The singularity characteristics of the gradient of the wire kernel are different than those of the wire kernel, and the axial and radial components have different singularities. The characteristics of the gradient of the wire kernel are discussed in [2]. To evaluate the near electric and magnetic fields of a wire, the integration of the gradient of the wire kernel needs to be calculated over the source wire. Since the vector bases for current have constant direction on linear wire segments, these integrals reduce to integrals of the form

Using a Gradient Vector to Find Multiple Periodic Oscillations in Suspension Bridge Models

ERIC Educational Resources Information Center

Humphreys, L. D.; McKenna, P. J.

2005-01-01

This paper describes how the method of steepest descent can be used to find periodic solutions of differential equations. Applications to two suspension bridge models are discussed, and the method is used to find non-obvious large-amplitude solutions.

A High Order Finite Difference Scheme with Sharp Shock Resolution for the Euler Equations

NASA Technical Reports Server (NTRS)

Gerritsen, Margot; Olsson, Pelle

1996-01-01

We derive a high-order finite difference scheme for the Euler equations that satisfies a semi-discrete energy estimate, and present an efficient strategy for the treatment of discontinuities that leads to sharp shock resolution. The formulation of the semi-discrete energy estimate is based on a symmetrization of the Euler equations that preserves the homogeneity of the flux vector, a canonical splitting of the flux derivative vector, and the use of difference operators that satisfy a discrete analogue to the integration by parts procedure used in the continuous energy estimate. Around discontinuities or sharp gradients, refined grids are created on which the discrete equations are solved after adding a newly constructed artificial viscosity. The positioning of the sub-grids and computation of the viscosity are aided by a detection algorithm which is based on a multi-scale wavelet analysis of the pressure grid function. The wavelet theory provides easy to implement mathematical criteria to detect discontinuities, sharp gradients and spurious oscillations quickly and efficiently.

MATLAB Simulation of Gradient-Based Neural Network for Online Matrix Inversion

NASA Astrophysics Data System (ADS)

Zhang, Yunong; Chen, Ke; Ma, Weimu; Li, Xiao-Dong

This paper investigates the simulation of a gradient-based recurrent neural network for online solution of the matrix-inverse problem. Several important techniques are employed as follows to simulate such a neural system. 1) Kronecker product of matrices is introduced to transform a matrix-differential-equation (MDE) to a vector-differential-equation (VDE); i.e., finally, a standard ordinary-differential-equation (ODE) is obtained. 2) MATLAB routine "ode45" is introduced to solve the transformed initial-value ODE problem. 3) In addition to various implementation errors, different kinds of activation functions are simulated to show the characteristics of such a neural network. Simulation results substantiate the theoretical analysis and efficacy of the gradient-based neural network for online constant matrix inversion.

Bidirectional composition on lie groups for gradient-based image alignment.

PubMed

Mégret, Rémi; Authesserre, Jean-Baptiste; Berthoumieu, Yannick

2010-09-01

In this paper, a new formulation based on bidirectional composition on Lie groups (BCL) for parametric gradient-based image alignment is presented. Contrary to the conventional approaches, the BCL method takes advantage of the gradients of both template and current image without combining them a priori. Based on this bidirectional formulation, two methods are proposed and their relationship with state-of-the-art gradient based approaches is fully discussed. The first one, i.e., the BCL method, relies on the compositional framework to provide the minimization of the compensated error with respect to an augmented parameter vector. The second one, the projected BCL (PBCL), corresponds to a close approximation of the BCL approach. A comparative study is carried out dealing with computational complexity, convergence rate and frequence of convergence. Numerical experiments using a conventional benchmark show the performance improvement especially for asymmetric levels of noise, which is also discussed from a theoretical point of view.

Sewall Wright's equation Deltaq=(q(1-q) partial differentialw/ partial differentialq)/2w.

PubMed

Edwards, A W

2000-02-01

An equation of Sewall Wright's expresses the change in the frequency of an allele under selection at a multiallelic locus as a function of the gradient of the mean fitness "surface" in the direction in which the relative proportions of the other alleles do not change. An attempt to derive this equation using conventional vector calculus shows that this description leads to a different equation and that the purported gradient in Wright's equation is not a gradient of the mean fitness surface except in the diallelic case, where the two equations are the same. It is further shown that if Fisher's angular transformation is applied to the diallelic case the genic variance is exactly equal to one-eighth of the square of the gradient of the mean fitness with respect to the transformed gene frequency. Copyright 2000 Academic Press.

Production of SV40-derived vectors.

PubMed

Strayer, David S; Mitchell, Christine; Maier, Dawn A; Nichols, Carmen N

2010-06-01

Recombinant simian virus 40 (rSV40)-derived vectors are particularly useful for gene delivery to bone marrow progenitor cells and their differentiated derivatives, certain types of epithelial cells (e.g., hepatocytes), and central nervous system neurons and microglia. They integrate rapidly into cellular DNA to provide long-term gene expression in vitro and in vivo in both resting and dividing cells. Here we describe a protocol for production and purification of these vectors. These procedures require only packaging cells (e.g., COS-7) and circular vector genome DNA. Amplification involves repeated infection of packaging cells with vector produced by transfection. Cotransfection is not required in any step. Viruses are purified by centrifugation using discontinuous sucrose or cesium chloride (CsCl) gradients and resulting vectors are replication-incompetent and contain no detectable wild-type SV40 revertants. These approaches are simple, give reproducible results, and may be used to generate vectors that are deleted only for large T antigen (Tag), or for all SV40-coding sequences capable of carrying up to 5 kb of foreign DNA. These vectors are best applied to long-term expression of proteins normally encoded by mammalian cells or by viruses that infect mammalian cells, or of untranslated RNAs (e.g., RNA interference). The preparative approaches described facilitate application of these vectors and allow almost any laboratory to exploit their strengths for diverse gene delivery applications.

The ecological foundations of transmission potential and vector-borne disease in urban landscapes.

PubMed

LaDeau, Shannon L; Allan, Brian F; Leisnham, Paul T; Levy, Michael Z

2015-07-01

Urban transmission of arthropod-vectored disease has increased in recent decades. Understanding and managing transmission potential in urban landscapes requires integration of sociological and ecological processes that regulate vector population dynamics, feeding behavior, and vector-pathogen interactions in these unique ecosystems. Vectorial capacity is a key metric for generating predictive understanding about transmission potential in systems with obligate vector transmission. This review evaluates how urban conditions, specifically habitat suitability and local temperature regimes, and the heterogeneity of urban landscapes can influence the biologically-relevant parameters that define vectorial capacity: vector density, survivorship, biting rate, extrinsic incubation period, and vector competence.Urban landscapes represent unique mosaics of habitat. Incidence of vector-borne disease in urban host populations is rarely, if ever, evenly distributed across an urban area. The persistence and quality of vector habitat can vary significantly across socio-economic boundaries to influence vector species composition and abundance, often generating socio-economically distinct gradients of transmission potential across neighborhoods.Urban regions often experience unique temperature regimes, broadly termed urban heat islands (UHI). Arthropod vectors are ectothermic organisms and their growth, survival, and behavior are highly sensitive to environmental temperatures. Vector response to UHI conditions is dependent on regional temperature profiles relative to the vector's thermal performance range. In temperate climates UHI can facilitate increased vector development rates while having countervailing influence on survival and feeding behavior. Understanding how urban heat island (UHI) conditions alter thermal and moisture constraints across the vector life cycle to influence transmission processes is an important direction for both empirical and modeling research.There remain persistent gaps in understanding of vital rates and drivers in mosquito-vectored disease systems, and vast holes in understanding for other arthropod vectored diseases. Empirical studies are needed to better understand the physiological constraints and socio-ecological processes that generate heterogeneity in critical transmission parameters, including vector survival and fitness. Likewise, laboratory experiments and transmission models must evaluate vector response to realistic field conditions, including variability in sociological and environmental conditions.

Orientational Order on Surfaces: The Coupling of Topology, Geometry, and Dynamics

NASA Astrophysics Data System (ADS)

Nestler, M.; Nitschke, I.; Praetorius, S.; Voigt, A.

2018-02-01

We consider the numerical investigation of surface bound orientational order using unit tangential vector fields by means of a gradient flow equation of a weak surface Frank-Oseen energy. The energy is composed of intrinsic and extrinsic contributions, as well as a penalization term to enforce the unity of the vector field. Four different numerical discretizations, namely a discrete exterior calculus approach, a method based on vector spherical harmonics, a surface finite element method, and an approach utilizing an implicit surface description, the diffuse interface method, are described and compared with each other for surfaces with Euler characteristic 2. We demonstrate the influence of geometric properties on realizations of the Poincaré-Hopf theorem and show examples where the energy is decreased by introducing additional orientational defects.

The effect of wind and moisture gradients on the arbitrary assignment of cloud motions to a vertical coordinate system in two Sesame cases

NASA Technical Reports Server (NTRS)

Peslen, C. A.; Koch, S. E.; Uccellini, L. W.

1984-01-01

Satellite-derived cloud motion 'wind' vectors (CMV) are increasingly used in mesoscale and in global analyses, and questions have been raised regarding the uncertainty of the level assignment for the CMV. One of two major problems in selecting a level for the CMV is related to uncertainties in assigning the motion vector to either the cloud top or base. The second problem is related to the inability to transfer the 'wind' derived from the CMV at individually specified heights to a standard coordinated surface. The present investigation has the objective to determine if the arbitrary level assignment represents a serious obstacle to the use of cloud motion wind vectors in the mesoscale analysis of a severe storm environment.

A Guided Tour of Mathematical Methods

NASA Astrophysics Data System (ADS)

Snieder, Roel

2009-04-01

1. Introduction; 2. Dimensional analysis; 3. Power series; 4. Spherical and cylindrical co-ordinates; 5. The gradient; 6. The divergence of a vector field; 7. The curl of a vector field; 8. The theorem of Gauss; 9. The theorem of Stokes; 10. The Laplacian; 11. Conservation laws; 12. Scale analysis; 13. Linear algebra; 14. The Dirac delta function; 15. Fourier analysis; 16. Analytic functions; 17. Complex integration; 18. Green's functions: principles; 19. Green's functions: examples; 20. Normal modes; 21. Potential theory; 22. Cartesian tensors; 23. Perturbation theory; 24. Asymptotic evaluation of integrals; 25. Variational calculus; 26. Epilogue, on power and knowledge; References.

A Study of Shaped-Charge Collapse and Jet Formation Using the HEMP (hydrodynamic, Elastic, Magneto, and Plastic) Code and a Comparison with Experimental Observations

DTIC Science & Technology

1984-12-01

Octol** explosive. The experimental charges were lightly confined with aluminum bodies and had cone diameters of 84mm. The charges modelled using HEMP...solved using the following relationships: .Final Final V 0 1 IV sin 9, where Voz aj i teailcmpnn fj~Fnl n Final F where V is the adial component of Fnan h...velocity vector is equal to the vector addition of the flow and 8. MMiles L. Lampson, "The Influence of Convergence - Velocity Gradients on the Formation

Lefschetz thimbles in fermionic effective models with repulsive vector-field

NASA Astrophysics Data System (ADS)

Mori, Yuto; Kashiwa, Kouji; Ohnishi, Akira

2018-06-01

We discuss two problems in complexified auxiliary fields in fermionic effective models, the auxiliary sign problem associated with the repulsive vector-field and the choice of the cut for the scalar field appearing from the logarithmic function. In the fermionic effective models with attractive scalar and repulsive vector-type interaction, the auxiliary scalar and vector fields appear in the path integral after the bosonization of fermion bilinears. When we make the path integral well-defined by the Wick rotation of the vector field, the oscillating Boltzmann weight appears in the partition function. This "auxiliary" sign problem can be solved by using the Lefschetz-thimble path-integral method, where the integration path is constructed in the complex plane. Another serious obstacle in the numerical construction of Lefschetz thimbles is caused by singular points and cuts induced by multivalued functions of the complexified scalar field in the momentum integration. We propose a new prescription which fixes gradient flow trajectories on the same Riemann sheet in the flow evolution by performing the momentum integration in the complex domain.

A contemporary view of the ventricular gradient of Wilson.

PubMed

Plonsey, R

1979-10-01

We have derived quantitative expressions for QRS, T, and QRST areas of the scalar electrocardiogram. The QRST area, or ventricular gradient, is seen to be essentially independent of the activation sequence and to reflect recovery properties of the tissue as weighted by the vector lead field of a given lead. The results are derived for uniform isotropic conditions and under the assumption that the temporal waveforms everywhere are identical except for possible variations in the duration of the plateau. However, it is noted that the results are, probably, valid under anisotropic conditions as well. The examination of ventricular gradients from epicardial and intramural leads should reflect local recovery properties and be a useful tool in study of the physiology of recovery, as well as the study of arrhythmias.

Evaluation of Temperature Gradient in Advanced Automated Directional Solidification Furnace (AADSF) by Numerical Simulation

NASA Technical Reports Server (NTRS)

Bune, Andris V.; Gillies, Donald C.; Lehoczky, Sandor L.

1996-01-01

A numerical model of heat transfer using combined conduction, radiation and convection in AADSF was used to evaluate temperature gradients in the vicinity of the crystal/melt interface for variety of hot and cold zone set point temperatures specifically for the growth of mercury cadmium telluride (MCT). Reverse usage of hot and cold zones was simulated to aid the choice of proper orientation of crystal/melt interface regarding residual acceleration vector without actual change of furnace location on board the orbiter. It appears that an additional booster heater will be extremely helpful to ensure desired temperature gradient when hot and cold zones are reversed. Further efforts are required to investigate advantages/disadvantages of symmetrical furnace design (i.e. with similar length of hot and cold zones).

Stabilization of lower hybrid drift modes by finite parallel wavenumber and electron temperature gradients in field-reversed configurations

NASA Astrophysics Data System (ADS)

Farengo, R.; Guzdar, P. N.; Lee, Y. C.

1989-08-01

The effect of finite parallel wavenumber and electron temperature gradients on the lower hybrid drift instability is studied in the parameter regime corresponding to the TRX-2 device [Fusion Technol. 9, 48 (1986)]. Perturbations in the electrostatic potential and all three components of the vector potential are considered and finite beta electron orbit modifications are included. The electron temperature gradient decreases the growth rate of the instability but, for kz=0, unstable modes exist for ηe(=T'en0/Ten0)>6. Since finite kz effects completely stabilize the mode at small values of kz/ky(≂5×10-3), magnetic shear could be responsible for stabilizing the lower hybrid drift instability in field-reversed configurations.

Conjugate gradient based projection - A new explicit methodology for frictional contact

NASA Technical Reports Server (NTRS)

Tamma, Kumar K.; Li, Maocheng; Sha, Desong

1993-01-01

With special attention towards the applicability to parallel computation or vectorization, a new and effective explicit approach for linear complementary formulations involving a conjugate gradient based projection methodology is proposed in this study for contact problems with Coulomb friction. The overall objectives are focussed towards providing an explicit methodology of computation for the complete contact problem with friction. In this regard, the primary idea for solving the linear complementary formulations stems from an established search direction which is projected to a feasible region determined by the non-negative constraint condition; this direction is then applied to the Fletcher-Reeves conjugate gradient method resulting in a powerful explicit methodology which possesses high accuracy, excellent convergence characteristics, fast computational speed and is relatively simple to implement for contact problems involving Coulomb friction.

A fast pulse design for parallel excitation with gridding conjugate gradient.

PubMed

Feng, Shuo; Ji, Jim

2013-01-01

Parallel excitation (pTx) is recognized as a crucial technique in high field MRI to address the transmit field inhomogeneity problem. However, it can be time consuming to design pTx pulses which is not desirable. In this work, we propose a pulse design with gridding conjugate gradient (CG) based on the small-tip-angle approximation. The two major time consuming matrix-vector multiplications are substituted by two operators which involves with FFT and gridding only. Simulation results have shown that the proposed method is 3 times faster than conventional method and the memory cost is reduced by 1000 times.

Global Solutions to Repulsive Hookean Elastodynamics

NASA Astrophysics Data System (ADS)

Hu, Xianpeng; Masmoudi, Nader

2017-01-01

The global existence of classical solutions to the three dimensional repulsive Hookean elastodynamics around an equilibrium is considered. By linearization and Hodge's decomposition, the compressible part of the velocity, the density, and the compressible part of the transpose of the deformation gradient satisfy Klein-Gordon equations with speed {√{2}}, while the incompressible parts of the velocity and of the transpose of the deformation gradient satisfy wave equations with speed one. The space-time resonance method combined with the vector field method is used in a novel way to obtain the decay of the solution and hence global existence.

Aztec user`s guide. Version 1

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hutchinson, S.A.; Shadid, J.N.; Tuminaro, R.S.

1995-10-01

Aztec is an iterative library that greatly simplifies the parallelization process when solving the linear systems of equations Ax = b where A is a user supplied n x n sparse matrix, b is a user supplied vector of length n and x is a vector of length n to be computed. Aztec is intended as a software tool for users who want to avoid cumbersome parallel programming details but who have large sparse linear systems which require an efficiently utilized parallel processing system. A collection of data transformation tools are provided that allow for easy creation of distributed sparsemore » unstructured matrices for parallel solution. Once the distributed matrix is created, computation can be performed on any of the parallel machines running Aztec: nCUBE 2, IBM SP2 and Intel Paragon, MPI platforms as well as standard serial and vector platforms. Aztec includes a number of Krylov iterative methods such as conjugate gradient (CG), generalized minimum residual (GMRES) and stabilized biconjugate gradient (BICGSTAB) to solve systems of equations. These Krylov methods are used in conjunction with various preconditioners such as polynomial or domain decomposition methods using LU or incomplete LU factorizations within subdomains. Although the matrix A can be general, the package has been designed for matrices arising from the approximation of partial differential equations (PDEs). In particular, the Aztec package is oriented toward systems arising from PDE applications.« less

Iterative free-energy optimization for recurrent neural networks (INFERNO).

PubMed

Pitti, Alexandre; Gaussier, Philippe; Quoy, Mathias

2017-01-01

The intra-parietal lobe coupled with the Basal Ganglia forms a working memory that demonstrates strong planning capabilities for generating robust yet flexible neuronal sequences. Neurocomputational models however, often fails to control long range neural synchrony in recurrent spiking networks due to spontaneous activity. As a novel framework based on the free-energy principle, we propose to see the problem of spikes' synchrony as an optimization problem of the neurons sub-threshold activity for the generation of long neuronal chains. Using a stochastic gradient descent, a reinforcement signal (presumably dopaminergic) evaluates the quality of one input vector to move the recurrent neural network to a desired activity; depending on the error made, this input vector is strengthened to hill-climb the gradient or elicited to search for another solution. This vector can be learned then by one associative memory as a model of the basal-ganglia to control the recurrent neural network. Experiments on habit learning and on sequence retrieving demonstrate the capabilities of the dual system to generate very long and precise spatio-temporal sequences, above two hundred iterations. Its features are applied then to the sequential planning of arm movements. In line with neurobiological theories, we discuss its relevance for modeling the cortico-basal working memory to initiate flexible goal-directed neuronal chains of causation and its relation to novel architectures such as Deep Networks, Neural Turing Machines and the Free-Energy Principle.

Characterization of a Dynamic String Method for the Construction of Transition Pathways in Molecular Reactions

PubMed Central

Johnson, Margaret E.; Hummer, Gerhard

2012-01-01

We explore the theoretical foundation of different string methods used to find dominant reaction pathways in high-dimensional configuration spaces. Pathways are assessed by the amount of reactive flux they carry and by their orientation relative to the committor function. By examining the effects of transforming between different collective coordinates that span the same underlying space, we unmask artificial coordinate dependences in strings optimized to follow the free energy gradient. In contrast, strings optimized to follow the drift vector produce reaction pathways that are significantly less sensitive to reparameterizations of the collective coordinates. The differences in these paths arise because the drift vector depends on both the free energy gradient and the diffusion tensor of the coarse collective variables. Anisotropy and position dependence of diffusion tensors arise commonly in spaces of coarse variables, whose generally slow dynamics are obtained by nonlinear projections of the strongly coupled atomic motions. We show here that transition paths constructed to account for dynamics by following the drift vector will (to a close approximation) carry the maximum reactive flux both in systems with isotropic position dependent diffusion, and in systems with constant but anisotropic diffusion. We derive a simple method for calculating the committor function along paths that follow the reactive flux. Lastly, we provide guidance for the practical implementation of the dynamic string method. PMID:22616575

SAGRAD: A Program for Neural Network Training with Simulated Annealing and the Conjugate Gradient Method

PubMed Central

Bernal, Javier; Torres-Jimenez, Jose

2015-01-01

SAGRAD (Simulated Annealing GRADient), a Fortran 77 program for computing neural networks for classification using batch learning, is discussed. Neural network training in SAGRAD is based on a combination of simulated annealing and Møller’s scaled conjugate gradient algorithm, the latter a variation of the traditional conjugate gradient method, better suited for the nonquadratic nature of neural networks. Different aspects of the implementation of the training process in SAGRAD are discussed, such as the efficient computation of gradients and multiplication of vectors by Hessian matrices that are required by Møller’s algorithm; the (re)initialization of weights with simulated annealing required to (re)start Møller’s algorithm the first time and each time thereafter that it shows insufficient progress in reaching a possibly local minimum; and the use of simulated annealing when Møller’s algorithm, after possibly making considerable progress, becomes stuck at a local minimum or flat area of weight space. Outlines of the scaled conjugate gradient algorithm, the simulated annealing procedure and the training process used in SAGRAD are presented together with results from running SAGRAD on two examples of training data. PMID:26958442

Numerical study of comparison of vorticity and passive vectors in turbulence and inviscid flows

NASA Astrophysics Data System (ADS)

Ohkitani, Koji

2002-04-01

The nonlinear vortex stretching in incompressible Navier-Stokes turbulence is compared with a linear stretching process of passive vectors (PVs). In particular, we pay special attention to the difference of these processes under long and short time evolutions. For finite time evolution, we confirm our previous finding that the stretching effect of vorticity is weaker than that of general passive vectors for a majority of the initial conditions with the same energy spectra. The above difference can be explained qualitatively by examining the Biot-Savart formula. In order to see to what extent infinitesimal time development explains the above difference, we examine the probability density functions (PDFs) of the stretching rates of the passive vectors in the vicinity of a solution of Navier-Stokes equations. It is found that the PDFs are found to have a Gaussian distribution, suggesting that there are equally many PVs that stretched less and more than the vorticity. This suggests the importance of the vorticity-strain correlation built up over finite time in turbulence. We also discuss the case of Euler equations, where the dynamics of the Jacobian matrix relating the physical and material coordinates is examined numerically. A kind of alignment problem associated with the Cauchy-Green tensor is proposed and studied using the results of numerical simulations. It is found that vorticity tends to align itself with the most compressing eigenvector of the Cauchy-Green tensor. A two-dimensional counterpart of active-passive comparison is briefly studied. There is no essential difference between stretching of vorticity gradients and that of passive scalar gradients and a physical interpretation is given to it.

Effect of handling characteristics on minimum time cornering with torque vectoring

NASA Astrophysics Data System (ADS)

Smith, E. N.; Velenis, E.; Tavernini, D.; Cao, D.

2018-02-01

In this paper, the effect of both passive and actively-modified vehicle handling characteristics on minimum time manoeuvring for vehicles with 4-wheel torque vectoring (TV) capability is studied. First, a baseline optimal TV strategy is sought, independent of any causal control law. An optimal control problem (OCP) is initially formulated considering 4 independent wheel torque inputs, together with the steering angle rate, as the control variables. Using this formulation, the performance benefit using TV against an electric drive train with a fixed torque distribution, is demonstrated. The sensitivity of TV-controlled manoeuvre time to the passive understeer gradient of the vehicle is then studied. A second formulation of the OCP is introduced where a closed-loop TV controller is incorporated into the system dynamics of the OCP. This formulation allows the effect of actively modifying a vehicle's handling characteristic via TV on its minimum time cornering performance of the vehicle to be assessed. In particular, the effect of the target understeer gradient as the key tuning parameter of the literature-standard steady-state linear single-track model yaw rate reference is analysed.

Angle-domain common imaging gather extraction via Kirchhoff prestack depth migration based on a traveltime table in transversely isotropic media

NASA Astrophysics Data System (ADS)

Liu, Shaoyong; Gu, Hanming; Tang, Yongjie; Bingkai, Han; Wang, Huazhong; Liu, Dingjin

2018-04-01

Angle-domain common image-point gathers (ADCIGs) can alleviate the limitations of common image-point gathers in an offset domain, and have been widely used for velocity inversion and amplitude variation with angle (AVA) analysis. We propose an effective algorithm for generating ADCIGs in transversely isotropic (TI) media based on the gradient of traveltime by Kirchhoff pre-stack depth migration (KPSDM), as the dynamic programming method for computing the traveltime in TI media would not suffer from the limitation of shadow zones and traveltime interpolation. Meanwhile, we present a specific implementation strategy for ADCIG extraction via KPSDM. Three major steps are included in the presented strategy: (1) traveltime computation using a dynamic programming approach in TI media; (2) slowness vector calculation by the gradient of a traveltime table calculated previously; (3) construction of illumination vectors and subsurface angles in the migration process. Numerical examples are included to demonstrate the effectiveness of our approach, which henceforce shows its potential application for subsequent tomographic velocity inversion and AVA.

An efficient parallel algorithm for matrix-vector multiplication

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hendrickson, B.; Leland, R.; Plimpton, S.

The multiplication of a vector by a matrix is the kernel computation of many algorithms in scientific computation. A fast parallel algorithm for this calculation is therefore necessary if one is to make full use of the new generation of parallel supercomputers. This paper presents a high performance, parallel matrix-vector multiplication algorithm that is particularly well suited to hypercube multiprocessors. For an n x n matrix on p processors, the communication cost of this algorithm is O(n/[radical]p + log(p)), independent of the matrix sparsity pattern. The performance of the algorithm is demonstrated by employing it as the kernel in themore » well-known NAS conjugate gradient benchmark, where a run time of 6.09 seconds was observed. This is the best published performance on this benchmark achieved to date using a massively parallel supercomputer.« less

An efficient formulation and implementation of the analytic energy gradient method to the single and double excitation coupled-cluster wave function - Application to Cl2O2

NASA Technical Reports Server (NTRS)

Rendell, Alistair P.; Lee, Timothy J.

1991-01-01

The analytic energy gradient for the single and double excitation coupled-cluster (CCSD) wave function has been reformulated and implemented in a new set of programs. The reformulated set of gradient equations have a smaller computational cost than any previously published. The iterative solution of the linear equations and the construction of the effective density matrices are fully vectorized, being based on matrix multiplications. The new method has been used to investigate the Cl2O2 molecule, which has recently been postulated as an important intermediate in the destruction of ozone in the stratosphere. In addition to reporting computational timings, the CCSD equilibrium geometries, harmonic vibrational frequencies, infrared intensities, and relative energetics of three isomers of Cl2O2 are presented.

A Guided Tour of Mathematical Methods for the Physical Sciences

NASA Astrophysics Data System (ADS)

Snieder, Roel; van Wijk, Kasper

2015-05-01

1. Introduction; 2. Dimensional analysis; 3. Power series; 4. Spherical and cylindrical coordinates; 5. Gradient; 6. Divergence of a vector field; 7. Curl of a vector field; 8. Theorem of Gauss; 9. Theorem of Stokes; 10. The Laplacian; 11. Scale analysis; 12. Linear algebra; 13. Dirac delta function; 14. Fourier analysis; 15. Analytic functions; 16. Complex integration; 17. Green's functions: principles; 18. Green's functions: examples; 19. Normal modes; 20. Potential-field theory; 21. Probability and statistics; 22. Inverse problems; 23. Perturbation theory; 24. Asymptotic evaluation of integrals; 25. Conservation laws; 26. Cartesian tensors; 27. Variational calculus; 28. Epilogue on power and knowledge.

A parallel Jacobson-Oksman optimization algorithm. [parallel processing (computers)

NASA Technical Reports Server (NTRS)

Straeter, T. A.; Markos, A. T.

1975-01-01

A gradient-dependent optimization technique which exploits the vector-streaming or parallel-computing capabilities of some modern computers is presented. The algorithm, derived by assuming that the function to be minimized is homogeneous, is a modification of the Jacobson-Oksman serial minimization method. In addition to describing the algorithm, conditions insuring the convergence of the iterates of the algorithm and the results of numerical experiments on a group of sample test functions are presented. The results of these experiments indicate that this algorithm will solve optimization problems in less computing time than conventional serial methods on machines having vector-streaming or parallel-computing capabilities.

Genetic algorithm and graph theory based matrix factorization method for online friend recommendation.

PubMed

Li, Qu; Yao, Min; Yang, Jianhua; Xu, Ning

2014-01-01

Online friend recommendation is a fast developing topic in web mining. In this paper, we used SVD matrix factorization to model user and item feature vector and used stochastic gradient descent to amend parameter and improve accuracy. To tackle cold start problem and data sparsity, we used KNN model to influence user feature vector. At the same time, we used graph theory to partition communities with fairly low time and space complexity. What is more, matrix factorization can combine online and offline recommendation. Experiments showed that the hybrid recommendation algorithm is able to recommend online friends with good accuracy.

Removal of envelope protein-free retroviral vectors by anion-exchange chromatography to improve product quality.

PubMed

Rodrigues, Teresa; Alves, Ana; Lopes, António; Carrondo, Manuel J T; Alves, Paula M; Cruz, Pedro E

2008-10-01

We have investigated the role of the retroviral lipid bilayer and envelope proteins in the adsorption of retroviral vectors (RVs) to a Fractogel DEAE matrix. Intact RVs and their degradation components (envelope protein-free vectors and solubilized vector components) were adsorbed to this matrix and eluted using a linear gradient. Envelope protein-free RVs (Env(-)) and soluble envelope proteins (gp70) eluted in a significantly lower range of conductivities than intact RVs (Env(+)) (13.7-30 mS/cm for Env(-) and gp70 proteins vs. 47-80 mS/cm for Env(+)). The zeta (zeta)-potential of Env(+) and Env(-) vectors was evaluated showing that envelope proteins define the pI of the viral particles (pI (Env(+)) < 2 versus 3 < pI (Env(-)) < 4) and that Env(+) and Env(-) vectors have similar zeta-potentials within pH 5 and 8. The results presented herein indicate that the adsorption of retroviral particles occurs through multi-point interaction of the envelope proteins with the cationic groups on the chromatographic matrix. The strength of this adsorption is thus dependent on the amount of envelope protein present in the viral lipid bilayer. In conclusion, AEXc enables the separation of gp70 proteins as well as envelope protein-free vectors constituting a significant improvement to the quality of retroviral preparations for gene therapy applications.

Observable cosmological vector mode in the dark ages

NASA Astrophysics Data System (ADS)

Saga, Shohei

2016-09-01

The second-order vector mode is inevitably induced from the coupling of first-order scalar modes in cosmological perturbation theory and might hinder a possible detection of primordial gravitational waves from inflation through 21 cm lensing observations. Here, we investigate the weak lensing signal in 21 cm photons emitted by neutral hydrogen atoms in the dark ages induced by the second-order vector mode by decomposing the deflection angle of the 21 cm lensing signal into the gradient and curl modes. The curl mode is a good tracer of the cosmological vector and tensor modes since the scalar mode does not induce the curl one. By comparing angular power spectra of the 21 cm lensing curl mode induced by the second-order vector mode and primordial gravitational waves whose amplitude is parametrized by the tensor-to-scalar ratio r , we find that the 21 cm curl mode from the second-order vector mode dominates over that from primordial gravitational waves on almost all scales if r ≲10-5. If we use the multipoles of the power spectrum up to ℓmax=1 05 and 1 06 in reconstructing the curl mode from 21 cm temperature maps, the signal-to-noise ratios of the 21 cm curl mode from the second-order vector mode achieve S /N ≈0.46 and 73, respectively. Observation of 21 cm radiation is, in principle, a powerful tool to explore not only the tensor mode but also the cosmological vector mode.

Fast gravity, gravity partials, normalized gravity, gravity gradient torque and magnetic field: Derivation, code and data

NASA Technical Reports Server (NTRS)

Gottlieb, Robert G.

1993-01-01

Derivation of first and second partials of the gravitational potential is given in both normalized and unnormalized form. Two different recursion formulas are considered. Derivation of a general gravity gradient torque algorithm which uses the second partial of the gravitational potential is given. Derivation of the geomagnetic field vector is given in a form that closely mimics the gravitational algorithm. Ada code for all algorithms that precomputes all possible data is given. Test cases comparing the new algorithms with previous data are given, as well as speed comparisons showing the relative efficiencies of the new algorithms.

Modification of the magnetization dynamics of a NiFe nanodot due to thermal spin injection

NASA Astrophysics Data System (ADS)

Asam, Nagarjuna; Yamanoi, Kazuto; Kimura, Takashi

2018-06-01

An array of NiFe nanodots has been prepared on a Cu/CoFeAl film. Since a thermal spin current is expected to be excited owing to a large spin-dependent Seebeck coefficient for the CoFeAl, we investigate the magnetization dynamics of the NiFe dots under the temperature gradient along the vertical direction. By using vector network analyzer measurements, we have demonstrated that the temperature gradient produces modulations of the frequency of ferromagnetic resonance and the linewidth of the resonance spectra. The observed parabolic dependences are well explained by the damping-like and field-like components of spin transfer torque.

Cellular Mechanisms of Gravitropic Response in Higher Plants

NASA Astrophysics Data System (ADS)

Medvedev, Sergei; Smolikova, Galina; Pozhvanov, Gregory; Suslov, Dmitry

The evolutionary success of land plants in adaptation to the vectorial environmental factors was based mainly on the development of polarity systems. In result, normal plant ontogenesis is based on the positional information. Polarity is a tool by which the developing plant organs and tissues are mapped and the specific three-dimensional structure of the organism is created. It is due to their polar organization plants are able to orient themselves relative to the gravity vector and different vectorial cues, and to respond adequately to various stimuli. Gravitation is one of the most important polarized environmental factor that guides the development of plant organisms in space. Every plant can "estimate" its position relative to the gravity vector and correct it, if necessary, by means of polarized growth. The direction and the magnitude of gravitational stimulus are constant during the whole plant ontogenesis. The key plant response to the action of gravity is gravitropism, i.e. the directed growth of organs with respect to the gravity vector. This response is a very convenient model to study the mechanisms of plant orientation in space. The present report is focused on the main cellular mechanisms responsible for graviropic bending in higher plants. These mechanisms and structures include electric polarization of plant cells, Ca ({2+) }gradients, cytoskeleton, G-proteins, phosphoinositides and the machinery responsible for asymmetric auxin distribution. Those mechanisms tightly interact demonstrating some hierarchy and multiple feedbacks. The Ca (2+) gradients provide the primary physiological basis of polarity in plant cells. Calcium ions influence on the bioelectric potentials, the organization of actin cytoskeleton, the activity of Ca (2+) -binding proteins and Ca (2+) -dependent protein kinases. Protein kinases modulate transcription factors activity thereby regulating the gene expression and switching the developmental programs. Actin cytoskeleton affects the molecular machinery of polar auxin transport. It results in the changes of auxin gradients in plant organs and tissues, which modulate all cellular mechanisms of polarity via multiple feedback loops. The understanding of the mechanisms of plant organism orientation relative to the gravity vector will allow us to develop efficient technologies for plant growing in microgravity conditions at orbital space stations and during long piloted space flights. This work was supported by the grant of Russian Foundation for Basic Research (N 14-04-01-624) and by the grant of St.-Petersburg State University (N 1.38.233.2014).

[Factors for Degaussing of a Cochlear Implant Magnet in the MR Scanner].

PubMed

Koganezawa, Takumi; Uchiyama, Naoko; Teshigawara, Mai; Ogura, Akio

This study examined the conditions influencing degauss of the magnet using magnetic resonance imaging (MRI). Poly methyl methacrylate (PMMA) was used to fix the measurement magnets to the MRI bed at angles from 0° to 180° for the magnetic flux vector of static magnetic field. The PMMA was moved in the MRI magnetic field. Magnetic flux density was measured before and after bed movement, and the rate of degauss was calculated. The contents examined are as follows: (1) the angle of the magnetic flux vector of the measurement magnets for the magnetic flux vector of the static magnetic field, (2) the number of movements, (3) moving velocity, and (4) the movement on the spatial gradient of magnetic field. Mann-Whitney U test was used for statistical analysis of the data. In conclusion, the effect of the angle of the magnetic flux vector of the implant magnet was high under the conditions of degauss in this study. Therefore, during the MRI examination of a patient with a cochlear implant magnet, the operators identified the directions of the magnetic flux vector and static magnetic field of the implant magnet.

Fast software-based volume rendering using multimedia instructions on PC platforms and its application to virtual endoscopy

NASA Astrophysics Data System (ADS)

Mori, Kensaku; Suenaga, Yasuhito; Toriwaki, Jun-ichiro

2003-05-01

This paper describes a software-based fast volume rendering (VolR) method on a PC platform by using multimedia instructions, such as SIMD instructions, which are currently available in PCs' CPUs. This method achieves fast rendering speed through highly optimizing software rather than an improved rendering algorithm. In volume rendering using a ray casting method, the system requires fast execution of the following processes: (a) interpolation of voxel or color values at sample points, (b) computation of normal vectors (gray-level gradient vectors), (c) calculation of shaded values obtained by dot-products of normal vectors and light source direction vectors, (d) memory access to a huge area, and (e) efficient ray skipping at translucent regions. The proposed software implements these fundamental processes in volume rending by using special instruction sets for multimedia processing. The proposed software can generate virtual endoscopic images of a 3-D volume of 512x512x489 voxel size by volume rendering with perspective projection, specular reflection, and on-the-fly normal vector computation on a conventional PC without any special hardware at thirteen frames per second. Semi-translucent display is also possible.

Spherical Harmonic Inductive Detection Coils and their use In Dynamic Pre-emphasis for Magnetic Resonance Imaging

NASA Astrophysics Data System (ADS)

Edler, Karl T.

The issue of eddy currents induced by the rapid switching of magnetic field gradients is a long-standing problem in magnetic resonance imaging. A new method for dealing with this problem is presented whereby spatial harmonic components of the magnetic field are continuously sensed, through their temporal rates of change, and corrected. In this way, the effects of the eddy currents on multiple spatial harmonic components of the magnetic field can be detected and corrections applied during the rise time of the gradients. Sensing the temporal changes in each spatial harmonic is made possible with specially designed detection coils. However to make the design of these coils possible, general relationships between the spatial harmonics of the field, scalar potential, and vector potential are found within the quasi-static approximation. These relationships allow the vector potential to be found from the field -- an inverse curl operation -- and may be of use beyond the specific problem of detection coil design. Using the detection coils as sensors, methods are developed for designing a negative feedback system to control the eddy current effects and optimizing that system with respect to image noise and distortion. The design methods are successfully tested in a series of proof-of-principle experiments which lead to a discussion of how to incorporate similar designs into an operational MRI. Keywords: magnetic resonance imaging, eddy currents, dynamic shimming, negative feedback, quasi-static fields, vector potential, inverse curl

Research in computer science

NASA Technical Reports Server (NTRS)

Ortega, J. M.

1984-01-01

The research efforts of University of Virginia students under a NASA sponsored program are summarized and the status of the program is reported. The research includes: testing method evaluations for N version programming; a representation scheme for modeling three dimensional objects; fault tolerant protocols for real time local area networks; performance investigation of Cyber network; XFEM implementation; and vectorizing incomplete Cholesky conjugate gradients.

Wavelet Transforms in Parallel Image Processing

DTIC Science & Technology

1994-01-27

NUMBER OF PAGES Object Segmentation, Texture Segmentation, Image Compression, Image 137 Halftoning , Neural Network, Parallel Algorithms, 2D and 3D...Vector Quantization of Wavelet Transform Coefficients ........ ............................. 57 B.1.f Adaptive Image Halftoning based on Wavelet...application has been directed to the adaptive image halftoning . The gray information at a pixel, including its gray value and gradient, is represented by

Developing operation algorithms for vision subsystems in autonomous mobile robots

NASA Astrophysics Data System (ADS)

Shikhman, M. V.; Shidlovskiy, S. V.

2018-05-01

The paper analyzes algorithms for selecting keypoints on the image for the subsequent automatic detection of people and obstacles. The algorithm is based on the histogram of oriented gradients and the support vector method. The combination of these methods allows successful selection of dynamic and static objects. The algorithm can be applied in various autonomous mobile robots.

Toward an integrated view of ionospheric plasma instabilities: Altitudinal transitions and strong gradient case

NASA Astrophysics Data System (ADS)

Makarevich, Roman A.

2016-04-01

A general dispersion relation is derived that integrates the Farley-Buneman, gradient-drift, and current-convective plasma instabilities (FBI, GDI, and CCI) within the same formalism for an arbitrary altitude, wave propagation vector, and background density gradient. The limiting cases of the FBI/GDI in the E region for nearly field-aligned irregularities, GDI/CCI in the main F region at long wavelengths, and GDI at high altitudes are successfully recovered using analytic analysis. Numerical solutions are found for more general representative cases spanning the entire ionosphere. It is demonstrated that the results are consistent with those obtained using a general FBI/GDI/CCI theory developed previously at and near E region altitudes under most conditions. The most significant differences are obtained for strong gradients (scale lengths of 100 m) at high altitudes such as those that may occur during highly structured soft particle precipitation events. It is shown that the strong gradient case is dominated by inertial effects and, for some scales, surprisingly strong additional damping due to higher-order gradient terms. The growth rate behavior is examined with a particular focus on the range of wave propagations with positive growth (instability cone) and its transitions between altitudinal regions. It is shown that these transitions are largely controlled by the plasma density gradients even when FBI is operational.

Krylov subspace methods on supercomputers

NASA Technical Reports Server (NTRS)

Saad, Youcef

1988-01-01

A short survey of recent research on Krylov subspace methods with emphasis on implementation on vector and parallel computers is presented. Conjugate gradient methods have proven very useful on traditional scalar computers, and their popularity is likely to increase as three-dimensional models gain importance. A conservative approach to derive effective iterative techniques for supercomputers has been to find efficient parallel/vector implementations of the standard algorithms. The main source of difficulty in the incomplete factorization preconditionings is in the solution of the triangular systems at each step. A few approaches consisting of implementing efficient forward and backward triangular solutions are described in detail. Polynomial preconditioning as an alternative to standard incomplete factorization techniques is also discussed. Another efficient approach is to reorder the equations so as to improve the structure of the matrix to achieve better parallelism or vectorization. An overview of these and other ideas and their effectiveness or potential for different types of architectures is given.

PRECONDITIONED CONJUGATE-GRADIENT 2 (PCG2), a computer program for solving ground-water flow equations

USGS Publications Warehouse

Hill, Mary C.

1990-01-01

This report documents PCG2 : a numerical code to be used with the U.S. Geological Survey modular three-dimensional, finite-difference, ground-water flow model . PCG2 uses the preconditioned conjugate-gradient method to solve the equations produced by the model for hydraulic head. Linear or nonlinear flow conditions may be simulated. PCG2 includes two reconditioning options : modified incomplete Cholesky preconditioning, which is efficient on scalar computers; and polynomial preconditioning, which requires less computer storage and, with modifications that depend on the computer used, is most efficient on vector computers . Convergence of the solver is determined using both head-change and residual criteria. Nonlinear problems are solved using Picard iterations. This documentation provides a description of the preconditioned conjugate gradient method and the two preconditioners, detailed instructions for linking PCG2 to the modular model, sample data inputs, a brief description of PCG2, and a FORTRAN listing.

Four-body trajectory optimization

NASA Technical Reports Server (NTRS)

Pu, C. L.; Edelbaum, T. N.

1974-01-01

A comprehensive optimization program has been developed for computing fuel-optimal trajectories between the earth and a point in the sun-earth-moon system. It presents methods for generating fuel optimal two-impulse trajectories which may originate at the earth or a point in space and fuel optimal three-impulse trajectories between two points in space. The extrapolation of the state vector and the computation of the state transition matrix are accomplished by the Stumpff-Weiss method. The cost and constraint gradients are computed analytically in terms of the terminal state and the state transition matrix. The 4-body Lambert problem is solved by using the Newton-Raphson method. An accelerated gradient projection method is used to optimize a 2-impulse trajectory with terminal constraint. The Davidon's Variance Method is used both in the accelerated gradient projection method and the outer loop of a 3-impulse trajectory optimization problem.

Content-based unconstrained color logo and trademark retrieval with color edge gradient co-occurrence histograms

NASA Astrophysics Data System (ADS)

Phan, Raymond; Androutsos, Dimitrios

2008-01-01

In this paper, we present a logo and trademark retrieval system for unconstrained color image databases that extends the Color Edge Co-occurrence Histogram (CECH) object detection scheme. We introduce more accurate information to the CECH, by virtue of incorporating color edge detection using vector order statistics. This produces a more accurate representation of edges in color images, in comparison to the simple color pixel difference classification of edges as seen in the CECH. Our proposed method is thus reliant on edge gradient information, and as such, we call this the Color Edge Gradient Co-occurrence Histogram (CEGCH). We use this as the main mechanism for our unconstrained color logo and trademark retrieval scheme. Results illustrate that the proposed retrieval system retrieves logos and trademarks with good accuracy, and outperforms the CECH object detection scheme with higher precision and recall.

Association of Anthropogenic Land Use Change and Increased Abundance of the Chagas Disease Vector Rhodnius pallescens in a Rural Landscape of Panama

PubMed Central

Gottdenker, Nicole L.; Calzada, José E.; Saldaña, Azäel; Carroll, C. Ronald

2011-01-01

Anthropogenic disturbance is associated with increased vector-borne infectious disease transmission in wildlife, domestic animals, and humans. The objective of this study was to evaluate how disturbance of a tropical forest landscape impacts abundance of the triatomine bug Rhodnius pallescens, a vector of Chagas disease, in the region of the Panama Canal in Panama. Rhodnius pallescens was collected (n = 1,186) from its primary habitat, the palm Attalea butyracea, in five habitat types reflecting a gradient of anthropogenic disturbance. There was a high proportion of palms infested with R. pallescens across all habitat types (range = 77.1–91.4%). Results show that disturbed habitats are associated with increased vector abundance compared with relatively undisturbed habitats. Bugs collected in disturbed sites, although in higher abundance, tended to be in poor body condition compared with bugs captured in protected forest sites. Abundance data suggests that forest remnants may be sources for R. pallescens populations within highly disturbed areas of the landscape. PMID:21212205

Online Sequential Projection Vector Machine with Adaptive Data Mean Update

PubMed Central

Chen, Lin; Jia, Ji-Ting; Zhang, Qiong; Deng, Wan-Yu; Wei, Wei

2016-01-01

We propose a simple online learning algorithm especial for high-dimensional data. The algorithm is referred to as online sequential projection vector machine (OSPVM) which derives from projection vector machine and can learn from data in one-by-one or chunk-by-chunk mode. In OSPVM, data centering, dimension reduction, and neural network training are integrated seamlessly. In particular, the model parameters including (1) the projection vectors for dimension reduction, (2) the input weights, biases, and output weights, and (3) the number of hidden nodes can be updated simultaneously. Moreover, only one parameter, the number of hidden nodes, needs to be determined manually, and this makes it easy for use in real applications. Performance comparison was made on various high-dimensional classification problems for OSPVM against other fast online algorithms including budgeted stochastic gradient descent (BSGD) approach, adaptive multihyperplane machine (AMM), primal estimated subgradient solver (Pegasos), online sequential extreme learning machine (OSELM), and SVD + OSELM (feature selection based on SVD is performed before OSELM). The results obtained demonstrated the superior generalization performance and efficiency of the OSPVM. PMID:27143958

Online Sequential Projection Vector Machine with Adaptive Data Mean Update.

PubMed

Chen, Lin; Jia, Ji-Ting; Zhang, Qiong; Deng, Wan-Yu; Wei, Wei

2016-01-01

We propose a simple online learning algorithm especial for high-dimensional data. The algorithm is referred to as online sequential projection vector machine (OSPVM) which derives from projection vector machine and can learn from data in one-by-one or chunk-by-chunk mode. In OSPVM, data centering, dimension reduction, and neural network training are integrated seamlessly. In particular, the model parameters including (1) the projection vectors for dimension reduction, (2) the input weights, biases, and output weights, and (3) the number of hidden nodes can be updated simultaneously. Moreover, only one parameter, the number of hidden nodes, needs to be determined manually, and this makes it easy for use in real applications. Performance comparison was made on various high-dimensional classification problems for OSPVM against other fast online algorithms including budgeted stochastic gradient descent (BSGD) approach, adaptive multihyperplane machine (AMM), primal estimated subgradient solver (Pegasos), online sequential extreme learning machine (OSELM), and SVD + OSELM (feature selection based on SVD is performed before OSELM). The results obtained demonstrated the superior generalization performance and efficiency of the OSPVM.

Internal and external potential-field estimation from regional vector data at varying satellite altitude

NASA Astrophysics Data System (ADS)

Plattner, Alain; Simons, Frederik J.

2017-10-01

When modelling satellite data to recover a global planetary magnetic or gravitational potential field, the method of choice remains their analysis in terms of spherical harmonics. When only regional data are available, or when data quality varies strongly with geographic location, the inversion problem becomes severely ill-posed. In those cases, adopting explicitly local methods is to be preferred over adapting global ones (e.g. by regularization). Here, we develop the theory behind a procedure to invert for planetary potential fields from vector observations collected within a spatially bounded region at varying satellite altitude. Our method relies on the construction of spatiospectrally localized bases of functions that mitigate the noise amplification caused by downward continuation (from the satellite altitude to the source) while balancing the conflicting demands for spatial concentration and spectral limitation. The `altitude-cognizant' gradient vector Slepian functions (AC-GVSF) enjoy a noise tolerance under downward continuation that is much improved relative to the `classical' gradient vector Slepian functions (CL-GVSF), which do not factor satellite altitude into their construction. Furthermore, venturing beyond the realm of their first application, published in a preceding paper, in the present article we extend the theory to being able to handle both internal and external potential-field estimation. Solving simultaneously for internal and external fields under the limitation of regional data availability reduces internal-field artefacts introduced by downward-continuing unmodelled external fields, as we show with numerical examples. We explain our solution strategies on the basis of analytic expressions for the behaviour of the estimation bias and variance of models for which signal and noise are uncorrelated, (essentially) space- and band-limited, and spectrally (almost) white. The AC-GVSF are optimal linear combinations of vector spherical harmonics. Their construction is not altogether very computationally demanding when the concentration domains (the regions of spatial concentration) have circular symmetry, for example, on spherical caps or rings—even when the spherical-harmonic bandwidth is large. Data inversion proceeds by solving for the expansion coefficients of truncated function sequences, by least-squares analysis in a reduced-dimensional space. Hence, our method brings high-resolution regional potential-field modelling from incomplete and noisy vector-valued satellite data within reach of contemporary desktop machines.

Host Life History Strategy, Species Diversity, and Habitat Influence Trypanosoma cruzi Vector Infection in Changing Landscapes

PubMed Central

Gottdenker, Nicole L.; Chaves, Luis Fernando; Calzada, José E.; Saldaña, Azael; Carroll, C. Ronald

2012-01-01

Background Anthropogenic land use may influence transmission of multi-host vector-borne pathogens by changing diversity, relative abundance, and community composition of reservoir hosts. These reservoir hosts may have varying competence for vector-borne pathogens depending on species-specific characteristics, such as life history strategy. The objective of this study is to evaluate how anthropogenic land use change influences blood meal species composition and the effects of changing blood meal species composition on the parasite infection rate of the Chagas disease vector Rhodnius pallescens in Panama. Methodology/Principal Findings R. pallescens vectors (N = 643) were collected in different habitat types across a gradient of anthropogenic disturbance. Blood meal species in DNA extracted from these vectors was identified in 243 (40.3%) vectors by amplification and sequencing of a vertebrate-specific fragment of the 12SrRNA gene, and T. cruzi vector infection was determined by pcr. Vector infection rate was significantly greater in deforested habitats as compared to contiguous forests. Forty-two different species of blood meal were identified in R. pallescens, and species composition of blood meals varied across habitat types. Mammals (88.3%) dominated R. pallescens blood meals. Xenarthrans (sloths and tamanduas) were the most frequently identified species in blood meals across all habitat types. A regression tree analysis indicated that blood meal species diversity, host life history strategy (measured as rmax, the maximum intrinsic rate of population increase), and habitat type (forest fragments and peridomiciliary sites) were important determinants of vector infection with T. cruzi. The mean intrinsic rate of increase and the skewness and variability of rmax were positively associated with higher vector infection rate at a site. Conclusions/Significance In this study, anthropogenic landscape disturbance increased vector infection with T. cruzi, potentially by changing host community structure to favor hosts that are short-lived with high reproductive rates. Study results apply to potential environmental management strategies for Chagas disease. PMID:23166846

A vectorized algorithm for 3D dynamics of a tethered satellite

NASA Technical Reports Server (NTRS)

Wilson, Howard B.

1989-01-01

Equations of motion characterizing the three dimensional motion of a tethered satellite during the retrieval phase are studied. The mathematical model involves an arbitrary number of point masses connected by weightless cords. Motion occurs in a gravity gradient field. The formulation presented accounts for general functions describing support point motion, rate of tether retrieval, and arbitrary forces applied to the point masses. The matrix oriented program language MATLAB is used to produce an efficient vectorized formulation for computing natural frequencies and mode shapes for small oscillations about the static equilibrium configuration; and for integrating the nonlinear differential equations governing large amplitude motions. An example of time response pertaining to the skip rope effect is investigated.

Three axis vector atomic magnetometer utilizing polarimetric technique

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pradhan, Swarupananda, E-mail: spradhan@barc.gov.in, E-mail: pradhans75@gmail.com

2016-09-15

The three axis vector magnetic field measurement based on the interaction of a single elliptically polarized light beam with an atomic system is described. The magnetic field direction dependent atomic responses are extracted by the polarimetric detection in combination with laser frequency modulation and magnetic field modulation techniques. The magnetometer geometry offers additional critical requirements like compact size and large dynamic range for space application. Further, the three axis magnetic field is measured using only the reflected signal (one polarization component) from the polarimeter and thus can be easily expanded to make spatial array of detectors and/or high sensitivity fieldmore » gradient measurement as required for biomedical application.« less

Computational Approaches to Image Understanding.

DTIC Science & Technology

1981-10-01

represnting points, edges, surfaces, and volumes to facilitate display. The geometry or perspective and parailcl (or orthographic) projection has...of making the image forming process explicit. This in turn leads to a concern with geometry , such as the properties f the gradient, stereographic, and...dual spaces. Combining geometry and smoothness leads naturally to multi-variate vector analysis, and to differential geometry . For the most part, a

Nonlinear optimization with linear constraints using a projection method

NASA Technical Reports Server (NTRS)

Fox, T.

1982-01-01

Nonlinear optimization problems that are encountered in science and industry are examined. A method of projecting the gradient vector onto a set of linear contraints is developed, and a program that uses this method is presented. The algorithm that generates this projection matrix is based on the Gram-Schmidt method and overcomes some of the objections to the Rosen projection method.

Simulating Flaring Events via an Intelligent Cellular Automata Mechanism

NASA Astrophysics Data System (ADS)

Dimitropoulou, M.; Vlahos, L.; Isliker, H.; Georgoulis, M.

2010-07-01

We simulate flaring events through a Cellular Automaton (CA) model, in which, for the first time, we use observed vector magnetograms as initial conditions. After non-linear force free extrapolation of the magnetic field from the vector magnetograms, we identify magnetic discontinuities, using two alternative criteria: (1) the average magnetic field gradient, or (2) the normalized magnetic field curl (i.e. the current). Magnetic discontinuities are identified at the grid-sites where the magnetic field gradient or curl exceeds a specified threshold. We then relax the magnetic discontinuities according to the rules of Lu and Hamilton (1991) or Lu et al. (1993), i.e. we redistribute the magnetic field locally so that the discontinuities disappear. In order to simulate the flaring events, we consider several alternative scenarios with regard to: (1) The threshold above which magnetic discontinuities are identified (applying low, high, and height-dependent threshold values); (2) The driving process that occasionally causes new discontinuities (at randomly chosen grid sites, magnetic field increments are added that are perpendicular (or may-be also parallel) to the existing magnetic field). We address the question whether the coronal active region magnetic fields can indeed be considered to be in the state of self-organized criticality (SOC).

Delta-configurations - Flare activity and magnetic-field structure

NASA Technical Reports Server (NTRS)

Patty, S. R.; Hagyard, M. J.

1986-01-01

Complex sunspots in four active regions of April and May 1980, all exhibiting regions of magnetic classification delta, were studied using data from the NASA Marshall Space Flight Center vector magnetograph. The vector magnetic field structure in the vicinity of each delta was determined, and the location of the deltas in each active region was correlated with the locations and types of flare activity for the regions. Two types of delta-configuration were found to exist, active and inactive, as defined by the relationships between magnetic field structure and activity. The active delta exhibited high flare activity, strong horizontal gradients of the longitudinal (line-of-sight) magnetic field component, a strong transverse (perpendicular to line-of-sight) component, and a highly nonpotential orientation of the photospheric magnetic field, all indications of a highly sheared magnetic field. The inactive delta, on the other hand, exhibited little or no flare production, weaker horizontal gradients of the longitudinal component, weaker transverse components, and a nearly potential, nonsheared orientation of the magnetic field. It is concluded that the presence of such sheared fields is the primary signature by which the active delta may be distinguished, and that it is this shear which produces the flare activity of the active delta.

Analysis of nonlocal phonon thermal conductivity simulations showing the ballistic to diffusive crossover

NASA Astrophysics Data System (ADS)

Allen, Philip B.

2018-04-01

Simulations [e.g., X. W. Zhou et al., Phys. Rev. B 79, 115201 (2009), 10.1103/PhysRevB.79.115201] show nonlocal effects of the ballistic/diffusive crossover. The local temperature has nonlinear spatial variation not contained in the local Fourier law j ⃗(r ⃗) =-κ ∇ ⃗T (r ⃗) . The heat current j ⃗(r ⃗) depends not just on the local temperature gradient ∇ ⃗T (r ⃗) but also on temperatures at points r⃗' within phonon mean free paths, which can be micrometers long. This paper uses the Peierls-Boltzmann transport theory in nonlocal form to analyze the spatial variation Δ T (r ⃗) . The relaxation-time approximation (RTA) is used because the full solution is very challenging. Improved methods of extrapolation to obtain the bulk thermal conductivity κ are proposed. Callaway invented an approximate method of correcting RTA for the q ⃗ (phonon wave vector or crystal momentum) conservation of N (Normal as opposed to Umklapp) anharmonic collisions. This method is generalized to the nonlocal case where κ (k ⃗) depends on the wave vector of the current j ⃗(k ⃗) and temperature gradient i k ⃗Δ T (k ⃗) .

Effects of space weather on GOCE electrostatic gravity gradiometer measurements

NASA Astrophysics Data System (ADS)

Ince, E. Sinem; Pagiatakis, Spiros D.

2016-12-01

We examine the presence of residual nongravitational signatures in gravitational gradients measured by GOCE electrostatic gravity gradiometer. These signatures are observed over the magnetic poles during geomagnetically active days and can contaminate the trace of the gravitational gradient tensor by up to three to five times the expected noise level of the instrument (˜ 11 mE). We investigate these anomalies in the gradiometer measurements along many satellite tracks and examine possible causes using external datasets, such as interplanetary electric field measurements from the ACE (advanced composition explorer) and WIND spacecraft, and Poynting vector (flux) estimated from equivalent ionospheric currents derived from spherical elementary current systems over North America and Greenland. We show that the variations in the east-west and vertical electrical currents and Poynting vector components at the satellite position are highly correlated with the disturbances observed in the gradiometer measurements. The results presented in this paper reveal that the disturbances are due to intense ionospheric current variations that are enhanced by increased solar activity that causes a very dynamic drag environment. Moreover, successful modelling and removal of a high percentage of these disturbances are possible using external geomagnetic field observations.

Forms of null Lagrangians in field theories of continuum mechanics

NASA Astrophysics Data System (ADS)

Kovalev, V. A.; Radaev, Yu. N.

2012-02-01

The divergence representation of a null Lagrangian that is regular in a star-shaped domain is used to obtain its general expression containing field gradients of order ≤ 1 in the case of spacetime of arbitrary dimension. It is shown that for a static three-component field in the three-dimensional space, a null Lagrangian can contain up to 15 independent elements in total. The general form of a null Lagrangian in the four-dimensional Minkowski spacetime is obtained (the number of physical field variables is assumed arbitrary). A complete theory of the null Lagrangian for the n-dimensional spacetime manifold (including the four-dimensional Minkowski spacetime as a special case) is given. Null Lagrangians are then used as a basis for solving an important variational problem of an integrating factor. This problem involves searching for factors that depend on the spacetime variables, field variables, and their gradients and, for a given system of partial differential equations, ensure the equality between the scalar product of a vector multiplier by the system vector and some divergence expression for arbitrary field variables and, hence, allow one to formulate a divergence conservation law on solutions to the system.

Using Chebyshev polynomials and approximate inverse triangular factorizations for preconditioning the conjugate gradient method

NASA Astrophysics Data System (ADS)

Kaporin, I. E.

2012-02-01

In order to precondition a sparse symmetric positive definite matrix, its approximate inverse is examined, which is represented as the product of two sparse mutually adjoint triangular matrices. In this way, the solution of the corresponding system of linear algebraic equations (SLAE) by applying the preconditioned conjugate gradient method (CGM) is reduced to performing only elementary vector operations and calculating sparse matrix-vector products. A method for constructing the above preconditioner is described and analyzed. The triangular factor has a fixed sparsity pattern and is optimal in the sense that the preconditioned matrix has a minimum K-condition number. The use of polynomial preconditioning based on Chebyshev polynomials makes it possible to considerably reduce the amount of scalar product operations (at the cost of an insignificant increase in the total number of arithmetic operations). The possibility of an efficient massively parallel implementation of the resulting method for solving SLAEs is discussed. For a sequential version of this method, the results obtained by solving 56 test problems from the Florida sparse matrix collection (which are large-scale and ill-conditioned) are presented. These results show that the method is highly reliable and has low computational costs.

Data Images and Other Graphical Displays for Directional Data

NASA Technical Reports Server (NTRS)

Morphet, Bill; Symanzik, Juergen

2005-01-01

Vectors, axes, and periodic phenomena have direction. Directional variation can be expressed as points on a unit circle and is the subject of circular statistics, a relatively new application of statistics. An overview of existing methods for the display of directional data is given. The data image for linear variables is reviewed, then extended to directional variables by displaying direction using a color scale composed of a sequence of four or more color gradients with continuity between sequences and ordered intuitively in a color wheel such that the color of the 0deg angle is the same as the color of the 360deg angle. Cross over, which arose in automating the summarization of historical wind data, and color discontinuity resulting from the use a single color gradient in computational fluid dynamics visualization are eliminated. The new method provides for simultaneous resolution of detail on a small scale and overall structure on a large scale. Example circular data images are given of a global view of average wind direction of El Nino periods, computed rocket motor internal combustion flow, a global view of direction of the horizontal component of earth's main magnetic field on 9/15/2004, and Space Shuttle solid rocket motor nozzle vectoring.

Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large-Gap Magnetic Suspension System

NASA Technical Reports Server (NTRS)

Groom, Nelson J.

1997-01-01

The expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system are presented. The core is assumed to be uniformly magnetized, and equations are developed for two orientations of the magnetization vector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment. Fields at a given point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately defined by expansion up to second-order terms. Examination of the expanded equations for the case where the magnetization vector is perpendicular to the axis of symmetry reveals that some of the second-order gradient terms provide a method of generating torque about the axis of magnetization and therefore provide the ability to produce six-degree-of-freedom control.

Fine-scale features in the far-field of a turbulent jet

NASA Astrophysics Data System (ADS)

Buxton, Oliver; Ganapathisubramani, Bharathram

2008-11-01

The structure of a fully turbulent axisymmetric jet, at Reynolds number based on jet exit conditions of 5000, is investigated with cinematographic (1 kHz) stereoscopic PIV in a plane normal to the jet axis. Taylor's hypothesis is employed to calculate all three velocity gradients in the axial direction. The technique's resolution allows all terms of the velocity gradient tensor, hence strain rate tensor and kinetic energy dissipation, to be computed at each point within the plane. The data reveals that the vorticity field is dominated by high enstrophy tube-like structures. Conversely, the dissipation field appears to consist of sheet-like structures. Several criteria for isolating these strongly swirling vortical structures from the background turbulence were employed. One such technique involves isolating points in which the velocity gradient tensor has a real and a pair of complex conjugate eigenvectors. Once identified, the alignment of the various structures with relation to the vorticity vector and the real velocity gradient tensor eigenvector is investigated. The effect of the strain field on the geometry of the structures is also examined.

AUTOMATED CELL SEGMENTATION WITH 3D FLUORESCENCE MICROSCOPY IMAGES.

PubMed

Kong, Jun; Wang, Fusheng; Teodoro, George; Liang, Yanhui; Zhu, Yangyang; Tucker-Burden, Carol; Brat, Daniel J

2015-04-01

A large number of cell-oriented cancer investigations require an effective and reliable cell segmentation method on three dimensional (3D) fluorescence microscopic images for quantitative analysis of cell biological properties. In this paper, we present a fully automated cell segmentation method that can detect cells from 3D fluorescence microscopic images. Enlightened by fluorescence imaging techniques, we regulated the image gradient field by gradient vector flow (GVF) with interpolated and smoothed data volume, and grouped voxels based on gradient modes identified by tracking GVF field. Adaptive thresholding was then applied to voxels associated with the same gradient mode where voxel intensities were enhanced by a multiscale cell filter. We applied the method to a large volume of 3D fluorescence imaging data of human brain tumor cells with (1) small cell false detection and missing rates for individual cells; and (2) trivial over and under segmentation incidences for clustered cells. Additionally, the concordance of cell morphometry structure between automated and manual segmentation was encouraging. These results suggest a promising 3D cell segmentation method applicable to cancer studies.

Parallel Implementations of Gradient Based Iterative Algorithms for a Class of Discrete Optimal Control Problems.

DTIC Science & Technology

1987-02-28

g Ru + Ft~FzQZ...0 otherwise, and let X be the nN costate vector X =(Xt, X, , t,)t )X in E defined by X = Fz-tQz. Then, given u and z 0, the gradient g may be...34 ’ ’-" ’.’.’.’ ’ ". ."’.’,’,’ ./, ’ ’- ," ," .- ’- ". %.’’.’ "- ’" " ,,’’ " "- . "-" ’ ,r .’’"" ’ "." .r -".’ ". ." ." ." " . ." " . €"". ’.- -8- g Ru + FIX. (6) With the notation

Analytic energy gradients for the coupled-cluster singles and doubles method with the density-fitting approximation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bozkaya, Uğur, E-mail: ugur.bozkaya@hacettepe.edu.tr; Department of Chemistry, Atatürk University, Erzurum 25240; Sherrill, C. David

2016-05-07

An efficient implementation is presented for analytic gradients of the coupled-cluster singles and doubles (CCSD) method with the density-fitting approximation, denoted DF-CCSD. Frozen core terms are also included. When applied to a set of alkanes, the DF-CCSD analytic gradients are significantly accelerated compared to conventional CCSD for larger molecules. The efficiency of our DF-CCSD algorithm arises from the acceleration of several different terms, which are designated as the “gradient terms”: computation of particle density matrices (PDMs), generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the relaxed PDMs and GFM, back-transformation of PDMs and GFM to the atomic orbitalmore » (AO) basis, and evaluation of gradients in the AO basis. For the largest member of the alkane set (C{sub 10}H{sub 22}), the computational times for the gradient terms (with the cc-pVTZ basis set) are 2582.6 (CCSD) and 310.7 (DF-CCSD) min, respectively, a speed up of more than 8-folds. For gradient related terms, the DF approach avoids the usage of four-index electron repulsion integrals. Based on our previous study [U. Bozkaya, J. Chem. Phys. 141, 124108 (2014)], our formalism completely avoids construction or storage of the 4-index two-particle density matrix (TPDM), using instead 2- and 3-index TPDMs. The DF approach introduces negligible errors for equilibrium bond lengths and harmonic vibrational frequencies.« less

Leishmania amazonensis chemotaxis under glucose gradient studied by the strength and directionality of forces measured with optical tweezers

NASA Astrophysics Data System (ADS)

de Ysasa Pozzo, Liliana; Fontes, Adriana; de Thomaz, André A.; Barbosa, Luiz Carlos; Ayres, Diana Copi; Giorgio, Selma; Cesar, Carlos Lenz

2007-02-01

Chemotaxis is the mechanism microorganisms use to sense the environment surrounding them and to direct their movement towards attractive, or away from the repellent, chemicals. The biochemical sensing is almost the only way for communication between unicellular organisms. Prokaryote and Eukaryote chemotaxis has been mechanically studied mainly by observing the directionality and timing of the microorganisms movements subjected to a chemical gradient, but not through the directionality and strength of the forces it generates. To observe the vector force of microorganisms under a chemical gradient we developed a system composed of two large chambers connected by a tiny duct capable to keep the chemical gradient constant for more than ten hours. We also used the displacements of a microsphere trapped in an Optical Tweezers as the force transducer to measure the direction and the strength of the propulsion forces of flagellum of the microorganism under several gradient conditions. A 9μm diameter microsphere particle was trapped with a Nd:YAG laser and its movement was measured through the light scattered focused on a quadrant detector. We observed the behavior of the protozoa Leishmania amazonensis (eukaryote) under several glucose gradients. This protozoa senses the gradient around it by swimming in circles for three to five times following by tumbling, and not by the typical straight swimming/tumbling of bacteria. Our results also suggest that force direction and strength are also used to control its movement, not only the timing of swimming/tumbling, because we observed a higher force strength clearly directed towards the glucose gradient.

Global gyrokinetic simulations of intrinsic rotation in ASDEX Upgrade Ohmic L-mode plasmas

NASA Astrophysics Data System (ADS)

Hornsby, W. A.; Angioni, C.; Lu, Z. X.; Fable, E.; Erofeev, I.; McDermott, R.; Medvedeva, A.; Lebschy, A.; Peeters, A. G.; The ASDEX Upgrade Team

2018-05-01

Non-linear, radially global, turbulence simulations of ASDEX Upgrade (AUG) plasmas are performed and the nonlinear generated intrinsic flow shows agreement with the intrinsic flow gradients measured in the core of Ohmic L-mode plasmas at nominal parameters. Simulations utilising the kinetic electron model show hollow intrinsic flow profiles as seen in a predominant number of experiments performed at similar plasma parameters. In addition, significantly larger flow gradients are seen than in a previous flux-tube analysis (Hornsby et al 2017 Nucl. Fusion 57 046008). Adiabatic electron model simulations can show a flow profile with opposing sign in the gradient with respect to a kinetic electron simulation, implying a reversal in the sign of the residual stress due to kinetic electrons. The shaping of the intrinsic flow is strongly determined by the density gradient profile. The sensitivity of the residual stress to variations in density profile curvature is calculated and seen to be significantly stronger than to neoclassical flows (Hornsby et al 2017 Nucl. Fusion 57 046008). This variation is strong enough on its own to explain the large variations in the intrinsic flow gradients seen in some AUG experiments. Analysis of the symmetry breaking properties of the turbulence shows that profile shearing is the dominant mechanism in producing a finite parallel wave-number, with turbulence gradient effects contributing a smaller portion of the parallel wave-vector.

Novel technique for ST-T interval characterization in patients with acute myocardial ischemia.

PubMed

Correa, Raúl; Arini, Pedro David; Correa, Lorena Sabrina; Valentinuzzi, Max; Laciar, Eric

2014-07-01

The novel signal processing techniques have allowed and improved the use of vectorcardiography (VCG) to diagnose and characterize myocardial ischemia. Herein, we studied vectorcardiographic dynamic changes of ventricular repolarization in 80 patients before (control) and during Percutaneous Transluminal Coronary Angioplasty (PTCA). We propose four vectorcardiographic ST-T parameters, i.e., (a) ST Vector Magnitude Area (aSTVM); (b) T-wave Vector Magnitude Area (aTVM); (c) ST-T Vector Magnitude Difference (ST-TVD), and (d) T-wave Vector Magnitude Difference (TVD). For comparison, the conventional ST-Change Vector Magnitude (STCVM) and Spatial Ventricular Gradient (SVG) were also calculated. Our results indicate that several vectorcardiographic parameters show significant differences (p-value<0.05) before starting and during PTCA. Statistical minute-by-minute PTCA comparison against the control situation showed that ischemic monitoring reached a sensitivity=90.5% and a specificity=92.6% at the 5th minute of the PTCA, when aSTVM and ST-TVD were used as classifiers. We conclude that the sensitivity and specificity for acute ischemia monitoring could be increased with the use of only two vectorcardiographic parameters. Hence, the proposed technique based on vectorcardiography could be used in addition to the conventional ST-T analysis for better monitoring of ischemic patients. Copyright © 2014 Elsevier Ltd. All rights reserved.

Shadowgraph Study of Gradient Driven Fluctuations

NASA Technical Reports Server (NTRS)

Cannell, David; Nikolaenko, Gennady; Giglio, Marzio; Vailati, Alberto; Croccolo, Fabrizio; Meyer, William

2002-01-01

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Research in computer science

NASA Technical Reports Server (NTRS)

Ortega, J. M.

1986-01-01

Various graduate research activities in the field of computer science are reported. Among the topics discussed are: (1) failure probabilities in multi-version software; (2) Gaussian Elimination on parallel computers; (3) three dimensional Poisson solvers on parallel/vector computers; (4) automated task decomposition for multiple robot arms; (5) multi-color incomplete cholesky conjugate gradient methods on the Cyber 205; and (6) parallel implementation of iterative methods for solving linear equations.

Laser cooling of rubidium atoms in a 2D optical lattice

NASA Astrophysics Data System (ADS)

Wei, Chunhua; Kuhn, Carlos C. N.

2018-06-01

Lossless polarization gradient cooling of ?? atoms in a far-detuned 2D optical lattice is demonstrated. Temperatures down to ?K and phase space densities as high as 1 / 1000 are achieved in a total duty cycle of ?. It is shown that utilizing the vector component of the optical lattice allows lower temperatures to be achieved when compared with pure scalar lattices.

DensToolKit: A comprehensive open-source package for analyzing the electron density and its derivative scalar and vector fields

NASA Astrophysics Data System (ADS)

Solano-Altamirano, J. M.; Hernández-Pérez, Julio M.

2015-11-01

DensToolKit is a suite of cross-platform, optionally parallelized, programs for analyzing the molecular electron density (ρ) and several fields derived from it. Scalar and vector fields, such as the gradient of the electron density (∇ρ), electron localization function (ELF) and its gradient, localized orbital locator (LOL), region of slow electrons (RoSE), reduced density gradient, localized electrons detector (LED), information entropy, molecular electrostatic potential, kinetic energy densities K and G, among others, can be evaluated on zero, one, two, and three dimensional grids. The suite includes a program for searching critical points and bond paths of the electron density, under the framework of Quantum Theory of Atoms in Molecules. DensToolKit also evaluates the momentum space electron density on spatial grids, and the reduced density matrix of order one along lines joining two arbitrary atoms of a molecule. The source code is distributed under the GNU-GPLv3 license, and we release the code with the intent of establishing an open-source collaborative project. The style of DensToolKit's code follows some of the guidelines of an object-oriented program. This allows us to supply the user with a simple manner for easily implement new scalar or vector fields, provided they are derived from any of the fields already implemented in the code. In this paper, we present some of the most salient features of the programs contained in the suite, some examples of how to run them, and the mathematical definitions of the implemented fields along with hints of how we optimized their evaluation. We benchmarked our suite against both a freely-available program and a commercial package. Speed-ups of ˜2×, and up to 12× were obtained using a non-parallel compilation of DensToolKit for the evaluation of fields. DensToolKit takes similar times for finding critical points, compared to a commercial package. Finally, we present some perspectives for the future development and growth of the suite.

Can invertebrates see the e-vector of polarization as a separate modality of light?

PubMed

Labhart, Thomas

2016-12-15

The visual world is rich in linearly polarized light stimuli, which are hidden from the human eye. But many invertebrate species make use of polarized light as a source of valuable visual information. However, exploiting light polarization does not necessarily imply that the electric (e)-vector orientation of polarized light can be perceived as a separate modality of light. In this Review, I address the question of whether invertebrates can detect specific e-vector orientations in a manner similar to that of humans perceiving spectral stimuli as specific hues. To analyze e-vector orientation, the signals of at least three polarization-sensitive sensors (analyzer channels) with different e-vector tuning axes must be compared. The object-based, imaging polarization vision systems of cephalopods and crustaceans, as well as the water-surface detectors of flying backswimmers, use just two analyzer channels. Although this excludes the perception of specific e-vector orientations, a two-channel system does provide a coarse, categoric analysis of polarized light stimuli, comparable to the limited color sense of dichromatic, 'color-blind' humans. The celestial compass of insects employs three or more analyzer channels. However, that compass is multimodal, i.e. e-vector information merges with directional information from other celestial cues, such as the solar azimuth and the spectral gradient in the sky, masking e-vector information. It seems that invertebrate organisms take no interest in the polarization details of visual stimuli, but polarization vision grants more practical benefits, such as improved object detection and visual communication for cephalopods and crustaceans, compass readings to traveling insects, or the alert 'water below!' to water-seeking bugs. © 2016. Published by The Company of Biologists Ltd.

Can invertebrates see the e-vector of polarization as a separate modality of light?

PubMed Central

2016-01-01

ABSTRACT The visual world is rich in linearly polarized light stimuli, which are hidden from the human eye. But many invertebrate species make use of polarized light as a source of valuable visual information. However, exploiting light polarization does not necessarily imply that the electric (e)-vector orientation of polarized light can be perceived as a separate modality of light. In this Review, I address the question of whether invertebrates can detect specific e-vector orientations in a manner similar to that of humans perceiving spectral stimuli as specific hues. To analyze e-vector orientation, the signals of at least three polarization-sensitive sensors (analyzer channels) with different e-vector tuning axes must be compared. The object-based, imaging polarization vision systems of cephalopods and crustaceans, as well as the water-surface detectors of flying backswimmers, use just two analyzer channels. Although this excludes the perception of specific e-vector orientations, a two-channel system does provide a coarse, categoric analysis of polarized light stimuli, comparable to the limited color sense of dichromatic, ‘color-blind’ humans. The celestial compass of insects employs three or more analyzer channels. However, that compass is multimodal, i.e. e-vector information merges with directional information from other celestial cues, such as the solar azimuth and the spectral gradient in the sky, masking e-vector information. It seems that invertebrate organisms take no interest in the polarization details of visual stimuli, but polarization vision grants more practical benefits, such as improved object detection and visual communication for cephalopods and crustaceans, compass readings to traveling insects, or the alert ‘water below!’ to water-seeking bugs. PMID:27974532

Leishmania sand fly interaction: progress and challenges.

PubMed

Bates, Paul A

2008-08-01

Complex interactions occurs between Leishmania parasites and their sand fly vectors. Promastigotes of Leishmania live exclusively within the gut, possess flagella and are motile, and kinesins, kinases and G proteins have been described that play a role in regulating flagellar assembly. Movement within the gut is not random: promastigotes can detect gradients of solutes via chemotaxis and osmotaxis. Further they use their flagella to attach to the fly midgut using surface glyconconjugates, a key step in establishment of the infection. Differentiation of mammal-infective stages is characterised by significant biochemical and cellular remodelling. Further, the parasites can manipulate the behaviour of the vector to maximise their transmission, and flies may even deliver altruistic apoptotic forms to aid transmission of infective stages.

Application of Support Vector Machine to Forex Monitoring

NASA Astrophysics Data System (ADS)

Kamruzzaman, Joarder; Sarker, Ruhul A.

Previous studies have demonstrated superior performance of artificial neural network (ANN) based forex forecasting models over traditional regression models. This paper applies support vector machines to build a forecasting model from the historical data using six simple technical indicators and presents a comparison with an ANN based model trained by scaled conjugate gradient (SCG) learning algorithm. The models are evaluated and compared on the basis of five commonly used performance metrics that measure closeness of prediction as well as correctness in directional change. Forecasting results of six different currencies against Australian dollar reveal superior performance of SVM model using simple linear kernel over ANN-SCG model in terms of all the evaluation metrics. The effect of SVM parameter selection on prediction performance is also investigated and analyzed.

Alternative expression of the Bloch wave group velocity in loss-less periodic media using the electromagnetic field energy

NASA Astrophysics Data System (ADS)

Deparis, Olivier; Lambin, Philippe

2018-01-01

In periodic optical media, the group velocity is defined as the gradient with respect to wave-vector of the corresponding Bloch mode frequency dispersion curve, forming the photonic band structure. Instead of deducing it from the numerically computed photonic crystal band structure, the group velocity can be calculated directly from the integral of the Poynting vector over the crystal unit cell, the physical meaning of which is immediately perceivable. The related formula, which can be regarded as the application of Hellmann-Feynman theorem to electromagnetism, has been reported previously though without proof. We provide hereafter a full derivation of that formula starting from Maxwell's equations and we discuss its usefulness in photonics.

Fully automated detection of diabetic macular edema and dry age-related macular degeneration from optical coherence tomography images

PubMed Central

Srinivasan, Pratul P.; Kim, Leo A.; Mettu, Priyatham S.; Cousins, Scott W.; Comer, Grant M.; Izatt, Joseph A.; Farsiu, Sina

2014-01-01

We present a novel fully automated algorithm for the detection of retinal diseases via optical coherence tomography (OCT) imaging. Our algorithm utilizes multiscale histograms of oriented gradient descriptors as feature vectors of a support vector machine based classifier. The spectral domain OCT data sets used for cross-validation consisted of volumetric scans acquired from 45 subjects: 15 normal subjects, 15 patients with dry age-related macular degeneration (AMD), and 15 patients with diabetic macular edema (DME). Our classifier correctly identified 100% of cases with AMD, 100% cases with DME, and 86.67% cases of normal subjects. This algorithm is a potentially impactful tool for the remote diagnosis of ophthalmic diseases. PMID:25360373

Computation of output feedback gains for linear stochastic systems using the Zangwill-Powell method

NASA Technical Reports Server (NTRS)

Kaufman, H.

1977-01-01

Because conventional optimal linear regulator theory results in a controller which requires the capability of measuring and/or estimating the entire state vector, it is of interest to consider procedures for computing controls which are restricted to be linear feedback functions of a lower dimensional output vector and which take into account the presence of measurement noise and process uncertainty. To this effect a stochastic linear model has been developed that accounts for process parameter and initial uncertainty, measurement noise, and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was then performed for both finite and infinite time performance indices without gradient computation by using Zangwill's modification of a procedure originally proposed by Powell.

Soft learning vector quantization and clustering algorithms based on ordered weighted aggregation operators.

PubMed

Karayiannis, N B

2000-01-01

This paper presents the development and investigates the properties of ordered weighted learning vector quantization (LVQ) and clustering algorithms. These algorithms are developed by using gradient descent to minimize reformulation functions based on aggregation operators. An axiomatic approach provides conditions for selecting aggregation operators that lead to admissible reformulation functions. Minimization of admissible reformulation functions based on ordered weighted aggregation operators produces a family of soft LVQ and clustering algorithms, which includes fuzzy LVQ and clustering algorithms as special cases. The proposed LVQ and clustering algorithms are used to perform segmentation of magnetic resonance (MR) images of the brain. The diagnostic value of the segmented MR images provides the basis for evaluating a variety of ordered weighted LVQ and clustering algorithms.

Principal fiber bundle description of number scaling for scalars and vectors: application to gauge theory

NASA Astrophysics Data System (ADS)

Benioff, Paul

2015-05-01

The purpose of this paper is to put the description of number scaling and its effects on physics and geometry on a firmer foundation, and to make it more understandable. A main point is that two different concepts, number and number value are combined in the usual representations of number structures. This is valid as long as just one structure of each number type is being considered. It is not valid when different structures of each number type are being considered. Elements of base sets of number structures, considered by themselves, have no meaning. They acquire meaning or value as elements of a number structure. Fiber bundles over a space or space time manifold, M, are described. The fiber consists of a collection of many real or complex number structures and vector space structures. The structures are parameterized by a real or complex scaling factor, s. A vector space at a fiber level, s, has, as scalars, real or complex number structures at the same level. Connections are described that relate scalar and vector space structures at both neighbor M locations and at neighbor scaling levels. Scalar and vector structure valued fields are described and covariant derivatives of these fields are obtained. Two complex vector fields, each with one real and one imaginary field, appear, with one complex field associated with positions in M and the other with position dependent scaling factors. A derivation of the covariant derivative for scalar and vector valued fields gives the same vector fields. The derivation shows that the complex vector field associated with scaling fiber levels is the gradient of a complex scalar field. Use of these results in gauge theory shows that the imaginary part of the vector field associated with M positions acts like the electromagnetic field. The physical relevance of the other three fields, if any, is not known.

Implicit continuum mechanics approach to heat conduction in granular materials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Massoudi, M.; Mehrabadi, M.

In this paper, we derive a properly frame-invariant implicit constitutive relationship for the heat flux vector for a granular medium (or a density-gradient-type fluid). The heat flux vector is commonly modeled by Fourier’s law of heat conduction, and for complex materials such as nonlinear fluids, porous media, or granular materials, the coefficient of thermal conductivity is generalized by assuming that it would depend on a host of material and kinematic parameters such as temperature, shear rate, porosity, concentration, etc. In this paper, we extend the approach of Massoudi [Massoudi, M. Math. Methods Appl. Sci. 2006, 29, 1585; Massoudi, M. Math.more » Methods Appl. Sci. 2006, 29, 1599], who provided explicit constitutive relations for the heat flux vector for flowing granular materials; in order to do so, we use the implicit scheme suggested by Fox [Fox, N. Int. J. Eng. Sci. 1969, 7, 437], who obtained implicit relations in thermoelasticity.« less

Support Vector Machines Trained with Evolutionary Algorithms Employing Kernel Adatron for Large Scale Classification of Protein Structures.

PubMed

Arana-Daniel, Nancy; Gallegos, Alberto A; López-Franco, Carlos; Alanís, Alma Y; Morales, Jacob; López-Franco, Adriana

2016-01-01

With the increasing power of computers, the amount of data that can be processed in small periods of time has grown exponentially, as has the importance of classifying large-scale data efficiently. Support vector machines have shown good results classifying large amounts of high-dimensional data, such as data generated by protein structure prediction, spam recognition, medical diagnosis, optical character recognition and text classification, etc. Most state of the art approaches for large-scale learning use traditional optimization methods, such as quadratic programming or gradient descent, which makes the use of evolutionary algorithms for training support vector machines an area to be explored. The present paper proposes an approach that is simple to implement based on evolutionary algorithms and Kernel-Adatron for solving large-scale classification problems, focusing on protein structure prediction. The functional properties of proteins depend upon their three-dimensional structures. Knowing the structures of proteins is crucial for biology and can lead to improvements in areas such as medicine, agriculture and biofuels.

Real-Time Visual Tracking through Fusion Features

PubMed Central

Ruan, Yang; Wei, Zhenzhong

2016-01-01

Due to their high-speed, correlation filters for object tracking have begun to receive increasing attention. Traditional object trackers based on correlation filters typically use a single type of feature. In this paper, we attempt to integrate multiple feature types to improve the performance, and we propose a new DD-HOG fusion feature that consists of discriminative descriptors (DDs) and histograms of oriented gradients (HOG). However, fusion features as multi-vector descriptors cannot be directly used in prior correlation filters. To overcome this difficulty, we propose a multi-vector correlation filter (MVCF) that can directly convolve with a multi-vector descriptor to obtain a single-channel response that indicates the location of an object. Experiments on the CVPR2013 tracking benchmark with the evaluation of state-of-the-art trackers show the effectiveness and speed of the proposed method. Moreover, we show that our MVCF tracker, which uses the DD-HOG descriptor, outperforms the structure-preserving object tracker (SPOT) in multi-object tracking because of its high-speed and ability to address heavy occlusion. PMID:27347951

The Ozobranchus leech is a candidate mechanical vector for the fibropapilloma-associated turtle herpesvirus found latently infecting skin tumors on Hawaiian green turtles (Chelonia mydas)

USGS Publications Warehouse

Greenblatt, R.J.; Work, Thierry M.; Balazs, G.; Sutton, C.A.; Casey, R.N.; Casey, J.W.

2004-01-01

Fibropapillomatosis (FP) of marine turtles is a neoplastic disease of ecological concern. A fibropapilloma-associated turtle herpesvirus (FPTHV) is consistently present, usually at loads exceeding one virus copy per tumor cell. DNA from an array of parasites of green turtles (Chelonia mydas) was examined with quantitative PCR (qPCR) to determine whether any carried viral loads are sufficient to implicate them as vectors for FPTHV. Marine leeches (Ozobranchus spp.) were found to carry high viral DNA loads; some samples approached 10 million copies per leech. Isopycnic sucrose density gradient/qPCR analysis confirmed that some of these copies were associated with particles of the density of enveloped viruses. The data implicate the marine leech Ozobranchus as a mechanical vector for FPTHV. Quantitative RT-PCR analysis of FPTHV gene expression indicated that most of the FPTHV copies in a fibropapilloma have restricted DNA polymerase expression, suggestive of latent infection.

Analysis of the vector magnetic fields of complex sunspots

NASA Technical Reports Server (NTRS)

Patty, S. R.

1981-01-01

An analysis of the vector magnetic field in the delta-configurations of two complex sunspot groups is presented, noting several characteristics identified in the delta-configurations. The observations of regions 2469 (S12E80) and 2470 (S21E83) took place in May, 1980 with a vector magnetograph, verified by optical viewing. Longitudinal magnetic field plots located the delta-configurations in relation to the transverse field neutral line. It is shown that data on the polarization yields qualitative information on the magnetic field strengths, while the azimuth of the transverse field can be obtained from the relative intensities of linear polarization measurements aligned with respect to the magnetograph analyses axis at 0 and 90 deg, and at the plus and minus 45 deg positions. Details of the longitudinal fields are discussed. A strong, sheared transverse field component is found to be a signature of strong delta. A weak delta is accompanied by a weak longitudinal gradient with an unsheared transverse component of variable strength.

The use of computer-generated color graphic images for transient thermal analysis. [for hypersonic aircraft

NASA Technical Reports Server (NTRS)

Edwards, C. L. W.; Meissner, F. T.; Hall, J. B.

1979-01-01

Color computer graphics techniques were investigated as a means of rapidly scanning and interpreting large sets of transient heating data. The data presented were generated to support the conceptual design of a heat-sink thermal protection system (TPS) for a hypersonic research airplane. Color-coded vector and raster displays of the numerical geometry used in the heating calculations were employed to analyze skin thicknesses and surface temperatures of the heat-sink TPS under a variety of trajectory flight profiles. Both vector and raster displays proved to be effective means for rapidly identifying heat-sink mass concentrations, regions of high heating, and potentially adverse thermal gradients. The color-coded (raster) surface displays are a very efficient means for displaying surface-temperature and heating histories, and thereby the more stringent design requirements can quickly be identified. The related hardware and software developments required to implement both the vector and the raster displays for this application are also discussed.

Flies dynamically anti-track, rather than ballistically escape, aversive odor during flight.

PubMed

Wasserman, Sara; Lu, Patrick; Aptekar, Jacob W; Frye, Mark A

2012-08-15

Tracking distant odor sources is crucial to foraging, courtship and reproductive success for many animals including fish, flies and birds. Upon encountering a chemical plume in flight, Drosophila melanogaster integrates the spatial intensity gradient and temporal fluctuations over the two antennae, while simultaneously reducing the amplitude and frequency of rapid steering maneuvers, stabilizing the flight vector. There are infinite escape vectors away from a noxious source, in contrast to a single best tracking vector towards an attractive source. Attractive and aversive odors are segregated into parallel neuronal pathways in flies; therefore, the behavioral algorithms for avoidance may be categorically different from tracking. Do flies plot random ballistic or otherwise variable escape vectors? Or do they instead make use of temporally dynamic mechanisms for continuously and directly avoiding noxious odors in a manner similar to tracking appetitive ones? We examine this question using a magnetic tether flight simulator that permits free yaw movements, such that flies can actively orient within spatially defined odor plumes. We show that in-flight aversive flight behavior shares all of the key features of attraction such that flies continuously 'anti-track' the noxious source.

Flies dynamically anti-track, rather than ballistically escape, aversive odor during flight

PubMed Central

Wasserman, Sara; Lu, Patrick; Aptekar, Jacob W.; Frye, Mark A.

2012-01-01

SUMMARY Tracking distant odor sources is crucial to foraging, courtship and reproductive success for many animals including fish, flies and birds. Upon encountering a chemical plume in flight, Drosophila melanogaster integrates the spatial intensity gradient and temporal fluctuations over the two antennae, while simultaneously reducing the amplitude and frequency of rapid steering maneuvers, stabilizing the flight vector. There are infinite escape vectors away from a noxious source, in contrast to a single best tracking vector towards an attractive source. Attractive and aversive odors are segregated into parallel neuronal pathways in flies; therefore, the behavioral algorithms for avoidance may be categorically different from tracking. Do flies plot random ballistic or otherwise variable escape vectors? Or do they instead make use of temporally dynamic mechanisms for continuously and directly avoiding noxious odors in a manner similar to tracking appetitive ones? We examine this question using a magnetic tether flight simulator that permits free yaw movements, such that flies can actively orient within spatially defined odor plumes. We show that in-flight aversive flight behavior shares all of the key features of attraction such that flies continuously ‘anti-track’ the noxious source. PMID:22837456

Regional and seasonal response of a West Nile virus vector to climate change.

PubMed

Morin, Cory W; Comrie, Andrew C

2013-09-24

Climate change will affect the abundance and seasonality of West Nile virus (WNV) vectors, altering the risk of virus transmission to humans. Using downscaled general circulation model output, we calculate a WNV vector's response to climate change across the southern United States using process-based modeling. In the eastern United States, Culex quinquefasciatus response to projected climate change displays a latitudinal and elevational gradient. Projected summer population depressions as a result of increased immature mortality and habitat drying are most severe in the south and almost absent further north; extended spring and fall survival is ubiquitous. Much of California also exhibits a bimodal pattern. Projected onset of mosquito season is delayed in the southwestern United States because of extremely dry and hot spring and summers; however, increased temperature and late summer and fall rains extend the mosquito season. These results are unique in being a broad-scale calculation of the projected impacts of climate change on a WNV vector. The results show that, despite projected widespread future warming, the future seasonal response of C. quinquefasciatus populations across the southern United States will not be homogeneous, and will depend on specific combinations of local and regional conditions.

Vectorial mask optimization methods for robust optical lithography

NASA Astrophysics Data System (ADS)

Ma, Xu; Li, Yanqiu; Guo, Xuejia; Dong, Lisong; Arce, Gonzalo R.

2012-10-01

Continuous shrinkage of critical dimension in an integrated circuit impels the development of resolution enhancement techniques for low k1 lithography. Recently, several pixelated optical proximity correction (OPC) and phase-shifting mask (PSM) approaches were developed under scalar imaging models to account for the process variations. However, the lithography systems with larger-NA (NA>0.6) are predominant for current technology nodes, rendering the scalar models inadequate to describe the vector nature of the electromagnetic field that propagates through the optical lithography system. In addition, OPC and PSM algorithms based on scalar models can compensate for wavefront aberrations, but are incapable of mitigating polarization aberrations in practical lithography systems, which can only be dealt with under the vector model. To this end, we focus on developing robust pixelated gradient-based OPC and PSM optimization algorithms aimed at canceling defocus, dose variation, wavefront and polarization aberrations under a vector model. First, an integrative and analytic vector imaging model is applied to formulate the optimization problem, where the effects of process variations are explicitly incorporated in the optimization framework. A steepest descent algorithm is then used to iteratively optimize the mask patterns. Simulations show that the proposed algorithms can effectively improve the process windows of the optical lithography systems.

Joint Smoothed l₀-Norm DOA Estimation Algorithm for Multiple Measurement Vectors in MIMO Radar.

PubMed

Liu, Jing; Zhou, Weidong; Juwono, Filbert H

2017-05-08

Direction-of-arrival (DOA) estimation is usually confronted with a multiple measurement vector (MMV) case. In this paper, a novel fast sparse DOA estimation algorithm, named the joint smoothed l 0 -norm algorithm, is proposed for multiple measurement vectors in multiple-input multiple-output (MIMO) radar. To eliminate the white or colored Gaussian noises, the new method first obtains a low-complexity high-order cumulants based data matrix. Then, the proposed algorithm designs a joint smoothed function tailored for the MMV case, based on which joint smoothed l 0 -norm sparse representation framework is constructed. Finally, for the MMV-based joint smoothed function, the corresponding gradient-based sparse signal reconstruction is designed, thus the DOA estimation can be achieved. The proposed method is a fast sparse representation algorithm, which can solve the MMV problem and perform well for both white and colored Gaussian noises. The proposed joint algorithm is about two orders of magnitude faster than the l 1 -norm minimization based methods, such as l 1 -SVD (singular value decomposition), RV (real-valued) l 1 -SVD and RV l 1 -SRACV (sparse representation array covariance vectors), and achieves better DOA estimation performance.

Mass flow and velocity profiles in Neurospora hyphae: partial plug flow dominates intra-hyphal transport.

PubMed

Abadeh, Aryan; Lew, Roger R

2013-11-01

Movement of nuclei, mitochondria and vacuoles through hyphal trunks of Neurospora crassa were vector-mapped using fluorescent markers and green fluorescent protein tags. The vectorial movements of all three were strongly correlated, indicating the central role of mass (bulk) flow in cytoplasm movements in N. crassa. Profiles of velocity versus distance from the hyphal wall did not match the parabolic shape predicted by the ideal Hagen-Poiseuille model of flow at low Reynolds number. Instead, the profiles were flat, consistent with a model of partial plug flow due to the high concentration of organelles in the flowing cytosol. The intra-hyphal pressure gradients were manipulated by localized external osmotic treatments to demonstrate the dependence of velocity (and direction) on pressure gradients within the hyphae. The data support the concept that mass transport, driven by pressure gradients, dominates intra-hyphal transport. The transport occurs by partial plug flow due to the organelles in the cytosol.

Propagation properties of the chirped Airy beams through the gradient-index medium

NASA Astrophysics Data System (ADS)

Feng, Liyan; Zhang, Jianbin; Pang, Zihao; Wang, Linyi; Zhong, Tianfen; Yang, Xiangbo; Deng, Dongmei

2017-11-01

Through analytical derivation and numerical analysis, the propagation properties of the chirped Airy(CAi) beams in the gradient-index medium are investigated. The intensity and the phase distributions, the propagation trajectory and the Poynting vector of the CAi beams are demonstrated to investigate the propagation properties. Owing to the special and symmetrical refractive index profile of the gradient-index medium, the CAi beams propagate periodically. The effects of the distribution factor and the chirped parameter on the propagation of the CAi beams are analyzed. As the increasing of the distribution factor, the intensity distribution of the CAi beams is more scattering. However, with the chirped parameter increasing, the focusing property of the CAi beams strengthens. The variation of the chirped parameter can change the position of the peak intensity maximum, but it cannot alter the period of the peak intensity. The variations of the initial phase and the energy of the beams in the transverse plane expedite accordingly.

How Artificial Should the Treatment of a Plasma's Viscosity Be?

NASA Astrophysics Data System (ADS)

Whitney, K. G.; Velikovich, A. L.; Thornhill, J. W.; Davis, J.

1999-11-01

Electron viscosity dominates over ion viscosity and is important in describing the generation of shock fronts in highly ionizable plasmas. The sizes of shock front jumps in electron and ion temperature are determined from the magnitudes of the heat flow vector and pressure tensor, which, in turn, acquire non-negligible nonlinear contributions from the temperature and density gradients when these gradients are large. Thus, a consistent treatment of steep gradient formation in plasmas must come from investigations that include the effects of these nonlinear contributions to heat and momentum transport. Coefficients for each of five nonlinear contributions to the pressure tensor for an (r,z) Z-pinch geometry are presented and discussed in this talk. Hydrodynamic code calculations generally are not designed to provide a testbed for directly evaluating the kinetic energy dissipation that occurs at shock fronts; therefore, the strength of these nonlinear pressure tensor terms will be estimated by post-processing a Z-pinch hydrodynamics calculation and a steady-state planar shock wave calculation.

No-reference image quality assessment based on natural scene statistics and gradient magnitude similarity

NASA Astrophysics Data System (ADS)

Jia, Huizhen; Sun, Quansen; Ji, Zexuan; Wang, Tonghan; Chen, Qiang

2014-11-01

The goal of no-reference/blind image quality assessment (NR-IQA) is to devise a perceptual model that can accurately predict the quality of a distorted image as human opinions, in which feature extraction is an important issue. However, the features used in the state-of-the-art "general purpose" NR-IQA algorithms are usually natural scene statistics (NSS) based or are perceptually relevant; therefore, the performance of these models is limited. To further improve the performance of NR-IQA, we propose a general purpose NR-IQA algorithm which combines NSS-based features with perceptually relevant features. The new method extracts features in both the spatial and gradient domains. In the spatial domain, we extract the point-wise statistics for single pixel values which are characterized by a generalized Gaussian distribution model to form the underlying features. In the gradient domain, statistical features based on neighboring gradient magnitude similarity are extracted. Then a mapping is learned to predict quality scores using a support vector regression. The experimental results on the benchmark image databases demonstrate that the proposed algorithm correlates highly with human judgments of quality and leads to significant performance improvements over state-of-the-art methods.

A technique for accelerating the convergence of restarted GMRES

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baker, A H; Jessup, E R; Manteuffel, T

2004-03-09

We have observed that the residual vectors at the end of each restart cycle of restarted GMRES often alternate direction in a cyclic fashion, thereby slowing convergence. We present a new technique for accelerating the convergence of restarted GMRES by disrupting this alternating pattern. The new algorithm resembles a full conjugate gradient method with polynomial preconditioning, and its implementation requires minimal changes to the standard restarted GMRES algorithm.

Direct Volume Rendering with Shading via Three-Dimensional Textures

NASA Technical Reports Server (NTRS)

VanGelder, Allen; Kim, Kwansik

1996-01-01

A new and easy-to-implement method for direct volume rendering that uses 3D texture maps for acceleration, and incorporates directional lighting, is described. The implementation, called Voltx, produces high-quality images at nearly interactive speeds on workstations with hardware support for three-dimensional texture maps. Previously reported methods did not incorporate a light model, and did not address issues of multiple texture maps for large volumes. Our research shows that these extensions impact performance by about a factor of ten. Voltx supports orthographic, perspective, and stereo views. This paper describes the theory and implementation of this technique, and compares it to the shear-warp factorization approach. A rectilinear data set is converted into a three-dimensional texture map containing color and opacity information. Quantized normal vectors and a lookup table provide efficiency. A new tesselation of the sphere is described, which serves as the basis for normal-vector quantization. A new gradient-based shading criterion is described, in which the gradient magnitude is interpreted in the context of the field-data value and the material classification parameters, and not in isolation. In the rendering phase, the texture map is applied to a stack of parallel planes, which effectively cut the texture into many slabs. The slabs are composited to form an image.

Correlation between solar flare productivity and photospheric vector magnetic fields

NASA Astrophysics Data System (ADS)

Cui, Yanmei; Wang, Huaning

2008-11-01

Studying the statistical correlation between the solar flare productivity and photospheric magnetic fields is very important and necessary. It is helpful to set up a practical flare forecast model based on magnetic properties and improve the physical understanding of solar flare eruptions. In the previous study ([Cui, Y.M., Li, R., Zhang, L.Y., He, Y.L., Wang, H.N. Correlation between solar flare productivity and photospheric magnetic field properties 1. Maximum horizontal gradient, length of neutral line, number of singular points. Sol. Phys. 237, 45 59, 2006]; from now on we refer to this paper as ‘Paper I’), three measures of the maximum horizontal gradient, the length of the neutral line, and the number of singular points are computed from 23990 SOHO/MDI longitudinal magnetograms. The statistical relationship between the solar flare productivity and these three measures is well fitted with sigmoid functions. In the current work, the three measures of the length of strong-shear neutral line, total unsigned current, and total unsigned current helicity are computed from 1353 vector magnetograms observed at Huairou Solar Observing Station. The relationship between the solar flare productivity and the current three measures can also be well fitted with sigmoid functions. These results are expected to be beneficial to future operational flare forecasting models.

Infrared modified gravity with propagating torsion: Instability of torsionfull de Sitter-like solutions

NASA Astrophysics Data System (ADS)

Nikiforova, Vasilisa; Damour, Thibault

2018-06-01

We continue the exploration of the consistency of a modified-gravity theory that generalizes general relativity by including a dynamical torsion in addition to the dynamical metric. The six-parameter theory we consider was found to be consistent around arbitrary torsionless Einstein backgrounds, in spite of its containing a (notoriously delicate) massive spin-2 excitation. At zero bare cosmological constant, this theory was found to admit a self-accelerating solution whose exponential expansion is sustained by a nonzero torsion background. The scalar-type perturbations of the latter torsionfull self-accelerating solution were recently studied and were found to preserve the number of propagating scalar degrees of freedom, but to exhibit, for some values of the torsion background, some exponential instabilities (of a rather mild type). Here, we study the tensor-type and vector-type perturbations of the torsionfull self-accelerating solution, and of its deformation by a nonzero bare cosmological constant. We find strong, "gradient" instabilities in the vector sector. No tuning of the parameters of the theory can kill these instabilities without creating instabilities in the other sectors. Further work is needed to see whether generic torsionfull backgrounds are prone to containing gradient instabilities, or if the instabilities we found are mainly due to the (generalized) self-accelerating nature of the special de Sitter backgrounds we considered.

Evaluating the statistical performance of less applied algorithms in classification of worldview-3 imagery data in an urbanized landscape

NASA Astrophysics Data System (ADS)

Ranaie, Mehrdad; Soffianian, Alireza; Pourmanafi, Saeid; Mirghaffari, Noorollah; Tarkesh, Mostafa

2018-03-01

In recent decade, analyzing the remotely sensed imagery is considered as one of the most common and widely used procedures in the environmental studies. In this case, supervised image classification techniques play a central role. Hence, taking a high resolution Worldview-3 over a mixed urbanized landscape in Iran, three less applied image classification methods including Bagged CART, Stochastic gradient boosting model and Neural network with feature extraction were tested and compared with two prevalent methods: random forest and support vector machine with linear kernel. To do so, each method was run ten time and three validation techniques was used to estimate the accuracy statistics consist of cross validation, independent validation and validation with total of train data. Moreover, using ANOVA and Tukey test, statistical difference significance between the classification methods was significantly surveyed. In general, the results showed that random forest with marginal difference compared to Bagged CART and stochastic gradient boosting model is the best performing method whilst based on independent validation there was no significant difference between the performances of classification methods. It should be finally noted that neural network with feature extraction and linear support vector machine had better processing speed than other.

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

NASA Technical Reports Server (NTRS)

Bune, Andris V.; Sen, Subhayu; Mukherjee, Sundeep; Catalina, Adrian; Stefanescu, Doru M.

1999-01-01

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Tire-breeding mosquitoes of public health importance along an urbanisation gradient in Buenos Aires, Argentina.

PubMed

Rubio, Alejandra; Cardo, María Victoria; Vezzani, Darío

2011-09-01

Used vehicle tires are a source of mosquito vectors and a means of their introduction and expansion. With the aim of assessing the effects of urbanisation on the main mosquito vectors in temperate Argentina, the infestation levels of Aedes aegypti (L.) and Culex pipiens L. were studied in used tires from highly urbanised cities to low-urbanised small towns in Buenos Aires. Immatures of both species accounted for 96% of the 9,722 individuals collected; the total individuals collected represented seven species. The percentage of water-filled tires containing mosquitoes [container index (CI)] was 33% and the percentage of infested sites [site index (SI)] was 65.2%. These indexes decreased significantly from low to high urbanisation levels for both mosquito species. The relative abundance (RA) of Ae. aegypti immatures was slightly higher toward large cities, but showed no difference for Cx. pipiens. The CI of shaded tires was significantly higher than the CI of exposed tires for both mosquito species. There was no difference in RA values between shaded and sunlit tires. The CI and the SI were highest during the summer across the urbanisation levels, except for Cx. pipiens, which continued to increase during the autumn in small towns. Results related to urbanisation gradient, sunlit exposure and seasonality are discussed.

Automatic lesion boundary detection in dermoscopy images using gradient vector flow snakes

PubMed Central

Erkol, Bulent; Moss, Randy H.; Stanley, R. Joe; Stoecker, William V.; Hvatum, Erik

2011-01-01

Background Malignant melanoma has a good prognosis if treated early. Dermoscopy images of pigmented lesions are most commonly taken at × 10 magnification under lighting at a low angle of incidence while the skin is immersed in oil under a glass plate. Accurate skin lesion segmentation from the background skin is important because some of the features anticipated to be used for diagnosis deal with shape of the lesion and others deal with the color of the lesion compared with the color of the surrounding skin. Methods In this research, gradient vector flow (GVF) snakes are investigated to find the border of skin lesions in dermoscopy images. An automatic initialization method is introduced to make the skin lesion border determination process fully automated. Results Skin lesion segmentation results are presented for 70 benign and 30 melanoma skin lesion images for the GVF-based method and a color histogram analysis technique. The average errors obtained by the GVF-based method are lower for both the benign and melanoma image sets than for the color histogram analysis technique based on comparison with manually segmented lesions determined by a dermatologist. Conclusions The experimental results for the GVF-based method demonstrate promise as an automated technique for skin lesion segmentation in dermoscopy images. PMID:15691255

15N backbone dynamics of the S-peptide from ribonuclease A in its free and S-protein bound forms: toward a site-specific analysis of entropy changes upon folding.

PubMed Central

Alexandrescu, A. T.; Rathgeb-Szabo, K.; Rumpel, K.; Jahnke, W.; Schulthess, T.; Kammerer, R. A.

1998-01-01

Backbone 15N relaxation parameters (R1, R2, 1H-15N NOE) have been measured for a 22-residue recombinant variant of the S-peptide in its free and S-protein bound forms. NMR relaxation data were analyzed using the "model-free" approach (Lipari & Szabo, 1982). Order parameters obtained from "model-free" simulations were used to calculate 1H-15N bond vector entropies using a recently described method (Yang & Kay, 1996), in which the form of the probability density function for bond vector fluctuations is derived from a diffusion-in-a-cone motional model. The average change in 1H-15N bond vector entropies for residues T3-S15, which become ordered upon binding of the S-peptide to the S-protein, is -12.6+/-1.4 J/mol.residue.K. 15N relaxation data suggest a gradient of decreasing entropy values moving from the termini toward the center of the free peptide. The difference between the entropies of the terminal and central residues is about -12 J/mol residue K, a value comparable to that of the average entropy change per residue upon complex formation. Similar entropy gradients are evident in NMR relaxation studies of other denatured proteins. Taken together, these observations suggest denatured proteins may contain entropic contributions from non-local interactions. Consequently, calculations that model the entropy of a residue in a denatured protein as that of a residue in a di- or tri-peptide, might over-estimate the magnitude of entropy changes upon folding. PMID:9521116

Computer-aided diagnosis of lung nodule using gradient tree boosting and Bayesian optimization.

PubMed

Nishio, Mizuho; Nishizawa, Mitsuo; Sugiyama, Osamu; Kojima, Ryosuke; Yakami, Masahiro; Kuroda, Tomohiro; Togashi, Kaori

2018-01-01

We aimed to evaluate a computer-aided diagnosis (CADx) system for lung nodule classification focussing on (i) usefulness of the conventional CADx system (hand-crafted imaging feature + machine learning algorithm), (ii) comparison between support vector machine (SVM) and gradient tree boosting (XGBoost) as machine learning algorithms, and (iii) effectiveness of parameter optimization using Bayesian optimization and random search. Data on 99 lung nodules (62 lung cancers and 37 benign lung nodules) were included from public databases of CT images. A variant of the local binary pattern was used for calculating a feature vector. SVM or XGBoost was trained using the feature vector and its corresponding label. Tree Parzen Estimator (TPE) was used as Bayesian optimization for parameters of SVM and XGBoost. Random search was done for comparison with TPE. Leave-one-out cross-validation was used for optimizing and evaluating the performance of our CADx system. Performance was evaluated using area under the curve (AUC) of receiver operating characteristic analysis. AUC was calculated 10 times, and its average was obtained. The best averaged AUC of SVM and XGBoost was 0.850 and 0.896, respectively; both were obtained using TPE. XGBoost was generally superior to SVM. Optimal parameters for achieving high AUC were obtained with fewer numbers of trials when using TPE, compared with random search. Bayesian optimization of SVM and XGBoost parameters was more efficient than random search. Based on observer study, AUC values of two board-certified radiologists were 0.898 and 0.822. The results show that diagnostic accuracy of our CADx system was comparable to that of radiologists with respect to classifying lung nodules.

Integration of celestial compass cues in the central complex of the locust brain.

PubMed

Pegel, Uta; Pfeiffer, Keram; Homberg, Uwe

2018-01-29

Many insects rely on celestial compass cues such as the polarization pattern of the sky for spatial orientation. In the desert locust, the central complex (CX) houses multiple sets of neurons, sensitive to the oscillation plane of polarized light and thus probably acts as an internal polarization compass. We investigated whether other sky compass cues like direct sunlight or the chromatic gradient of the sky might contribute to this compass. We recorded from polarization-sensitive CX neurons while an unpolarized green or ultraviolet light spot was moved around the head of the animal. All types of neuron that were sensitive to the plane of polarization ( E -vector) above the animal also responded to the unpolarized light spots in an azimuth-dependent way. The tuning to the unpolarized light spots was independent of wavelength, suggesting that the neurons encode solar azimuth based on direct sunlight and not on the sky chromatic gradient. Two cell types represented the natural 90 deg relationship between solar azimuth and zenithal E -vector orientation, providing evidence to suggest that solar azimuth information supports the internal polarization compass. Most neurons showed advances in their tuning to the E -vector and the unpolarized light spots dependent on rotation direction, consistent with anticipatory signaling. The amplitude of responses and its variability were dependent on the level of background firing, possibly indicating different internal states. The integration of polarization and solar azimuth information strongly suggests that besides the polarization pattern of the sky, direct sunlight might be an important cue for sky compass navigation in the locust. © 2018. Published by The Company of Biologists Ltd.

On efficient randomized algorithms for finding the PageRank vector

NASA Astrophysics Data System (ADS)

Gasnikov, A. V.; Dmitriev, D. Yu.

2015-03-01

Two randomized methods are considered for finding the PageRank vector; in other words, the solution of the system p T = p T P with a stochastic n × n matrix P, where n ˜ 107-109, is sought (in the class of probability distributions) with accuracy ɛ: ɛ ≫ n -1. Thus, the possibility of brute-force multiplication of P by the column is ruled out in the case of dense objects. The first method is based on the idea of Markov chain Monte Carlo algorithms. This approach is efficient when the iterative process p {/t+1 T} = p {/t T} P quickly reaches a steady state. Additionally, it takes into account another specific feature of P, namely, the nonzero off-diagonal elements of P are equal in rows (this property is used to organize a random walk over the graph with the matrix P). Based on modern concentration-of-measure inequalities, new bounds for the running time of this method are presented that take into account the specific features of P. In the second method, the search for a ranking vector is reduced to finding the equilibrium in the antagonistic matrix game where S n (1) is a unit simplex in ℝ n and I is the identity matrix. The arising problem is solved by applying a slightly modified Grigoriadis-Khachiyan algorithm (1995). This technique, like the Nazin-Polyak method (2009), is a randomized version of Nemirovski's mirror descent method. The difference is that randomization in the Grigoriadis-Khachiyan algorithm is used when the gradient is projected onto the simplex rather than when the stochastic gradient is computed. For sparse matrices P, the method proposed yields noticeably better results.

Monthly evaporation forecasting using artificial neural networks and support vector machines

NASA Astrophysics Data System (ADS)

Tezel, Gulay; Buyukyildiz, Meral

2016-04-01

Evaporation is one of the most important components of the hydrological cycle, but is relatively difficult to estimate, due to its complexity, as it can be influenced by numerous factors. Estimation of evaporation is important for the design of reservoirs, especially in arid and semi-arid areas. Artificial neural network methods and support vector machines (SVM) are frequently utilized to estimate evaporation and other hydrological variables. In this study, usability of artificial neural networks (ANNs) (multilayer perceptron (MLP) and radial basis function network (RBFN)) and ɛ-support vector regression (SVR) artificial intelligence methods was investigated to estimate monthly pan evaporation. For this aim, temperature, relative humidity, wind speed, and precipitation data for the period 1972 to 2005 from Beysehir meteorology station were used as input variables while pan evaporation values were used as output. The Romanenko and Meyer method was also considered for the comparison. The results were compared with observed class A pan evaporation data. In MLP method, four different training algorithms, gradient descent with momentum and adaptive learning rule backpropagation (GDX), Levenberg-Marquardt (LVM), scaled conjugate gradient (SCG), and resilient backpropagation (RBP), were used. Also, ɛ-SVR model was used as SVR model. The models were designed via 10-fold cross-validation (CV); algorithm performance was assessed via mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R 2). According to the performance criteria, the ANN algorithms and ɛ-SVR had similar results. The ANNs and ɛ-SVR methods were found to perform better than the Romanenko and Meyer methods. Consequently, the best performance using the test data was obtained using SCG(4,2,2,1) with R 2 = 0.905.

Production and purification of non replicative canine adenovirus type 2 derived vectors.

PubMed

Szelechowski, Marion; Bergeron, Corinne; Gonzalez-Dunia, Daniel; Klonjkowski, Bernard

2013-12-03

Adenovirus (Ad) derived vectors have been widely used for short or long-term gene transfer, both for gene therapy and vaccine applications. Because of the frequent pre-existing immunity against the classically used human adenovirus type 5, canine adenovirus type 2 (CAV2) has been proposed as an alternative vector for human gene transfer. The well-characterized biology of CAV2, together with its ease of genetic manipulation, offer major advantages, notably for gene transfer into the central nervous system, or for inducing a wide range of protective immune responses, from humoral to cellular immunity. Nowadays, CAV2 represents one of the most appealing nonhuman adenovirus for use as a vaccine vector. This protocol describes a simple method to construct, produce and titer recombinant CAV2 vectors. After cloning the expression cassette of the gene of interest into a shuttle plasmid, the recombinant genomic plasmid is obtained by homologous recombination in the E. coli BJ5183 bacterial strain. The resulting genomic plasmid is then transfected into canine kidney cells expressing the complementing CAV2-E1 genes (DK-E1). A viral amplification enables the production of a large viral stock, which is purified by ultracentrifugation through cesium chloride gradients and desalted by dialysis. The resulting viral suspension routinely has a titer of over 10(10) infectious particles per ml and can be directly administrated in vivo.

Hybrid simulations of radial transport driven by the Rayleigh-Taylor instability

NASA Astrophysics Data System (ADS)

Delamere, P. A.; Stauffer, B. H.; Ma, X.

2017-12-01

Plasma transport in the rapidly rotating giant magnetospheres is thought to involve a centrifugally-driven flux tube interchange instability, similar to the Rayleigh-Taylor (RT) instability. In three dimensions, the convective flow patterns associated with the RT instability can produce strong guide field reconnection, allowing plasma mass to move radially outward while conserving magnetic flux (Ma et al., 2016). We present a set of hybrid (kinetic ion / fluid electron) plasma simulations of the RT instability using high plasma beta conditions appropriate for Jupiter's inner and middle magnetosphere. A density gradient, combined with a centrifugal force, provide appropriate RT onset conditions. Pressure balance is achieved by initializing two ion populations: one with fixed temperature, but varying density, and the other with fixed density, but a temperature gradient that offsets the density gradient from the first population and the centrifugal force (effective gravity). We first analyze two-dimensional results for the plane perpendicular to the magnetic field by comparing growth rates as a function of wave vector following Huba et al. (1998). Prescribed perpendicular wave modes are seeded with an initial velocity perturbation. We then extend the model to three dimensions, introducing a stabilizing parallel wave vector. Boundary conditions in the parallel direction prohibit motion of the magnetic field line footprints to model the eigenmodes of the magnetodisc's resonant cavity. We again compare growth rates based on perpendicular wave number, but also on the parallel extent of the resonant cavity, which fixes the size of the largest parallel wavelength. Finally, we search for evidence of strong guide field magnetic reconnection within the domain by identifying areas with large parallel electric fields or changes in magnetic field topology.

Subtle but ubiquitous selection on body size in a natural population of collared flycatchers over 33 years.

PubMed

Björklund, M; Gustafsson, L

2017-07-01

Understanding the magnitude and long-term patterns of selection in natural populations is of importance, for example, when analysing the evolutionary impact of climate change. We estimated univariate and multivariate directional, quadratic and correlational selection on four morphological traits (adult wing, tarsus and tail length, body mass) over a time period of 33 years (≈ 19 000 observations) in a nest-box breeding population of collared flycatchers (Ficedula albicollis). In general, selection was weak in both males and females over the years regardless of fitness measure (fledged young, recruits and survival) with only few cases with statistically significant selection. When data were analysed in a multivariate context and as time series, a number of patterns emerged; there was a consistent, but weak, selection for longer wings in both sexes, selection was stronger on females when the number of fledged young was used as a fitness measure, there were no indications of sexually antagonistic selection, and we found a negative correlation between selection on tarsus and wing length in both sexes but using different fitness measures. Uni- and multivariate selection gradients were correlated only for wing length and mass. Multivariate selection gradient vectors were longer than corresponding vector of univariate gradients and had more constrained direction. Correlational selection had little importance. Overall, the fitness surface was more or less flat with few cases of significant curvature, indicating that the adaptive peak with regard to body size in this species is broader than the phenotypic distribution, which has resulted in weak estimates of selection. © 2017 European Society For Evolutionary Biology. Journal of Evolutionary Biology © 2017 European Society For Evolutionary Biology.

Vaccination of plasmid DNA encoding ORF81 gene of CJ strains of KHV provides protection to immunized carp.

PubMed

Zhou, Jingxiang; Xue, Jiangdong; Wang, Qiuju; Zhu, Xia; Li, Xingwei; Lv, Wenliang; Zhang, Dongming

2014-06-01

In order to construct the recombinant plasmid of pIRES-ORF81, the nucleic acid isolated from Koi herpes virus-CJ (KHV-CJ) strains was used as a template to insert the ORF81 gene fragments amplified by PCR into the pIRES-neo, a kind of eukaryotic expression vector. Using Western blotting analysis, it was verified that ORF81 gene protein can be expressed correctly by pIRES-ORF81, after MFC cells were transfected. The recombinant plasmid pIRES-ORF81 was set into three immunization dose gradients: 1, 10, and 50 μg/carp. Empty plasmid group, PBS group, and blank control group were set simultaneously. Giving intramuscular injections to healthy carps with an average body mass of 246 ± 20 g, indirect ELISA was used to regularly determine antibody levels after three times immunization injection. Neutralizing antibodies were detected by neutralization assay. The results of inoculation tests showed that the pIRES-ORF81 recombinant plasmid can induce the production of carp-specific antibodies. The differences of immune effect between the three different doses of immune gradients were not significant (P > 0.05), but they can induce the production of neutralizing antibodies. After 25 d of inoculation, carp mortality of pIRES-neo empty vector treatment groups was 85%, while the carp mortality of eukaryotic expression recombinant plasmid pIRES-ORF81 injected with three different doses of immune gradients was 20, 17.5, and 12.5%, respectively. Differences in comparison to the control group were highly significant (P < 0.01). However, histopathological section of immunohistochemistry organization revealed no significant changes. It demonstrated that the DNA vaccine pIRES-ORF81 constructed in the experiment displayed a good protective effect against KHV, which had the potential to industrial applications.

Orientationally invariant metrics of apparent compartment eccentricity from double pulsed field gradient diffusion experiments.

PubMed

Jespersen, Sune Nørhøj; Lundell, Henrik; Sønderby, Casper Kaae; Dyrby, Tim B

2013-12-01

Pulsed field gradient diffusion sequences (PFG) with multiple diffusion encoding blocks have been indicated to offer new microstructural tissue information, such as the ability to detect nonspherical compartment shapes in macroscopically isotropic samples, i.e. samples with negligible directional signal dependence on diffusion gradients in standard diffusion experiments. However, current acquisition schemes are not rotationally invariant in the sense that the derived metrics depend on the orientation of the sample, and are affected by the interplay of sampling directions and compartment orientation dispersion when applied to macroscopically anisotropic systems. Here we propose a new framework, the d-PFG 5-design, to enable rotationally invariant estimation of double wave vector diffusion metrics (d-PFG). The method is based on the idea that an appropriate orientational average of the signal emulates the signal from a powder preparation of the same sample, where macroscopic anisotropy is absent by construction. Our approach exploits the theory of exact numerical integration (quadrature) of polynomials on the rotation group, and we exemplify the general procedure with a set consisting of 60 pairs of diffusion wave vectors (the d-PFG 5-design) facilitating a theoretically exact determination of the fourth order Taylor or cumulant expansion of the orientationally averaged signal. The d-PFG 5-design is evaluated with numerical simulations and ex vivo high field diffusion MRI experiments in a nonhuman primate brain. Specifically, we demonstrate rotational invariance when estimating compartment eccentricity, which we show offers new microstructural information, complementary to that of fractional anisotropy (FA) from diffusion tensor imaging (DTI). The imaging observations are supported by a new theoretical result, directly relating compartment eccentricity to FA of individual pores. Copyright © 2013 John Wiley & Sons, Ltd.

Turbulence-flame interactions in DNS of a laboratory high Karlovitz premixed turbulent jet flame

NASA Astrophysics Data System (ADS)

Wang, Haiou; Hawkes, Evatt R.; Chen, Jacqueline H.

2016-09-01

In the present work, direct numerical simulation (DNS) of a laboratory premixed turbulent jet flame was performed to study turbulence-flame interactions. The turbulent flame features moderate Reynolds number and high Karlovitz number (Ka). The orientations of the flame normal vector n, the vorticity vector ω and the principal strain rate eigenvectors ei are examined. The in-plane and out-of-plane angles are introduced to quantify the vector orientations, which also measure the flame geometry and the vortical structures. A general observation is that the distributions of these angles are more isotropic downstream as the flame and the flow become more developed. The out-of-plane angle of the flame normal vector, β, is a key parameter in developing the correction of 2D measurements to estimate the corresponding 3D quantities. The DNS results show that the correction factor is unity at the inlet and approaches its theoretical value of an isotropic distribution downstream. The alignment characteristics of n, ω and ei, which reflect the interactions of turbulence and flame, are also studied. Similar to a passive scalar gradient in non-reacting flows, the flame normal has a tendency to align with the most compressive strain rate, e3, in the flame, indicating that turbulence contributes to the production of scalar gradient. The vorticity dynamics are examined via the vortex stretching term, which was found to be the predominant source of vorticity generation balanced by dissipation, in the enstrophy transport equation. It is found that although the vorticity preferentially aligns with the intermediate strain rate, e2, the contribution of the most extensive strain rate, e1, to vortex stretching is comparable with that of the intermediate strain rate, e2. This is because the eigenvalue of the most extensive strain rate, λ1, is always large and positive. It is confirmed that the vorticity vector is preferentially positioned along the flame tangential plane, contributing to the dominance of cylindrical curvature of the flame front. Finally, the effect of heat release on the turbulence-flame interactions is examined. It is found that heat release has only limited impact on the statistics due to the minor role played by the strain rate induced by heat release rate in the current high Ka flame.

Epidemiological analysis of bovine ephemeral fever in 2012-2013 in the subtropical islands of Japan.

PubMed

Hayama, Yoko; Moriguchi, Sachiko; Yanase, Tohru; Suzuki, Moemi; Niwa, Tsuyoshi; Ikemiyagi, Kazufumi; Nitta, Yoshiki; Yamamoto, Takehisa; Kobayashi, Sota; Murai, Kiyokazu; Tsutsui, Toshiyuki

2016-03-09

Bovine ephemeral fever (BEF) is a febrile disease of cattle that is transmitted by arthropod vectors such as mosquitoes and Culicoides biting midges. An outbreak of BEF recently occurred in Ishigaki Island and surrounding islands that are located southwest of Japan. In this study, an epidemiological analysis was conducted to understand the temporal and spatial characteristics of the outbreak. Factors associated with the disease spread within Ishigaki Island were investigated by hierarchical Bayesian models. The possibility of between-island transmission by windborne vectors and transmission by long-distance migration of infected vectors were examined using atmospheric dispersion models. In September 2012, the first case of the disease was detected in the western part of Ishigaki Island. In 1 month, it had rapidly spread to the southern part of the island and to surrounding islands, and led to 225 suspected cases of BEF during the outbreak. The dispersion model demonstrated the high possibility of between-island transmission by wind. Spatial analysis showed that paddy fields, farmlands, and slope gradients had a significant impact on the 1-km cell-level incidence risk. These factors may have influenced the habitats and movements of the vectors with regard to the spread of BEF. A plausible incursion event of infected vectors from Southeast Asia to Ishigaki Island was estimated to have occurred at the end of August. This study revealed that the condition of a terrain and land use significantly influenced disease transmission. These factors are important in assessing favorable environments for related vectors. The results of the dispersion model indicated the likely transmission of the infected vectors by wind on the local scale and on the long-distance scale. These findings would be helpful for developing a surveillance program and developing preventive measures against BEF.

Using geometric morphometric visualizations of directional selection gradients to investigate morphological differentiation.

PubMed

Weaver, Timothy D; Gunz, Philipp

2018-04-01

Researchers studying extant and extinct taxa are often interested in identifying the evolutionary processes that have lead to the morphological differences among the taxa. Ideally, one could distinguish the influences of neutral evolutionary processes (genetic drift, mutation) from natural selection, and in situations for which selection is implicated, identify the targets of selection. The directional selection gradient is an effective tool for investigating evolutionary process, because it can relate form (size and shape) differences between taxa to the variation and covariation found within taxa. However, although most modern morphometric analyses use the tools of geometric morphometrics (GM) to analyze landmark data, to date, selection gradients have mainly been calculated from linear measurements. To address this methodological gap, here we present a GM approach for visualizing and comparing between-taxon selection gradients with each other, associated difference vectors, and "selection" gradients from neutral simulations. To exemplify our approach, we use a dataset of 347 three-dimensional landmarks and semilandmarks recorded on the crania of 260 primate specimens (112 humans, 67 common chimpanzees, 36 bonobos, 45 gorillas). Results on this example dataset show how incorporating geometric information can provide important insights into the evolution of the human braincase, and serve to demonstrate the utility of our approach for understanding morphological evolution. © 2018 The Author(s). Evolution © 2018 The Society for the Study of Evolution.

A Multi-Gradient Generator in a Single Microfluidic Device for Optical Microscopy and Interferometry

NASA Astrophysics Data System (ADS)

Bedrossian, Manuel; Nadeau, Jay; Lindensmith, Chris

2016-11-01

The goal of this work was to create a single microfluidic device capable of establishing multiple types of gradients in a quantifiable manner. Many microbial species are known to exhibit directed motility in the presence of stimuli. This phenomenon, known as taxis, can be used as a bio-signature and a means of identifying microorganisms. Directed microbial motility has been seen as a response to the presence of certain chemicals, light, heat, magnetic fields, and other stimuli. Microbial movement along the gradient vector, that cannot be explained by passive hydrodynamics or Brownian motion, can shed light on whether the sample contains living microbes or not. The ability to create multiple types of gradients in a single microfluidic device allows for high throughput testing of heterogeneous samples to detect taxis. There has been increased interest in the search for life within our solar system where liquid water is known to exist. Induced directional motility can serve as a viable method for detecting living organisms that actively respond to their environment. The device developed here includes a chemical, photonic, thermal, and magnetic gradient generator, while maintaining high optical quality in order to be used for microscopy as well as quantitative phase imaging This work was funded by the Gordon and Betty Moore Foundation, who the authors wish to thank for their generosity.

Wave propagation through an inhomogeneous slab sandwiched by the piezoelectric and the piezomagnetic half spaces.

PubMed

Jiao, Fengyu; Wei, Peijun; Li, Li

2017-01-01

Wave propagation through a gradient slab sandwiched by the piezoelectric and the piezomagnetic half spaces are studied in this paper. First, the secular equations in the transverse isotropic piezoelectric/piezomagnetic half spaces are derived from the general dynamic equation. Then, the state vectors at piezoelectric and piezomagnetic half spaces are related to the amplitudes of various possible waves. The state transfer equation of the functionally graded slab is derived from the equations of motion by the reduction of order, and the transfer matrix of the functionally gradient slab is obtained by solving the state transfer equation with the spatial-varying coefficient. Finally, the continuous interface conditions are used to lead to the resultant algebraic equations. The algebraic equations are solved to obtain the amplitude ratios of various waves which are further used to obtain the energy reflection and transmission coefficients of various waves. The numerical results are shown graphically and are validated by the energy conservation law. Based on the numerical results on the fives of gradient profiles, the influences of the graded slab on the wave propagation are discussed. It is found that the reflection and transmission coefficients are obviously dependent upon the gradient profile. The various surface waves are more sensitive to the gradient profile than the bulk waves. Copyright © 2016 Elsevier B.V. All rights reserved.

Parallel and fault-tolerant algorithms for hypercube multiprocessors

DOE Office of Scientific and Technical Information (OSTI.GOV)

Aykanat, C.

1988-01-01

Several techniques for increasing the performance of parallel algorithms on distributed-memory message-passing multi-processor systems are investigated. These techniques are effectively implemented for the parallelization of the Scaled Conjugate Gradient (SCG) algorithm on a hypercube connected message-passing multi-processor. Significant performance improvement is achieved by using these techniques. The SCG algorithm is used for the solution phase of an FE modeling system. Almost linear speed-up is achieved, and it is shown that hypercube topology is scalable for an FE class of problem. The SCG algorithm is also shown to be suitable for vectorization, and near supercomputer performance is achieved on a vectormore » hypercube multiprocessor by exploiting both parallelization and vectorization. Fault-tolerance issues for the parallel SCG algorithm and for the hypercube topology are also addressed.« less

Spatial Resolution Effect on Forest Road Gradient Calculation and Erosion Modelling

NASA Astrophysics Data System (ADS)

Cao, L.; Elliot, W.

2017-12-01

Road erosion is one of the main sediment sources in a forest watershed and should be properly evaluated. With the help of GIS technology, road topography can be determined and soil loss can be predicted at a watershed scale. As a vector geographical feature, the road gradient should be calculated following road direction rather than hillslope direction. This calculation might be difficult with a coarse (30-m) DEM which only provides the underlying topography information. This study was designed to explore the effect of road segmentation and DEM resolution on the road gradient calculation and erosion prediction at a watershed scale. The Water Erosion Prediction Project (WEPP) model was run on road segments of 9 lengths ranging from 40m to 200m. Road gradient was calculated from three DEM data sets: 1m LiDAR, and 10m and 30m USGS DEMs. The 1m LiDAR DEM calculated gradients were very close to the field observed road gradients, so we assumed the 1m LiDAR DEM predicted the true road gradient. The results revealed that longer road segments skipped detail topographical undulations and resulted in lower road gradients. Coarser DEMs computed steeper road gradients as larger grid cells covered more adjacent areas outside road resulting in larger elevation differences. Field surveyed results also revealed that coarser DEM might result in more gradient deviation in a curved road segment when it passes through a convex or concave slope. As road segment length increased, the gradient difference between three DEMs was reduced. There were no significant differences between road gradients of different segment lengths and DEM resolution when segments were longer than 100m. For long segments, the 10m DEM calculated road gradient was similar to the 1m LiDAR gradient. When evaluating the effects of road segment length, the predicted erosion rate decreased with increasing length when road gradient was less than 3%. In cases where the road gradients exceed 3% and rill erosion dominates, predicted erosion rates exponentially increased with segment length. At the watershed scale, most of the predicted soil loss occurred on segments with gradients ranging from 3% to 9%. Based on the road gradient calculated with the 10-m and 30-m DEMs, soil loss was overestimated when compared to the 1m LiDAR DEM. Both the 10m and 30m DEM result in similar total road soil loss.

A framework for developing a mimetic tensor artificial viscosity for Lagrangian hydrocodes on arbitrary polygonal and polyhedral meshes (u)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lipnikov, Konstantin; Shashkov, Mikhail

2011-01-11

We construct a new mimetic tensor artificial viscosity on general polygonal and polyhedral meshes. The tensor artificial viscosity is based on a mimetic discretization of coordinate invariant operators, divergence of a tensor and gradient of a vector. The focus of this paper is on the symmetric form, div ({mu},{var_epsilon}(u)), of the tensor artificial viscosity where {var_epsilon}(u) is the symmetrized gradient of u and {mu}, is a tensor. The mimetic discretizations of this operator is derived for the case of a full tensor coefficient {mu}, that may reflect a shock direction. We demonstrate performance of the new viscosity for the Nohmore » implosion, Sedov explosion and Saltzman piston problems in both Cartesian and axisymmetric coordinate systems.« less

Three-dimensional Gravity Inversion with a New Gradient Scheme on Unstructured Grids

NASA Astrophysics Data System (ADS)

Sun, S.; Yin, C.; Gao, X.; Liu, Y.; Zhang, B.

2017-12-01

Stabilized gradient-based methods have been proved to be efficient for inverse problems. Based on these methods, setting gradient close to zero can effectively minimize the objective function. Thus the gradient of objective function determines the inversion results. By analyzing the cause of poor resolution on depth in gradient-based gravity inversion methods, we find that imposing depth weighting functional in conventional gradient can improve the depth resolution to some extent. However, the improvement is affected by the regularization parameter and the effect of the regularization term becomes smaller with increasing depth (shown as Figure 1 (a)). In this paper, we propose a new gradient scheme for gravity inversion by introducing a weighted model vector. The new gradient can improve the depth resolution more efficiently, which is independent of the regularization parameter, and the effect of regularization term will not be weakened when depth increases. Besides, fuzzy c-means clustering method and smooth operator are both used as regularization terms to yield an internal consecutive inverse model with sharp boundaries (Sun and Li, 2015). We have tested our new gradient scheme with unstructured grids on synthetic data to illustrate the effectiveness of the algorithm. Gravity forward modeling with unstructured grids is based on the algorithm proposed by Okbe (1979). We use a linear conjugate gradient inversion scheme to solve the inversion problem. The numerical experiments show a great improvement in depth resolution compared with regular gradient scheme, and the inverse model is compact at all depths (shown as Figure 1 (b)). AcknowledgeThis research is supported by Key Program of National Natural Science Foundation of China (41530320), China Natural Science Foundation for Young Scientists (41404093), and Key National Research Project of China (2016YFC0303100, 2017YFC0601900). ReferencesSun J, Li Y. 2015. Multidomain petrophysically constrained inversion and geology differentiation using guided fuzzy c-means clustering. Geophysics, 80(4): ID1-ID18. Okabe M. 1979. Analytical expressions for gravity anomalies due to homogeneous polyhedral bodies and translations into magnetic anomalies. Geophysics, 44(4), 730-741.

Wavefront analysis from its slope data

NASA Astrophysics Data System (ADS)

Mahajan, Virendra N.; Acosta, Eva

2017-08-01

In the aberration analysis of a wavefront over a certain domain, the polynomials that are orthogonal over and represent balanced wave aberrations for this domain are used. For example, Zernike circle polynomials are used for the analysis of a circular wavefront. Similarly, the annular polynomials are used to analyze the annular wavefronts for systems with annular pupils, as in a rotationally symmetric two-mirror system, such as the Hubble space telescope. However, when the data available for analysis are the slopes of a wavefront, as, for example, in a Shack- Hartmann sensor, we can integrate the slope data to obtain the wavefront data, and then use the orthogonal polynomials to obtain the aberration coefficients. An alternative is to find vector functions that are orthogonal to the gradients of the wavefront polynomials, and obtain the aberration coefficients directly as the inner products of these functions with the slope data. In this paper, we show that an infinite number of vector functions can be obtained in this manner. We show further that the vector functions that are irrotational are unique and propagate minimum uncorrelated additive random noise from the slope data to the aberration coefficients.

Vertical amplitude phase structure of a low-frequency acoustic field in shallow water

NASA Astrophysics Data System (ADS)

Kuznetsov, G. N.; Lebedev, O. V.; Stepanov, A. N.

2016-11-01

We obtain in integral and analytic form the relations for calculating the amplitude and phase characteristics of an interference structure of orthogonal projections of the oscillation velocity vector in shallow water. For different frequencies and receiver depths, we numerically study the source depth dependences of the effective phase velocities of an equivalent plane wave, the orthogonal projections of the sound pressure phase gradient, and the projections of the oscillation velocity vector. We establish that at low frequencies in zones of interference maxima, independently of source depth, weakly varying effective phase velocity values are observed, which exceed the sound velocity in water by 5-12%. We show that the angles of arrival of the equivalent plane wave and the oscillation velocity vector in the general case differ; however, they virtually coincide in the zone of the interference maximum of the sound pressure under the condition that the horizontal projections of the oscillation velocity appreciably exceed the value of the vertical projection. We give recommendations on using the sound field characteristics in zones with maximum values for solving rangefinding and signal-detection problems.

The effect of the arbitrary level assignment of satellite cloud motion wind vectors on wind analyses in the pre-thunderstorm environment

NASA Technical Reports Server (NTRS)

Peslen, C. A.; Koch, S. E.; Uccellini, L. W.

1985-01-01

The impact of satellite-derived cloud motion vectors on SESAME rawinsonde wind fields was studied in two separate cases. The effect of wind and moisture gradients on the arbitrary assignment of the satellite data is assessed to coordinate surfaces in a severe storm environment marked by strong vertical wind shear. Objective analyses of SESAME rawinsonde winds and combined winds are produced and differences between these two analyzed fields are used to make an assessment of coordinate level choice. It is shown that the standard method of arbitrarily assigning wind vectors to a low level coordinate surface yields systematic differences between the rawinsonde and combined wind analyses. Arbitrary assignment of cloud motions to the 0.9 sigma surface produces smaller differences than assignment to the 825 mb pressure surface. Systematic differences occur near moisture discontinuities and in regions of horizontal and vertical wind shears. The differences between the combined and SESAME wind fields are made smallest by vertically interpolating cloud motions to either a pressure or sigma surface.

Molecular Cardiac Surgery with Recirculating Delivery (MCARD): Procedure and Vector Transfer.

PubMed

Katz, Michael G; Fargnoli, Anthony S; Kendle, Andrew P; Bridges, Charles R

2017-01-01

Despite progress in clinical treatment, cardiovascular diseases are still the leading cause of morbidity and mortality worldwide. Therefore, novel therapeutic approaches are needed, targeting the underlying molecular mechanisms of disease with improved outcomes for patients. Gene therapy is one of the most promising fields for the development of new treatments for the advanced stages of cardiovascular diseases. The establishment of clinically relevant methods of gene transfer remains one of the principal limitations on the effectiveness of gene therapy. Recently, there have been significant advances in direct and transvascular gene delivery methods. The ideal gene transfer method should be explored in clinically relevant large animal models of heart disease to evaluate the roles of specific molecular pathways in disease pathogenesis. Characteristics of the optimal technique for gene delivery include low morbidity, an increased myocardial transcapillary gradient, esxtended vector residence time in the myocytes, and the exclusion of residual vector from the systemic circulation after delivery to minimize collateral expression and immune response. Here we describe myocardial gene transfer techniques with molecular cardiac surgery with recirculating delivery in a large animal model of post ischemic heart failure.

Observed Structure and Characteristics of Cold Pools over Tropical Oceans using Vector Wind Retrievals and WRF simulations

NASA Astrophysics Data System (ADS)

Garg, P.; Nesbitt, S. W.; Lang, T. J.; Chronis, T.; Thayer, J. D.; Hence, D. A.

2017-12-01

Cold pools generated in the wake of convective activity can enhance the surface sensible heat flux, latent heat flux, and also changes in evaporation out of, and fresh water flux into, the ocean. Recent studies have shown that over the open ocean, cold pool outflow boundaries and their intersections can organize and initiate a spectrum of deep convective clouds, which is a key driver of shallow and deep convection over conditionally-unstable tropical oceans. The primary goal of this study is to understand the structure and characteristics of cold pools over the tropical oceans using observations. With the idea that cold pools will have strong wind gradients at their boundaries, we use ASCAT vector wind retrievals. We identify regions of steep gradients in wind vectors as gradient features (GFs), akin to cold pools. Corresponding to these GFs, sensible and latent heat fluxes were calculated using the observed winds and background temperatures from MERRA-2 reanalysis. To evaluate the proposed technique, cold pools were observed using S-PolKa radar from the DYNAMO/AMIE field campaign in the Indian Ocean for the period of 1 October 2011 to 31 March 2012 and were compared with ASCAT GFs. To relate the thermodynamic and kinematic characteristics of observed and simulated cold pools, simulations were carried out on WRF on a 3-km domain explicitly. The areas of cold pools were identified in the models using virtual temperature (Tv), which is a direct measure of air density, while GFs were identified using model simulated winds. Quantitative measures indicate that GFs are highly correspondent with model-simulated cold pools. In global measurements of cold pools from 2007-2015, it is possible to examine the characteristics of GFs across all tropical ocean basins, and relate them to meteorological conditions, as well as the characteristics of the parent precipitation systems. Our results indicate that while there is a general relationship between the amount of precipitation and the number of cold pools, the largest cold pools exist over the Eastern Pacific basin, where the most stratiform rain is produced from oceanic MCSs. It is anticipated that improved understanding of cold pools, which are a primary triggering mechanism of oceanic shallow and deep convection, will improve prediction of this important component of the climate system.

Evaluation of Spatially Targeted Strategies to Control Non-Domiciliated Triatoma dimidiata Vector of Chagas Disease

PubMed Central

Barbu, Corentin; Dumonteil, Eric; Gourbière, Sébastien

2011-01-01

Background Chagas disease is a major neglected tropical disease with deep socio-economical effects throughout Central and South America. Vector control programs have consistently reduced domestic populations of triatomine vectors, but non-domiciliated vectors still have to be controlled efficiently. Designing control strategies targeting these vectors is challenging, as it requires a quantitative description of the spatio-temporal dynamics of village infestation, which can only be gained from combinations of extensive field studies and spatial population dynamic modelling. Methodology/Principal Findings A spatially explicit population dynamic model was combined with a two-year field study of T. dimidiata infestation dynamics in the village of Teya, Mexico. The parameterized model fitted and predicted accurately both intra-annual variation and the spatial gradient in vector abundance. Five different control strategies were then applied in concentric rings to mimic spatial design targeting the periphery of the village, where vectors were most abundant. Indoor insecticide spraying and insect screens reduced vector abundance by up to 80% (when applied to the whole village), and half of this effect was obtained when control was applied only to the 33% of households closest to the village periphery. Peri-domicile cleaning was able to eliminate up to 60% of the vectors, but at the periphery of the village it has a low effect, as it is ineffective against sylvatic insects. The use of lethal traps and the management of house attractiveness provided similar levels of control. However this required either house attractiveness to be null, or ≥5 lethal traps, at least as attractive as houses, to be installed in each household. Conclusion/Significance Insecticide and insect screens used in houses at the periphery of the village can contribute to reduce house infestation in more central untreated zones. However, this beneficial effect remains insufficient to allow for a unique spatially targeted strategy to offer protection to all households. Most efficiently, control should combine the use of insect screens in outer zones to reduce infestation by both sylvatic and peri-domiciliated vectors, and cleaning of peri-domicile in the centre of the village where sylvatic vectors are absent. The design of such spatially mixed strategies of control offers a promising avenue to reduce the economic cost associated with the control of non-domiciliated vectors. PMID:21610862

Concepts and procedures required for successful reduction of tensor magnetic gradiometer data obtained from an unexploded ordnance detection demonstration at Yuma Proving Grounds, Arizona

USGS Publications Warehouse

Bracken, Robert E.; Brown, Philip J.

2006-01-01

On March 12, 2003, data were gathered at Yuma Proving Grounds, in Arizona, using a Tensor Magnetic Gradiometer System (TMGS). This report shows how these data were processed and explains concepts required for successful TMGS data reduction. Important concepts discussed include extreme attitudinal sensitivity of vector measurements, low attitudinal sensitivity of gradient measurements, leakage of the common-mode field into gradient measurements, consequences of thermal drift, and effects of field curvature. Spatial-data collection procedures and a spin-calibration method are addressed. Discussions of data-reduction procedures include tracking of axial data by mathematically matching transfer functions among the axes, derivation and application of calibration coefficients, calculation of sensor-pair gradients, thermal-drift corrections, and gradient collocation. For presentation, the magnetic tensor at each data station is converted to a scalar quantity, the I2 tensor invariant, which is easily found by calculating the determinant of the tensor. At important processing junctures, the determinants for all stations in the mapped area are shown in shaded relief map-view. Final processed results are compared to a mathematical model to show the validity of the assumptions made during processing and the reasonableness of the ultimate answer obtained.

Ship detection based on rotation-invariant HOG descriptors for airborne infrared images

NASA Astrophysics Data System (ADS)

Xu, Guojing; Wang, Jinyan; Qi, Shengxiang

2018-03-01

Infrared thermal imagery is widely used in various kinds of aircraft because of its all-time application. Meanwhile, detecting ships from infrared images attract lots of research interests in recent years. In the case of downward-looking infrared imagery, in order to overcome the uncertainty of target imaging attitude due to the unknown position relationship between the aircraft and the target, we propose a new infrared ship detection method which integrates rotation invariant gradient direction histogram (Circle Histogram of Oriented Gradient, C-HOG) descriptors and the support vector machine (SVM) classifier. In details, the proposed method uses HOG descriptors to express the local feature of infrared images to adapt to changes in illumination and to overcome sea clutter effects. Different from traditional computation of HOG descriptor, we subdivide the image into annular spatial bins instead of rectangle sub-regions, and then Radial Gradient Transform (RGT) on the gradient is applied to achieve rotation invariant histogram information. Considering the engineering application of airborne and real-time requirements, we use SVM for training ship target and non-target background infrared sample images to discriminate real ships from false targets. Experimental results show that the proposed method has good performance in both the robustness and run-time for infrared ship target detection with different rotation angles.

Nuclear Energy Gradients for Internally Contracted Complete Active Space Second-Order Perturbation Theory: Multistate Extensions.

PubMed

Vlaisavljevich, Bess; Shiozaki, Toru

2016-08-09

We report the development of the theory and computer program for analytical nuclear energy gradients for (extended) multistate complete active space perturbation theory (CASPT2) with full internal contraction. The vertical shifts are also considered in this work. This is an extension of the fully internally contracted CASPT2 nuclear gradient program recently developed for a state-specific variant by us [MacLeod and Shiozaki, J. Chem. Phys. 2015, 142, 051103]; in this extension, the so-called λ equation is solved to account for the variation of the multistate CASPT2 energies with respect to the change in the amplitudes obtained in the preceding state-specific CASPT2 calculations, and the Z vector equations are modified accordingly. The program is parallelized using the MPI3 remote memory access protocol that allows us to perform efficient one-sided communication. The optimized geometries of the ground and excited states of a copper corrole and benzophenone are presented as numerical examples. The code is publicly available under the GNU General Public License.

Research in computer science

NASA Technical Reports Server (NTRS)

Ortega, J. M.

1985-01-01

Synopses are given for NASA supported work in computer science at the University of Virginia. Some areas of research include: error seeding as a testing method; knowledge representation for engineering design; analysis of faults in a multi-version software experiment; implementation of a parallel programming environment; two computer graphics systems for visualization of pressure distribution and convective density particles; task decomposition for multiple robot arms; vectorized incomplete conjugate gradient; and iterative methods for solving linear equations on the Flex/32.

Cloud-to-ground lightning activity in Colombia and the influence of topography

NASA Astrophysics Data System (ADS)

Aranguren, D.; López, J.; Inampués, J.; Torres, H.; Betz, H.

2017-02-01

Lightning activity on the Colombian mountains, where the altitude varies from 0 to more than 5000 MSL, is studied based on VLF/LF lightning detection data and using a 2012-2013 dataset. The influence of altitude is observed by evaluating cloud-to-ground lightning incidence at different altitude intervals. The relationship between ground flash density and altitude gradient vectors is studied. Results show a clear dependence of the flash density on elevation.

An Alternative Treatment of Heat Flow for Charge Transport in Semiconductor Devices (Postprint)

DTIC Science & Technology

2010-07-01

is tantamount to treating them as ideal gases. A three-dimensional ideal Fermi gas is spherically symmetric in momentum space, and its distribution in...the first mo- ment of the Boltzmann equation using the momentum relax- ation time and effective mass approximations.13 Neglecting any magnetic field and...where the integral is over all momentum vectors k, v is electron velocity, k is the momentum relaxation time, and kf denotes the gradient in momentum

Dynamics of the EEG of human brain in the gradient magnetic fields of geological faults in different geographical and climatic zones

NASA Astrophysics Data System (ADS)

Pobachenko, S. V.; Sokolov, M. V.; Grigoriev, P. E.; Vasilieva, I. V.

2017-11-01

There are presented the results of experimental studies of the dynamics of indices of the functional state of a person located within the zones characterized by anomalous parameters of spatial distribution of magnetic field vector values. It is shown that these geophysical modifications have a pronounced effect on the dynamics of electrical activity indices of the human brain, regardless of geographic and climatic conditions.

Parametric FEM for geometric biomembranes

NASA Astrophysics Data System (ADS)

Bonito, Andrea; Nochetto, Ricardo H.; Sebastian Pauletti, M.

2010-05-01

We consider geometric biomembranes governed by an L2-gradient flow for bending energy subject to area and volume constraints (Helfrich model). We give a concise derivation of a novel vector formulation, based on shape differential calculus, and corresponding discretization via parametric FEM using quadratic isoparametric elements and a semi-implicit Euler method. We document the performance of the new parametric FEM with a number of simulations leading to dumbbell, red blood cell and toroidal equilibrium shapes while exhibiting large deformations.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Benioff, Paul

This paper describes the effects of a complex scalar scaling field on quantum mechanics. The field origin is an extension of the gauge freedom for basis choice in gauge theories to the underlying scalar field. The extension is based on the idea that the value of a number at one space time point does not determine the value at another point. This, combined with the description of mathematical systems as structures of different types, results in the presence of separate number fields and vector spaces as structures, at different space time locations. Complex number structures and vector spaces at eachmore » location are scaled by a complex space time dependent scaling factor. The effect of this scaling factor on several physical and geometric quantities has been described in other work. Here the emphasis is on quantum mechanics of one and two particles, their states and properties. Multiparticle states are also briefly described. The effect shows as a complex, nonunitary, scalar field connection on a fiber bundle description of nonrelativistic quantum mechanics. Here, the lack of physical evidence for the presence of this field so far means that the coupling constant of this field to fermions is very small. It also means that the gradient of the field must be very small in a local region of cosmological space and time. Outside this region, there are no restrictions on the field gradient.« less

Effects of a scalar scaling field on quantum mechanics

DOE PAGES

Benioff, Paul

2016-04-18

This paper describes the effects of a complex scalar scaling field on quantum mechanics. The field origin is an extension of the gauge freedom for basis choice in gauge theories to the underlying scalar field. The extension is based on the idea that the value of a number at one space time point does not determine the value at another point. This, combined with the description of mathematical systems as structures of different types, results in the presence of separate number fields and vector spaces as structures, at different space time locations. Complex number structures and vector spaces at eachmore » location are scaled by a complex space time dependent scaling factor. The effect of this scaling factor on several physical and geometric quantities has been described in other work. Here the emphasis is on quantum mechanics of one and two particles, their states and properties. Multiparticle states are also briefly described. The effect shows as a complex, nonunitary, scalar field connection on a fiber bundle description of nonrelativistic quantum mechanics. Here, the lack of physical evidence for the presence of this field so far means that the coupling constant of this field to fermions is very small. It also means that the gradient of the field must be very small in a local region of cosmological space and time. Outside this region, there are no restrictions on the field gradient.« less

Diagnosis of Tempromandibular Disorders Using Local Binary Patterns.

PubMed

Haghnegahdar, A A; Kolahi, S; Khojastepour, L; Tajeripour, F

2018-03-01

Temporomandibular joint disorder (TMD) might be manifested as structural changes in bone through modification, adaptation or direct destruction. We propose to use Local Binary Pattern (LBP) characteristics and histogram-oriented gradients on the recorded images as a diagnostic tool in TMD assessment. CBCT images of 66 patients (132 joints) with TMD and 66 normal cases (132 joints) were collected and 2 coronal cut prepared from each condyle, although images were limited to head of mandibular condyle. In order to extract features of images, first we use LBP and then histogram of oriented gradients. To reduce dimensionality, the linear algebra Singular Value Decomposition (SVD) is applied to the feature vectors matrix of all images. For evaluation, we used K nearest neighbor (K-NN), Support Vector Machine, Naïve Bayesian and Random Forest classifiers. We used Receiver Operating Characteristic (ROC) to evaluate the hypothesis. K nearest neighbor classifier achieves a very good accuracy (0.9242), moreover, it has desirable sensitivity (0.9470) and specificity (0.9015) results, when other classifiers have lower accuracy, sensitivity and specificity. We proposed a fully automatic approach to detect TMD using image processing techniques based on local binary patterns and feature extraction. K-NN has been the best classifier for our experiments in detecting patients from healthy individuals, by 92.42% accuracy, 94.70% sensitivity and 90.15% specificity. The proposed method can help automatically diagnose TMD at its initial stages.

Spinor Geometry and Signal Transmission in Three-Space

NASA Astrophysics Data System (ADS)

Binz, Ernst; Pods, Sonja; Schempp, Walter

2002-09-01

For a singularity free gradient field in an open set of an oriented Euclidean space of dimension three we define a natural principal bundle out of an immanent complex line bundle. The elements of both bundles are called internal variables. Several other natural bundles are associated with the principal bundle and, in turn, determine the vector field. Two examples are given and it is shown that for a constant vector field circular polarized waves travelling along a field line can be considered as waves of internal variables. Einstein's equation epsilon = m [middle dot] c2 is derived from the geometry of the principal bundle. On SU(2) a relation between spin representations and Schrodinger representations is established. The link between the spin 1/2-model and the Schrodinger representations yields a connection between a microscopic and a macroscopic viewpoint.

A real-time neutron-gamma discriminator based on the support vector machine method for the time-of-flight neutron spectrometer

NASA Astrophysics Data System (ADS)

Wei, ZHANG; Tongyu, WU; Bowen, ZHENG; Shiping, LI; Yipo, ZHANG; Zejie, YIN

2018-04-01

A new neutron-gamma discriminator based on the support vector machine (SVM) method is proposed to improve the performance of the time-of-flight neutron spectrometer. The neutron detector is an EJ-299-33 plastic scintillator with pulse-shape discrimination (PSD) property. The SVM algorithm is implemented in field programmable gate array (FPGA) to carry out the real-time sifting of neutrons in neutron-gamma mixed radiation fields. This study compares the ability of the pulse gradient analysis method and the SVM method. The results show that this SVM discriminator can provide a better discrimination accuracy of 99.1%. The accuracy and performance of the SVM discriminator based on FPGA have been evaluated in the experiments. It can get a figure of merit of 1.30.

Reduced state feedback gain computation. [optimization and control theory for aircraft control

NASA Technical Reports Server (NTRS)

Kaufman, H.

1976-01-01

Because application of conventional optimal linear regulator theory to flight controller design requires the capability of measuring and/or estimating the entire state vector, it is of interest to consider procedures for computing controls which are restricted to be linear feedback functions of a lower dimensional output vector and which take into account the presence of measurement noise and process uncertainty. Therefore, a stochastic linear model that was developed is presented which accounts for aircraft parameter and initial uncertainty, measurement noise, turbulence, pilot command and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was performed for both finite and infinite time performance indices without gradient computation by using Zangwill's modification of a procedure originally proposed by Powell. Results using a seventh order process show the proposed procedures to be very effective.

Computation of output feedback gains for linear stochastic systems using the Zangnill-Powell Method

NASA Technical Reports Server (NTRS)

Kaufman, H.

1975-01-01

Because conventional optimal linear regulator theory results in a controller which requires the capability of measuring and/or estimating the entire state vector, it is of interest to consider procedures for computing controls which are restricted to be linear feedback functions of a lower dimensional output vector and which take into account the presence of measurement noise and process uncertainty. To this effect a stochastic linear model has been developed that accounts for process parameter and initial uncertainty, measurement noise, and a restricted number of measurable outputs. Optimization with respect to the corresponding output feedback gains was then performed for both finite and infinite time performance indices without gradient computation by using Zangwill's modification of a procedure originally proposed by Powell. Results using a seventh order process show the proposed procedures to be very effective.

A Thermo-Optic Propagation Modeling Capability.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schrader, Karl; Akau, Ron

2014-10-01

A new theoretical basis is derived for tracing optical rays within a finite-element (FE) volume. The ray-trajectory equations are cast into the local element coordinate frame and the full finite-element interpolation is used to determine instantaneous index gradient for the ray-path integral equation. The FE methodology (FEM) is also used to interpolate local surface deformations and the surface normal vector for computing the refraction angle when launching rays into the volume, and again when rays exit the medium. The method is implemented in the Matlab(TM) environment and compared to closed- form gradient index models. A software architecture is also developedmore » for implementing the algorithms in the Zemax(TM) commercial ray-trace application. A controlled thermal environment was constructed in the laboratory, and measured data was collected to validate the structural, thermal, and optical modeling methods.« less

Development and implementation of (Q)SAR modeling within the CHARMMing web-user interface.

PubMed

Weidlich, Iwona E; Pevzner, Yuri; Miller, Benjamin T; Filippov, Igor V; Woodcock, H Lee; Brooks, Bernard R

2015-01-05

Recent availability of large publicly accessible databases of chemical compounds and their biological activities (PubChem, ChEMBL) has inspired us to develop a web-based tool for structure activity relationship and quantitative structure activity relationship modeling to add to the services provided by CHARMMing (www.charmming.org). This new module implements some of the most recent advances in modern machine learning algorithms-Random Forest, Support Vector Machine, Stochastic Gradient Descent, Gradient Tree Boosting, so forth. A user can import training data from Pubchem Bioassay data collections directly from our interface or upload his or her own SD files which contain structures and activity information to create new models (either categorical or numerical). A user can then track the model generation process and run models on new data to predict activity. © 2014 Wiley Periodicals, Inc.

Parallel Conjugate Gradient: Effects of Ordering Strategies, Programming Paradigms, and Architectural Platforms

NASA Technical Reports Server (NTRS)

Oliker, Leonid; Heber, Gerd; Biswas, Rupak

2000-01-01

The Conjugate Gradient (CG) algorithm is perhaps the best-known iterative technique to solve sparse linear systems that are symmetric and positive definite. A sparse matrix-vector multiply (SPMV) usually accounts for most of the floating-point operations within a CG iteration. In this paper, we investigate the effects of various ordering and partitioning strategies on the performance of parallel CG and SPMV using different programming paradigms and architectures. Results show that for this class of applications, ordering significantly improves overall performance, that cache reuse may be more important than reducing communication, and that it is possible to achieve message passing performance using shared memory constructs through careful data ordering and distribution. However, a multi-threaded implementation of CG on the Tera MTA does not require special ordering or partitioning to obtain high efficiency and scalability.

Solving large test-day models by iteration on data and preconditioned conjugate gradient.

PubMed

Lidauer, M; Strandén, I; Mäntysaari, E A; Pösö, J; Kettunen, A

1999-12-01

A preconditioned conjugate gradient method was implemented into an iteration on a program for data estimation of breeding values, and its convergence characteristics were studied. An algorithm was used as a reference in which one fixed effect was solved by Gauss-Seidel method, and other effects were solved by a second-order Jacobi method. Implementation of the preconditioned conjugate gradient required storing four vectors (size equal to number of unknowns in the mixed model equations) in random access memory and reading the data at each round of iteration. The preconditioner comprised diagonal blocks of the coefficient matrix. Comparison of algorithms was based on solutions of mixed model equations obtained by a single-trait animal model and a single-trait, random regression test-day model. Data sets for both models used milk yield records of primiparous Finnish dairy cows. Animal model data comprised 665,629 lactation milk yields and random regression test-day model data of 6,732,765 test-day milk yields. Both models included pedigree information of 1,099,622 animals. The animal model ¿random regression test-day model¿ required 122 ¿305¿ rounds of iteration to converge with the reference algorithm, but only 88 ¿149¿ were required with the preconditioned conjugate gradient. To solve the random regression test-day model with the preconditioned conjugate gradient required 237 megabytes of random access memory and took 14% of the computation time needed by the reference algorithm.

Standard Free Droplet Digital Polymerase Chain Reaction as a New Tool for the Quality Control of High-Capacity Adenoviral Vectors in Small-Scale Preparations

PubMed Central

Boehme, Philip; Stellberger, Thorsten; Solanki, Manish; Zhang, Wenli; Schulz, Eric; Bergmann, Thorsten; Liu, Jing; Doerner, Johannes; Baiker, Armin E.

2015-01-01

Abstract High-capacity adenoviral vectors (HCAdVs) are promising tools for gene therapy as well as for genetic engineering. However, one limitation of the HCAdV vector system is the complex, time-consuming, and labor-intensive production process and the following quality control procedure. Since HCAdVs are deleted for all viral coding sequences, a helper virus (HV) is needed in the production process to provide the sequences for all viral proteins in trans. For the purification procedure of HCAdV, cesium chloride density gradient centrifugation is usually performed followed by buffer exchange using dialysis or comparable methods. However, performing these steps is technically difficult, potentially error-prone, and not scalable. Here, we establish a new protocol for small-scale production of HCAdV based on commercially available adenovirus purification systems and a standard method for the quality control of final HCAdV preparations. For titration of final vector preparations, we established a droplet digital polymerase chain reaction (ddPCR) that uses a standard free-end-point PCR in small droplets of defined volume. By using different probes, this method is capable of detecting and quantifying HCAdV and HV in one reaction independent of reference material, rendering this method attractive for accurately comparing viral titers between different laboratories. In summary, we demonstrate that it is possible to produce HCAdV in a small scale of sufficient quality and quantity to perform experiments in cell culture, and we established a reliable protocol for vector titration based on ddPCR. Our method significantly reduces time and required equipment to perform HCAdV production. In the future the ddPCR technology could be advantageous for titration of other viral vectors commonly used in gene therapy. PMID:25640117

Purification of adenoviral vectors by combined anion exchange and gel filtration chromatography.

PubMed

Eglon, Marc N; Duffy, Aoife M; O'Brien, Timothy; Strappe, Padraig M

2009-11-01

Adenoviral vectors are used extensively in human gene therapy trials and in vaccine development. Large-scale GMP production requires a downstream purification process, and liquid chromatography is emerging as the most powerful mode of purification, enabling the production of vectors at a clinically relevant scale and quality. The present study describes the development of a two-step high-performance liquid chromatography (HPLC) process combining anion exchange (AIEX) and gel filtration (GF) in comparison with the caesium chloride density gradient method. HEK-293 cells were cultured in ten-layer CellStacks() and infected with 10 pfu/cell of adenoviral vector expressing green fluorescent protein (Ad5-GFP). Cell-bound virus was harvested and benzonase added to digest DNA, crude lysate was clarified by centrifugation and filtration prior to HPLC. Chromatography fractions were added to HEK-293 cells and GFP expression measured using a fluorescent plate reader. Using AIEX then GF resulted in an adenoviral vector with purity comparable to Ad5-GFP purified by CsCl, whereas the reverse process (GF-AIEX) showed a reduced purity by electrophoresis and required further buffer exchange of the product. The optimal process (AIEX-GF) resulted in a vector yield of 2.3 x 10(7) pfu/cm(2) of cell culture harvested compared to 3.3 x 10(7) pfu/cm(2) for CsCl. The process recovery for the HPLC process was 36% compared to 27.5% for CsCl and total virion to infectious particle ratios of 18 and 11, respectively, were measured. We present a simple two-step chromatography process that is capable of producing high-quality adenovirus at a titre suitable for scale-up and clinical translation.

Mapping the Conjugate Gradient Algorithm onto High Performance Heterogeneous Computers

DTIC Science & Technology

2014-05-01

Matrix Storage Formats According to J . Dongarra (Dongerra 2000), the efficiency of most iterative methods, such as CG, can be attributed to the...valh = aij) ⇒ (colh = j ). The ptr integer vector is of length n + 1 and contains the index in val where each matrix row starts. For example, the...first nonzero element of matrix rowm is found at index ptrm of val. By convention, ptrn+1 ≡ nz + 1. Notice that (aij) ⇒ (ptri ≤ j < ptri+1) for all i. An

Time-Parallel Solutions to Ordinary Differential Equations on GPUs with a New Functional Optimization Approach Related to the Sobolev Gradient Method

DTIC Science & Technology

2012-10-01

black and approximations in cyan and magenta. The second ODE is the pendulum equation, given by: This ODE was also implemented using Crank...The drawback of approaches like the one proposed can be observed with a very simple example. Suppose vector is found by applying 4 linear...public release; distribution unlimited Figure 2. A phase space plot of the Pendulum example. Fine solution (black) contains 32768 time steps

Direct Visualization of Surfaces from Computed Tomography Data,

DTIC Science & Technology

1988-01-01

dominant method - slice-by-slice - makes comprehension of convoluted, small, or faint structures difficult. From a densitometric point of view, the human ...Di t SPeclal 4 N(xI) = Vf(lx) IVf (XI)I* There are many ways to estimate the gradient vector Vf(xi). The selection of an operator depends on the...the data is relatively noise-free or has been pre-smoothed. A graph of a(xi) as a function of f(xi) and IVf (xl)l for three tissue types A, B, and C

Generation of Human Induced Pluripotent Stem Cells from Peripheral Blood Mononuclear Cells Using Sendai Virus.

PubMed

Soares, Filipa A C; Pedersen, Roger A; Vallier, Ludovic

2016-01-01

This protocol describes the efficient isolation of peripheral blood mononuclear cells from circulating blood via density gradient centrifugation and subsequent generation of integration-free human induced pluripotent stem cells. Peripheral blood mononuclear cells are cultured for 9 days to allow expansion of the erythroblast population. The erythroblasts are then used to derive human induced pluripotent stem cells using Sendai viral vectors, each expressing one of the four reprogramming factors Oct4, Sox2, Klf4, and c-Myc.

On-Time 3D Time-Domain EMI and Tensor Magnetic Gradiometry for UXO Detection and Discrimination

DTIC Science & Technology

2008-06-04

three-axis fluxgate magnetometers mounted on a tetrahedral structure as shown in figure 5.1.3.1. TMGS was intended to measure gradients of vector...manufacturer of the fluxgate magnetometers supplied specifications for each of the triaxial sensors (figure 5.3.2.25). Figure 5.3.2.25...Manufacturer’s specifications for the four triaxial fluxgate magnetometers used in the TMGS planar array. As shown previously in figure 5.3.2.4

The effect of axial crucible rotation on microstructural uniformity during horizontal directional solidification

NASA Technical Reports Server (NTRS)

Grugel, R. N.; Kim, Shinwoo; Woodward, Tracey; Wang, T. G.

1992-01-01

The effects on microstructure of crucible orientation with respect to the earth's gravitational vector, g, during directional solidification of low-volume fraction copper and aluminum, Pb-Cu, and Sn-Al alloys are examined. It is demonstrated that horizontal alignment (i.e. perpendicular to g) in combination with axial rotation of the crucible during growth is sufficient to negate factors which initiate macrosegregation, e.g. density gradients attributed to temperature and/or compositional differences, and promotes a uniform microstructure.

Deployment and Intelligent Nanosatellite Operations Colorado Final Technical Report

DTIC Science & Technology

2006-09-28

environmental factors will cause disturbance torques during orbit around the Earth . These factors are solar radiation pressure from the sun , aerodynamic...software. The 3- axis sensing of the magnetometer allows a vector the B- field of the Earth to be sensed. Geopack 2003 then can be utilized with the orbit ...gradient torque can be represented as the following: g, ; 3wo21 Eq. 2-11 where ow is the angular velocity of the spacecraft as it orbits the earth . DINO’s

Observations of wave-induced pore pressure gradients and bed level response on a surf zone sandbar

NASA Astrophysics Data System (ADS)

Anderson, Dylan; Cox, Dan; Mieras, Ryan; Puleo, Jack A.; Hsu, Tian-Jian

2017-06-01

Horizontal and vertical pressure gradients may be important physical mechanisms contributing to onshore sediment transport beneath steep, near-breaking waves in the surf zone. A barred beach was constructed in a large-scale laboratory wave flume with a fixed profile containing a mobile sediment layer on the crest of the sandbar. Horizontal and vertical pore pressure gradients were obtained by finite differences of measurements from an array of pressure transducers buried within the upper several centimeters of the bed. Colocated observations of erosion depth were made during asymmetric wave trials with wave heights between 0.10 and 0.98 m, consistently resulting in onshore sheet flow sediment transport. The pore pressure gradient vector within the bed exhibited temporal rotations during each wave cycle, directed predominantly upward under the trough and then rapidly rotating onshore and downward as the wavefront passed. The magnitude of the pore pressure gradient during each phase of rotation was correlated with local wave steepness and relative depth. Momentary bed failures as deep as 20 grain diameters were coincident with sharp increases in the onshore-directed pore pressure gradients, but occurred at horizontal pressure gradients less than theoretical critical values for initiation of the motion for compact beds. An expression combining the effects of both horizontal and vertical pore pressure gradients with bed shear stress and soil stability is used to determine that failure of the bed is initiated at nonnegligible values of both forces.