需要像 Maple 一样的 Mathematica 短代码

Question

我对 Mathematica 的导出结果有一个问题。两个矩阵 A 和 B 必须以特殊形式导出。 这两个代码列出了从 Maple 导出的数据。 重要的是，用写字板打开的导出文件看起来像列（附加文件）。

如果您已经检查过它是否正常工作，请给我写答案，谢谢！您可以检查您的答案与文件的比较。

代码在这里

矩阵 A 和 B 以及 Maple 中的代码和导出文件

http://www.2shared.com/file/49wW8Z0-/EXAMPLE_EXPORT_MAPLE_FINAL.html

我也会在这里展示给大家看的容易

代码 1)

A := Matrix(2, 2, {(1, 1) = (455200000000/6133413)*w(1), (1, 2) = -(1792000000000/116534847)*w(1), (2, 1) = (455200000000/6133413)*w(2), (2, 2) = -(1792000000000/116534847)*w(2)})


precision := double: writeto(`Aexport.for`):
for i from 1 to 2 do:for j from 1 to 2 do:
if  A[i,j]<>0  then codegen[fortran]([A00[i,j]=A[i,j]],optimized):
fi:od:od:writeto(terminal):

代码 2)

B := Matrix(2, 2, {(1, 1) = 6436781.609, (1, 2) = 0, (2, 1) = 0, (2, 2) = 3862068.966})

  writeto(Bexport);
    for i to 2 do 
    for j to 2 do 
     printf("%016.15E\n", B[i, j]) 
      end do:
        end do:
         writeto(terminal)

Answer 1

这只是（B）部分的翻译：

matrix = {{6436781.609, 0}, {0, 3862068.966}}

Export["Bexport", Map[FortranForm, N@Flatten[matrix]], "Table"]

请测试一下，让我知道它是否适合你。

与 Maple 版本相比的不同之处：E 写为小写，输出的位数不固定（但是，如您所见，所有有效数字都被保留） .这些差异会导致您的应用出现问题吗？

Answer 2

我相信这可以满足您对矩阵 B 的要求：

b = {{6436781.609, 0}, {0, 3862068.966}}

bformatted = 
  NumberForm[
    Flatten@b,
    {16, 15}, 
    NumberFormat -> (Row[{#, "E+", StringTake["00" <> #3, -2]}] &)
  ];

bstring = 
  StringReplace[
    ToString@bformatted,
    {"{"|"}"|" " -> "", "," -> "\n"}
  ];

WriteString["Bexport.dat", bstring, "\n"]

Close["Bexport.dat"]