Genetic relationships between reproductive and production traits in Jersey crossbred cattle.

The study aimed to estimate the genetic parameters of different reproductive traits namely age at first calving (AFC), calving interval (CI), days open (DO) and number of service per conception (NSPC) and their associations with productive traits including 305-day milk yield (305DMY), total lactation milk yield (TLMY) and lactation length (LL) of Jersey crossbred cattle maintained at Kalyani, Nadia, West Bengal, India. Genetic parameters of reproductive traits and their correlations with productive traits were estimated by Restricted Maximum Likelihood method and Bayesian approach. Using both analytical approaches, the estimates of heritability for AFC, CI, DO and NSPC ranged from 0.12 -0.15, 0.05-0.08, 0.08-0.09 and 0.04-0.06, respectively. Low proportion of variances associated with permanent environmental effect of animals (c2 effect) were detected for CI (0.08-0.10), DO (0.09-0.11) and NSPC (0.05-0.06) in both the methods. Repeatability measures for all the reproductive traits considered in this study were low to moderate in nature, which ranged from 0.09 to 0.17. Genetic correlations between different reproductive traits were positive and low (0.05) to high (0.98) in magnitude except AFC-NSPC. Low and negative genetic correlations of AFC with 305DMY and TLMY were favourable and indicated animals with high milk yield had early age of maturity. Positive genetic correlations between CI, DO and NSPC with all production traits implied the antagonism relationships among these traits, therefore in any breeding program for improvement of production traits via selection, the reproductive traits should be taken into account as well.

Toward GPGPU accelerated human electromechanical cardiac simulations.

In this paper, we look at the acceleration of weakly coupled electromechanics using the graphics processing unit (GPU). Specifically, we port to the GPU a number of components of CHeart--a CPU-based finite element code developed for simulating multi-physics problems. On the basis of a criterion of computational cost, we implemented on the GPU the ODE and PDE solution steps for the electrophysiology problem and the Jacobian and residual evaluation for the mechanics problem. Performance of the GPU implementation is then compared with single core CPU (SC) execution as well as multi-core CPU (MC) computations with equivalent theoretical performance. Results show that for a human scale left ventricle mesh, GPU acceleration of the electrophysiology problem provided speedups of 164 × compared with SC and 5.5 times compared with MC for the solution of the ODE model. Speedup of up to 72 × compared with SC and 2.6 × compared with MC was also observed for the PDE solve. Using the same human geometry, the GPU implementation of mechanics residual/Jacobian computation provided speedups of up to 44 × compared with SC and 2.0 × compared with MC.

Quality metrics for high order meshes: analysis of the mechanical simulation of the heart beat.

The quality of a computational mesh is an important characteristic for stable and accurate simulations. Quality depends on the regularity of the initial mesh, and in mechanical simulations it evolves in time, with deformations causing changes in volume and distortion of mesh elements. Mesh quality metrics are therefore relevant for both mesh personalization and the monitoring of the simulation process. This work evaluates the significance, in meshes with high order interpolation, of four quality metrics described in the literature, applying them to analyse the stability of the simulation of the heart beat. It also investigates how image registration and mesh warping parameters affect the quality and stability of meshes. Jacobian-based metrics outperformed or matched the results of coarse geometrical metrics of aspect ratio or orthogonality, although they are more expensive computationally. The stability of simulations of a complete heart cycle was best predicted with a specificity of 61%, sensitivity of 85%, and only nominal differences were found changing the intra-element and per-element combination of quality values. A compromise between fitting accuracy and mesh stability and quality was found. Generic geometrical quality metrics have a limited success predicting stability, and an analysis of the simulation problem may be required for an optimal definition of quality.

Verification of cardiac tissue electrophysiology simulators using an N-version benchmark.

Ongoing developments in cardiac modelling have resulted, in particular, in the development of advanced and increasingly complex computational frameworks for simulating cardiac tissue electrophysiology. The goal of these simulations is often to represent the detailed physiology and pathologies of the heart using codes that exploit the computational potential of high-performance computing architectures. These developments have rapidly progressed the simulation capacity of cardiac virtual physiological human style models; however, they have also made it increasingly challenging to verify that a given code provides a faithful representation of the purported governing equations and corresponding solution techniques. This study provides the first cardiac tissue electrophysiology simulation benchmark to allow these codes to be verified. The benchmark was successfully evaluated on 11 simulation platforms to generate a consensus gold-standard converged solution. The benchmark definition in combination with the gold-standard solution can now be used to verify new simulation codes and numerical methods in the future.

An accurate, fast and robust method to generate patient-specific cubic Hermite meshes.

In-silico continuum simulations of organ and tissue scale physiology often require a discretisation or mesh of the solution domain. Cubic Hermite meshes provide a smooth representation of anatomy that is well-suited for simulating large deformation mechanics. Models of organ mechanics and deformation have demonstrated significant potential for clinical application. However, the production of a personalised mesh from patient's anatomy using medical images remains a major bottleneck in simulation workflows. To address this issue, we have developed an accurate, fast and automatic method for deriving patient-specific cubic Hermite meshes. The proposed solution customises a predefined template with a fast binary image registration step and a novel cubic Hermite mesh warping constructed using a variational technique. Image registration is used to retrieve the mapping field between the template mesh and the patient images. The variational warping technique then finds a smooth and accurate projection of this field into the basis functions of the mesh. Applying this methodology, cubic Hermite meshes are fitted to the binary description of shape with sub-voxel accuracy and within a few minutes, which is a significant advance over the existing state of the art methods. To demonstrate its clinical utility, a generic cubic Hermite heart biventricular model is personalised to the anatomy of four patients, and the resulting mechanical stability of these customised meshes is successfully demonstrated.

Comparative study of different discretizations of the phi(4) model.

We examine various recently proposed translationally invariant discretizations of the well-known phi(4) field theory. We compare and contrast the properties of their fundamental solutions including the nature of their kink-type solitary waves and the spectral properties of the linearization around such waves. We study these features as a function of the lattice spacing h , as one deviates from the continuum limit of h --> 0. We then proceed to a more "stringent" comparison of the models, by discussing the scattering properties of a kink-antikink pair for the different discretizations. These collisions are well known to possess properties that quite sensitively depend on the initial speed even at the continuum limit. We examine how typical model behaviors are modified in the presence (and as a function) of discreteness. One of the surprising trends that we observe is the increasing elasticity of kink collisions with deviation from the continuum limit. Another general feature is that the most inelastic kink collisions are observed in the classical discrete phi(4) model, while they are more elastic in the four studied translationally invariant models.