On the Convergence of Two-Stage Iterative Processes for Solving Linear Equations

This paper considers two-stage iterative processes for solving the linear system $Af = b$. The outer iteration is defined by $Mf^{k + 1} = Nf^k + b$, where $M$ is a nonsingular matrix such that $M - N = A$. At each stage $f^{k + 1} $ is computed approximately using an inner iteration process to solve $Mv = Nf^k + b$ for $v$. At the $k$th outer iteration, $p_k $ inner iterations are performed. It is shown that this procedure converges if $p_k \geqq P$ for some $P$ provided that the inner iteration is convergent and that the outer process would converge if $f^{k + 1} $ were determined exactly at every step. Convergence is also proved under more specialized conditions, and for the procedure where $p_k = p$ for all $k$, an estimate for $p$ is obtained which optimizes the convergence rate. Examples are given for systems arising from the numerical solution of elliptic partial differential equations and numerical results are presented.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Asymptotic behavior at infinity of solutions of multidimensional second kind integral equations

We consider second kind integral equations of the form x(s) - (abbreviated x - K x = y ), in which Ω is some unbounded subset of Rn. Let Xp denote the weighted space of functions x continuous on Ω and satisfying x (s) = O(|s|-p ),s → ∞We show that if the kernel k(s,t) decays like |s — t|-q as |s — t| → ∞ for some sufficiently large q (and some other mild conditions on k are satisfied), then K ∈ B(XP) (the set of bounded linear operators on Xp), for 0 ≤ p ≤ q. If also (I - K)-1 ∈ B(X0) then (I - K)-1 ∈ B(XP) for 0 < p < q, and (I- K)-1∈ B(Xq) if further conditions on k hold. Thus, if k(s, t) = O(|s — t|-q). |s — t| → ∞, and y(s)=O(|s|-p), s → ∞, the asymptotic behaviour of the solution x may be estimated as x (s) = O(|s|-r), |s| → ∞, r := min(p, q). The case when k(s,t) = к(s — t), so that the equation is of Wiener-Hopf type, receives especial attention. Conditions, in terms of the symbol of I — K, for I — K to be invertible or Fredholm on Xp are established for certain cases (Ω a half-space or cone). A boundary integral equation, which models three-dimensional acoustic propaga-tion above flat ground, absorbing apart from an infinite rigid strip, illustrates the practical application and sharpness of the above results. This integral equation mod-els, in particular, road traffic noise propagation along an infinite road surface sur-rounded by absorbing ground. We prove that the sound propagating along the rigid road surface eventually decays with distance at the same rate as sound propagating above the absorbing ground.

Aggregate vs disaggregate data analysis: a paradox in the estimation of money demand function of Japan under the low interest rate policy

We investigate the issue of whether there was a stable money demand function for Japan in 1990's using both aggregate and disaggregate time series data. The aggregate data appears to support the contention that there was no stable money demand function. The disaggregate data shows that there was a stable money demand function. Neither was there any indication of the presence of liquidity trapo Possible sources of discrepancy are explored and the diametrically opposite results between the aggregate and disaggregate analysis are attributed to the neglected heterogeneity among micro units. We also conduct simulation analysis to show that when heterogeneity among micro units is present. The prediction of aggregate outcomes, using aggregate data is less accurate than the prediction based on micro equations. Moreover. policy evaluation based on aggregate data can be grossly misleading.

Development of equations to estimate microbial contamination in ruminal incubation residues of forage produced under tropical conditions using 15N as a label1

The objective of this study was to use 15N to label microbial cells to allow development of equations for estimating the microbial contamination in ruminal in situ incubation residues of forage produced under tropical conditions. A total of 24 tropical forages were ruminal incubated in 3 steers at 3 separate times. To determine microbial contamination of the incubated residues, ruminal bacteria were labeled with 15N by continuous intraruminal infusion 60 h before the first incubation and continued until the last day of incubation. Ruminal digesta was collected for the isolation of bacteria before the first infusion of 15N on adaptation period and after the infusion of 15N on collection period. To determine the microbial contamination of CP fractions, restricted models were compared with the full model using the model identity test. A value of the corrected fraction A was estimated from the corresponding noncorrected fraction by this equation: Corrected A fraction (ACPC) = 1.99286 + 0.98256 × A fraction without correction (ACPWC). The corrected fraction B was estimated from the corresponding noncorrected fraction and from CP, NDF, neutral detergent insoluble protein (NDIP), and indigestible NDF (iNDF) using the equation corrected B fraction (BCPC) = -17.2181 - 0.0344 × fraction B without correction (BCPWC) + 0.65433 × CP + 1.03787 × NDF + 2.66010 × NDIP - 0.85979 × iNDF. The corrected degradation rate of B fraction (kd)was estimated using the equation corrected degradation rate of B fraction (kdCPC) = 0.04667 + 0.35139 × degradation rate of B fraction without correction (kdCPWC) + 0.0020 × CP - 0.00055839 × NDF - 0.00336 × NDIP + 0.00075089 × iNDF. This equation was obtained to estimate the contamination using CP of the feeds: %C = 79.21 × (1 - e-0.0555t) × e-0.0874CP. It was concluded that A and B fractions and kd of CP could be highly biased by microbial CP contamination, and therefore these corrected values could be obtained mathematically, replacing the use of microbial markers. The percentage of contamination and the corrected apparent degradability of CP could be obtained from values of CP and time of incubation for each feed, which could reduce cost and labor involved when using 15N. © 2013 American Society of Animal Science. All rights reserved.

