gnuplot and usleep on C

Question

我正在尝试在Linux上进行实时绘图。一切运行正常，只是当程序结束时，会有一段时间占用大量CPU。此外，因为我在usleep上使用的时间很短（比如0.3秒），我不知道为什么bash shell由gnuplot接口主导，即使对于我试图打开的其他bash窗口，我也不知道，所以我不能使用ctrl+c来停止程序。

除了usleep，我还有其他功能可以用吗？或者我可以纠正一些事情来避免gnuplot对windows的统治？

提前感谢并为我糟糕的英语感到抱歉

我让这里的代码

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
float s=10.;
float r=28.;
float b=8.0/3.0;
/* Definimos las funciones */
float f(float x,float y,float z){
    return s*(y-x);
}
float g(float x,float y,float z){
    return x*(r-z)-y;
}
float h(float x,float y,float z){
    return x*y-b*z;
}
FILE *output;
FILE *gp;
int main(){
    gp = popen("gnuplot -","w");
    output = fopen("lorenzgplot.dat","w");
    float t=0.; 
    float dt=0.01;
    float tf=30;
    float x=3.;
    float y=2.;
    float z=0.;
    float k1x,k1y,k1z, k2x,k2y,k2z,k3x,k3y,k3z,k4x,k4y,k4z;
    fprintf(output,"%f %f %f n",x,y,z);
    fprintf(gp, "splot '/home/david/documents/lorenzgplot.dat' with lines n");
/* Ahora Runge Kutta de orden 4 */  
    while(t<tf){
        /* RK4 paso 1 */
        k1x = f(x,y,z)*dt;
        k1y = g(x,y,z)*dt;
        k1z = h(x,y,z)*dt;
        /* RK4 paso 2 */
        k2x = f(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        k2y = g(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        k2z = h(x+0.5*k1x,y+0.5*k1y,z+0.5*k1z)*dt;
        /* RK4 paso 3 */
        k3x = f(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        k3y = g(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        k3z = h(x+0.5*k2x,y+0.5*k2y,z+0.5*k2z)*dt;
        /* RK4 paso 4 */
        k4x = f(x+k3x,y+k3y,z+k3z)*dt;
        k4y = g(x+k3x,y+k3y,z+k3z)*dt;
        k4z = h(x+k3x,y+k3y,z+k3z)*dt;
        /* Actualizamos las variables y el tiempo */
        x += (k1x/6.0 + k2x/3.0 + k3x/3.0 + k4x/6.0);
        y += (k1y/6.0 + k2y/3.0 + k3y/3.0 + k4y/6.0);
        z += (k1z/6.0 + k2z/3.0 + k3z/3.0 + k4z/6.0);
        /* finalmente escribimos sobre el archivo */
        fprintf(output,"%f %f %f n",x,y,z);
        fflush(output); 
        usleep(100000);
        fprintf(gp, "replot n");
        fflush(gp);
        t += dt;
    }
    fclose(gp);
    fclose(output);
    return 0;
}

Answer 1

考虑：

gp = popen("nice gnuplot -","w");

这将在一定程度上防止gnuplot破坏您的系统，但您是否想过在启动gnuplot之前运行runge-kutta代码？将输出写入一个文件，然后在后台作业中将该文件馈送到gnuplot。

例如，这一切都集成在一个脚本中，如果你想的话，你可以很好地和/或在一起使用。

runge-kutta > rkfile.dat
cat rkfile.dat | gnuplot - &
wait

例如：

at -k now <<!
 cd /path/to/rk
 nice ./my_rkscript.sh
!

所有这些都是防止工作侵蚀当前互动环境的一些方法。可以说，这样一来，"窗口"就不会被束缚住。示例脚本需要重定向stderr，并根据需要添加一些其他改进。