On acceleration of Krylov-subspace-based Newton and Arnoldi iterations for incompressible CFD: replacing time steppers and generation of initial guess

We propose two techniques aimed at improving the convergence rate of steady state and eigenvalue solvers preconditioned by the inverse Stokes operator and realized via time-stepping. First, we suggest a generalization of the Stokes operator so that the resulting preconditioner operator depends on several parameters and whose action preserves zero divergence and boundary conditions. The parameters can be tuned for each problem to speed up the convergence of a Krylov-subspace-based linear algebra solver. This operator can be inverted by the Uzawa-like algorithm, and does not need a time-stepping. Second, we propose to generate an initial guess of steady flow, leading eigenvalue and eigenvector using orthogonal projection on a divergence-free basis satisfying all boundary conditions. The approach, including the two proposed techniques, is illustrated on the solution of the linear stability problem for laterally heated square and cubic cavities

An in vitro evaluation of standard rotational thromboelastography in monitoring of effects of recombinant factor VIIa on coagulopathy induced by hydroxy ethyl starch

BACKGROUND: Rotational thromboelastography (ROTEG) has been proposed as a monitoring tool that can be used to monitor treatment of hemophilia with recombinant factor VIIa (rFVIIa). In these studies special non-standard reagents were used as activators of the coagulation. The aim of this study was to evaluate if standard ROTEG analysis could be used for monitoring of effects of recombinant factor VIIa (rFVIIa) on Hydroxy Ethyl Starch-induced dilutional coagulopathy. METHODS: The study was performed in vitro on healthy volunteers. Prothrombin time (PT) and ROTEG analysis were performed after dilution with 33% hydroxy ethyl starch and also after addition of rFVIIa to the diluted blood. RESULTS: PT was impaired with INR changing from 0.9 before dilution to 1.2 after dilution while addition of rFVIIa to diluted blood lead to an overcorrection of the PT to an International Normalized Ratio (INR) value of 0.6 (p = 0.01). ROTEG activated with the contact activator ellagic acid was impaired by hemodilution (p = 0.01) while addition of rFVIIa had no further effects. ROTEG activated with tissue factor (TF) was also impaired by hemodilution (p = 0.01) while addition of rFVIIa lead to further impairment of the coagulation (p = 0.01). CONCLUSIONS: The parameters affected in the ROTEG analysis were Clot Formation Time and Amplitude after 15 minutes while the Clotting Time was unaffected. We believe these effects to be due to methodological problems when using standard activators of the coagulation in the ROTEG analysis in combination with rFVIIa

Effect of non-random mating on genomic and BLUP selection schemes

<p>Abstract</p> <p>Background</p> <p>The risk of long-term unequal contribution of mating pairs to the gene pool is that deleterious recessive genes can be expressed. Such consequences could be alleviated by appropriately designing and optimizing breeding schemes i.e. by improving selection and mating procedures.</p> <p>Methods</p> <p>We studied the effect of mating designs, random, minimum coancestry and minimum covariance of ancestral contributions on rate of inbreeding and genetic gain for schemes with different information sources, i.e. sib test or own performance records, different genetic evaluation methods, i.e. BLUP or genomic selection, and different family structures, i.e. factorial or pair-wise.</p> <p>Results</p> <p>Results showed that substantial differences in rates of inbreeding due to mating design were present under schemes with a pair-wise family structure, for which minimum coancestry turned out to be more effective to generate lower rates of inbreeding. Specifically, substantial reductions in rates of inbreeding were observed in schemes using sib test records and BLUP evaluation. However, with a factorial family structure, differences in rates of inbreeding due mating designs were minor. Moreover, non-random mating had only a small effect in breeding schemes that used genomic evaluation, regardless of the information source.</p> <p>Conclusions</p> <p>It was concluded that minimum coancestry remains an efficient mating design when BLUP is used for genetic evaluation or when the size of the population is small, whereas the effect of non-random mating is smaller in schemes using genomic evaluation.</p

Application of the Jacobi Davidson method for spectral low-rank preconditioning in computational electromagnetics problems

