OpenAD/code_heat_transport_MPI/timeRatio.c

#include <sys/time.h>
#include <stdio.h>

void timeratio_() { 
  /* 
     three invocations produce
     the ratio of the timedifference of 
     the 3rd vs. the 2nd invocation over
     the difference of the 2nd vs. the 1st invocation
  */
  static short stage=0;
#define tvCount 3
  static struct timeval tv_a[tvCount];
  double deltaThis, deltaPrevious;
  if (stage<tvCount) 
    gettimeofday(&(tv_a[stage]),0);
  printf("JU: TIMING: stamp %i: %i.%06i\n",
         stage,
         tv_a[stage].tv_sec,
         tv_a[stage].tv_usec);
  if (stage>0) { 
    printf("JU: TIMING: delta stamps %i-%i: %i.%06i\n",
           stage,
           stage-1,
           (tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?tv_a[stage].tv_sec-tv_a[stage-1].tv_sec-1:tv_a[stage].tv_sec-tv_a[stage-1].tv_sec,
           (tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?1000000-(tv_a[stage-1].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[stage-1].tv_usec);
  }
  if (stage>1) { 
    printf("JU: TIMING: delta stamps %i-%i: %i.%06i\n",
           stage,
           0,
           (tv_a[stage].tv_usec<tv_a[0].tv_usec)?tv_a[stage].tv_sec-tv_a[0].tv_sec-1:tv_a[stage].tv_sec-tv_a[0].tv_sec,
           (tv_a[stage].tv_usec<tv_a[0].tv_usec)?1000000-(tv_a[0].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[0].tv_usec);
  }
  if (stage==tvCount-1 && stage>1) { 
    deltaThis=(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?tv_a[stage].tv_sec-tv_a[stage-1].tv_sec-1:tv_a[stage].tv_sec-tv_a[stage-1].tv_sec;
    deltaThis*=1.0e6;
    deltaThis+=(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?1000000-(tv_a[stage-1].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[stage-1].tv_usec;
    deltaPrevious=(tv_a[stage-1].tv_usec<tv_a[stage-2].tv_usec)?tv_a[stage-1].tv_sec-tv_a[stage-2].tv_sec-1:tv_a[stage-1].tv_sec-tv_a[stage-2].tv_sec;
    deltaPrevious*=1.0e6;
    deltaPrevious+=(tv_a[stage-1].tv_usec<tv_a[stage-2].tv_usec)?1000000-(tv_a[stage-2].tv_usec-tv_a[stage-1].tv_usec):tv_a[stage-1].tv_usec-tv_a[stage-2].tv_usec;
    printf("JU: TIMING: ratio stamps (%i-%i)/(%i-%i): %e/%e=%e\n",
           stage,
           stage-1,
           stage-1,
           stage-2,
           deltaThis,
           deltaPrevious,
           deltaThis/deltaPrevious);
    stage=0;  
  }
  else
    stage++;
}
1	#include <sys/time.h>
2	#include <stdio.h>
3
4	void timeratio_() {
5	/*
6	three invocations produce
7	the ratio of the timedifference of
8	the 3rd vs. the 2nd invocation over
9	the difference of the 2nd vs. the 1st invocation
10	*/
11	static short stage=0;
12	#define tvCount 3
13	static struct timeval tv_a[tvCount];
14	double deltaThis, deltaPrevious;
15	if (stage<tvCount)
16	gettimeofday(&(tv_a[stage]),0);
17	printf("JU: TIMING: stamp %i: %i.%06i\n",
18	stage,
19	tv_a[stage].tv_sec,
20	tv_a[stage].tv_usec);
21	if (stage>0) {
22	printf("JU: TIMING: delta stamps %i-%i: %i.%06i\n",
23	stage,
24	stage-1,
25	(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?tv_a[stage].tv_sec-tv_a[stage-1].tv_sec-1:tv_a[stage].tv_sec-tv_a[stage-1].tv_sec,
26	(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?1000000-(tv_a[stage-1].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[stage-1].tv_usec);
27	}
28	if (stage>1) {
29	printf("JU: TIMING: delta stamps %i-%i: %i.%06i\n",
30	stage,
31	0,
32	(tv_a[stage].tv_usec<tv_a[0].tv_usec)?tv_a[stage].tv_sec-tv_a[0].tv_sec-1:tv_a[stage].tv_sec-tv_a[0].tv_sec,
33	(tv_a[stage].tv_usec<tv_a[0].tv_usec)?1000000-(tv_a[0].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[0].tv_usec);
34	}
35	if (stage==tvCount-1 && stage>1) {
36	deltaThis=(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?tv_a[stage].tv_sec-tv_a[stage-1].tv_sec-1:tv_a[stage].tv_sec-tv_a[stage-1].tv_sec;
37	deltaThis*=1.0e6;
38	deltaThis+=(tv_a[stage].tv_usec<tv_a[stage-1].tv_usec)?1000000-(tv_a[stage-1].tv_usec-tv_a[stage].tv_usec):tv_a[stage].tv_usec-tv_a[stage-1].tv_usec;
39	deltaPrevious=(tv_a[stage-1].tv_usec<tv_a[stage-2].tv_usec)?tv_a[stage-1].tv_sec-tv_a[stage-2].tv_sec-1:tv_a[stage-1].tv_sec-tv_a[stage-2].tv_sec;
40	deltaPrevious*=1.0e6;
41	deltaPrevious+=(tv_a[stage-1].tv_usec<tv_a[stage-2].tv_usec)?1000000-(tv_a[stage-2].tv_usec-tv_a[stage-1].tv_usec):tv_a[stage-1].tv_usec-tv_a[stage-2].tv_usec;
42	printf("JU: TIMING: ratio stamps (%i-%i)/(%i-%i): %e/%e=%e\n",
43	stage,
44	stage-1,
45	stage-1,
46	stage-2,
47	deltaThis,
48	deltaPrevious,
49	deltaThis/deltaPrevious);
50	stage=0;
51	}
52	else
53	stage++;
54	}