myokit 1.35.4__py3-none-any.whl → 1.36.1__py3-none-any.whl

myokit/tests/data/decker.model

@@ -113,7 +113,7 @@ calcium_alpha = 0.05;
 calcium_beta = 0.00068;
 csqnBar = 10.0;
 diff_fTrap = calcium_alpha * active * (active - fTrap) - calcium_beta * fTrap;
-diff_uCa_cal = (-ICaL) * AF / (v_cal * 2.0) + CaL * v_sr / v_cal;
+diff_uCa_cal = -ICaL * AF / (v_cal * 2.0) + CaL * v_sr / v_cal;
 diff_uCa_i = Ca * v_sr / v_cyt - (ICab + IpCa - 2.0 * INaCa) * AF / (v_cyt * 2.0) - (Iup - Ileak) * v_nsr / v_cyt;
 diff_uCa_jsr = Itr - Irel;
 diff_uCa_sr = INaCaSS * AF / v_sr + Irel * v_jsr / v_sr - (Ca + CaL);
@@ -122,11 +122,11 @@ diff_uCa_sr = INaCaSS * AF / v_sr + Irel * v_jsr / v_sr - (Ca + CaL);
 BSLBar = 1.124;
 BSRBar = 0.047;
 BS_sum = BSRBar + BSLBar;
-Ca_cal = (-b_2) / 3.0 + 2.0 / 3.0 * sqrt(b_2 * b_2 - 3.0 * c_2) * cos(acos((9.0 * b_2 * c_2 - 2.0 * pow(b_2, 3.0) - 27.0 * d_1) / (2.0 * pow(b_2 * b_2 - 3.0 * c_2, 1.5))) / 3.0);
+Ca_cal = -b_2 / 3.0 + 2.0 / 3.0 * sqrt(b_2 * b_2 - 3.0 * c_2) * cos(acos((9.0 * b_2 * c_2 - 2.0 * pow(b_2, 3.0) - 27.0 * d_1) / (2.0 * pow(b_2 * b_2 - 3.0 * c_2, 1.5))) / 3.0);
 b_2 = sum2 - uCa_cal;
 c_2 = pro_sum2 - Km_sum2 * uCa_cal;
-d_1 = (-Km_pro2) * uCa_cal;
-Ca_i = (-b_3) / 3.0 + 2.0 / 3.0 * sqrt(b_3 * b_3 - 3.0 * c_3) * cos(acos((9.0 * b_3 * c_3 - 2.0 * pow(b_3, 3.0) - 27.0 * d_2) / (2.0 * pow(b_3 * b_3 - 3.0 * c_3, 1.5))) / 3.0);
+d_1 = -Km_pro2 * uCa_cal;
+Ca_i = -b_3 / 3.0 + 2.0 / 3.0 * sqrt(b_3 * b_3 - 3.0 * c_3) * cos(acos((9.0 * b_3 * c_3 - 2.0 * pow(b_3, 3.0) - 27.0 * d_2) / (2.0 * pow(b_3 * b_3 - 3.0 * c_3, 1.5))) / 3.0);
 Km_cmdn = 0.00238;
 Km_pro = Km_cmdn * Km_trpn;
 Km_sum = Km_trpn + Km_cmdn;
@@ -134,14 +134,14 @@ Km_trpn = 0.0005;
 b_3 = sum - uCa_i;
 c_3 = pro - uCa_i * Km_sum;
 cmdnBar = 0.05;
-d_2 = (-Km_pro) * uCa_i;
+d_2 = -Km_pro * uCa_i;
 pro = Km_pro + trpnBar * Km_cmdn + cmdnBar * Km_trpn;
 sum = cmdnBar + trpnBar + Km_trpn + Km_cmdn;
 trpnBar = 0.07;
