1

私は非常に重くて遅い科学シミュレーションを行っていますが、ユーザーがキーを押したときにプログラムにそのステータスに関する情報を表示させることができるかどうか疑問に思いました。

Linuxの場合だけ、移植可能である必要はありません。

現在、私はこの種のコードを持っています：

main()
{
    //SOME GLOBAL INITIALIZATIONS

    //this loop is parallelized
    for(i=0;i<max;i++)
    {
        //Some particular to i initializations

        for(j=0;max2;j++)
        {
            // Here are all the slow and ugly calculations
            // the status should be able to be shown when
            // my program is here
        }
    }
    // Final calculations and final data
}

それが可能かどうかはわかりませんが、コードが停止しているか、何もしていないと思うことがあります。これは非常に役立つ場合があります。

また、これを実行すると同時にompを使用できない場合は、より「手動」の並列化を使用できます。

前もって感謝します！

4

4 に答える 4

1

ncurses の getch() を使用できます

#include <ctype.h>
#include <ncurses.h>
#include <stdio.h>

int keywaiting(int *key) {
  int ch;
  ch = getch();
  if ((ch != ERR) && key) *key = ch;
  return (ch != ERR);
}

int yourprogram(void) {
  int max = 1000, max2 = 1000;
  int i, j;
  int ch;
  /* SOME GLOBAL INITIALIZATIONS */
  /* this loop is parallelized */
  for (i = 0; i < max; i++) {
    /* Some particular to i initializations */
    for (j = 0; j < max2; j++) {
    /* Here are all the slow and ugly calculations */
    /* the status should be able to be shown when */
    /* my program is here */
      printf("%03d%03d\r", i, j);
      if (keywaiting(&ch)) {
        printf("\r\n\r\nAborted @ i = %d; j = %d\r\n", i, j);
        printf("Keypress %d ('%c')\r\n\n", ch, isprint((unsigned char)ch) ? ch : '.');
        printf("Press any key ...");
        refresh();
        i = max;
        j = max2;
      }
    }
  }
  /* Final calculations and final data */
  return 0;
}
int main(void) {
  /* initialize ncurses */
  initscr();
  nodelay(stdscr, TRUE);
  noecho();
  keypad(stdscr, TRUE);
  curs_set(0);

  yourprogram();

  getchar();
  /* done with ncurses */
  clear();
  refresh();
  endwin();

  return 0;
}

curses ライブラリとのリンクを忘れないでください。

于 2013-02-18T19:17:05.983 に答える
1

ライブラリを使用できますpthreads。1 つのスレッドがキーまたは任意のシグナルを待機し、それがトリガーされたときを待ちます。出力変数を表示するために高価な計算を行う必要がある場合は、ロックを使用してシミュレーションの実行を停止し、変数を出力してロックを解放します。それがopenmpの考えでうまくいくかどうかはわかりません。

変数を出力するだけであれば、pthreads と openmp を混在させても問題はないはずです。

于 2013-02-18T19:00:48.873 に答える
0

わかりました、すべての提案に感謝します。欲しかったものに似たものを手に入れることができました。同じではありませんが、役に立ちます (^_^)

コードは次のとおりです。文字を待機するスレッドを生成し、レポート配列を0から1に「切り替え」ます。各ompスレッドが各「j」反復を開始すると、レポートをチェックし、必要に応じてレポートします.

#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
#include <pthread.h>
#include <ncurses.h>

unsigned char *report; // this will store whether a thread has to report or not
int max_threads;

void *check_keyboard(void *whatever)
{
    int c;
    int i;

    newterm(NULL,stdin,stdout); // I want to still see the terminal :)
    noecho();
    cbreak();

    while(1)
    {
        if(getch())
        {
            for(i = 0; i < max_threads; i++)
                report[i]=1;
        }
    }

    return NULL;
}

int main(int argc, const char *argv[])
{
    int i,j;
    int tid;
    pthread_t thread0;

    max_threads = omp_get_max_threads();
    report = malloc(max_threads*sizeof(unsigned char));

    pthread_create(&thread0,NULL,check_keyboard,NULL);

#pragma omp parallel private(i,j,tid) shared(report)
    {
    tid = omp_get_thread_num();

#pragma omp for
    for(i = 0; i < 10; i++)
    {
        for(j = 0; j < 10; j++)
        {
#pragma omp critical
            {
            if(report[tid])
            {
                fprintf(stderr,"Hello World from thread %d\n\r",tid);
                fprintf(stderr,"This is iteration i: %d, j: %d\n\r",i,j);
                report[tid] = 0;
            }
            }
            sleep(1); // This 'is' a test for the long and ugly code
        }
    }
    }

    pthread_cancel(thread0);
    pthread_exit(NULL);
    return 0;
}
于 2013-02-19T19:36:27.437 に答える
0

シグナルを使用する例を次に示します (ただし、部分的に C++ で記述されていますが、シグナル処理はそうではありません)。

//ctrl+c to output status, ctrl+z to quit

//#include <iostream> c++ syntax
#include <stdio.h>
#include <signal.h>
volatile sig_atomic_t showStatus = 0;

void signalHandler(int sig){
    if (sig==SIGINT){
        showStatus=1;
    }
    if (sig==SIGTSTP){
        showStatus=2;
    }
}


int main(int argc, char** argv){
    signal(SIGINT,signalHandler);
    signal(SIGTSTP,signalHandler);
    while(true){
        if (showStatus==1){
            showStatus=0;
            //std::cout << "status" << std::endl; c++ syntax
                    printf("status\n");
        }else if(showStatus>1){
            break;
        }

    }

    return 0;

}
于 2013-02-18T18:18:38.053 に答える