//Created by libAntimony v2.4
model guccione_mcculloch_waldman_1991_version03__equations(E11, E22, E33, E12, E13, E23, aff, afn, afs, ann, ans, ass, bff, bfn, bfs, bnn, bns, bss, kff, kfn, kfs, knn, kns, kss, Tdev11, Tdev22, Tdev33, Tdev12, Tdev13, Tdev23)

  // Assignment Rules:
  Tdev11 := piecewise( power(aMeabs11, -1.0 * bff) * kff * E11 * (2.0 + E11 * bff * signE11 * power(aMeabs11, -1.0)) , eabs11 <= elim11 , p11 * E11 + q11 );
  Tdev22 := piecewise( power(aMeabs22, -1.0 * bss) * kss * E22 * (2.0 + E22 * bss * signE22 * power(aMeabs22, -1.0)) , eabs22 <= elim22 , p22 * E22 + q22 );
  Tdev33 := piecewise( power(aMeabs33, -1.0 * bnn) * knn * E33 * (2.0 + E33 * bnn * signE33 * power(aMeabs33, -1.0)) , eabs33 <= elim33 , p33 * E33 + q33 );
  Tdev12 := piecewise( power(aMeabs12, -1.0 * bfs) * kfs * E12 * (2.0 + E12 * bfs * signE12 * power(aMeabs12, -1.0)) , eabs12 <= elim12 , p12 * E12 + q12 );
  Tdev13 := piecewise( power(aMeabs13, -1.0 * bfn) * kfn * E13 * (2.0 + E13 * bfn * signE13 * power(aMeabs13, -1.0)) , eabs13 <= elim13 , p13 * E13 + q13 );
  Tdev23 := piecewise( power(aMeabs23, -1.0 * bns) * kns * E23 * (2.0 + E23 * bns * signE23 * power(aMeabs23, -1.0)) , eabs23 <= elim23 , p23 * E23 + q23 );
  elim11 := tol * aff;
  elim22 := tol * ass;
  elim33 := tol * ann;
  elim12 := tol * afs;
  elim13 := tol * afn;
  elim23 := tol * ans;
  etol11 := aff - elim11;
  etol22 := ass - elim22;
  etol33 := ann - elim33;
  etol12 := afs - elim12;
  etol13 := afn - elim13;
  etol23 := ans - elim23;
  signE11 := piecewise( 1.0 , geq(E11, 0.0) , -1.0 );
  signE22 := piecewise( 1.0 , geq(E22, 0.0) , -1.0 );
  signE33 := piecewise( 1.0 , geq(E33, 0.0) , -1.0 );
  signE12 := piecewise( 1.0 , geq(E12, 0.0) , -1.0 );
  signE13 := piecewise( 1.0 , geq(E13, 0.0) , -1.0 );
  signE23 := piecewise( 1.0 , geq(E23, 0.0) , -1.0 );
  eabs11 := signE11 * E11;
  eabs22 := signE22 * E22;
  eabs33 := signE33 * E33;
  eabs12 := signE12 * E12;
  eabs13 := signE13 * E13;
  eabs23 := signE23 * E23;
  aMeabs11 := aff - eabs11;
  aMeabs22 := ass - eabs22;
  aMeabs33 := ann - eabs33;
  aMeabs12 := afs - eabs12;
  aMeabs13 := afn - eabs13;
  aMeabs23 := ans - eabs23;
  mbffM1 := -1 - bff;
  mbfnM1 := -1 - bfn;
  mbfsM1 := -1 - bfs;
  mbnnM1 := -1 - bnn;
  mbnsM1 := -1 - bns;
  mbssM1 := -1 - bss;
  p11 := power(etol11, -1.0 * bff) * kff * (2 + 4 * bff * elim11 * signE11 * power(etol11, -1.0) + -1.0 * bff * power(elim11, 2.0) * mbffM1 * power(etol11, -2.0));
  p12 := power(etol12, -1.0 * bfs) * kfs * (2 + 4 * bfs * elim12 * signE12 * power(etol12, -1.0) + -1.0 * bfs * power(elim12, 2.0) * mbfsM1 * power(etol12, -2.0));
  p13 := power(etol13, -1.0 * bfn) * kfn * (2 + 4 * bfn * elim13 * signE13 * power(etol13, -1.0) + -1.0 * bfn * power(elim13, 2.0) * mbfnM1 * power(etol13, -2.0));
  p22 := power(etol22, -1.0 * bss) * kss * (2 + 4 * bss * elim22 * signE22 * power(etol22, -1.0) + -1.0 * bss * power(elim22, 2.0) * mbssM1 * power(etol22, -2.0));
  p23 := power(etol23, -1.0 * bns) * kns * (2 + 4 * bns * elim23 * signE23 * power(etol23, -1.0) + -1.0 * bns * power(elim23, 2.0) * mbnsM1 * power(etol23, -2.0));
  p33 := power(etol33, -1.0 * bnn) * knn * (2 + 4 * bnn * elim33 * signE33 * power(etol33, -1.0) + -1.0 * bnn * power(elim33, 2.0) * mbnnM1 * power(etol33, -2.0));
  q11 := power(etol11, -1.0 * bff) * kff * elim11 * (2.0 + elim11 * bff * signE11 * power(etol11, -1.0)) - p11 * elim11;
  q12 := power(etol12, -1.0 * bfs) * kfs * elim12 * (2.0 + elim12 * bfs * signE12 * power(etol12, -1.0)) - p12 * elim12;
  q13 := power(etol13, -1.0 * bfn) * kfn * elim13 * (2.0 + elim13 * bfn * signE13 * power(etol13, -1.0)) - p13 * elim13;
  q22 := power(etol22, -1.0 * bss) * kss * elim22 * (2.0 + elim22 * bss * signE22 * power(etol22, -1.0)) - p22 * elim22;
  q23 := power(etol23, -1.0 * bns) * kns * elim23 * (2.0 + elim23 * bns * signE23 * power(etol23, -1.0)) - p23 * elim23;
  q33 := power(etol33, -1.0 * bnn) * knn * elim33 * (2.0 + elim33 * bnn * signE33 * power(etol33, -1.0)) - p33 * elim33;

  // Variable initializations:
  E11 = ;
  E22 = ;
  E33 = ;
  E12 = ;
  E13 = ;
  E23 = ;
  aff = ;
  afn = ;
  afs = ;
  ann = ;
  ans = ;
  ass = ;
  bff = ;
  bfn = ;
  bfs = ;
  bnn = ;
  bns = ;
  bss = ;
  kff = ;
  kfn = ;
  kfs = ;
  knn = ;
  kns = ;
  kss = ;
  tol = 0.9;
