如何在 Maxima 中格式化 Fortran 输出

Question

我正在使用 Maxima（又名 Macsyma）生成一些系数，而 Maxima 必须在外部文件中写入表达式的 fortran 表示，以便 python 程序读取。

我正在使用以下代码写下系数：

with_stdout ("An.txt",
    for n : 0 thru N do
        fortran(A[n]));

但不幸的是，它们被写在文件中 An.txt 完全像这样：

      exp(beta*t)*(-alpha*(beta**2*t**2+2*beta*t+2)*exp(-beta*t)*D/beta*
     1   *3+alpha*(beta*t+1)*exp(-beta*t)*D/beta**2-alpha**2*(beta**4*t*
     2   *4+4*beta**3*t**3+12*beta**2*t**2+24*beta*t+24)*exp(-beta*t)/be
     3   ta**5+2*alpha**2*(beta**3*t**3+3*beta**2*t**2+6*beta*t+6)*exp(-
     4   beta*t)/beta**4-alpha**2*(beta**2*t**2+2*beta*t+2)*exp(-beta*t)
     5   /beta**3+%c)
      alpha*(t-1)*t
      beta
      (2*(beta*D+3*alpha*beta*(t-1)*t)+alpha*t+alpha*(t-1))/A_0(t)/6.0E+
     1   0

是的，在每个条目之前有一个五（或六个）-space 缩进，数字表示表达式的延续。

我有没有机会得到一个更容易阅读的版本？一个完整行中的每个条目都是理想的。我已经尝试过函数 f90 但它什么也没打印出来，我也尝试过使用 fortindent 和 fortspaces 但没有任何变化（不过我可能对它们使用了错误的语法） .

提前致谢。

Answer 1

我认为你可以在告诉它允许很长的行之后使用 f90。请注意，默认情况下不加载 f90。

(%i1) load (f90) $
(%i2) :lisp (setq *f90-output-line-length-max* 1000000000)
1000000000
(%i2) A : makelist (expand ((foo + exp(bar*t))^n), n, [5, 10, 20]) $
(%i3) with_stdout ("A.txt", for e in A do f90(e));
(%o3)                                done

请注意，您必须通过 Lisp 设置行长；那是个疣。这是我得到的输出：

$ cat A.txt
exp(5*bar*t)+5*foo*exp(4*bar*t)+10*foo**2*exp(3*bar*t)+10*foo**3*exp(2*bar*t)+5*foo**4*exp(bar*t)+foo**5
exp(10*bar*t)+10*foo*exp(9*bar*t)+45*foo**2*exp(8*bar*t)+120*foo**3*exp(7*bar*t)+210*foo**4*exp(6*bar*t)+252*foo**5*exp(5*bar*t)+210*foo**6*exp(4*bar*t)+120*foo**7*exp(3*bar*t)+45*foo**8*exp(2*bar*t)+10*foo**9*exp(bar*t)+foo**10
exp(20*bar*t)+20*foo*exp(19*bar*t)+190*foo**2*exp(18*bar*t)+1140*foo**3*exp(17*bar*t)+4845*foo**4*exp(16*bar*t)+15504*foo**5*exp(15*bar*t)+38760*foo**6*exp(14*bar*t)+77520*foo**7*exp(13*bar*t)+125970*foo**8*exp(12*bar*t)+167960*foo**9*exp(11*bar*t)+184756*foo**10*exp(10*bar*t)+167960*foo**11*exp(9*bar*t)+125970*foo**12*exp(8*bar*t)+77520*foo**13*exp(7*bar*t)+38760*foo**14*exp(6*bar*t)+15504*foo**15*exp(5*bar*t)+4845*foo**16*exp(4*bar*t)+1140*foo**17*exp(3*bar*t)+190*foo**18*exp(2*bar*t)+20*foo**19*exp(bar*t)+foo**20
$ wc A.txt
  3   3 856 A.txt

我看到你的输出中有一个 %c；也许你需要用一些东西来代替它以使输出可以接受 Python?

对了，定宽连字符是打卡时代遗留下来的神器....

如何在 Maxima 中格式化 Fortran 输出

How to format fortran output in Maxima

formatting

indentation

maxima

output