Гнуплот и усли на С - PullRequest
Новая
       25

Гнуплот и усли на С

2 голосов
/ 23 февраля 2012

Я пытаюсь сделать сюжет в реальном времени на Linux.Все работает хорошо, просто когда программа заканчивается, некоторое время загружается процессор.Кроме того, из-за того, что я использую очень короткое время (например, 0,3 с), я не знаю, почему в оболочке bash преобладает интерфейс gnuplot, даже для других окон bash я пытаюсь открыть, поэтому я не могу использовать ctrl + cостановите программу.

Есть ли какая-либо другая функция, которую я мог бы использовать вместо usleep?Или что-то, что я мог бы исправить, чтобы избежать доминирования gnuplot над окнами?

Заранее спасибо и извините за мой плохой английский

Я даю код здесь 

#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;
}

1 Ответ

4 голосов
/ 23 февраля 2012

учтите: 

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

не позволит gnuplot немного съесть вашу систему, но задумывались ли вы о запуске кода runge-kutta перед тем, как даже запустить gnuplot?Записать вывод в файл, а затем передать его в gnuplot в фоновом режиме.

Например, все это происходит в одном скрипте, который вы можете использовать и / или использовать в котором, если хотите.*

в примере: 

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

Все это - некоторые из способов не дать работе съесть вашу текущую интерактивную среду.Таким образом, «окно», так сказать, не связано.Скрипт примера требует перенаправления stderr и некоторых других улучшений, которые вы добавляете, если хотите.

...