2

我有一个粒子系统,我基本上跟踪它的坐标。到目前为止,我正在输出文件中的所有坐标,然后使用 gnuplot 脚本绘制它们并使用 ffmpeg 从 PNG 中创建视频。为了跳过所有的文件 I/O,我首先尝试使用 OpenGL,但无济于事,所以我想为 gnuplot 设置一个管道。

我在 MSVC 2012 中这样做:

FILE *gnuplotpipe = _popen(GNUPLOT_NAME, "w");
setvbuf(gnuplotpipe, NULL, _IONBF, NULL);
//setbuf(gnuplotpipe, NULL);
fprintf(gnuplotpipe,"set term png\n");
fprintf(gnuplotpipe,"set size 720, 480\n");
fprintf(gnuplotpipe,"set size ratio -1\n");
fprintf(gnuplotpipe,"unset key\n");
fprintf(gnuplotpipe,"set xrange [-%d:%d]\n", RANGE+100,RANGE+100);
fprintf(gnuplotpipe,"set yrange [-%d:%d]\n", RANGE+100,RANGE+100);
for(time-steps)
    //do stuff
    fprintf(gnuplotpipe,"plot '-'\n");
    for(all particles)
          // calculate coordinates
          fprintf(gnuplotpipe,"%f %f\n", particle[i]->x, particle[i]->y);
    end for
    fprintf(gnuplotpipe, "e\n");
    //fflush(gnuplotpipe);
end for
_pcloce(gnuplotpipe);

并相应地在 linux 中使用gnuplot -persist而不是pgnuplot -persist,popenpclose, 而setbuf不是setvbuf. 这已经有很多试验和错误。

在linux中它是相当一致的。我可以plot '-'在控制台中看到每个时间步的命令,然后是难以理解的文本。什么都没有打开。在 Windows 中这很神奇,因为有时它会起作用,我可以在 gnuplot“MS Windows”窗口上看到数据流(我的意思不是控制台输出,这本身就很神奇,因为我没有使用 wgnuplot_pipes.exe)。没有绘制任何内容。窗口关闭。其他时间会绘制几个时间步,但每个时间步在不同的窗口中。有时,它会崩溃或它会工作但发送难以理解的文本。

我大部分时间都在使用 100 个粒子。我可以有多达 10,000 个时间步长。不知道是不是很多。尽管我将缓冲区设置为 NULL,但我怀疑 Windows 中存在中断或其他问题。我在 VirtualBox 中的 Win7 x64 中的 MSVC 和 Linux Mint x64 中的 Qt 中工作。我正在使用 gnuplot 4.7 (win) 和 4.6.0.8 (linux)。我在 SO 中进行了相当广泛的研究,我发现的唯一问题是它-persist在 Windows 中并没有真正持续存在。任何帮助,将不胜感激。

4

1 回答 1

1

只是这样问题就不会没有答案..这是您的程序的 python 版本,带有 x11 终端。它建立了一个带有愚蠢动画的窗口。

不过请注意,如果其他一些驱动程序尝试为每一帧打开一个窗口,我一点也不感到惊讶

import subprocess
f=subprocess.Popen(['/usr/bin/gnuplot','-persist'],
                      stdin=subprocess.PIPE,shell=False)
f.stdin.write('set term x11 \n')
f.stdin.write('set yrange [-2:2]\n')
import math
i=0
d=1
count=0
while count<1000 :
 count+=1
 if abs(i)==50 : d *= -1
 i+=d 
 f.stdin.write('plot \'-\' with linespoints pt 4 ps 2 \n')
 for j in range(0,500):
  f.stdin.write('%f %f\n'%(float(j),math.sin(float(j)/float(abs(i)+1))))
 f.stdin.write('e\n')
于 2013-03-29T20:33:19.043 回答