OpenCMISS / iron

Source code repository for OpenCMISS-Iron
9 stars 62 forks source link

Implemented analytic Jacobian evaluation for Guccione active contraction #163

Open MarioHabenbacher opened 6 years ago

MarioHabenbacher commented 6 years ago

Updated analytic evaluation of the Jacobian for the Guccione active contraction to incorporate length dependence:

  lambda(1) = 1/(SQRT(2.0_DP*E(1)+1)) !deriv of lambda wrt E
  lambda(2) = 1/(SQRT(2.0_DP*E(2)+1)) !deriv of lambda wrt E
  lambda(3) = 1/(SQRT(2.0_DP*E(3)+1)) !deriv of lambda wrt E

...

      STRESS_TENSOR(component_idx)=STRESS_TENSOR(component_idx)+VALUE*(1.0_DP+1.45_DP*(lambda(component_idx)-1.0_DP))