int main (int argc, char *argv[])
{
int i, j, k;
long s1, s2;
int nodes, cells, cell_size;
if ( argc != 2 ){
// Print usage
printf("usage: %s filename \n", argv[0]);
printf("The filename must be of the form: /Path/To/Meshes/Folder/meshname_without_extension\n");
exit(0);
}
utils_init(CORES);
/*
* Read in 2D mesh informations, coordinates of the vertices and the
* map from trinagles to vertices.
*/
char * node_path = str_cat(argv[1], ".node");
char * cell_path = str_cat(argv[1], ".ele");
double * coords_2D = read_coords_2D(node_path, &nodes);
int * map_2D = read_cell_node_map_2D(cell_path, &cells, &cell_size);
free(node_path);
free(cell_path);
/*
* 3D coordinate field.
*/
double * coords_3D = extrude_coords(coords_2D, nodes, LAYERS, LAYER_HEIGHT);
free(coords_2D);
/*
* 3D map from cells to vertices.
*/
int * map_3D = extrude_map(map_2D, cells, cell_size, LAYERS);
int off_3D[6] = {1, 1, 1, 1, 1, 1};
free(map_2D);
/*
* Helmholtz Assembly
*
* Assembly of the LHS and RHS of a Helmholtz Equation.
*/
wrap_init(CORES);
/*
* Evaluate an expression over the mesh.
*
*/
double *expr1 = (double*)malloc(sizeof(double) * nodes * LAYERS);
printf(" Evaluating expression... ");
s1 = stamp();
wrap_expression_1(0, cells,
expr1, map_3D,
coords_3D, map_3D,
off_3D, off_3D, LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
//fprint(expr1, 150, 1);
/*
* Zero an array
*/
double *expr2 = (double*)malloc(sizeof(double) * nodes * LAYERS);
printf(" Set array to zero... ");
s1 = stamp();
wrap_zero_1(0, nodes * LAYERS,
expr2,
LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
/*
* Interpolation operation.
*/
printf(" Interpolate expression... ");
s1 = stamp();
wrap_rhs_1(0, cells,
expr2, map_3D,
coords_3D, map_3D,
expr1, map_3D,
off_3D, off_3D, off_3D, LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
/*
* Another expression kernel
*/
double *expr3 = (double*)malloc(sizeof(double) * nodes * LAYERS);
printf(" Evaluating expression... ");
s1 = stamp();
wrap_expression_2(0, nodes * LAYERS,
expr2,
expr3,
LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
/*
* RHS assembly loop
*/
double *expr4 = (double*)malloc(sizeof(double) * nodes * LAYERS);
printf(" Assembling right-hand side... ");
s1 = stamp();
wrap_rhs(0, cells,
expr4, map_3D,
coords_3D, map_3D,
expr2, map_3D,
expr3, map_3D,
off_3D, off_3D, off_3D, off_3D, LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
/*
* Matrix assembly loop
*/
double *expr5 = (double*)malloc(sizeof(double) * nodes * LAYERS);
printf(" Assembling left-hand side... ");
s1 = stamp();
wrap_lhs(0, cells,
expr5, map_3D, map_3D,
coords_3D, map_3D,
off_3D, off_3D, off_3D, LAYERS);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
/*
* RHS and LHS output
*/
write: s1 = stamp();
printf(" Numerical results written to output files... ");
output(FILE_RHS, expr4, nodes * LAYERS, 1);
output(FILE_LHS, expr5, nodes * LAYERS, 1);
s2 = stamp();
printf("%g s\n", (s2 - s1)/1e9);
free(coords_3D);
free(map_3D);
free(expr1);
free(expr2);
free(expr3);
free(expr4);
free(expr5);
return 0;
}
int __wrap_main(int argc, char *argv[], char *envp[])
{
int c;
char *config_file = NULL;
char *log_file = NULL;
bool daemon = false;
while ((c = getopt(argc, argv, "c:dhl:vV")) != -1) {
switch (c) {
case 'c':
config_file = strdup(optarg);
break;
case 'd':
daemon = true;
break;
case 'h':
usage();
_exit(EXIT_SUCCESS);
case 'l':
log_file = strdup(optarg);
break;
case 'v':
++verbose;
break;
case 'V':
printf("%s\n", openrelease_version());
_exit(EXIT_SUCCESS);
case '?':
usage();
_exit(EXIT_FAILURE);
default:
say_error("unmatched option: -%c", c);
_exit(EXIT_FAILURE);
}
}
create_log(log_file);
config_init(config_file);
if (getenv("OPENRELEASE_STAGE2") == NULL) {
if (daemon) {
if (daemonize(1, 0) == -1) {
say_error("can't daemonize");
_exit(EXIT_FAILURE);
}
}
if (stdio_wrap(config.input, config.output) == -1) {
say_error("stdio_wrap failed");
_exit(EXIT_FAILURE);
}
putenv("OPENRELEASE_STAGE2=1");
execvp(argv[0], argv);
}
mmaps_init();
symfile_load(config.symfile);
wrap_init();
debug_init();
key_action_init();
say_info("dive into RELEASE");
unsetenv("LD_PRELOAD");
argc = 1;
if (optind != argc) {
memmove(&argv[1], &argv[optind], (optind - argc) * sizeof(argv[0]));
argc += optind - argc;
}
argv[argc] = NULL;
__real_main(argc, argv, envp);
_exit(EXIT_SUCCESS);
return 0;
}