Modeling the growth rate and lag time of different strains of Salmonella enterica and Listeria monocytogenes in ready-to-eat lettuce

The growth parameters (growth rate, mu and lag time, lambda) of three different strains each of Salmonella enterica and Listeria monocytogenes in minimally processed lettuce (MPL) and their changes as a function of temperature were modeled. MPL were packed under modified atmosphere (5% O-2, 15% CO2 and 80% N-2), stored at 7-30 degrees C and samples collected at different time intervals were enumerated for S. enterica and L monocytogenes. Growth curves and equations describing the relationship between mu and lambda as a function of temperature were constructed using the DMFit Excel add-in and through linear regression, respectively. The predicted growth parameters for the pathogens observed in this study were compared to ComBase, Pathogen modeling program (PMP) and data from the literature. High R-2 values (0.97 and 0.93) were observed for average growth curves of different strains of pathogens grown on MPL Secondary models of mu and lambda for both pathogens followed a linear trend with high R2 values (>0.90). Root mean square error (RMSE) showed that the models obtained are accurate and suitable for modeling the growth of S. enterica and L monocytogenes in MP lettuce. The current study provides growth models for these foodborne pathogens that can be used in microbial risk assessment. (C) 2011 Elsevier Ltd. 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.

Irreversible spherical model and its stationary entropy production rate

The nonequilibrium stationary state of an irreversible spherical model is investigated on hypercubic lattices. The model is defined by Langevin equations similar to the reversible case, but with asymmetric transition rates. In spite of being irreversible, we have succeeded in finding an explicit form for the stationary probability distribution, which turns out to be of the Boltzmann-Gibbs type. This enables one to evaluate the exact form of the entropy production rate at the stationary state, which is non-zero if the dynamical rules of the transition rates are asymmetric.

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)

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.

[Calculating the basal metabolic rate and severe and morbid obesity]

The aim of this study was to evaluate the currently available predictive equations for basal metabolic rate (BMR) in subjects with obesity class II and III, and to assess the contribution by the components of a two-compartment model of body composition, namely the lean body mass (LBM) and the fat mass (FM) to the prediction. A second objective was to examine the reliability of the Harris Benedict equation in obese subjects, especially with a weight > or = 120 kg.

A comparison of estimated glomerular filtration rates using Cockcroft−Gault and the Chronic Kidney Disease Epidemiology Collaboration estimating equations in HIV infection

OBJECTIVES: The aim of this study was to determine whether the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI)- or Cockcroft-Gault (CG)-based estimated glomerular filtration rates (eGFRs) performs better in the cohort setting for predicting moderate/advanced chronic kidney disease (CKD) or end-stage renal disease (ESRD). METHODS: A total of 9521 persons in the EuroSIDA study contributed 133 873 eGFRs. Poisson regression was used to model the incidence of moderate and advanced CKD (confirmed eGFR < 60 and < 30 mL/min/1.73 m(2) , respectively) or ESRD (fatal/nonfatal) using CG and CKD-EPI eGFRs. RESULTS: Of 133 873 eGFR values, the ratio of CG to CKD-EPI was ≥ 1.1 in 22 092 (16.5%) and the difference between them (CG minus CKD-EPI) was ≥ 10 mL/min/1.73 m(2) in 20 867 (15.6%). Differences between CKD-EPI and CG were much greater when CG was not standardized for body surface area (BSA). A total of 403 persons developed moderate CKD using CG [incidence 8.9/1000 person-years of follow-up (PYFU); 95% confidence interval (CI) 8.0-9.8] and 364 using CKD-EPI (incidence 7.3/1000 PYFU; 95% CI 6.5-8.0). CG-derived eGFRs were equal to CKD-EPI-derived eGFRs at predicting ESRD (n = 36) and death (n = 565), as measured by the Akaike information criterion. CG-based moderate and advanced CKDs were associated with ESRD [adjusted incidence rate ratio (aIRR) 7.17; 95% CI 2.65-19.36 and aIRR 23.46; 95% CI 8.54-64.48, respectively], as were CKD-EPI-based moderate and advanced CKDs (aIRR 12.41; 95% CI 4.74-32.51 and aIRR 12.44; 95% CI 4.83-32.03, respectively). CONCLUSIONS: Differences between eGFRs using CG adjusted for BSA or CKD-EPI were modest. In the absence of a gold standard, the two formulae predicted clinical outcomes with equal precision and can be used to estimate GFR in HIV-positive persons.