[EN] We consider the numerical solution of linear systems arising from computational electromagnetics applications. For large scale problems the solution is usually obtained iteratively with a Krylov subspace method. It is well known that for ill conditioned problems the convergence of these methods can be very slow or even it may be impossible to obtain a satisfactory solution. To improve the convergence a preconditioner can be used, but in some cases additional strategies are needed. In this work we study the application of spectral lowrank updates (SLRU) to a previously computed sparse approximate inverse preconditioner.The updates are based on the computation of a small subset of the eigenpairs closest to the origin. Thus, the performance of the SLRU technique depends on the method available to compute the eigenpairs of interest. The SLRU method was first used using the IRA s method implemented in ARPACK. In this work we investigate the use of a Jacobi Davidson method, in particular its JDQR variant. The results of the numerical experiments show that the application of the JDQR method to obtain the spectral low-rank updates can be quite competitive compared with the IRA s method.Mas Marí, J.; Cerdán Soriano, JM.; Malla Martínez, N.; Marín Mateos-Aparicio, J. (2015). Application of the Jacobi Davidson method for spectral low-rank preconditioning in computational electromagnetics problems. Journal of the Spanish Society of Applied Mathematics. 67:39-50. doi:10.1007/s40324-014-0025-6S395067Bergamaschi, L., Pini, G., Sartoretto, F.: Computational experience with sequential, and parallel, preconditioned Jacobi–Davidson for large sparse symmetric matrices. J. Comput. Phys. 188(1), 318–331 (2003)Carpentieri, B.: Sparse preconditioners for dense linear systems from electromagnetics applications. PhD thesis, Institut National Polytechnique de Toulouse, CERFACS (2002)Carpentieri, B., Duff, I.S., Giraud, L.: Sparse pattern selection strategies for robust Frobenius-norm minimization preconditioners in electromagnetism. Numer. Linear Algebr. Appl. 7(7–8), 667–685 (2000)Carpentieri, B., Duff, I.S., Giraud, L.: A class of spectral two-level preconditioners. SIAM J. Sci. Comput. 25(2), 749–765 (2003)Carpentieri, B., Duff, I.S., Giraud, L., Magolu monga Made, M.: Sparse symmetric preconditioners for dense linear systems in electromagnetism. Numer. Linear Algebr. Appl. 11(8–9), 753–771 (2004)Carpentieri, B., Duff, I.S., Giraud, L., Sylvand, G.: Combining fast multipole techniques and an approximate inverse preconditioner for large electromagnetism calculations. SIAM J. Sci. Comput. 27(3), 774–792 (2005)Darve, E.: The fast multipole method I: error analysis and asymptotic complexity. SIAM J. Numer. Anal. 38(1), 98–128 (2000)Fokkema, D.R., Sleijpen, G.L., Van der Vorst, H.A.: Jacobi–Davidson style QR and QZ algorithms for the reduction of matrix pencils. SIAM J. Sci. Comput. 20(1), 94–125 (1998)Greengard, L., Rokhlin, V.: A fast algorithm for particle simulations. J. Comput. Phys. 73(3), 325–348 (1987)Grote, M., Huckle, T.: Parallel preconditioning with sparse approximate inverses. SIAM J. Sci. Comput. 18(3), 838–853 (1997)Harrington, R.: Origin and development of the method of moments for field computation. IEEE Antenna Propag. Mag. (1990)Kunz, K.S., Luebbers, R.J.: The finite difference time domain method for electromagnetics. SIAM J. Sci. Comput. 18(3), 838–853 (1997)Maxwell, J.C.: A dynamical theory of the electromagnetic field. Roy. S. Trans. CLV, (1864). Reprinted in Tricker, R. A. R. The Contributions of Faraday and Maxwell to Electrial Science, Pergamon Press (1966)Marín, J., Malla M.: Some experiments preconditioning via spectral low rank updates for electromagnetism applications. In: Proceedings of the international conference on preconditioning techniques for large sparse matrix problems in scientific and industrial applications (Preconditioning 2007), Toulouse (2007)Meijerink, J.A., van der Vorst, H.A.: An iterative solution method for linear systems of which the coefficient matrix is a symmetric M-matrix. Math. Comput. 31, 148–162 (1977)Sorensen, D.C., Lehoucq, R.B., Yang, C.: ARPACK users’ guide: solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods. SIAM, Philadelphia (1998)Rao, S.M., Wilton, D.R., Glisson, A.W.: Electromagnetic scattering by surfaces of arbitrary shape. IEEE Trans. Antenna Propag. 30, 409–418 (1982)Saad, Y.: Iterative methods for sparse linear systems. PWS Publishing Company, Boston (1996)Silvester, P.P., Ferrari, R.L.: Finite elements for electrical engineers. Cambridge University Press, Cambridge (1990)Sleijpen, S.L., van der Vorst, H.A.: A Jacobi–Davidson iteration method for linear eigenvalue problems. SIAM J. Matrix Anal. Appl. 17, 401–425 (1996)van der Vorst, H.A.: Bi-CGSTAB: a fast and smoothly converging variant of Bi-CG for the solution of non-symmetric linear systems. SIAM J. Sci. Stat. Comput. 12(6), 631–644 (1992

Selective advantage of implementing optimal contributions selection and timescales for the convergence of long-term genetic contributions

Additional file 1: Figure S1. Accumulation of long-term genetic contributions over time with a line of best fit

A survey of failed post-retained restorations

Survival of endodontically treated, post-restored teeth depends on a multitude of factors, all of which are practically impossible to include in a randomized, controlled clinical study. The purpose of this survey was to characterize and analyze reported failures of post-retained restorations to identify factors critical to failure and to type of failure. A questionnaire was mailed to private practitioners in Denmark with a request to complete the questionnaire whenever a patient presented with a failed post-retained restoration. Information was gathered on factors related to the patient, the tooth, the restorative materials, and the techniques. Two-hundred and sixty questionnaires were collected from 171 practitioners over a 3-year period. Functioning time until failure varied between 3 months and 38 years. Mean survival time until failure was 11 years. Of the failed restorations, 61% had functioned for 10 years or less. Fracture of the tooth was the most common type of failure reported, followed by loosening of the post and fracture of the post. Tapered posts implied an increased risk of tooth fracture compared to loosening or fracture of the post, and the relative risk of tooth fracture increased with the functioning time until failure. Fracture of the post was more common among male than female patients. On the basis of this survey of failed post-retained restorations, it was concluded that tapered posts were associated with a higher risk of tooth fracture than were parallel-sided posts

Prognosis of ovarian cancer subsequent to venous thromboembolism: a nationwide Danish cohort study

BACKGROUND: Venous thromboembolism (VTE) is associated with ovarian cancer and may impact the prognosis of ovarian cancer. Our aims were to examine the extent of disease at the time of the diagnosis of ovarian cancer and to estimate the impact of VTE on survival of ovarian cancer. METHODS: We identified 12,835 ovarian cancer patients diagnosed from 1980 to 2003 in the Danish Cancer Registry and obtained information on previous primary VTE diagnosis from the Danish National Hospital Discharge Registry. Ovarian cancer patients with previous VTE related to other cancers, surgery, or pregnancy were excluded. The vital status was determined by linking data to the Civil Registration System. RESULTS: We identified 50 ovarian cancer patients diagnosed less than 4 months after the VTE and 78 ovarian cancer patients diagnosed more than 4 months after the VTE diagnosis. Advanced stages tended to be more common among patients with VTE. One-year survivals were 44% and 54% among the two VTE groups, compared with 63% among patients without VTE. Adjusted (for age, calendar time, comorbidity, and FIGO-stage) mortality ratios were 1.7 (95% CI = 1.2–2.5) and 1.2 (95% CI = 0.8–1.7), respectively. CONCLUSION: Ovarian cancer diagnosed less than four months before VTE is associated with an advanced stage and a poorer prognosis

Inexperienced clinicians can extract pathoanatomic information from MRI narrative reports with high reproducibility for use in research/quality assurance

<p>Abstract</p> <p>Background</p> <p>Although reproducibility in reading MRI images amongst radiologists and clinicians has been studied previously, no studies have examined the reproducibility of inexperienced clinicians in extracting pathoanatomic information from magnetic resonance imaging (MRI) narrative reports and transforming that information into quantitative data. However, this process is frequently required in research and quality assurance contexts. The purpose of this study was to examine inter-rater reproducibility (agreement and reliability) among an inexperienced group of clinicians in extracting spinal pathoanatomic information from radiologist-generated MRI narrative reports.</p> <p>Methods</p> <p>Twenty MRI narrative reports were randomly extracted from an institutional database. A group of three physiotherapy students independently reviewed the reports and coded the presence of 14 common pathoanatomic findings using a categorical electronic coding matrix. Decision rules were developed after initial coding in an effort to resolve ambiguities in narrative reports. This process was repeated a further three times using separate samples of 20 MRI reports until no further ambiguities were identified (total n = 80). Reproducibility between trainee clinicians and two highly trained raters was examined in an arbitrary coding round, with agreement measured using percentage agreement and reliability measured using unweighted Kappa (<it>k</it>). Reproducibility was then examined in another group of three trainee clinicians who had not participated in the production of the decision rules, using another sample of 20 MRI reports.</p> <p>Results</p> <p>The mean percentage agreement for paired comparisons between the initial trainee clinicians improved over the four coding rounds (97.9-99.4%), although the greatest improvement was observed after the first introduction of coding rules. High inter-rater reproducibility was observed between trainee clinicians across 14 pathoanatomic categories over the four coding rounds (agreement range: 80.8-100%; reliability range <it>k </it>= 0.63-1.00). Concurrent validity was high in paired comparisons between trainee clinicians and highly trained raters (agreement 97.8-98.1%, reliability <it>k </it>= 0.83-0.91). Reproducibility was also high in the second sample of trainee clinicians (inter-rater agreement 96.7-100.0% and reliability <it>k </it>= 0.76-1.00; intra-rater agreement 94.3-100.0% and reliability <it>k </it>= 0.61-1.00).</p> <p>Conclusions</p> <p>A high level of radiological training is not required in order to transform MRI-derived pathoanatomic information from a narrative format to a quantitative format with high reproducibility for research or quality assurance purposes.</p

Integrated Polygenic Tool Substantially Enhances Coronary Artery Disease Prediction

BACKGROUND: There is considerable interest in whether genetic data can be used to improve standard cardiovascular disease risk calculators, as the latter are routinely used in clinical practice to manage preventative treatment. METHODS: Using the UK Biobank resource, we developed our own polygenic risk score for coronary artery disease (CAD). We used an additional 60 000 UK Biobank individuals to develop an integrated risk tool (IRT) that combined our polygenic risk score with established risk tools (either the American Heart Association/American College of Cardiology pooled cohort equations [PCE] or UK QRISK3), and we tested our IRT in an additional, independent set of 186 451 UK Biobank individuals. RESULTS: The novel CAD polygenic risk score shows superior predictive power for CAD events, compared with other published polygenic risk scores, and is largely uncorrelated with PCE and QRISK3. When combined with PCE into an IRT, it has superior predictive accuracy. Overall, 10.4% of incident CAD cases were misclassified as low risk by PCE and correctly classified as high risk by the IRT, compared with 4.4% misclassified by the IRT and correctly classified by PCE. The overall net reclassification improvement for the IRT was 5.9% (95% CI, 4.7–7.0). When individuals were stratified into age-by-sex subgroups, the improvement was larger for all subgroups (range, 8.3%–15.4%), with the best performance in 40- to 54-year-old men (15.4% [95% CI, 11.6–19.3]). Comparable results were found using a different risk tool (QRISK3) and also a broader definition of cardiovascular disease. Use of the IRT is estimated to avoid up to 12 000 deaths in the United States over a 5-year period. CONCLUSIONS: An IRT that includes polygenic risk outperforms current risk stratification tools and offers greater opportunity for early interventions. Given the plummeting costs of genetic tests, future iterations of CAD risk tools would be enhanced with the addition of a person’s polygenic risk