end

model guccione_mcculloch_waldman_1991_version03__interface(E11, E22, E33, E12, E13, E23, aff, afn, afs, ann, ans, ass, bff, bfn, bfs, bnn, bns, bss, kff, kfn, kfs, knn, kns, kss, Tdev11, Tdev22, Tdev33, Tdev12, Tdev13, Tdev23)

  // Sub-modules, and any changes to those submodules:
  equations: guccione_mcculloch_waldman_1991_version03__equations(E11, E22, E33, E12, E13, E23, aff, afn, afs, ann, ans, ass, bff, bfn, bfs, bnn, bns, bss, kff, kfn, kfs, knn, kns, kss, Tdev11, Tdev22, Tdev33, Tdev12, Tdev13, Tdev23);

  // Variable initializations:
  aff = 0;
  afn = 0;
  afs = 0;
  ann = 0;
  ans = 0;
  ass = 0;
  bff = 0;
  bfn = 0;
  bfs = 0;
  bnn = 0;
  bns = 0;
  bss = 0;
  kff = 0;
  kfn = 0;
  kfs = 0;
  knn = 0;
  kns = 0;
  kss = 0;
end

model *guccione_mcculloch_waldman_1991_version03____main()

  // Sub-modules, and any changes to those submodules:
  interface: guccione_mcculloch_waldman_1991_version03__interface(E11, E22, E33, E12, E13, E23, aff, afn, afs, ann, ans, ass, bff, bfn, bfs, bnn, bns, bss, kff, kfn, kfs, knn, kns, kss, Tdev11, Tdev22, Tdev33, Tdev12, Tdev13, Tdev23);
end