Hi all,
I am trying to write the expanded form of a matrix determinant to file.
I have successfully done this for a number of matrices, however I am
encountering 'stringout errors' with certain matrices.
The following matrix creates the error:
E:matrix([1,0,-ecvATPase_Pc,0,0,0,0,0,0,0,0,0,0,0],
[0,1,0,-ecvGlyK_Pg,-ecvGlyK_Glc6P,-ecvGlyK_Fru6P,-ecvGlyK_Fru16BP,0,-ecvGlyK_GAP,0,-ecvGlyK_DHAP,-ecvGlyK_BPGA13,0,0],
[1.0,0,-ecvPK_Pc,0,0,0,0,0,0,0,0,0,-ecvPK_Pyr,-ecvPK_Nb],
[JvATPase/(2*JvHK),JvGlyK/(2*JvHK),0,-ecvHK_Pg,-ecvHK_Glc6P,0,0,-ecvHK_GlcI,0,0,0,0,0,0],
[JvATPase/(2*JvHK),JvGlyK/(2*JvHK),0,0,-ecvPGI_Glc6P,-ecvPGI_Fru6P,0,0,0,0,0,0,0,0],
[JvATPase/(2*JvHK),JvGlyK/(2*JvHK),0,-ecvPFK_Pg,0,-ecvPFK_Fru6P,-ecvPFK_Fru16BP,0,0,0,0,0,0,0],
[JvATPase/(2*JvHK),JvGlyK/(2*JvHK),0,-ecvALD_Pg,-ecvALD_Glc6P,-ecvALD_Fru6P,-ecvALD_Fru16BP,0,-ecvALD_GAP,0,-ecvALD_DHAP,-ecvALD_BPGA13,0,0],
[JvATPase/(2*JvHK),JvGlyK/(2*JvHK),0,0,0,0,0,-ecvGlcTr_GlcI,0,0,0,0,0,0],
[1.0,0,0,0,0,0,0,0,-ecvGAPdh_GAP,ecvGAPdh_NADH*NAD/NADH-ecvGAPdh_NAD,0,-ecvGAPdh_BPGA13,0,0],
[1.0,0,0,-ecvGDH_Pg,-ecvGDH_Glc6P,-ecvGDH_Fru6P,-ecvGDH_Fru16BP,0,-ecvGDH_GAP,ecvGDH_NADH*NAD/NADH-ecvGDH_NAD,-ecvGDH_DHAP,-ecvGDH_BPGA13,0,0],
[JvATPase/JvGPO,-JvGlyK/JvGPO,0,-ecvGPO_Pg,-ecvGPO_Glc6P,-ecvGPO_Fru6P,-ecvGPO_Fru16BP,0,-ecvGPO_GAP,0,-ecvGPO_DHAP,-ecvGPO_BPGA13,0,0],
[1.0,0,0,0,0,0,0,0,0,0,0,0,-ecvPyrTr_Pyr,0],
[1.0,0,0,-ecvPGK_Pg,0,0,0,0,0,0,0,-ecvPGK_BPGA13,0,-ecvPGK_Nb],
[JvATPase/(2*JvTPI),-JvGlyK/(2*JvTPI),0,-ecvTPI_Pg,-ecvTPI_Glc6P,-ecvTPI_Fru6P,-ecvTPI_Fru16BP,0,-ecvTPI_GAP,0,-ecvTPI_DHAP,-ecvTPI_BPGA13,0,0]);
I extract the determinant as folows:
den:determinant(E)$
den:expand(den)$
stringout(outfile,den);
To which I get the following error message:
Error in `stringout' attempt -- an error. To debug this try
debugmode(true);
Is there a workaround for this, or an alternative to write data for
really large expressions to a file.
Thanks,
Tim