A new high-order stable numerical method for matrix inversion.
Ontology highlight
ABSTRACT:
SUBMITTER: Haghani FK
PROVIDER: S-EPMC3933421 | biostudies-other | 2014
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We develop a finite element scheme to approximate the dynamics of two and three dimensional fluidic membranes in Navier-Stokes flow. Local inextensibility of the membrane is ensured by solving a tangential Navier-Stokes equation, taking surface viscosity effects of Boussinesq-Scriven type into account. In our approach the bulk and surface degrees of freedom are discretized independently, which leads to an unfitted finite element approximation of the underlying free boundary problem. Bending elastic forces resulting from an elastic membrane energy are discretized using an approximation introduced by Dziuk (Numer Math 111:55-80, 2008). The obtained numerical scheme can be shown to be stable and to have good mesh properties. Finally, the evolution of membrane shapes is studied numerically in different flow situations in two and three space dimensions. The numerical results demonstrate the robustness of the method, and it is observed that the conservation properties are fulfilled to a high precision.
| S-EPMC5444514 | biostudies-literature
Project description:Genomics has revolutionized biology, enabling the interrogation of whole transcriptomes, genome-wide binding sites for proteins, and many other molecular processes. However, individual genomic assays measure elements that interact in vivo as components of larger molecular machines. Understanding how these high-order interactions drive gene expression presents a substantial statistical challenge. Building on random forests (RFs) and random intersection trees (RITs) and through extensive, biologically inspired simulations, we developed the iterative random forest algorithm (iRF). iRF trains a feature-weighted ensemble of decision trees to detect stable, high-order interactions with the same order of computational cost as the RF. We demonstrate the utility of iRF for high-order interaction discovery in two prediction problems: enhancer activity in the early Drosophila embryo and alternative splicing of primary transcripts in human-derived cell lines. In Drosophila, among the 20 pairwise transcription factor interactions iRF identifies as stable (returned in more than half of bootstrap replicates), 80% have been previously reported as physical interactions. Moreover, third-order interactions, e.g., between Zelda (Zld), Giant (Gt), and Twist (Twi), suggest high-order relationships that are candidates for follow-up experiments. In human-derived cells, iRF rediscovered a central role of H3K36me3 in chromatin-mediated splicing regulation and identified interesting fifth- and sixth-order interactions, indicative of multivalent nucleosomes with specific roles in splicing regulation. By decoupling the order of interactions from the computational cost of identification, iRF opens additional avenues of inquiry into the molecular mechanisms underlying genome biology.
| S-EPMC5828575 | biostudies-literature
Project description:A new iterative scheme has been constructed for finding minimal solution of a rational matrix equation of the form X + A*X (-1) A = I. The new method is inversion-free per computing step. The convergence of the method has been studied and tested via numerical experiments.
| S-EPMC4151580 | biostudies-other
Project description:A high-order finite difference scheme is proposed for solving time fractional heat equations. The time fractional derivative is described in the Riemann-Liouville sense. In the proposed scheme a new second-order discretization, which is based on Crank-Nicholson method, is applied for the time fractional part and fourth-order accuracy compact approximation is applied for the second-order space derivative. The spectral stability and the Fourier stability analysis of the difference scheme are shown. Finally a detailed numerical analysis, including tables, figures, and error comparison, is given to demonstrate the theoretical results and high accuracy of the proposed scheme.
| S-EPMC3947870 | biostudies-other
Project description: Not available
| S-EPMC6492311 | biostudies-literature
Project description:Epistatic interactions are referred to as SNPs (single nucleotide polymorphisms) that affect disease development and trait expression nonlinearly, and hence identifying epistatic interactions plays a great role in explaining the pathogenesis and genetic heterogeneity of complex diseases. Many methods have been proposed for epistasis detection; nevertheless, they mainly focus on low-order epistatic interactions, two-order or three-order for instance, and often ignore high-order interactions due to computational burden. In this paper, a module detection method called MDSN is proposed for identifying high-order epistatic interactions. First, an SNP network is constructed by a construction strategy of interaction complementary, which consists of low-order SNP interactions that can be obtained from fast computations. Then, a node evaluation measure that integrates multi-topological features is proposed to improve the node expansion algorithm, where the importance of a node is comprehensively evaluated by the topological characteristics of the neighborhood. Finally, modules are detected in the constructed SNP network, which have high-order epistatic interactions associated with the disease. The MDSN was compared with four state-of-the-art methods on simulation datasets and a real Age-related Macular Degeneration dataset. The results demonstrate that MDSN has higher performance on detecting high-order interactions.
| S-EPMC9778340 | biostudies-literature
Project description:Infectious diseases pose a serious threat to human life, the Genome Wide Association Studies (GWAS) can analyze susceptibility genes of infectious diseases from the genetic level and carry out targeted prevention and treatment. The susceptibility genes for infectious diseases often act in combination with multiple susceptibility sites; therefore, high-order epistasis detection has become an important means. However, due to intensive computational burden and diversity of disease models, existing methods have drawbacks on low detection power, high computation cost, and preference for some types of disease models. Furthermore, these methods are exposed to repeated query and model inversion attacks in the process of iterative optimization, which may disclose Single Nucleotide Polymorphism (SNP) information associated with individual privacy. Therefore, in order to solve these problems, this paper proposed a safe harmony search algorithm for high-order gene interaction detection, termed as HS-DP. Firstly, the linear weighting method was used to integrate 5 objective functions to screen out high-order SNP sets with high correlation, including K2-Score, JS divergence, logistic regression, mutual information, and Gini. Then, based on the Differential Privacy (DP) theory, the function disturbance mechanism was introduced to protect the security of individual privacy information associated with the objective function, and we proved the rationality of the disturbance mechanism theoretically. Finally, the practicability and superiority of the algorithm were verified by experiments. Experimental results showed that the algorithm proposed in this paper could improve the detection accuracy to the greatest extent while guaranteeing privacy.
| S-EPMC9050263 | biostudies-literature
Project description:As nanoscale photonic devices are densely integrated, multiple near-field optical eigenmodes take part in their functionalization. Inevitably, these eigenmodes are highly multiplexed in their spectra and superposed in their spatial distributions, making it extremely difficult for conventional near-field scanning optical microscopy (NSOM) to address individual eigenmodes. Here, we develop a near-field transmission matrix microscopy for mapping the high-order eigenmodes of nanostructures, which are invisible with conventional NSOM. At an excitation wavelength where multiple modes are superposed, we measure the near-field amplitude and phase maps for various far-field illumination angles, from which we construct a fully phase-referenced far- to near-field transmission matrix. By performing the singular value decomposition, we extract orthogonal near-field eigenmodes such as anti-symmetric mode and quadruple mode of multiple nano-slits whose gap size (50?nm) is smaller than the probe aperture (150?nm). Analytic model and numerical mode analysis validated the experimentally observed modes.
| S-EPMC7244505 | biostudies-literature
Project description:A new synthetic method has been developed to prepare fluorescent silica nanoparticles without employing isothiocyanated dye molecules and (3-aminopropyl)triethoxysilane (APS) for the thiourea linkage formation; the resulting fluorescent silica nanoparticles show excellent photochemical, thermal and pH stabilities and a good biocompatibility with over 85% viability from various cell types.
| S-EPMC2879016 | biostudies-literature
Project description:Stable and conservative numerical boundary schemes, for both compact and explicit (central) finite differences require a number of parameters that must be tuned for stability. Values of these coefficients for 4th, 6th, and 8th boundary schemes are given in this article. The stability of the schemes is demonstrated through a series of numerical tests in "High-Order, Stable, and Conservative Boundary Schemes for Central and Compact Finite Differences" Brady and Livescu, 2019. These tests include: a neutrally stable constant coefficient hyperbolic system, a two-dimensional varying coefficient hyperbolic scalar equation and, examining the transport of an inviscid vortex using the compressible Euler equations. The error norms for the variety of tests associated with different the schemes for different grid resolutions and time-step constraints are given in the accompanying databases.
| S-EPMC6595279 | biostudies-literature