-Ca_sr = (-b_4) / 3.0 + 2.0 / 3.0 * sqrt(b_4 * b_4 - 3.0 * c_4) * cos(acos((9.0 * b_4 * c_4 - 2.0 * pow(b_4, 3.0) - 27.0 * d_3) / (2.0 * pow(b_4 * b_4 - 3.0 * c_4, 1.5))) / 3.0);
+Ca_sr = -b_4 / 3.0 + 2.0 / 3.0 * sqrt(b_4 * b_4 - 3.0 * c_4) * cos(acos((9.0 * b_4 * c_4 - 2.0 * pow(b_4, 3.0) - 27.0 * d_3) / (2.0 * pow(b_4 * b_4 - 3.0 * c_4, 1.5))) / 3.0);
 b_4 = sum2 - uCa_sr;
 c_4 = pro_sum2 - Km_sum2 * uCa_sr;
-d_3 = (-Km_pro2) * uCa_sr;
+d_3 = -Km_pro2 * uCa_sr;
 Km_BSL = 0.0087;
 Km_BSR = 0.00087;
 Km_pro2 = Km_BSR * Km_BSL;
@@ -152,7 +152,7 @@ sum2 = BS_sum + Km_sum2;
 // cell
 A = 2.0 * geoArea;
 AF = A / F;
-geoArea = 3.14 * 2.0 * radius * radius + 2.0 *  3.14159265358979312e+00 * radius * length;
+geoArea = 3.14 * 2.0 * radius * radius + 2.0 * 3.14159265358979312e+00 * radius * length;
 length = 0.01;
 radius = 0.0011;
 v_cal = volume * 0.002;
@@ -181,9 +181,9 @@ z1_1 = EK - ECl;
 z2_1 = 87.8251;
 CTNaCl = NaClBar * z1_2 / (z1_2 + z2_2);
 z1_2 = pow(ENa - ECl, 4.0);
-z2_2 =  5.94941976406925395e+07;
+z2_2 = 5.94941976406925395e+07;
 KClBar = 1.77e-05;
-NaClBar =  2.46108000000000016e-05;
+NaClBar = 2.46108000000000016e-05;
 
 // diff
 Ca = (Ca_sr - Ca_i) / var_diff_tau_diff;
@@ -206,10 +206,10 @@ e2v_1 = exp(2.0 * vf);
 // ical
 ACT_ss = 1.0 / (1.0 + exp((p05_1 - V) / p06_1));
 ACT_tau = p00_1 + p01_1 * exp(p02_1 * (V + p03_1)) / (1.0 + exp(p04_1 * (V + p03_1)));
-diff_C = (-(ical_alpha + delta_1 + yy)) * C + ical_beta * O + theta * Cs + xx * CI;
-diff_CI = (-(ical_alpha + delta_I + xx)) * CI + yy * C + theta_I * CIs + ical_beta * OI;
-diff_CIs = (-(ical_alpha + theta_I + xs)) * CIs + ys * Cs + delta_I * CI + ical_beta * OIs;
-diff_Cs = (-(ical_alpha + theta + ys)) * Cs + delta_1 * C + ical_beta * Os + xs * CIs;
+diff_C = -(ical_alpha + delta_1 + yy) * C + ical_beta * O + theta * Cs + xx * CI;
+diff_CI = -(ical_alpha + delta_I + xx) * CI + yy * C + theta_I * CIs + ical_beta * OI;
+diff_CIs = -(ical_alpha + theta_I + xs) * CIs + ys * Cs + delta_I * CI + ical_beta * OIs;
+diff_Cs = -(ical_alpha + theta + ys) * Cs + delta_1 * C + ical_beta * Os + xs * CIs;
 ICaL = ICaLBar * g_1;
 ICaLBar = ((e2v_2 == 1.0) ? p21 * 4.0 * F * FRT * (gCai * Ca_cal - gCao * Cao) / (2.0 * FRT) : p21 * 4.0 * vff * (gCai * Ca_cal * e2v_2 - gCao * Cao) / (e2v_2 - 1.0));
 e2v_2 = exp(2.0 * vf);
@@ -217,16 +217,16 @@ g_1 = O + Os;
 gCai = 1.0;
 gCao = 0.341;
 I_V1 = p07_1 * (1.0 + exp((V + p08_1) / p09_1));
