| Exec Time (%) | Codelet Name | Error (%) |
|---|---|---|
| 93.2 | main_main_39 | 1.7 |
| 92.6 | lbm_LBM_performStreamCollide_186 | 8.8 |
| 2.1 | lbm_LBM_showGridStatistics_459 | 3.1 |
| 2.2 | lbm_LBM_initializeGrid_62 | 8.7 |
Exec Time : 6.408211e+09 cycles
Matching : 96.9%
__invivo__lbm_LBM_initializeGrid_62
Real time: 1.461459e+08 -- Predicted time: 1.601652e+08
Callcount: 2
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 1 | 8.014368e+07 | 7.197408e+07 | 10.1 |
| 2 | ### 2 | 1 | 8.002151e+07 | 7.347708e+07 | 8.1 |
1
/* $Id: main.c,v 1.4 2004/04/21 04:23:43 pohlt Exp $ */2
​3
/*############################################################################*/4
​5
6
7
8
9
​10
11
12
13
14
15
16
​17
18
​19
/*############################################################################*/20
​21
static LBM_GridPtr srcGrid, dstGrid;22
​23
/*############################################################################*/24
​25
int main( int nArgs, char* arg[] ) {26
MAIN_Param param;27
28
MAIN_Time time;29
30
int t;31
​32
MAIN_parseCommandLine( nArgs, arg, ¶m );33
MAIN_printInfo( ¶m );34
MAIN_initialize( ¶m );35
36
MAIN_startClock( &time );37
38
​39
for( t = 1; t <= param.nTimeSteps; t++ ) {40
if( param.simType == CHANNEL ) {41
LBM_handleInOutFlow( *srcGrid );42
}43
​44
LBM_performStreamCollide( *srcGrid, *dstGrid );45
LBM_swapGrids( &srcGrid, &dstGrid );46
​47
if( (t & 63) == 0 ) {48
printf( "timestep: %i\n", t );49
LBM_showGridStatistics( *srcGrid );50
}51
}52
​53
54
MAIN_stopClock( &time, ¶m );55
56
MAIN_finalize( ¶m );57
​58
return 0;59
}60
​61
/*############################################################################*/62
​63
void MAIN_parseCommandLine( int nArgs, char* arg[], MAIN_Param* param ) {64
struct stat fileStat;65
66
if( nArgs < 5 || nArgs > 6 ) {67
printf( "syntax: lbm <time steps> <result file> <0: nil, 1: cmp, 2: str> <0: ldc, 1: channel flow> [<obstacle file>]\n" );68
exit( 1 );69
}70
​71
param->nTimeSteps = atoi( arg[1] );72
param->resultFilename = arg[2];73
param->action = (MAIN_Action) atoi( arg[3] );74
param->simType = (MAIN_SimType) atoi( arg[4] );75
​76
if( nArgs == 6 ) {77
param->obstacleFilename = arg[5];78
​79
if( stat( param->obstacleFilename, &fileStat ) != 0 ) {80
printf( "MAIN_parseCommandLine: cannot stat obstacle file '%s'\n",81
param->obstacleFilename );82
exit( 1 );83
}84
if( fileStat.st_size != SIZE_X*SIZE_Y*SIZE_Z+(SIZE_Y+1)*SIZE_Z ) {85
printf( "MAIN_parseCommandLine:\n"86
"\tsize of file '%s' is %i bytes\n"87
"\texpected size is %i bytes\n",88
param->obstacleFilename, (int) fileStat.st_size,89
SIZE_X*SIZE_Y*SIZE_Z+(SIZE_Y+1)*SIZE_Z );90
exit( 1 );91
}92
}93
else param->obstacleFilename = NULL;94
​95
if( param->action == COMPARE &&96
stat( param->resultFilename, &fileStat ) != 0 ) {97
printf( "MAIN_parseCommandLine: cannot stat result file '%s'\n",98
param->resultFilename );99
exit( 1 );100
}101
}102
​103
/*############################################################################*/104
​105
void MAIN_printInfo( const MAIN_Param* param ) {106
const char actionString[3][32] = {"nothing", "compare", "store"};107
const char simTypeString[3][32] = {"lid-driven cavity", "channel flow"};108
printf( "MAIN_printInfo:\n"109
"\tgrid size : %i x %i x %i = %.2f * 10^6 Cells\n"110
"\tnTimeSteps : %i\n"111
"\tresult file : %s\n"112
"\taction : %s\n"113
"\tsimulation type: %s\n"114
"\tobstacle file : %s\n\n",115
SIZE_X, SIZE_Y, SIZE_Z, 1e-6*SIZE_X*SIZE_Y*SIZE_Z,116
param->nTimeSteps, param->resultFilename, 117
actionString[param->action], simTypeString[param->simType],118
(param->obstacleFilename == NULL) ? "<none>" :119
param->obstacleFilename );120
}121
​122
/*############################################################################*/123
​124
void MAIN_initialize( const MAIN_Param* param ) {125
LBM_allocateGrid( (double**) &srcGrid );126
LBM_allocateGrid( (double**) &dstGrid );127
​128
LBM_initializeGrid( *srcGrid );129
LBM_initializeGrid( *dstGrid );130
​131
if( param->obstacleFilename != NULL ) {132
LBM_loadObstacleFile( *srcGrid, param->obstacleFilename );133
LBM_loadObstacleFile( *dstGrid, param->obstacleFilename );134
}135
​136
if( param->simType == CHANNEL ) {137
LBM_initializeSpecialCellsForChannel( *srcGrid );138
LBM_initializeSpecialCellsForChannel( *dstGrid );139
}140
else {141
LBM_initializeSpecialCellsForLDC( *srcGrid );142
LBM_initializeSpecialCellsForLDC( *dstGrid );143
}144
​145
LBM_showGridStatistics( *srcGrid );146
}147
​148
/*############################################################################*/149
​150
void MAIN_finalize( const MAIN_Param* param ) {151
LBM_showGridStatistics( *srcGrid );152
​153
if( param->action == COMPARE )154
LBM_compareVelocityField( *srcGrid, param->resultFilename, TRUE );155
if( param->action == STORE )156
LBM_storeVelocityField( *srcGrid, param->resultFilename, TRUE );157
​158
LBM_freeGrid( (double**) &srcGrid );159
LBM_freeGrid( (double**) &dstGrid );160
}161
​162
163
/*############################################################################*/164
​165
void MAIN_startClock( MAIN_Time* time ) {166
time->timeScale = 1.0 / sysconf( _SC_CLK_TCK );167
time->tickStart = times( &(time->timeStart) );168
}169
​170
​171
/*############################################################################*/172
​173
void MAIN_stopClock( MAIN_Time* time, const MAIN_Param* param ) {174
time->tickStop = times( &(time->timeStop) );175
​176
printf( "MAIN_stopClock:\n"177
"\tusr: %7.2f sys: %7.2f tot: %7.2f wct: %7.2f MLUPS: %5.2f\n\n",178
(time->timeStop.tms_utime - time->timeStart.tms_utime) * time->timeScale,179
(time->timeStop.tms_stime - time->timeStart.tms_stime) * time->timeScale,180
(time->timeStop.tms_utime - time->timeStart.tms_utime +181
time->timeStop.tms_stime - time->timeStart.tms_stime) * time->timeScale,182
(time->tickStop - time->tickStart ) * time->timeScale,183
1.0e-6 * SIZE_X * SIZE_Y * SIZE_Z * param->nTimeSteps /184
(time->tickStop - time->tickStart ) / time->timeScale );185
}186
One invocation only
__invivo__main_main_39
Real time: 5.974737e+09 -- Predicted time: 6.083370e+09
Callcount: 1
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 1 | 6.083370e+09 | 5.936650e+09 | 2.4 |
__invivo__lbm_LBM_performStreamCollide_186
Real time: 5.934739e+09 -- Predicted time: 5.408247e+09
Callcount: 20
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 7 | ### 1 | 11.0001831267787 | 2.723869e+08 | 2.905747e+08 | 6.2 |
| 14 | ### 2 | 8.98596177740608 | 2.684121e+08 | 2.984825e+08 | 10.0 |
__invivo__lbm_LBM_showGridStatistics_459
Real time: 1.357075e+08 -- Predicted time: 1.401514e+08
Callcount: 2
| Invocation | Cluster | Part | Invitro (cycles) | Invivo (cycles) | Error (%) |
|---|---|---|---|---|---|
| 1 | ### 1 | 1 | 7.025516e+07 | 6.769733e+07 | 3.6 |
| 2 | ### 2 | 1 | 6.989629e+07 | 6.787814e+07 | 2.8 |