Source author record

Tomas Vejchodsky

Tomas Vejchodsky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Molecular Networks Quantitative Methods math.DS Subcellular Processes

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Two-sided bounds of eigenvalues - local efficiency and convergence of adaptive algorithm

We generalize and analyse the method for computing lower bounds of the principal eigenvalue proposed in our previous paper (I. Sebestova, T. Vejchodsky, SIAM J. Numer. Anal. 2014). This method is suitable for symmetric elliptic eigenvalue problems with mixed boundary conditions of Dirichlet, Neumann, and Robin type and it is based on a posteriori error analysis using flux reconstructions. We improve the original result in several aspects. We show how to obtain lower bounds even for higher eigenvalues. We present a local approach for the flux reconstruction enabling efficient implementation. We prove the equivalence of the resulting estimator with the classical residual estimator and consequently its local efficiency. We also prove the convergence of the corresponding adaptive algorithm. Finally, we illustrate the practical performance of the method by numerical examples.

preprint2015arXiv

Chemical Reaction Systems with a Homoclinic Bifurcation: an Inverse Problem

An inverse problem framework for constructing reaction systems with prescribed properties is presented. Kinetic transformations are defined and analysed as a part of the framework, allowing an arbitrary polynomial ordinary differential equation to be mapped to the one that can be represented as a reaction network. The framework is used for construction of specific two- and three-dimensional bistable reaction systems undergoing a supercritical homoclinic bifurcation, and the topology of their phase spaces is discussed.

preprint2015arXiv

Parameter estimation and bifurcation analysis of stochastic models of gene regulatory networks: tensor-structured methods

Stochastic modelling provides an indispensable tool for understanding how random events at the molecular level influence cellular functions. In practice, the common challenge is to calibrate a large number of model parameters against the experimental data. A related problem is to efficiently study how the behaviour of a stochastic model depends on its parameters, i.e. whether a change in model parameters can lead to a significant qualitative change in model behaviour (bifurcation). In this paper, tensor-structured parametric analysis (TPA) is presented. It is based on recently proposed low-parametric tensor-structured representations of classical matrices and vectors. This approach enables simultaneous computation of the model properties for all parameter values within a parameter space. This methodology is exemplified to study the parameter estimation, robustness, sensitivity and bifurcation structure in stochastic models of biochemical networks. The TPA has been implemented in Matlab and the codes are available at http://www.stobifan.org .

preprint2014arXiv

Noise-induced multistability in chemical systems: Discrete vs Continuum modeling

The noisy dynamics of chemical systems is commonly studied using either the chemical master equation (CME) or the chemical Fokker-Planck equation (CFPE). The latter is a continuum approximation of the discrete CME approach. We here show that the CFPE may fail to capture the CME's prediction of noise-induced multistability. In particular we find a simple chemical system for which the CME's marginal probability distribution changes from unimodal to multimodal as the system-size decreases below a critical value, while the CFPE's marginal probability distribution is unimodal for all physically meaningful system sizes.

preprint2011arXiv

Adaptive Finite Element Method Assisted by Stochastic Simulation of Chemical Systems

Stochastic models of chemical systems are often analysed by solving the corresponding Fokker-Planck equation which is a drift-diffusion partial differential equation for the probability distribution function. Efficient numerical solution of the Fokker-Planck equation requires adaptive mesh refinements. In this paper, we present a mesh refinement approach which makes use of a stochastic simulation of the underlying chemical system. By observing the stochastic trajectory for a relatively short amount of time, the areas of the state space with non-negligible probability density are identified. By refining the finite element mesh in these areas, and coarsening elsewhere, a suitable mesh is constructed and used for the computation of the probability density.