1
|
Moreno JD, Silva JR. Emerging methods to model cardiac ion channel and myocyte electrophysiology. BIOPHYSICS REVIEWS 2023; 4:011315. [PMID: 37034130 PMCID: PMC10071990 DOI: 10.1063/5.0127713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 02/28/2023] [Indexed: 04/03/2023]
Abstract
In the field of cardiac electrophysiology, modeling has played a central role for many decades. However, even though the effort is well-established, it has recently seen a rapid and sustained evolution in the complexity and predictive power of the models being created. In particular, new approaches to modeling have allowed the tracking of parallel and interconnected processes that span from the nanometers and femtoseconds that determine ion channel gating to the centimeters and minutes needed to describe an arrhythmia. The connection between scales has brought unprecedented insight into cardiac arrhythmia mechanisms and drug therapies. This review focuses on the generation of these models from first principles, generation of detailed models to describe ion channel kinetics, algorithms to create and numerically solve kinetic models, and new approaches toward data gathering that parameterize these models. While we focus on application of these models for cardiac arrhythmia, these concepts are widely applicable to model the physiology and pathophysiology of any excitable cell.
Collapse
Affiliation(s)
- Jonathan D. Moreno
- Division of Cardiology, Department of Medicine, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| | - Jonathan R. Silva
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, Missouri 63130, USA
| |
Collapse
|
2
|
Laudenschlager S, Cai XC. An inner-outer subcycling algorithm for parallel cardiac electrophysiology simulations. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3677. [PMID: 36573938 DOI: 10.1002/cnm.3677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Revised: 11/16/2022] [Accepted: 12/21/2022] [Indexed: 06/17/2023]
Abstract
This paper explores cardiac electrophysiological simulations of the monodomain equations and introduces a novel subcycling time integration algorithm to exploit the structure of the ionic model. The aim of this work is to improve upon the efficiency of parallel cardiac monodomain simulations by using our subcycling algorithm in the computation of the ionic model to handle the local sharp changes of the solution. This will reduce the turnaround time for the simulation of basic cardiac electrical function on both idealized and patient-specific geometry. Numerical experiments show that the proposed approach is accurate and also has close to linear parallel scalability on a computer with more than 1000 processor cores. Ultimately, the reduction in simulation time can be beneficial in clinical applications, where multiple simulations are often required to tune a model to match clinical measurements.
Collapse
Affiliation(s)
| | - Xiao-Chuan Cai
- Department of Mathematics, University of Macau, Macau, China
| |
Collapse
|
3
|
Ogiermann D, Perotti LE, Balzani D. A simple and efficient adaptive time stepping technique for low-order operator splitting schemes applied to cardiac electrophysiology. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3670. [PMID: 36510350 DOI: 10.1002/cnm.3670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 10/25/2022] [Accepted: 12/04/2022] [Indexed: 06/17/2023]
Abstract
We present a simple, yet efficient adaptive time stepping scheme for cardiac electrophysiology (EP) simulations based on standard operator splitting techniques. The general idea is to exploit the relation between the splitting error and the reaction's magnitude-found in a previous one-dimensional analytical study by Spiteri and Ziaratgahi-to construct the new time step controller for three-dimensional problems. Accordingly, we propose to control the time step length of the operator splitting scheme as a function of the reaction magnitude, in addition to the common approach of adapting the reaction time step. This conforms with observations in numerical experiments supporting the need for a significantly smaller time step length during depolarization than during repolarization. The proposed scheme is compared with classical proportional-integral-differential controllers using state-of-the-art error estimators, which are also presented in details as they have not been previously applied in the context of cardiac EP with operator splitting techniques. Benchmarks show that choosing the time step as a sigmoidal function of the reaction magnitude is highly efficient and full cardiac cycles can be computed with precision even in a realistic biventricular setup. The proposed scheme outperforms common adaptive time stepping techniques, while depending on fewer tuning parameters.
Collapse
Affiliation(s)
- Dennis Ogiermann
- Chair of Continuum Mechanics, Ruhr University Bochum, Bochum, Germany
| | - Luigi E Perotti
- Mechanical and Aerospace Engineering Department, University of Central Florida, Orlando, Florida, USA
| | - Daniel Balzani
- Chair of Continuum Mechanics, Ruhr University Bochum, Bochum, Germany
| |
Collapse
|
4
|
Lindner LP, Gerach T, Jahnke T, Loewe A, Weiss D, Wieners C. Efficient time splitting schemes for the monodomain equation in cardiac electrophysiology. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3666. [PMID: 36562492 DOI: 10.1002/cnm.3666] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Revised: 10/19/2022] [Accepted: 11/26/2022] [Indexed: 06/17/2023]
Abstract
Approximating the fast dynamics of depolarization waves in the human heart described by the monodomain model is numerically challenging. Splitting methods for the PDE-ODE coupling enable the computation with very fine space and time discretizations. Here, we compare different splitting approaches regarding convergence, accuracy, and efficiency. Simulations were performed for a benchmark problem with the Beeler-Reuter cell model on a truncated ellipsoid approximating the left ventricle including a localized stimulation. For this configuration, we provide a reference solution for the transmembrane potential. We found a semi-implicit approach with state variable interpolation to be the most efficient scheme. The results are transferred to a more physiological setup using a bi-ventricular domain with a complex external stimulation pattern to evaluate the accuracy of the activation time for different resolutions in space and time.
Collapse
Affiliation(s)
- Laura P Lindner
- Institute of Applied and Numerical Mathematics, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Tobias Gerach
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Tobias Jahnke
- Institute of Applied and Numerical Mathematics, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Axel Loewe
- Institute of Biomedical Engineering, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Daniel Weiss
- Institute of Applied and Numerical Mathematics, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Christian Wieners
- Institute of Applied and Numerical Mathematics, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
5
|
Woodworth LA, Cansız B, Kaliske M. A numerical study on the effects of spatial and temporal discretization in cardiac electrophysiology. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3443. [PMID: 33522111 DOI: 10.1002/cnm.3443] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Revised: 12/16/2020] [Accepted: 01/23/2021] [Indexed: 06/12/2023]
Abstract
Millions of degrees of freedom are often required to accurately represent the electrophysiology of the myocardium due to the presence of discretization effects. This study seeks to explore the influence of temporal and spatial discretization on the simulation of cardiac electrophysiology in conjunction with changes in modeling choices. Several finite element analyses are performed to examine how discretization affects solution time, conduction velocity and electrical excitation. Discretization effects are considered along with changes in the electrophysiology model and solution approach. Two action potential models are considered: the Aliev-Panfilov model and the ten Tusscher-Noble-Noble-Panfilov model. The solution approaches consist of two time integration schemes and different treatments for solving the local system of ordinary differential equations. The efficiency and stability of the calculation approaches are demonstrated to be dependent on the action potential model. The dependency of the conduction velocity on the element size and time step is shown to be different for changes in material parameters. Finally, the discrepancies between the wave propagation in coarse and fine meshes are analyzed based on the temporal evolution of the transmembrane potential at a node and its neighboring Gauss points. Insight obtained from this study can be used to suggest new methods to improve the efficiency of simulations in cardiac electrophysiology.
Collapse
Affiliation(s)
- Lucas A Woodworth
- Institute for Structural Analysis, Technische Universität Dresden, Dresden, Germany
| | - Barış Cansız
- Institute for Structural Analysis, Technische Universität Dresden, Dresden, Germany
| | - Michael Kaliske
- Institute for Structural Analysis, Technische Universität Dresden, Dresden, Germany
| |
Collapse
|
6
|
Combination of "generalized Trotter operator splitting" and "quadratic adaptive algorithm" method for tradeoff among speedup, stability, and accuracy in the Markov chain model of sodium ion channels in the ventricular cell model. Med Biol Eng Comput 2020; 58:2131-2141. [PMID: 32676840 DOI: 10.1007/s11517-020-02220-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Accepted: 06/25/2020] [Indexed: 10/23/2022]
Abstract
The fast hybrid operator splitting (HOS) and stable uniformization (UNI) methods have been proposed to save computation cost and enhance stability for Markov chain model in cardiac cell simulations. Moreover, Chen-Chen-Luo's quadratic adaptive algorithm (CCL) combined with HOS or UNI was used to improve the tradeoff between speedup and stability, but without considering accuracy. To compromise among stability, acceleration, and accuracy, we propose a generalized Trotter operator splitting (GTOS) method combined with CCL independent of the asymptotic property of a particular ion-channel model. Due to the accuracy underestimation of the mixed root mean square error (MRMSE) method, threshold root mean square error (TRMSE) is proposed to evaluate computation accuracy. With the fixed time-step RK4 as a reference, the second-order GTOS combined with CCL (30.8-fold speedup) for the wild-type Markov chain model with nine states (WT-9 model) or (7.4-fold) for the wild-type Markov chain model with eight states (WT-8 model) is faster than UNI combined with CCL (15.6-fold) for WT-9 model or (1.2-fold) for WT-8 model, separately. Besides, the second-order GTOS combined with CCL has 3.81% TRMSE for WT-9 model or 4.32% TRMSE for WT-8 model more accurate than 72.43% TRMSE for WT-9 model or 136.17% TRMSE for WT-8 model of HOS combined with CCL. To compromise speedup and accuracy, low-order GTOS combined with CCL is suggested to have the advantages of high precision and low computation cost. For high-accuracy requirements, high-order GTOS combined with CCL is recommended. Graphical abstract.
Collapse
|
7
|
Combination of “quadratic adaptive algorithm” and “hybrid operator splitting” or uniformization algorithms for stability against acceleration in the Markov model of sodium ion channels in the ventricular cell model. Med Biol Eng Comput 2019; 57:1367-1379. [DOI: 10.1007/s11517-019-01956-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2018] [Accepted: 01/28/2019] [Indexed: 10/27/2022]
|
8
|
Sachetto Oliveira R, Martins Rocha B, Burgarelli D, Meira W, Constantinides C, Weber Dos Santos R. Performance evaluation of GPU parallelization, space-time adaptive algorithms, and their combination for simulating cardiac electrophysiology. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e2913. [PMID: 28636811 DOI: 10.1002/cnm.2913] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/10/2016] [Revised: 06/09/2017] [Accepted: 06/16/2017] [Indexed: 05/23/2023]
Abstract
The use of computer models as a tool for the study and understanding of the complex phenomena of cardiac electrophysiology has attained increased importance nowadays. At the same time, the increased complexity of the biophysical processes translates into complex computational and mathematical models. To speed up cardiac simulations and to allow more precise and realistic uses, 2 different techniques have been traditionally exploited: parallel computing and sophisticated numerical methods. In this work, we combine a modern parallel computing technique based on multicore and graphics processing units (GPUs) and a sophisticated numerical method based on a new space-time adaptive algorithm. We evaluate each technique alone and in different combinations: multicore and GPU, multicore and GPU and space adaptivity, multicore and GPU and space adaptivity and time adaptivity. All the techniques and combinations were evaluated under different scenarios: 3D simulations on slabs, 3D simulations on a ventricular mouse mesh, ie, complex geometry, sinus-rhythm, and arrhythmic conditions. Our results suggest that multicore and GPU accelerate the simulations by an approximate factor of 33×, whereas the speedups attained by the space-time adaptive algorithms were approximately 48. Nevertheless, by combining all the techniques, we obtained speedups that ranged between 165 and 498. The tested methods were able to reduce the execution time of a simulation by more than 498× for a complex cellular model in a slab geometry and by 165× in a realistic heart geometry simulating spiral waves. The proposed methods will allow faster and more realistic simulations in a feasible time with no significant loss of accuracy.
Collapse
Affiliation(s)
- Rafael Sachetto Oliveira
- Departamento de Ciência da Computação, Universidade Federal de São João de Rei, São João del-rei MG, Brazil
| | - Bernardo Martins Rocha
- Departamento de Ciência da Computação e Programa em Modelagem Computacional, Universidade Federal de Juiz de Fora, Juiz de Fora, MG, Brazil
| | - Denise Burgarelli
- Departamento de Matemática, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
| | - Wagner Meira
- Departamento de Ciência da Computação, Universidade Federal de Minas Gerais, Belo Horizonte, MG, Brazil
| | | | - Rodrigo Weber Dos Santos
- Departamento de Ciência da Computação e Programa em Modelagem Computacional, Universidade Federal de Juiz de Fora, Juiz de Fora, MG, Brazil
| |
Collapse
|
9
|
Starý T, Biktashev VN. Fast-slow asymptotics for a Markov chain model of fast sodium current. CHAOS (WOODBURY, N.Y.) 2017; 27:093937. [PMID: 28964114 DOI: 10.1063/1.5003013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We explore the feasibility of using fast-slow asymptotics to eliminate the computational stiffness of discrete-state, continuous-time deterministic Markov chain models of ionic channels underlying cardiac excitability. We focus on a Markov chain model of fast sodium current, and investigate its asymptotic behaviour with respect to small parameters identified in different ways.
Collapse
Affiliation(s)
- Tomáš Starý
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Harrison Building, North Park Road, Exeter EX4 4QF, United Kingdom
| | - Vadim N Biktashev
- College of Engineering, Mathematics and Physical Sciences, University of Exeter, Harrison Building, North Park Road, Exeter EX4 4QF, United Kingdom
| |
Collapse
|
10
|
Myokit: A simple interface to cardiac cellular electrophysiology. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2015; 120:100-14. [PMID: 26721671 DOI: 10.1016/j.pbiomolbio.2015.12.008] [Citation(s) in RCA: 64] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2015] [Revised: 11/07/2015] [Accepted: 12/16/2015] [Indexed: 11/24/2022]
Abstract
Myokit is a new powerful and versatile software tool for modeling and simulation of cardiac cellular electrophysiology. Myokit consists of an easy-to-read modeling language, a graphical user interface, single and multi-cell simulation engines and a library of advanced analysis tools accessible through a Python interface. Models can be loaded from Myokit's native file format or imported from CellML. Model export is provided to C, MATLAB, CellML, CUDA and OpenCL. Patch-clamp data can be imported and used to estimate model parameters. In this paper, we review existing tools to simulate the cardiac cellular action potential to find that current tools do not cater specifically to model development and that there is a gap between easy-to-use but limited software and powerful tools that require strong programming skills from their users. We then describe Myokit's capabilities, focusing on its model description language, simulation engines and import/export facilities in detail. Using three examples, we show how Myokit can be used for clinically relevant investigations, multi-model testing and parameter estimation in Markov models, all with minimal programming effort from the user. This way, Myokit bridges a gap between performance, versatility and user-friendliness.
Collapse
|
11
|
Cooper J, Spiteri RJ, Mirams GR. Cellular cardiac electrophysiology modeling with Chaste and CellML. Front Physiol 2015; 5:511. [PMID: 25610400 PMCID: PMC4285015 DOI: 10.3389/fphys.2014.00511] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2014] [Accepted: 12/09/2014] [Indexed: 11/18/2022] Open
Abstract
Chaste is an open-source C++ library for computational biology that has well-developed cardiac electrophysiology tissue simulation support. In this paper, we introduce the features available for performing cardiac electrophysiology action potential simulations using a wide range of models from the Physiome repository. The mathematics of the models are described in CellML, with units for all quantities. The primary idea is that the model is defined in one place (the CellML file), and all model code is auto-generated at compile or run time; it never has to be manually edited. We use ontological annotation to identify model variables describing certain biological quantities (membrane voltage, capacitance, etc.) to allow us to import any relevant CellML models into the Chaste framework in consistent units and to interact with them via consistent interfaces. This approach provides a great deal of flexibility for analysing different models of the same system. Chaste provides a wide choice of numerical methods for solving the ordinary differential equations that describe the models. Fixed-timestep explicit and implicit solvers are provided, as discussed in previous work. Here we introduce the Rush–Larsen and Generalized Rush–Larsen integration techniques, made available via symbolic manipulation of the model equations, which are automatically rearranged into the forms required by these approaches. We have also integrated the CVODE solvers, a ‘gold standard’ for stiff systems, and we have developed support for symbolic computation of the Jacobian matrix, yielding further increases in the performance and accuracy of CVODE. We discuss some of the technical details of this work and compare the performance of the available numerical methods. Finally, we discuss how this is generalized in our functional curation framework, which uses a domain-specific language for defining complex experiments as a basis for comparison of model behavior.
Collapse
Affiliation(s)
- Jonathan Cooper
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| | - Raymond J Spiteri
- Numerical Simulation Research Lab, Department of Computer Science, University of Saskatchewan Saskatoon, SK, Canada
| | - Gary R Mirams
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| |
Collapse
|
12
|
Starý T, Biktashev VN. Exponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes. IEEE Trans Biomed Eng 2014; 62:1070-6. [PMID: 25376030 DOI: 10.1109/tbme.2014.2366466] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The modern Markov chain models of ionic channels in excitable membranes are numerically stiff. The popular numerical methods for these models require very small time steps to ensure stability. Our objective is to formulate and test two methods addressing this issue, so that the timestep can be chosen based on accuracy rather than stability. Both proposed methods extend Rush-Larsen technique, which was originally developed to Hogdkin-Huxley type gate models. One method, "matrix Rush-Larsen" (MRL) uses a matrix reformulation of the Rush-Larsen scheme, where the matrix exponentials are calculated using precomputed tables of eigenvalues and eigenvectors. The other, "hybrid operator splitting" (HOS) method exploits asymptotic properties of a particular Markov chain model, allowing explicit analytical expressions for the substeps. We test both methods on the Clancy and Rudy (2002) I(Na)Markov chain model. With precomputed tables for functions of the transmembrane voltage, both methods are comparable to the forward Euler method in accuracy and computational cost, but allow longer time steps without numerical instability. We conclude that both methods are of practical interest. MRL requires more computations than HOS, but is formulated in general terms which can be readily extended to other Markov chain channel models, whereas the utility of HOS depends on the asymptotic properties of a particular model. The significance of the methods is that they allow a considerable speed-up of large-scale computations of cardiac excitation models by increasing the time step, while maintaining acceptable accuracy and preserving numerical stability.
Collapse
|
13
|
Gomes JM, Alvarenga A, Campos RS, Rocha BM, da Silva APC, dos Santos RW. Uniformization method for solving cardiac electrophysiology models based on the Markov-chain formulation. IEEE Trans Biomed Eng 2014; 62:600-8. [PMID: 25296402 DOI: 10.1109/tbme.2014.2361325] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
This paper compares different numerical methods for the solution of myocyte models of cardiac electrophysiology. In particular, it presents how the technique called uniformization method substantially increases the stability of simple first-order methods such as Euler explicit method and Rush-Larsen (RL) method, for the solution of modern electrophysiology models that are based on continuous-time Markov chains (MCs) for the description of subcellular structures, such as ion channels. The MCs are often associated with stiff ordinary differential equations that severely limit the time step used by these traditional methods. By using the uniformization method, we could significantly increase the time steps for the solution of different cardiac electrophysiology models and improve the computational performance up to 150 times compared to the performance of Euler's and RL's methods.
Collapse
|
14
|
Garcia-Molla VM, Liberos A, Vidal A, Guillem MS, Millet J, Gonzalez A, Martinez-Zaldivar FJ, Climent AM. Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units. Comput Biol Med 2013; 44:15-26. [PMID: 24377685 DOI: 10.1016/j.compbiomed.2013.10.023] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2013] [Revised: 10/21/2013] [Accepted: 10/24/2013] [Indexed: 10/26/2022]
Abstract
In this paper we studied the implementation and performance of adaptive step methods for large systems of ordinary differential equations systems in graphics processing units, focusing on the simulation of three-dimensional electric cardiac activity. The Rush-Larsen method was applied in all the implemented solvers to improve efficiency. We compared the adaptive methods with the fixed step methods, and we found that the fixed step methods can be faster while the adaptive step methods are better in terms of accuracy and robustness.
Collapse
Affiliation(s)
- V M Garcia-Molla
- Department of Information Systems and Computing, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - A Liberos
- ITACA Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - A Vidal
- Department of Information Systems and Computing, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - M S Guillem
- ITACA Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - J Millet
- ITACA Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - A Gonzalez
- ITEAM Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - F J Martinez-Zaldivar
- ITEAM Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| | - A M Climent
- ITACA Institute, Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain.
| |
Collapse
|