BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Stefan Engblom

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

ResearcherUnclaimed source record
math.NADistributed, Parallel, and Cluster ComputingNumerical AnalysisData Structures and Algorithmsmath.PRMolecular NetworksBiological PhysicsBiomoleculeseess.SYMachine Learningmath.APmath.DSMathematical SoftwareNeurons and Cognitionphysics.flu-dynPopulations and EvolutionQuantitative MethodsSystems and Control

Catalog footprint

What is connected

13works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

13 published item(s)

preprint2020arXiv

Distributed and Adaptive Fast Multipole Method In Three Dimensions

We develop a general distributed implementation of an adaptive fast multipole method in three space dimensions. We rely on a balanced type of adaptive space discretisation which supports a highly transparent and fully distributed implementation. A complexity analysis indicates favorable scaling properties and numerical experiments on up to 512 cores and 1 billion source points verify them. The parameters controlling the algorithm are subject to in-depth experiments and the performance response to the input parameters implies that the overall implementation is well-suited to automated tuning.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Initialization of a Disease Transmission Model

Approaches to the calculation of the full state vector of a larger epidemiological model for the spread of COVID-19 in Sweden at the initial time instant from available data and with a simplified dynamical model are proposed and evaluated. The larger epidemiological model is based on a continuous Markov chain and captures the demographic composition of and the transport flows between the counties of Sweden. Its intended use is to predict the outbreak development in temporal and spatial coordinates as well as across the demographic groups. It can also support evaluating and comparing of prospective intervention strategies in terms of e.g. lockdown in certain areas or isolation of specific age groups. The simplified model is a discrete time-invariant linear system that has cumulative infectious incidence, infected population, asymptomatic population, exposed population, and infectious pressure as the state variables. Since the system matrix of the model depends on a number transition rates, structural properties of the model are investigated for suitable parameter ranges. It is concluded that the model becomes unobservable for some parameter values. Two contrasting approaches to the initial state estimation are considered. One is a version of Rauch-Tung-Striebel smoother and another is based on solving a batch nonlinear optimization problem. The benefits and shortcomings of the considered estimation techniques are analyzed and compared on synthetic data for several Swedish counties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast Matlab compatible sparse assembly on multicore computers

We develop and implement in this paper a fast sparse assembly algorithm, the fundamental operation which creates a compressed matrix from raw index data. Since it is often a quite demanding and sometimes critical operation, it is of interest to design a highly efficient implementation. We show how to do this, and moreover, we show how our implementation can be parallelized to utilize the power of modern multicore computers. Our freely available code, fully Matlab compatible, achieves about a factor of 5 times in speedup on a typical 6-core machine and 10 times on a dual-socket 16 core machine compared to the built-in serial implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic focusing coupled with negative feedback enables robust regulation in biochemical reaction networks

Nature presents multiple intriguing examples of processes which proceed at high precision and regularity. This remarkable stability is frequently counter to modelers' experience with the inherent stochasticity of chemical reactions in the regime of low copy numbers. Moreover, the effects of noise and nonlinearities can lead to "counter-intuitive" behavior, as demonstrated for a basic enzymatic reaction scheme that can display stochastic focusing (SF). Under the assumption of rapid signal fluctuations, SF has been shown to convert a graded response into a threshold mechanism, thus attenuating the detrimental effects of signal noise. However, when the rapid fluctuation assumption is violated, this gain in sensitivity is generally obtained at the cost of very large product variance, and this unpredictable behavior may be one possible explanation of why, more than a decade after its introduction, SF has still not been observed in real biochemical systems. In this work we explore the noise properties of a simple enzymatic reaction mechanism with a small and fluctuating number of active enzymes that behaves as a high-gain, noisy amplifier due to SF caused by slow enzyme fluctuations. We then show that the inclusion of a plausible negative feedback mechanism turns the system from a noisy signal detector to a strong homeostatic mechanism by exchanging high gain with strong attenuation in output noise and robustness to parameter variations. Moreover, we observe that the discrepancy between deterministic and stochastic descriptions of stochastically focused systems in the evolution of the means almost completely disappears, despite very low molecule counts and the additional nonlinearity due to feedback. The reaction mechanism considered here can provide a possible resolution to the apparent conflict between intrinsic noise and high precision in critical intracellular processes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong convergence for split-step methods in stochastic jump kinetics

Mesoscopic models in the reaction-diffusion framework have gained recognition as a viable approach to describing chemical processes in cell biology. The resulting computational problem is a continuous-time Markov chain on a discrete and typically very large state space. Due to the many temporal and spatial scales involved many different types of computationally more effective multiscale models have been proposed, typically coupling different types of descriptions within the Markov chain framework. In this work we look at the strong convergence properties of the basic first order Strang, or Lie-Trotter, split-step method, which is formed by decoupling the dynamics in finite time-steps. Thanks to its simplicity and flexibility, this approach has been tried in many different combinations. We develop explicit sufficient conditions for path-wise well-posedness and convergence of the method, including error estimates, and we illustrate our findings with numerical examples. In doing so, we also suggest a certain partition of unity representation for the split-step method, which in turn implies a concrete simulation algorithm under which trajectories may be compared in a path-wise sense.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Approximations for the Moments of Nonstationary and State Dependent Birth-Death Queues