-I_V2 = p10_1 * (1.0 + exp((-(V + p11_1)) / p12_1));
+I_V2 = p10_1 * (1.0 + exp(-(V + p11_1) / p12_1));
 I_Vss = (p15_1 + 1.0 / (1.0 + exp((V + p13_1) / p14_1))) / (1.0 + p15_1);
 I_Vtau = 1.0 / (1.0 / I_V1 + 1.0 / I_V2);
 Is_Ca = p16 - p17 / (1.0 + zz_2);
-Is_V2 = Is_Ca * (1.0 + exp((-(V + p11_1)) / p12_1));
+Is_V2 = Is_Ca * (1.0 + exp(-(V + p11_1) / p12_1));
 Is_Vss = (p18 + 1.0 / (1.0 + exp((V + p13_1) / p14_1))) / (1.0 + p18);
 Is_Vtau = 1.0 / (1.0 / I_V1 + 1.0 / Is_V2);
-diff_O = (-(ical_beta + delta_1 + yy)) * O + ical_alpha * C + theta * Os + xx * OI;
-diff_OI = (-(ical_beta + delta_I + xx)) * OI + yy * O + theta_I * OIs + ical_alpha * CI;
-diff_OIs = (-(ical_beta + theta_I + xs)) * OIs + ys * Os + delta_I * OI + ical_alpha * CIs;
+diff_O = -(ical_beta + delta_1 + yy) * O + ical_alpha * C + theta * Os + xx * OI;
+diff_OI = -(ical_beta + delta_I + xx) * OI + yy * O + theta_I * OIs + ical_alpha * CI;
+diff_OIs = -(ical_beta + theta_I + xs) * OIs + ys * Os + delta_I * OI + ical_alpha * CIs;
 Os = 1.0 - C - CI - CIs - Cs - O - OI - OIs;
 ical_alpha = ACT_ss / ACT_tau;
 ical_beta = (1.0 - ACT_ss) / ACT_tau;
@@ -255,7 +255,7 @@ p17 = 17.5;
 p18 = 0.0001;
 p19 = 3.0;
 p20 = 1e-06;
-p21 =  1.55520000000000002e-04;
+p21 = 1.55520000000000002e-04;
 theta = 1.0;
 theta_I = p20;
 xs = Is_Vss / Is_Vtau;
