Matrix-vector multiplication and triangular linear solver using GPGPU for symmetric positive definite matrices derived from elliptic equations

The modern GPUs are well suited for intensive computational tasks and massive parallel computation. Sparse matrix multiplication and linear triangular solver are the most important and heavily used kernels in scientific computation, and several challenges in developing a high performance kernel with the two modules is investigated. The main interest it to solve linear systems derived from the elliptic equations with triangular elements. The resulting linear system has a symmetric positive definite matrix. The sparse matrix is stored in the compressed sparse row (CSR) format. It is proposed a CUDA algorithm to execute the matrix vector multiplication using directly the CSR format. A dependence tree algorithm is used to determine which variables the linear triangular solver can determine in parallel. To increase the number of the parallel threads, a coloring graph algorithm is implemented to reorder the mesh numbering in a pre-processing phase. The proposed method is compared with parallel and serial available libraries. The results show that the proposed method improves the computation cost of the matrix vector multiplication. The pre-processing associated with the triangular solver needs to be executed just once in the proposed method. The conjugate gradient method was implemented and showed similar convergence rate for all the compared methods. The proposed method showed significant smaller execution time.

Existence and concentration of solutions for a class of biharmonic equations

Some superlinear fourth order elliptic equations are considered. A family of solutions is proved to exist and to concentrate at a point in the limit. The proof relies on variational methods and makes use of a weak version of the Ambrosetti-Rabinowitz condition. The existence and concentration of solutions are related to a suitable truncated equation. (C) 2012 Elsevier Inc. All rights reserved.

MULTIPLICITY OF NONRADIAL SOLUTIONS FOR A CLASS OF QUASILINEAR EQUATIONS ON ANNULUS WITH EXPONENTIAL CRITICAL GROWTH

