- From: alan.garny at physiol.ox.ac.uk (Alan Garny)
- Subject: [cellml-discussion] pcenv development priorities
- Date: Mon, 30 Oct 2006 12:44:57 -0000
>
-----Original Message-----
>
From: cellml-discussion-bounces at cellml.org
>
[mailto:cellml-discussion-bounces at cellml.org] On Behalf Of
>
David Nickerson
>
Sent: 30 October 2006 10:50
>
To: For those interested in contributing to the development of CellML.
>
Subject: Re: [cellml-discussion] pcenv development priorities
>
>
Alan Garny wrote:
>
> As you know, there is a CellML 1.0 version of the
>
epicardial version
>
> of the ten Tusscher model, which we have. I computed that model for
>
> one-second worth of cardiac activity, plotting the trans-membrane
>
> potential. From there, we could extrapolate to 70 minutes by saying
>
> that the frequency of the stimulus is 1 Hz. Like David, I
>
have set the
>
> maximum time step to 0.1 ms. Here are some rough figures:
>
>
>
> Simulation time: 1037.4 s (i.e. ~17 min 17 sec) Computation time:
>
> 684.6 s (i.e. ~11 min 24 sec)
>
>
so the total (predicted) wall clock run time is 17minutes? or
>
28minutes?
Oops, sorry for not making that clear: it is 17 minutes.
>
any chance you could run something longer than a 1s
>
simulation and extrapolate from there, if not run a whole
>
70minutes worth? If not, I guess I can run it next time I'm
>
on a windows machine ;)
Ok, I have quickly set up a CellML 1.0 file with a 1 Hz stimulus protocol
(incidently, for those who are interested on how to do that in CellML, as
well as a few other things, have a look at
http://cor.physiol.ox.ac.uk/FAQ/,
towards the end of that page) and ran it for 70-minute worth of cardiac
activity on another computer (a Pentium 4 2.4 GHz machine with 1.5 GB of
RAM). Here are the figures (using CVODE: BDF+Newton(Dense) and 0.1 ms max
time step):
Simulation time (i.e. wall clock run time): 1107.193 (i.e. ~18 min 27 sec).
Computation time: 945.563 s (i.e. ~15 min 45 sec).
I think it's not bad for some 'manually' generated machine code, eh? :)
Alan.
Archive powered by MHonArc 2.6.18.