Ingestion rate and egg production of copepods collected in the upper 20m in the eastern Mediterranean Sea in April 2008 during SES_GR2

The ingestion on ciliates and phytoplankton dataset is based on samples taken during October 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod ingestion was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Copepods for the experiments were obtained with slow non-quantitative tows from the upper 20 m layer of the water column using 200 µm mesh size nets fitted with a large non-filtering cod end. For the grazing experiments we used the following copepod species: Clausocalanus furcatus, and Temoraa stylifera according to the relevant reference (Bamstedt et al. 2000). Copepod clearance rates on ciliates were calculated according to Frost equations (Frost 1972). Ingestion rates were calculated by multiplying clearance rates by the initial standing stocks (Bamstedt et al. 2000). The egg production dataset is based on samples taken during October 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod egg production was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Egg production rates of the dominant calanoid copepods were determined by incubation of fertilised females (eggs female/day) collected in the 0-20m layer. Copepod egg production was measured for the copepods Clausocalanus furcatus, Temora stylifera. On board experiments for the estimation of copepod egg production were taken place. For the estimation of copepod production (mgC/m**2/day), lengths (copepods and eggs) were converted to body carbon (Hopcroft et al., 1998) and production was estimated from biomass and weight-specific egg production rates, by assuming that those rates are representative for juvenile specific growth rates (Berggreen et al., 1988).

Ingestion rate and egg production of copepods collected in the Turkish Strait during Bilim 2 cruise in April 2008

The copepod Ingestion on ciliates, phytoplankton and the copepod production dataset is based on samples taken during April 2008 in Dardanelles Straits, Marmara Sea and Bosporus Straits at the third priority stations. These experiments were set up according to DoW of Sesame project. Copepods for the experiments were obtained with slow non-quantitative tows from the upper 50 m layer of the water column using 200 µm mesh size nets fitted with a large non-filtering cod end. For the grazing experiments we used the following copepod species: Centropages typicus and Acartia clausi according to the relevant reference (Bamstedt et al. 2000). Copepod clearance rates on ciliates were calculated according to Frost equations (Frost 1972). Ingestion rates were calculated by multiplying clearance rates by the initial standing stocks (Bamstedt et al. 2000). Egg production rates of the dominant calanoid copepods were determined by incubation of fertilised females (eggs/female/day) collected in the 0-20m layer. Copepod egg production was measured for the copepods Centropages typicus and Acartia clausi. On board experiments for the estimation of copepod egg production were taken place. For the estimation of copepod production (mg/m**2/day), lengths (copepods and eggs) were converted to body carbon (Hopcroft et al., 1998) and production was estimated from biomass and weight-specific egg production rates, by assuming that those rates are representative for juvenile specific growth rates (Berggreen et al., 1988).

Ingestion rate and egg production of copepods collected in the upper 20m in the eastern Mediterranean Sea in March and April 2008 during SES_GR1

The ingestion on ciliates and phytoplankton dataset is based on samples taken during April 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod ingestion was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Copepods for the experiments were obtained with slow non-quantitative tows from the upper 20 m layer of the water column using 200 µm mesh size nets fitted with a large non-filtering cod end. For the grazing experiments we used the following copepod species: Centropages typicus and Calanus helgolandicus according to the relevant reference (Bamstedt et al. 2000). Copepod clearance rates on ciliates were calculated according to Frost equations (Frost 1972). Ingestion rates were calculated by multiplying clearance rates by the initial standing stocks (Bamstedt et al. 2000). The egg production dataset is based on samples taken during April 2008 in Northern Aegean Sea, the area influenced by the Black Sea water outflow. A Lagrangian experiment was established and copepod egg production was estimated from experiments performed at stations according to the different positions of drifters during the cruise. Egg production rates of the dominant calanoid copepods were determined by incubation of fertilised females (eggs female/day) collected in the 0-20m layer. Copepod egg production was measured for the copepods Centropages typicus, Calanus helgolandicus. On board experiments for the estimation of copepod egg production were taken place. For the estimation of copepod production (mgC/ m**2 /day), lengths (copepods and eggs) were converted to body carbon (Hopcroft et al., 1998) and production was estimated from biomass and weight-specific egg production rates, by assuming that those rates are representative for juvenile specific growth rates (Berggreen et al., 1988).