In this paper, we establish the existence of many rotationally non-equivalent and nonradial solutions for the following class of quasilinear problems (p) {-Delta(N)u = lambda f(vertical bar x vertical bar, u) x is an element of Omega(r), u > 0 x is an element of Omega(r), u = 0 x is an element of Omega(r), where Omega(r) = {x is an element of R-N : r < vertical bar x vertical bar < r + 1}, N >= 2, N not equal 3, r >0, lambda > 0, Delta(N)u = div(vertical bar del u vertical bar(N-2)del u) is the N-Laplacian operator and f is a continuous function with exponential critical growth.

A nonlinear elliptic problem with terms concentrating in the boundary

In this paper, we investigate the behavior of a family of steady-state solutions of a nonlinear reaction diffusion equation when some reaction and potential terms are concentrated in a e-neighborhood of a portion G of the boundary. We assume that this e-neighborhood shrinks to G as the small parameter e goes to zero. Also, we suppose the upper boundary of this e-strip presents a highly oscillatory behavior. Our main goal here was to show that this family of solutions converges to the solutions of a limit problem, a nonlinear elliptic equation that captures the oscillatory behavior. Indeed, the reaction term and concentrating potential are transformed into a flux condition and a potential on G, which depends on the oscillating neighborhood. Copyright (C) 2012 John Wiley & Sons, Ltd.

MULTI-BUMP SOLUTIONS FOR A CLASS OF QUASILINEAR EQUATIONS ON R

This paper is concerned with the existence of multi-bump solutions to a class of quasilinear Schrodinger equations in R. The proof relies on variational methods and combines some arguments given by del Pino and Felmer, Ding and Tanaka, and Sere.

SUPERLINEAR ELLIPTIC PROBLEMS WITH SIGN CHANGING COEFFICIENTS

Via variational methods, we study multiplicity of solutions for the problem {-Delta u = lambda b(x)vertical bar u vertical bar(q-2)u + au + g(x, u) in Omega, u - 0 on partial derivative Omega, where a simple example for g( x, u) is |u|(p-2)u; here a, lambda are real parameters, 1 < q < 2 < p <= 2* and b(x) is a function in a suitable space L-sigma. We obtain a class of sign changing coefficients b(x) for which two non-negative solutions exist for any lambda > 0, and a total of five nontrivial solutions are obtained when lambda is small and a >= lambda(1). Note that this type of results are valid even in the critical case.

Continuity of Dynamical Structures for Nonautonomous Evolution Equations Under Singular Perturbations

In this paper we study the continuity of invariant sets for nonautonomous infinite-dimensional dynamical systems under singular perturbations. We extend the existing results on lower-semicontinuity of attractors of autonomous and nonautonomous dynamical systems. This is accomplished through a detailed analysis of the structure of the invariant sets and its behavior under perturbation. We prove that a bounded hyperbolic global solutions persists under singular perturbations and that their nonlinear unstable manifold behave continuously. To accomplish this, we need to establish results on roughness of exponential dichotomies under these singular perturbations. Our results imply that, if the limiting pullback attractor of a nonautonomous dynamical system is the closure of a countable union of unstable manifolds of global bounded hyperbolic solutions, then it behaves continuously (upper and lower) under singular perturbations.

On abstract differential equations with nonlocal conditions involving the temporal derivative of the solution

In this paper we discuss the existence of solutions for a class of abstract differential equations with nonlocal conditions for which the nonlocal term involves the temporal derivative of the solution. Some concrete applications to parabolic differential equations with nonlocal conditions are considered. (C) 2012 Royal Dutch Mathematical Society (KWG). Published by Elsevier B.V. All rights reserved.

Directed and elliptic flow of charged particles in Cu plus Cu collisions at root s(NN)=22.4 GeV

This paper reports results for directed flow v(1) and elliptic flow v(2) of charged particles in Cu + Cu collisions at root s(NN) = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4-GeV Cu + Cu collisions the prior observation that v1 is independent of the system size at 62.4 and 200 GeV and also extend the scaling of v(1) with eta/y(beam) to this system. The measured v(2)(p(T)) in Cu + Cu collisions is similar for root s(NN) throughout the range 22.4 to 200 GeV. We also report a comparison with results from transport model (ultrarelativistic quantum molecular dynamics and multiphase transport model) calculations. The model results do not agree quantitatively with the measured v(1)(eta), v(2)(p(T)), and v(2)(eta).

Measure functional differential equations and functional dynamic equations on time scales

We study measure functional differential equations and clarify their relation to generalized ordinary differential equations. We show that functional dynamic equations on time scales represent a special case of measure functional differential equations. For both types of equations, we obtain results on the existence and uniqueness of solutions, continuous dependence, and periodic averaging.

Boundedness of solutions of retarded functional differential equations with variable impulses via generalized ordinary differential equations

In this paper, we give sufficient conditions for the uniform boundedness and uniform ultimate boundedness of solutions of a class of retarded functional differential equations with impulse effects acting on variable times. We employ the theory of generalized ordinary differential equations to obtain our results. As an example, we investigate the boundedness of the solution of a circulating fuel nuclear reactor model.

GROUND STATE AND NON-GROUND STATE SOLUTIONS OF SOME STRONGLY COUPLED ELLIPTIC SYSTEMS

We study an elliptic system of the form Lu = vertical bar v vertical bar(p-1) v and Lv = vertical bar u vertical bar(q-1) u in Omega with homogeneous Dirichlet boundary condition, where Lu := -Delta u in the case of a bounded domain and Lu := -Delta u + u in the cases of an exterior domain or the whole space R-N. We analyze the existence, uniqueness, sign and radial symmetry of ground state solutions and also look for sign changing solutions of the system. More general non-linearities are also considered.

On the generalized canonical equations of Hamilton for a time-dependent mass particle

Fundamental principles of mechanics were primarily conceived for constant mass systems. Since the pioneering works of Meshcherskii (see historical review in Mikhailov (Mech. Solids 10(5):32-40, 1975), efforts have been made in order to elaborate an adequate mathematical formalism for variable mass systems. This is a current research field in theoretical mechanics. In this paper, attention is focused on the derivation of the so-called 'generalized canonical equations of Hamilton' for a variable mass particle. The applied technique consists in the consideration of the mass variation process as a dissipative phenomenon. Kozlov's (Stek. Inst. Math 223:178-184, 1998) method, originally devoted to the derivation of the generalized canonical equations of Hamilton for dissipative systems, is accordingly extended to the scenario of variable mass systems. This is done by conveniently writing the flux of kinetic energy from or into the variable mass particle as a 'Rayleigh-like dissipation function'. Cayley (Proc. R Soc. Lond. 8:506-511, 1857) was the first scholar to propose such an analogy. A deeper discussion on this particular subject will be left for a future paper.

Hyperbolic-like estimates for higher order equations

The main goal of this paper is to derive long time estimates of the energy for the higher order hyperbolic equations with time-dependent coefficients. in particular, we estimate the energy in the hyperbolic zone of the extended phase space by means of a function f (t) which depends on the principal part and on the coefficients of the terms of order m - 1. Then we look for sufficient conditions that guarantee the same energy estimate from above in all the extended phase space. We call this class of estimates hyperbolic-like since the energy behavior is deeply depending on the hyperbolic structure of the equation. In some cases, these estimates produce a dissipative effect on the energy. (C) 2012 Elsevier Inc. All rights reserved.

Vectorial capacity, basic reproduction number, force of infection and all that: formal notation to complete and adjust their classical concepts and equations

A dimensional analysis of the classical equations related to the dynamics of vector-borne infections is presented. It is provided a formal notation to complete the expressions for the Ross' threshold theorem, the Macdonald's basic reproduction "rate" and sporozoite "rate", Garret-Jones' vectorial capacity and Dietz-Molineaux-Thomas' force of infection. The analysis was intended to provide a formal notation that complete the classical equations proposed by these authors.