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.

The effect of flex fuel vehicles in the Brazilian light road transportation

The retaking of the ethanol program in the year 2003 as a fuel for light road transportation in Brazil through the introduction of flex fuel vehicles fleet was a good strategy to overcome the difficulties of the ethanol production sector and did work to increase its market share relative to gasoline. This process, however, may cause a future disequilibrium on the food production and on the refining oil derivates structure. In order to analyze the substitution process resultant of the competition between two opponents fighting for the same market, in this case the gasoline/ethanol substitution process, a method derived from the biomathematics based on the non-linear differential equations (NLDE) system is utilized. A brief description of the method is presented. Numerical adherence of the method to explain several substitution phenomena that occurred in the past is presented in the previous author`s paper, in which the urban gas pipeline system substitution of bottled LPG in the dwelling sector and the substitution of the urban diesel transportation fleet by compressed natural gas (CNG) buses is presented. The proposed method is particularly suitable for prospective analysis and scenarios assessment. (c) 2008 Elsevier Ltd. All rights reserved.

Age-Related Maximum Heart Rate Among Ischemic and Nonischemic Heart Failure Patients Receiving beta-Blockade Therapy

Background: Equations to predict maximum heart rate (HRmax) in heart failure (HF) patients receiving beta-adrenergic blocking (BB) agents do not consider the cause of HF. We determined equations to predict HRmax in patients with ischemic and nonischemic HF receiving BB therapy. Methods and Results: Using treadmill cardiopulmonary exercise testing, we studied HF patients receiving BB therapy being considered for transplantation from 1999 to 2010. Exclusions were pacemaker and/or implantable defibrillator, left ventricle ejection fraction (LVEF) >50%, peak respiratory exchange ratio (RER) <1.00, and Chagas disease. We used linear regression equations to predict HRmax based on age in ischemic and nonischemic patients. We analyzed 278 patients, aged 47 +/- 10 years, with ischemic (n = 75) and nonischemic (n = 203) HF. LVEF was 30.8 +/- 9.4% and 28.6 +/- 8.2% (P = .04), peak VO2 16.9 +/- 4.7 and 16.9 +/- 5.2 mL kg(-1) min(-1) (P = NS), and the HRmax 130.8 +/- 23.3 and 125.3 +/- 25.3 beats/min (P = .051) in ischemic and nonischemic patients, respectively. We devised the equation HRmax = 168 - 0.76 x age (R-2 = 0.095; P = .007) for ischemic HF patients, but there was no significant relationship between age and HRmax in nonischemic HF patients (R-2 = 0.006; P = NS). Conclusions: Our study suggests that equations to estimate HRmax should consider the cause of HF. (J Cardiac Fail 2012;18:831-836)

Maximal exercise test is a useful method for physical capacity and oxygen consumption determination in streptozotocin-diabetic rats

Abstract Background The aim of the present study was to investigate the relationship between speed during maximum exercise test (ET) and oxygen consumption (VO2) in control and STZ-diabetic rats, in order to provide a useful method to determine exercise capacity and prescription in researches involving STZ-diabetic rats. Methods Male Wistar rats were divided into two groups: control (CG, n = 10) and diabetic (DG, n = 8). The animals were submitted to ET on treadmill with simultaneous gas analysis through open respirometry system. ET and VO2 were assessed 60 days after diabetes induction (STZ, 50 mg/Kg). Results VO2 maximum was reduced in STZ-diabetic rats (72.5 ± 1 mL/Kg/min-1) compared to CG rats (81.1 ± 1 mL/Kg/min-1). There were positive correlations between ET speed and VO2 (r = 0.87 for CG and r = 0.8 for DG), as well as between ET speed and VO2 reserve (r = 0.77 for CG and r = 0.7 for DG). Positive correlations were also obtained between measured VO2 and VO2 predicted values (r = 0.81 for CG and r = 0.75 for DG) by linear regression equations to CG (VO2 = 1.54 * ET speed + 52.34) and DG (VO2 = 1.16 * ET speed + 51.99). Moreover, we observed that 60% of ET speed corresponded to 72 and 75% of VO2 reserve for CG and DG, respectively. The maximum ET speed was also correlated with VO2 maximum for both groups (CG: r = 0.7 and DG: r = 0.7). Conclusion These results suggest that: a) VO2 and VO2 reserve can be estimated using linear regression equations obtained from correlations with ET speed for each studied group; b) exercise training can be prescribed based on ET in control and diabetic-STZ rats; c) physical capacity can be determined by ET. Therefore, ET, which involves a relatively simple methodology and low cost, can be used as an indicator of cardio-respiratory capacity in future studies that investigate the physiological effect of acute or chronic exercise in control and STZ-diabetic male rats.

Error Autocorrelation and Linear Regression for Temperature-Based Evapotranspiration Estimates Improvement

Estimates of evapotranspiration on a local scale is important information for agricultural and hydrological practices. However, equations to estimate potential evapotranspiration based only on temperature data, which are simple to use, are usually less trustworthy than the Food and Agriculture Organization (FAO)Penman-Monteith standard method. The present work describes two correction procedures for potential evapotranspiration estimates by temperature, making the results more reliable. Initially, the standard FAO-Penman-Monteith method was evaluated with a complete climatologic data set for the period between 2002 and 2006. Then temperature-based estimates by Camargo and Jensen-Haise methods have been adjusted by error autocorrelation evaluated in biweekly and monthly periods. In a second adjustment, simple linear regression was applied. The adjusted equations have been validated with climatic data available for the Year 2001. Both proposed methodologies showed good agreement with the standard method indicating that the methodology can be used for local potential evapotranspiration estimates.

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.

Probabilistic crack growth analyses using a boundary element model: Applications in linear elastic fracture and fatigue problems

This paper addresses the numerical solution of random crack propagation problems using the coupling boundary element method (BEM) and reliability algorithms. Crack propagation phenomenon is efficiently modelled using BEM, due to its mesh reduction features. The BEM model is based on the dual BEM formulation, in which singular and hyper-singular integral equations are adopted to construct the system of algebraic equations. Two reliability algorithms are coupled with BEM model. The first is the well known response surface method, in which local, adaptive polynomial approximations of the mechanical response are constructed in search of the design point. Different experiment designs and adaptive schemes are considered. The alternative approach direct coupling, in which the limit state function remains implicit and its gradients are calculated directly from the numerical mechanical response, is also considered. The performance of both coupling methods is compared in application to some crack propagation problems. The investigation shows that direct coupling scheme converged for all problems studied, irrespective of the problem nonlinearity. The computational cost of direct coupling has shown to be a fraction of the cost of response surface solutions, regardless of experiment design or adaptive scheme considered. (C) 2012 Elsevier Ltd. All rights reserved.

The concept of quasi-integrability for modified non-linear Schrodinger models

We consider modifications of the nonlinear Schrodinger model (NLS) to look at the recently introduced concept of quasi-integrability. We show that such models possess an in finite number of quasi-conserved charges which present intriguing properties in relation to very specific space-time parity transformations. For the case of two-soliton solutions where the fields are eigenstates of this parity, those charges are asymptotically conserved in the scattering process of the solitons. Even though the charges vary in time their values in the far past and the far future are the same. Such results are obtained through analytical and numerical methods, and employ adaptations of algebraic techniques used in integrable field theories. Our findings may have important consequences on the applications of these models in several areas of non-linear science. We make a detailed numerical study of the modified NLS potential of the form V similar to (vertical bar psi vertical bar(2))(2+epsilon), with epsilon being a perturbation parameter. We perform numerical simulations of the scattering of solitons for this model and find a good agreement with the results predicted by the analytical considerations. Our paper shows that the quasi-integrability concepts recently proposed in the context of modifications of the sine-Gordon model remain valid for perturbations of the NLS model.

Optimal mean-variance control for discrete-time linear systems with Markovian jumps and multiplicative noises

In this paper, we consider the stochastic optimal control problem of discrete-time linear systems subject to Markov jumps and multiplicative noises under two criteria. The first one is an unconstrained mean-variance trade-off performance criterion along the time, and the second one is a minimum variance criterion along the time with constraints on the expected output. We present explicit conditions for the existence of an optimal control strategy for the problems, generalizing previous results in the literature. We conclude the paper by presenting a numerical example of a multi-period portfolio selection problem with regime switching in which it is desired to minimize the sum of the variances of the portfolio along the time under the restriction of keeping the expected value of the portfolio greater than some minimum values specified by the investor. (C) 2011 Elsevier Ltd. All rights reserved.