hecmw_result_copy_f2c.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 <string.h>
#include <errno.h>
#include "hecmw_config.h"
#include "hecmw_result.h"
#include "hecmw_util.h"
 
static struct hecmwST_result_data *Result;
static int NNode, NElem;
 
/*-----------------------------------------------------------------------------
 * SetFunc
 */
 
static int set_ng_component(void *src) {
  Result->ng_component = *((int *)src);
  return 0;
}
 
static int set_nn_component(void *src) {
  Result->nn_component = *((int *)src);
  return 0;
}
 
static int set_ne_component(void *src) {
  Result->ne_component = *((int *)src);
  return 0;
}
 
static int set_ng_dof(void *src) {
  int size;
 
  if (Result->ng_component <= 0) return 0;
  size           = sizeof(*Result->ng_dof) * Result->ng_component;
  Result->ng_dof = HECMW_malloc(size);
  if (Result->ng_dof == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->ng_dof, src, size);
  return 0;
}
 
static int set_nn_dof(void *src) {
  int size;
 
  if (Result->nn_component <= 0) return 0;
  size           = sizeof(*Result->nn_dof) * Result->nn_component;
  Result->nn_dof = HECMW_malloc(size);
  if (Result->nn_dof == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->nn_dof, src, size);
  return 0;
}
 
static int set_ne_dof(void *src) {
  int size;
 
  if (Result->ne_component <= 0) return 0;
  size           = sizeof(*Result->ne_dof) * Result->ne_component;
  Result->ne_dof = HECMW_malloc(size);
  if (Result->ne_dof == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->ne_dof, src, size);
  return 0;
}
 
