visualizer.c

Go to the documentation of this file.

/*****************************************************************************
 * Copyright (c) 2019 FrontISTR Commons
 * This software is released under the MIT License, see LICENSE.txt
 *****************************************************************************/
 
#include <stdio.h>
#include <stdlib.h>
#include "hecmw_struct.h"
#include "hecmw_result.h"
#include "hecmw_util.h"
#include "hecmw_io.h"
#include "hecmw_visualizer.h"
 
extern PSF_link *psf;
extern PVR_link *pvr;
 
int main(int argc, char **argv) {
  int i;
  struct hecmwST_local_mesh *mesh;
  struct hecmwST_result_data *data;
  char *resultfile, resultfile1[HECMW_FILENAME_LEN];
  int min_step, max_step, interval, timestep;
  PSF_link *tp1;
  PVR_link *tv1;
  int fg_text;
  int mynode;
 
  if (HECMW_init(&argc, &argv)) abort();
 
  mesh = HECMW_get_mesh("fstrMSH");
  if (mesh == NULL) HECMW_abort(HECMW_comm_get_comm());
 
  HECMW_Comm_rank(mesh->HECMW_COMM, &mynode);
 
  HECMW_visualize_init();
 
  min_step = 100000000;
  max_step = -10000000;
  tp1      = psf->next_psf;
  tv1      = pvr->next_pvr;
  for (i = 0; i < psf->num_of_psf; i++) {
    if ((tp1->visual_start_step != -1) && (tp1->visual_start_step < min_step))
      min_step = tp1->visual_start_step;
    if ((tp1->visual_end_step != -1) && (tp1->visual_end_step > max_step))
      max_step = tp1->visual_end_step;
    interval   = tp1->visual_interval_step;
    tp1        = tp1->next_psf;
  }
  for (i = 0; i < pvr->num_of_pvr; i++) {
    if ((tv1->visual_start_step != -1) && (tv1->visual_start_step < min_step))
      min_step = tv1->visual_start_step;
    if ((tv1->visual_end_step != -1) && (tv1->visual_end_step > max_step))
      max_step = tv1->visual_end_step;
    interval   = tp1->visual_interval_step;
    tv1        = tv1->next_pvr;
  }
  if ((min_step == 100000000) && (max_step == -10000000)) {
    min_step = 1;
    max_step = 1;
    interval = 1;
  }
 
  for (timestep = min_step; timestep <= max_step; timestep++) {
    if ((timestep % interval) != 0 && timestep != max_step) continue;
 
    resultfile = HECMW_ctrl_get_result_fileheader("fstrRES", timestep,
                                                  &fg_text);
    sprintf(resultfile1, "%s.%d.%d", resultfile, mynode, timestep);
    data = HECMW_result_read_by_fname(resultfile1);
    if (data == NULL) HECMW_abort(HECMW_comm_get_comm());
 
    HECMW_visualize(mesh, data, timestep);
 
    HECMW_result_free(data);
  }
 
  HECMW_visualize_finalize();
 
  HECMW_dist_free(mesh);
 
  HECMW_finalize();
 
  return 0;
}