Personal tools
You are here: Home cm Tracker Example 5d4(Active Contraction based on a user-defined CellML function) doesn't work

#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
Status: Pending
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:
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"