In this paper we propose a new method for approximating the nonstationary moment dynamics of one dimensional Markovian birth-death processes. By expanding the transition probabilities of the Markov process in terms of Poisson-Charlier polynomials, we are able to estimate any moment of the Markov process even though the system of moment equations may not be closed. Using new weighted discrete Sobolev spaces, we derive explicit error bounds of the transition probabilities and new weak a priori estimates for approximating the moments of the Markov processs using a truncated form of the expansion. Using our error bounds and estimates, we are able to show that our approximations converge to the true stochastic process as we add more terms to the expansion and give explicit bounds on the truncation error. As a result, we are the first paper in the queueing literature to provide error bounds and estimates on the performance of a moment closure approximation. Lastly, we perform several numerical experiments for some important models in the queueing theory literature and show that our expansion techniques are accurate at estimating the moment dynamics of these Markov process with only a few terms of the expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic autotuning of adaptive fast multipole methods on hybrid multicore CPU & GPU systems

We discuss an implementation of adaptive fast multipole methods targeting hybrid multicore CPU- and GPU-systems. From previous experiences with the computational profile of our version of the fast multipole algorithm, suitable parts are off-loaded to the GPU, while the remaining parts are threaded and executed concurrently by the CPU. The parameters defining the algorithm affects the performance and by measuring this effect we are able to dynamically balance the algorithm towards optimal performance. Our setup uses the dynamic nature of the computations and is therefore of general character.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Machine learning for ultrafast X-ray diffraction patterns on large-scale GPU clusters

The classical method of determining the atomic structure of complex molecules by analyzing diffraction patterns is currently undergoing drastic developments. Modern techniques for producing extremely bright and coherent X-ray lasers allow a beam of streaming particles to be intercepted and hit by an ultrashort high energy X-ray beam. Through machine learning methods the data thus collected can be transformed into a three-dimensional volumetric intensity map of the particle itself. The computational complexity associated with this problem is very high such that clusters of data parallel accelerators are required. We have implemented a distributed and highly efficient algorithm for inversion of large collections of diffraction patterns targeting clusters of hundreds of GPUs. With the expected enormous amount of diffraction data to be produced in the foreseeable future, this is the required scale to approach real time processing of data at the beam site. Using both real and synthetic data we look at the scaling properties of the application and discuss the overall computational viability of this exciting and novel imaging technique.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the stability of stochastic jump kinetics

Motivated by the lack of a suitable constructive framework for analyzing popular stochastic models of Systems Biology, we devise conditions for existence and uniqueness of solutions to certain jump stochastic differential equations (SDEs). Working from simple examples we find reasonable and explicit assumptions on the driving coefficients for the SDE representation to make sense. By `reasonable' we mean that stronger assumptions generally do not hold for systems of practical interest. In particular, we argue against the traditional use of global Lipschitz conditions and certain common growth restrictions. By `explicit', finally, we like to highlight the fact that the various constants occurring among our assumptions all can be determined once the model is fixed. We show how basic long time estimates and some limit results for perturbations can be derived in this setting such that these can be contrasted with the corresponding estimates from deterministic dynamics. The main complication is that the natural path-wise representation is generated by a counting measure with an intensity that depends nonlinearly on the state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite element analysis of neuronal electric fields: the effect of heterogeneous resistivity

Simulation of extracellular fields is one of the substantial methods used in the area of computational neuroscience. Its most common usage is validation of experimental methods as EEG and extracellular spike recordings or modeling of physiological phenomena which can not be easily determined empirically. Continuous experimental work has been re-raising the importance of polarization effects between neuronal structures to neuronal communication. As this effects relies on very small potential changes, better modeling methods are necessary to quantify the weak electrical fields in the microscopic scale in a more realistic way. An important factor of influence on local field effects in the hippocampal formation is the heterogeneous resistivity of extracellular tissue. The vast majority of modeling studies consider the extracellular space to be homogeneous while experimentally, it has been shown that the stratum pyramidale has two times higher resistivity then other hippocampal layers. Common simulation methods for extracellular electrical fields based on the point source approximation are bound to describe the resistance of the space with a single, linear factor. We propose that models should be based on the space- and time-dependent Maxwell equations in order to account for heterogeneous properties of the extracellular space and specific arrangements of neurons in dense hippocampal layers. To demonstrate the influence of heterogeneous extracellular resistivity and neuronal spatial orientation on modeling results, we combine solutions of classical compartment models with spatiotemporal PDEs solved by the FEM. With the help of these methods, we show that the inclusion of heterogeneous resistivity has a substantial impact on voltages in close proximity to emitting neurons, increasing the extracellular potentials substantially compared to the homogeneous variant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive fast multipole methods on the GPU

We present a highly general implementation of fast multipole methods on graphics processing units (GPUs). Our two-dimensional double precision code features an asymmetric type of adaptive space discretization leading to a particularly elegant and flexible implementation. All steps of the multipole algorithm are efficiently performed on the GPU, including the initial phase which assembles the topological information of the input data. Through careful timing experiments we investigate the effects of the various peculiarities of the GPU architecture.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On diffuse interface modeling and simulation of surfactants in two-phase fluid flow

An existing phase-field model of two immiscible fluids with a single soluble surfactant present is discussed in detail. We analyze the well-posedness of the model and provide strong evidence that it is mathematically ill-posed for a large set of physically relevant parameters. As a consequence, critical modifications to the model are suggested that substantially increase the domain of validity. Carefully designed numerical simulations offer informative demonstrations as to the sharpness of our theoretical results and the qualities of the physical model. A fully coupled hydrodynamic test-case demonstrates the potential to capture also non-trivial effects on the overall flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On well-separated sets and fast multipole methods

The notion of well-separated sets is crucial in fast multipole methods as the main idea is to approximate the interaction between such sets via cluster expansions. We revisit the one-parameter multipole acceptance criterion in a general setting and derive a relative error estimate. This analysis benefits asymmetric versions of the method, where the division of the multipole boxes is more liberal than in conventional codes. Such variants offer a particularly elegant implementation with a balanced multipole tree, a feature which might be very favorable on modern computer architectures.

0 citations0 reviewsNo code linkedOpen access