static int set_global_label(void *src) {
  int i;
 
  if (Result->ng_component <= 0) return 0;
 
  Result->global_label =
      HECMW_malloc(sizeof(*Result->global_label) * Result->ng_component);
  if (Result->global_label == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  for (i = 0; i < Result->ng_component; i++) {
    char *src_point       = (char *)src + HECMW_NAME_LEN * i;
    Result->global_label[i] = HECMW_strcpy_f2c(src_point, HECMW_NAME_LEN);
  }
 
  return 0;
}
 
static int set_node_label(void *src) {
  int i;
 
  if (Result->nn_component <= 0) return 0;
 
  Result->node_label =
      HECMW_malloc(sizeof(*Result->node_label) * Result->nn_component);
  if (Result->node_label == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  for (i = 0; i < Result->nn_component; i++) {
    char *src_point       = (char *)src + HECMW_NAME_LEN * i;
    Result->node_label[i] = HECMW_strcpy_f2c(src_point, HECMW_NAME_LEN);
  }
 
  return 0;
}
 
static int set_elem_label(void *src) {
  int i;
 
  if (Result->ne_component <= 0) return 0;
 
  Result->elem_label =
      HECMW_malloc(sizeof(*Result->elem_label) * Result->ne_component);
  if (Result->elem_label == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  for (i = 0; i < Result->ne_component; i++) {
    char *src_point       = (char *)src + HECMW_NAME_LEN * i;
    Result->elem_label[i] = HECMW_strcpy_f2c(src_point, HECMW_NAME_LEN);
  }
 
  return 0;
}
 
static int set_global_val_item(void *src) {
  int i, size;
  int n = 0;
 
  if (Result->ng_component <= 0) return 0;
  for (i = 0; i < Result->ng_component; i++) {
    n += Result->ng_dof[i];
  }
  size                  = sizeof(*Result->global_val_item) * n;
  Result->global_val_item = HECMW_malloc(size);
  if (Result->global_val_item == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->global_val_item, src, size);
  return 0;
}
 
static int set_node_val_item(void *src) {
  int i, size;
  int n = 0;
 
  if (Result->nn_component <= 0) return 0;
  for (i = 0; i < Result->nn_component; i++) {
    n += Result->nn_dof[i];
  }
  size                  = sizeof(*Result->node_val_item) * n * NNode;
  Result->node_val_item = HECMW_malloc(size);
  if (Result->node_val_item == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->node_val_item, src, size);
  return 0;
}
 
static int set_elem_val_item(void *src) {
  int i, size;
  int n = 0;
 
  if (Result->ne_component <= 0) return 0;
  for (i = 0; i < Result->ne_component; i++) {
    n += Result->ne_dof[i];
  }
  size                  = sizeof(*Result->elem_val_item) * n * NElem;
  Result->elem_val_item = HECMW_malloc(size);
  if (Result->elem_val_item == NULL) {
    HECMW_set_error(errno, "");
    return -1;
  }
  memcpy(Result->elem_val_item, src, size);
  return 0;
}
 
/*-----------------------------------------------------------------------------
 * SetFunc table
 */
 
typedef int (*SetFunc)(void *);
 
static struct func_table {
  char *struct_name;
  char *var_name;
  SetFunc set_func;
} functions[] = {
    /*  { Struct name, Variable name, memcpy function } */
    {"hecmwST_result_data", "ng_component", set_ng_component},
    {"hecmwST_result_data", "nn_component", set_nn_component},
    {"hecmwST_result_data", "ne_component", set_ne_component},
    {"hecmwST_result_data", "ng_dof", set_ng_dof},
    {"hecmwST_result_data", "nn_dof", set_nn_dof},
    {"hecmwST_result_data", "ne_dof", set_ne_dof},
    {"hecmwST_result_data", "global_label", set_global_label},
    {"hecmwST_result_data", "node_label", set_node_label},
    {"hecmwST_result_data", "elem_label", set_elem_label},
    {"hecmwST_result_data", "global_val_item", set_global_val_item},
    {"hecmwST_result_data", "node_val_item", set_node_val_item},
    {"hecmwST_result_data", "elem_val_item", set_elem_val_item},
};
 
static const int NFUNC = sizeof(functions) / sizeof(functions[0]);
 
static SetFunc get_set_func(char *struct_name, char *var_name) {
  int i;
 
  for (i = 0; i < NFUNC; i++) {
    if (strcmp(functions[i].struct_name, struct_name) == 0 &&
        strcmp(functions[i].var_name, var_name) == 0) {
      return functions[i].set_func;
    }
  }
  return NULL;
}
 
/*----------------------------------------------------------------------------*/
 
int HECMW_result_copy_f2c_init(struct hecmwST_result_data *result_data,
                               int n_node, int n_elem) {
  Result = result_data;
  NNode  = n_node;
  NElem  = n_elem;
  return 0;
}
 
int HECMW_result_copy_f2c_finalize(void) {
  Result = NULL;
  return 0;
}
 
/*----------------------------------------------------------------------------*/
 
void hecmw_result_copy_f2c_set_if(char *struct_name, char *var_name, void *src,
                                  int *err, int slen, int vlen) {
  SetFunc func;
  char sname[HECMW_NAME_LEN + 1];
  char vname[HECMW_NAME_LEN + 1];
 
  *err = 1;
 
  if (Result == NULL) {
    HECMW_set_error(
        HECMW_ALL_E0102,
        "hecmw_result_copy_f2c_set_if(): 'result' has not initialized yet");
    return;
  }
  if (struct_name == NULL) {
    HECMW_set_error(HECMW_ALL_E0101,
                    "hecmw_result_copy_f2c_set_if(): 'sname' is NULL");
    return;
  }
  if (var_name == NULL) {
    HECMW_set_error(HECMW_ALL_E0101,
                    "hecmw_result_copy_f2c_set_if(): 'vname' is NULL");
    return;
  }
 
  if (HECMW_strcpy_f2c_r(struct_name, slen, sname, sizeof(sname)) == NULL) {
    return;
  }
 
  if (HECMW_strcpy_f2c_r(var_name, vlen, vname, sizeof(vname)) == NULL) {
    return;
  }
 
  func = get_set_func(sname, vname);
  if (func == NULL) {
    HECMW_set_error(HECMW_ALL_E0102,
                    "hecmw_result_copy_f2c_set_if(): SetFunc not found");
    return;
  }
 
  if ((*func)(src)) {
    return;
  }
 
  *err = 0;
}
 
/*----------------------------------------------------------------------------*/
 
void hecmw_result_write_st_init_if(int *err) {
  *err   = 1;
  Result = HECMW_calloc(1, sizeof(*Result));
  if (Result == NULL) {
    HECMW_set_error(errno, "");
    return;
  }
 
  NNode = HECMW_result_get_nnode();
  NElem = HECMW_result_get_nelem();
  *err  = 0;
}
 
/*----------------------------------------------------------------------------*/
 
void hecmw_result_write_st_finalize_if(int *err) {
  HECMW_result_free(Result);
 
  *err = 0;
}
 
/*----------------------------------------------------------------------------*/
 
void hecmw_result_write_st_by_name_if(char *name_ID, int *err, int len) {
  char name_ID_str[HECMW_NAME_LEN + 1];
  char head[HECMW_HEADER_LEN + 1];
  char comment[HECMW_MSG_LEN + 1];
 
  *err = 1;
 
  if (HECMW_strcpy_f2c_r(name_ID, len, name_ID_str, sizeof(name_ID_str)) ==
      NULL)
    return;
 
  HECMW_result_get_header(head);
  HECMW_result_get_comment(comment);
  if (HECMW_result_write_ST_by_name(name_ID_str, Result, NNode, NElem, head, comment))
    return;
 
  *err = 0;
}