@@ -274,7 +274,7 @@ Ibar_1 = ((evf == 1.0) ? 9e-07 * vff * Clo : 9e-07 * vff * (Cl_i - Clo * evf) /
 KCaI = 1.0 - 1.0 / (1.0 + pow(Irel, 2.0) / 0.16);
 diff_i_fast_1 = (ss - i_fast_1) / 8.0;
 i_fast_alpha = 0.025 / (1.0 + exp((V + 58.0) / 5.0));
-i_fast_beta = 1.0 / (5.0 * (1.0 + exp((V + 19.0) / (-9.0))));
+i_fast_beta = 1.0 / (5.0 * (1.0 + exp((V + 19.0) / -9.0)));
 ss = i_fast_alpha / (i_fast_alpha + i_fast_beta);
 
 // ik1
@@ -283,17 +283,17 @@ Gbar_2 = 0.5;
 IK1 = GK1 * K1 * dv;
 K1 = ak / (ak + bk);
 ak = 1.02 / (1.0 + exp(0.2385 * (dv - 59.215)));
-bk = (0.49124 * exp(0.08032 * (dv + 5.476)) + exp(0.06175 * (dv - 594.31))) / (1.0 + exp((-0.5143) * (dv + 4.753)));
+bk = (0.49124 * exp(0.08032 * (dv + 5.476)) + exp(0.06175 * (dv - 594.31))) / (1.0 + exp(-0.5143 * (dv + 4.753)));
 dv = V - EK;
 
 // ikr
 G = Gbar_3 * sqrt(Ko / 5.4);
 Gbar_3 = 0.0138542;
 IKr = G * a_1 * i_1 * (V - EK);
-a_1_inf = 1.0 / (1.0 + exp((V + 10.085) / (-4.25)));
+a_1_inf = 1.0 / (1.0 + exp((V + 10.085) / -4.25));
 tau_a_1 = 1.0 / (0.0006 * ta - 0.0003 * tb);
 ta = ((dd_1 == 0.0) ? 1.0 / 0.136 : vv_1 / dd_1);
-dd_1 = -expm1((-0.136) * vv_1);
+dd_1 = -expm1(-0.136 * vv_1);
 vv_1 = V - 1.7384;
 tb = ((dd_2 == 0.0) ? 1.0 / 0.1522 : vv_2 / dd_2);
 dd_2 = -expm1(0.1522 * vv_2);
@@ -318,7 +318,7 @@ diff_C8 = 2.0 * a_2 * C7 + 3.0 * b_5 * C9 + 3.0 * g_2 * C4 + 2.0 * d_4 * C11 - C
 diff_C9 = a_2 * C8 + 4.0 * g_2 * C5 + 2.0 * d_4 * C12 - C9 * (3.0 * b_5 + 3.0 * g_2 + d_4);
 GKs = p15_2 * (1.0 + 0.6 / (1.0 + pow(3.8e-05 / Ca_i, 1.4)));
 IKs = GKs * (O1 + O2) * (V - EKs);
-diff_O1 = (-(e + p)) * O1 + o * O2 + t * C15;
+diff_O1 = -(e + p) * O1 + o * O2 + t * C15;
 O2 = 1.0 - C1 - C2 - C3 - C4 - C5 - C6 - C7 - C8 - C9 - C10 - C11 - C12 - C13 - C14 - C15 - O1;
 a_2 = p00_2 * exp(vf * p01_2);
 b_5 = p02_2 * exp(vf * p03_2);
@@ -327,21 +327,21 @@ e = p08_2 * exp(vf * p09_2);
 g_2 = p04_2 * exp(vf * p05_2);
 o = p11_2 * exp(vf * p12_2);
 p = p13_2 * exp(vf * p14_2);
-p00_2 =  1.48645979808599994e-02;
-p01_2 =  2.98773012358800015e-02;
-p02_2 =  8.39863121998299933e-02;
-p03_2 = (- 5.54610571266399982e-02);
-p04_2 =  1.46006611831600002e-02;
-p05_2 =  2.44649530996450010e-01;
-p06_2 =  3.11732688739999996e-03;
-p07_2 = (- 4.26254519443760027e-01);
-p08_2 =  7.73199009733099990e-02;
-p09_2 = (- 6.47261224887100001e-02);
-p10_2 =  8.95383064110199939e-02;
-p11_2 =  7.94054599586400012e-01;
-p12_2 = (- 8.01737819297699966e-02);
-p13_2 =  5.86382286630140026e-01;
-p14_2 =  2.82055543314960022e-01;
+p00_2 = 1.48645979808599994e-02;
+p01_2 = 2.98773012358800015e-02;
+p02_2 = 8.39863121998299933e-02;
+p03_2 = -5.54610571266399982e-02;
+p04_2 = 1.46006611831600002e-02;
+p05_2 = 2.44649530996450010e-01;
+p06_2 = 3.11732688739999996e-03;
+p07_2 = -4.26254519443760027e-01;
+p08_2 = 7.73199009733099990e-02;
+p09_2 = -6.47261224887100001e-02;
+p10_2 = 8.95383064110199939e-02;
+p11_2 = 7.94054599586400012e-01;
+p12_2 = -8.01737819297699966e-02;
+p13_2 = 5.86382286630140026e-01;
+p14_2 = 2.82055543314960022e-01;
 p15_2 = 0.0826;
 t = p10_2;
 
@@ -353,14 +353,14 @@ K = 1.0 / (1.0 + exp((7.488 - V) / 5.98));
 // ina
 Gbar_5 = 9.075;
 INa = Gbar_5 * pow(m, 3.0) * h * j * (V - ENa);
-alpha_h = ((V >= (-40.0)) ? 0.0 : 0.135 * exp((80.0 + V) / (-6.8)));
-beta_h = ((V >= (-40.0)) ? 1.0 / (0.13 * (1.0 + exp((V + 10.66) / (-11.1)))) : 3.56 * exp(0.079 * V) + 310000.0 * exp(0.35 * V));
-alpha_j = ((V >= (-40.0)) ? 0.0 : ((-127140.0) * exp(0.2444 * V) - 6.948e-05 * exp((-0.04391) * V)) * (V + 37.78) / (1.0 + exp(0.311 * (V + 79.23))));
-beta_j = ((V >= (-40.0)) ? 0.3 * exp((-2.535e-07) * V) / (1.0 + exp((-0.1) * (V + 32.0))) : 0.1212 * exp((-0.01052) * V) / (1.0 + exp((-0.1378) * (V + 40.14))));
+alpha_h = ((V >= -40.0) ? 0.0 : 0.135 * exp((80.0 + V) / -6.8));
+beta_h = ((V >= -40.0) ? 1.0 / (0.13 * (1.0 + exp((V + 10.66) / -11.1))) : 3.56 * exp(0.079 * V) + 310000.0 * exp(0.35 * V));
+alpha_j = ((V >= -40.0) ? 0.0 : (-127140.0 * exp(0.2444 * V) - 6.948e-05 * exp(-0.04391 * V)) * (V + 37.78) / (1.0 + exp(0.311 * (V + 79.23))));
+beta_j = ((V >= -40.0) ? 0.3 * exp(-2.535e-07 * V) / (1.0 + exp(-0.1 * (V + 32.0))) : 0.1212 * exp(-0.01052 * V) / (1.0 + exp(-0.1378 * (V + 40.14))));
 alpha_m = ((dd_3 == 0.0) ? 3.2 : 0.32 * vv_3 / dd_3);
-dd_3 = -expm1((-0.1) * vv_3);
+dd_3 = -expm1(-0.1 * vv_3);
 vv_3 = V + 47.13;
-beta_m = 0.08 * exp((-V) / 11.0);
+beta_m = 0.08 * exp(-V / 11.0);
 
 // inab
 INaB = ((evf == 1.0) ? P * F * (Nai - Nao) : P * vff * (Nai * evf - Nao) / (evf - 1.0));
@@ -398,21 +398,21 @@ vMax = 4.5;
 // inak
 INaK = Ibar_2 * fv * pow(Pna, 3.0) * Pk;
 Pna = Nai / (Nai + K_Nai);
-fv = 1.0 / (1.0 + exp((-phi)));
+fv = 1.0 / (1.0 + exp(-phi));
 phi = (V - VHalf) * FRT;
 Ibar_2 = 1.4;
 K_Ko = 1.5;
 K_Nai = 2.6;
 Pk = Ko / (Ko + K_Ko);
-VHalf = (-92.0);
+VHalf = -92.0;
 
 // inal
 Gbar_6 = 0.0065;
 INaL = Gbar_6 * pow(a_3, 3.0) * i_2 * (V - ENa);
 alpha_a_3 = ((dd_4 == 0.0) ? 3.2 : 0.32 * vv_4 / dd_4);
-dd_4 = -expm1((-0.1) * vv_4);
+dd_4 = -expm1(-0.1 * vv_4);
 vv_4 = V + 47.13;
-beta_a_3 = 0.08 * exp((-V) / 11.0);
+beta_a_3 = 0.08 * exp(-V / 11.0);
 diff_i_2 = (hLss - i_2) / 600.0;
 hLss = 1.0 / (1.0 + exp((V + 91.0) / 6.1));
 
@@ -422,7 +422,7 @@ IpCa_bar = 0.0575;
 Km_pCa = 0.0005;
 
 // irel
-diff_Irel = (-(irel_inf_var + Irel)) / irel_tau;
+diff_Irel = -(irel_inf_var + Irel) / irel_tau;
 irel_alpha = irel_beta * 0.1125;
 irel_beta = 4.75 * (1.0 + delta_2);
 delta_2 = 1.0 / (1.0 + pow(0.28 / active, 10.0));
@@ -433,11 +433,11 @@ irel_tau = irel_beta / (1.0 + 0.0123 / Ca_jsr);
 Gbar_7 = 0.4975;
 ITo = Gbar_7 * pow(a_4, 3.0) * i_fast_2 * i_slow * Rto1 * (V - EK);
 Rto1 = exp(V / 550.0);
-ito_a_alpha = 1.0 / (1.2089 * (1.0 + exp((V - 18.4099) / (-29.3814))));
+ito_a_alpha = 1.0 / (1.2089 * (1.0 + exp((V - 18.4099) / -29.3814)));
 ito_a_beta = 3.5 / (1.0 + exp((V + 100.0) / 29.3814));
-a_4_inf = 1.0 / (1.0 + exp((V + 9.437) / (-7.133)));
+a_4_inf = 1.0 / (1.0 + exp((V + 9.437) / -7.133));
 tau_a_4 = 1.0 / (ito_a_alpha + ito_a_beta);
-ito_beta_i = 1.0 / (9.7953 * (1.0 + exp((V + 19.0) / (-9.0))));
+ito_beta_i = 1.0 / (9.7953 * (1.0 + exp((V + 19.0) / -9.0)));
 beta_i_fast_2 = ito_beta_i;
 beta_i_slow = ito_beta_i;
 alpha_i_fast_2 = 0.025 / (1.0 + exp((V + 58.0) / 5.0));
@@ -446,10 +446,10 @@ alpha_i_slow = 1.0 / (250.0 * (1.0 + exp((V + 60.0) / 5.0)));
 // membrane
 i_ion = INa_tot + IK_tot + ICa_tot + ICl_tot;
 i_stim = 0.0 * stim_amplitude;
-stim_amplitude = (-80.0);
+stim_amplitude = -80.0;
 
 // nernst
-ECl = (-RTF) * log(Clo / Cl_i);
+ECl = -RTF * log(Clo / Cl_i);
 EK = RTF * log(Ko / K_i);
 EKs = RTF * log((Ko + P_NaK * Nao) / (K_i + P_NaK * Nai));
 ENa = RTF * log(Nao / Nai);
@@ -464,14 +464,14 @@ T = 310.0;
 
 // potassium
 IK_tot = IK1 + IKr + IKs + IKur + ITo - INaK * 2.0;
-diff_K_i = (-IK_tot) * (AF / v_cyt) + CTKCl;
+diff_K_i = -IK_tot * (AF / v_cyt) + CTKCl;
 
 // sodium
 INa_cyt = INa + INaL + INaB + INaK * 3.0 + INaCa * 3.0;
 INa_tot = INa_cyt + INaCaSS * 3.0;
-diff_Na_sr = (-(3.0 * INaCaSS * r1_3 + Na));
+diff_Na_sr = -(3.0 * INaCaSS * r1_3 + Na);
 r1_3 = AF / v_sr;
-diff_Nai = (-r1_4) * INa_cyt + r2_2 * Na + CTNaCl;
+diff_Nai = -r1_4 * INa_cyt + r2_2 * Na + CTNaCl;
 r1_4 = AF / v_cyt;
 r2_2 = v_sr / v_cyt;
 

myokit/tests/test_formats.py

@@ -5,6 +5,7 @@
 # This file is part of Myokit.
 # See http://myokit.org for copyright, sharing, and licensing details.
 #
+import os
 import unittest
 
 import myokit
@@ -14,9 +15,7 @@ from myokit.tests import TemporaryDirectory
 
 
 class FormatsTest(unittest.TestCase):
-    """
-    Test infrastructure from the formats module.
-    """
+    """ Test shared formats functionality. """
 
     def test_register_external_importer(self):
         # Tests registering an external importer.
@@ -101,10 +100,8 @@ class FormatsTest(unittest.TestCase):
 
 
 class ExporterTest(unittest.TestCase):
-    """
-    Tests some parts of :class:`Exporter` that are not tested when testing the
-    various exporter implementations.
-    """
+    """ Tests shared :class:`Exporter` functionality. """
+
     def test_writable_dir(self):
         # Test :meth:`_test_writable_dir` for existing paths.
 
@@ -126,6 +123,117 @@ class ExporterTest(unittest.TestCase):
             self.assertRaisesRegex(
                 myokit.ExportError, 'file exists', e.runnable, path, m, p)
 
+    def test_runnable_exporter(self):
+        # Test shared functionality of the TemplatedRunnableExporters.
+
+        e = myokit.formats.exporter('ansic')
+
+        # Load model, protocol
+        m, p, x = myokit.load('example')
+
+        # Create empty output directory as subdirectory of DIR_OUT
+        with TemporaryDirectory() as d:
+            path = d.path()
+
+            # Simple export
+            dpath = os.path.join(path, 'runnable1')
+            ret = e.runnable(dpath, m)
+            self.assertIsNone(ret)
+            self.assertTrue(os.path.isdir(dpath))
+            self.assertTrue(len(os.listdir(dpath)) > 0)
+
+            # Write to complex path
+            dpath = os.path.join(path, 'runnable2', 'nest', 'test')
+            ret = e.runnable(dpath, m, p)
+            self.assertIsNone(ret)
+            self.assertTrue(os.path.isdir(dpath))
+            self.assertTrue(len(os.listdir(dpath)) > 0)
+
+            # Overwrite existing path
+            ret = e.runnable(dpath, m, p)
+            self.assertIsNone(ret)
+            self.assertTrue(os.path.isdir(dpath))
+            self.assertTrue(len(os.listdir(dpath)) > 0)
+
+            # Path pointing to file
+            dpath = os.path.join(path, 'file')
+            with open(dpath, 'w') as f:
+                f.write('contents\n')
+            self.assertRaisesRegex(
+                myokit.ExportError, 'file exists', e.runnable, dpath, m, p)
+
+            # Directory exists where we're trying to write a file
+            dpath = os.path.join(path, 'runnable3')
+            fname = os.path.join(dpath, 'sim.c')
+            os.makedirs(fname)
+            self.assertRaisesRegex(
+                myokit.ExportError, 'Directory exists',
+                e.runnable, dpath, m, p)
+
+            # Directory embedded in the output file path
+            def embedded():
+                return {'sim.c': 'nested/sim.c'}
+
+            # 1. Normal operation
+            e._dict = embedded
+            dpath = os.path.join(path, 'runnable4')
+            ret = e.runnable(dpath, m, p)
+            self.assertIsNone(ret)
+            self.assertTrue(os.path.isdir(dpath))
+            self.assertTrue(len(os.listdir(dpath)) > 0)
+
+            # 2. Try to create directory where file exists
+            def embedded():
+                return {'sim.c': 'nested/sim.c/som.c'}
+
+            e._dict = embedded
+            dpath = os.path.join(path, 'runnable4')
+            self.assertRaisesRegex(
+                myokit.ExportError, 'file or link', e.runnable, dpath, m, p)
+
+
+class ImporterTest(unittest.TestCase):
+    """ Test shared importer functionality. """
+
+    def test_importer_interface(self):
+        # Test listing and creating importers.
+        ims = myokit.formats.importers()
+        self.assertTrue(len(ims) > 0)
+        for i in ims:
+            self.assertIsInstance(i, str)
+            i = myokit.formats.importer(i)
+            self.assertTrue(isinstance(i, myokit.formats.Importer))
+
+    def test_unknown(self):
+        # Test requesting an unknown importer.
+        # Test fetching using importer method
+        self.assertRaisesRegex(
+            KeyError, 'Importer not found', myokit.formats.importer, 'blip')
+
+
+class EWriterTest(unittest.TestCase):
+    """ Test shared ewriter functionality. """
+
+    def test_ewriter_interface(self):
+        # Test listing and creating expression writers.
+
+        # Test listing
+        es = myokit.formats.ewriters()
+        self.assertTrue(len(es) > 0)
+
+        # Create one of each
+        for e in es:
+            self.assertIsInstance(e, str)
+            e = myokit.formats.ewriter(e)
+            self.assertTrue(isinstance(e, myokit.formats.ExpressionWriter))
+
+    def test_unknown(self):
+        # Test requesting an unknown expression writer.
+        # Test fetching using ewriter method
+        self.assertRaisesRegex(
+            KeyError, 'Expression writer not found', myokit.formats.ewriter,
+            'dada')
+
 
 if __name__ == '__main__':
     unittest.main()