hecmw_io_get_mesh.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_util.h"
#include "hecmw_io_get_mesh.h"
#include "hecmw_io_mesh.h"
#include "hecmw_io_hec.h"
#include "hecmw_io_geofem.h"
#include "hecmw_io_inp.h"
#include "hecmw_io_dist.h"
#include "hecmw_dist_refine.h"
#include "hecmw_dist_free.h"
 
static struct hecmwST_local_mesh *get_entire_mesh(
    struct hecmw_ctrl_meshfiles *files) {
  int i;
  struct hecmwST_local_mesh *mesh;
 
  if (HECMW_io_init()) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "io_init done");
 
  if (HECMW_io_pre_process()) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "io_pre_process done");
 
  for (i = 0; i < files->n_mesh; i++) {
    struct hecmw_ctrl_meshfile *file = &files->meshfiles[i];
 
    switch (file->type) {
      case HECMW_CTRL_FTYPE_HECMW_ENTIRE:
        if (HECMW_read_entire_mesh(file->filename)) return NULL;
        break;
      case HECMW_CTRL_FTYPE_GEOFEM:
        if (HECMW_read_geofem_mesh(file->filename)) return NULL;
        break;
      case HECMW_CTRL_FTYPE_INP:
        if (HECMW_read_inp_mesh(file->filename)) return NULL;
        break;
      default:
        HECMW_assert(0);
    }
  }
  HECMW_log(HECMW_LOG_DEBUG, "reading mesh done\n");
 
  if (HECMW_io_post_process()) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "post_process done\n");
 
  if (HECMW_io_create_dummy_nodes()) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "creating rotational dummy nodes done\n");
 
  mesh = HECMW_io_make_local_mesh();
  if (mesh == NULL) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "converting mesh done\n");
 
  if (HECMW_io_finalize()) return NULL;
  HECMW_log(HECMW_LOG_DEBUG, "io_finalize done\n");
 
  return mesh;
}
 
struct hecmwST_local_mesh *HECMW_get_mesh(char *name_ID) {
  struct hecmw_ctrl_meshfiles *files;
  struct hecmwST_local_mesh *mesh;
  char filename[HECMW_FILENAME_LEN + 1];
  char *cad_filename;
  FILE *fp;
 
  files = HECMW_ctrl_get_meshfiles(name_ID);
  if (files == NULL) return NULL;
 
  if (files->n_mesh == 1 &&
      files->meshfiles[0].type == HECMW_CTRL_FTYPE_HECMW_DIST) {
    snprintf(filename, sizeof(filename), "%s", files->meshfiles[0].filename);
    mesh = HECMW_get_dist_mesh(filename);
  } else {
    mesh = get_entire_mesh(files);
  }
 
  char *dot;
  snprintf(filename, sizeof(filename), "%s", files->meshfiles[0].filename);
  dot = strchr(filename, '.');
  if (dot != NULL) *dot = '\0';
  strncat(filename, ".rnf", sizeof(filename) - strlen(filename) - 1);
  if ((fp = fopen(filename, "r")) == NULL) {
    cad_filename = NULL;
  } else {
    fclose(fp);
    cad_filename = filename;
  }
 
  if (HECMW_dist_refine(&mesh, files->meshfiles[0].refine, cad_filename,
                        NULL) != HECMW_SUCCESS) {
    HECMW_dist_free(mesh);
    return NULL;
  }
 
  HECMW_ctrl_free_meshfiles(files);
 
  return mesh;
}