#271: Example 5d4(Active Contraction based on a user-defined CellML function) doesn't work
View
(Anonymous)
issuedata
Classification:
Bug
Program:
cm
Categories:
CellML
Importance:
Medium
Assigned to:
Created by: rjyoung
Created at: 2007-12-07
progress
Deadline:
2008-02-05 02:16
Percent done:
0
contact
Name:
Robert Young
Description:
I downloaded CM from this website in order to do some finite elasticity calculations; the provided models don't quite fit my problem, so I'd like to use a CellML model, but I'm running into difficulties. Example 5d4 covers some features(CellML, spatially varying parameters) which I'd like to use, but it doesn't work on my machine. The CellML-defined material laws work fine, but CellML- defined active stress doesn't. At first, CM gave me some errors about ">>Incorrect #options" for the files Activation_cellml.ipsolv, Activation.ipsolv, and Activation.ipfit. Once I fixed those, the example ran, but gave incorrect results in both phases; in the first phase, one face of the cube is displaced by .35 instead of .1 as specified. The relevant portion of the .com file seems to be FEM change disp init FEM change disp translate nodes 2,4,6,8 by $disp,0,0 where $disp=.1 . The .ipinit file specifies a displacement of .25, and instead of resetting this to .1, it translates it by .1 In the second phase, when active stress is supposed to be turned on, nothing happens; the solver converges after 0 iterations. The problem seems to be that CellML-defined stress is missing from my build; if I run the example and then type fem def acti;p the first prompt is: Specify type of contraction model [2]: (1) SS tension-length-Ca relation (OLD) (2) Hunter/McCulloch/ter Keurs (coupled to cell/grid model) (3) Unused (4) Outer hair cell stiffness According to the Activation_cellml.ipacti file, the example chooses option 3, which should be (3) Enhance 2nd PKST by user-defined functions The logs on the example page don't seem to run into any of these problems, but they use a binary with a different date; mine is CMISS(cm) version 2.1 Thu Feb 23 09:27:49 NZDT 2006
Transcript
#1:
2007-12-07 02:52 (rjyoung)
categories: "()" ->
"['CellML', '']"
title: "" ->
"Example 5d4(Active Contraction based on a user-defined CellML function) doesn't work"
classification: "Feature" ->
"Bug"
description: "" ->
"I downloaded CM from this website in order to do some finite elasticity calculations; the provided models don't quite fit my problem, so I'd like to use a CellML model, but I'm running into difficulties. Example 5d4 covers some features(CellML, spatially varying parameters) which I'd like to use, but it doesn't work on my machine. The CellML-defined material laws work fine, but CellML-defined active stress doesn't.
At first, CM gave me some errors about ">>Incorrect #options" for the files Activation_cellml.ipsolv, Activation.ipsolv, and Activation.ipfit. Once I fixed those, the example ran, but gave incorrect results in both phases; in the first phase, one face of the cube is displaced by .35 instead of .1 as specified. The relevant portion of the .com file seems to be
FEM change disp init
FEM change disp translate nodes 2,4,6,8 by $disp,0,0
where $disp=.1 . The .ipinit file specifies a displacement of .25, and instead of resetting this to .1, it translates it by .1
In the second phase, when active stress is supposed to be turned on, nothing happens; the solver converges after 0 iterations. The problem seems to be that CellML-defined stress is missing from my build; if I run the example and then type
fem def acti;p
the first prompt is:
Specify type of contraction model [2]:
(1) SS tension-length-Ca relation (OLD)
(2) Hunter/McCulloch/ter Keurs (coupled to cell/grid model)
(3) Unused
(4) Outer hair cell stiffness
According to the Activation_cellml.ipacti file, the example chooses option 3, which should be
(3) Enhance 2nd PKST by user-defined functions
The logs on the example page don't seem to run into any of these problems, but they use a binary with a different date; mine is CMISS(cm) version 2.1 Thu Feb 23 09:27:49 NZDT 2006"
PloneCollectorNG (C) 2003-2004 by ZOPYX - Software Development and Consulting Andreas Jung