gmi.h

Go to the documentation of this file.

/****************************************************************************** 

  Copyright 2014 Scientific Computation Research Center, 
      Rensselaer Polytechnic Institute. All rights reserved.
  
  This work is open source software, licensed under the terms of the
  BSD license as described in the LICENSE file in the top-level directory.

*******************************************************************************/
#ifndef GMI_H
#define GMI_H

#include <stdio.h>

struct gmi_ent;

struct gmi_iter;

struct gmi_model;

struct gmi_set {
  int n;
  struct gmi_ent* e[1];
};

struct gmi_model_ops {
  struct gmi_iter* (*begin)(struct gmi_model* m, int dim);
  struct gmi_ent* (*next)(struct gmi_model* m, struct gmi_iter* i);
  void (*end)(struct gmi_model* m, struct gmi_iter* i);
  int (*dim)(struct gmi_model* m, struct gmi_ent* e);
  int (*tag)(struct gmi_model* m, struct gmi_ent* e);
  struct gmi_ent* (*find)(struct gmi_model* m, int dim, int tag);
  struct gmi_set* (*adjacent)(struct gmi_model* m, struct gmi_ent* e, int dim);
  void (*eval)(struct gmi_model* m, struct gmi_ent* e,
      double const p[2], double x[3]);
  void (*reparam)(struct gmi_model* m, struct gmi_ent* from,
      double const from_p[2], struct gmi_ent* to, double to_p[2]);
  int (*periodic)(struct gmi_model* m, struct gmi_ent* e, int dim);
  void (*range)(struct gmi_model* m, struct gmi_ent* e, int dim,
      double r[2]);
  void (*closest_point)(struct gmi_model* m, struct gmi_ent* e,
      double const from[3], double to[3], double to_p[2]);
  void (*normal)(struct gmi_model* m, struct gmi_ent* e,
    double const p[2], double n[3]);
  void (*first_derivative)(struct gmi_model* m, struct gmi_ent* e,
    double const p[2], double t0[3], double t1[3]);
  int (*is_point_in_region)(struct gmi_model* m, struct gmi_ent* e,
      double p[3]);
  void (*bbox)(struct gmi_model* m, struct gmi_ent* e,
      double bmin[3], double bmax[3]);
  int (*is_in_closure_of)(struct gmi_model* m, struct gmi_ent* e,
      struct gmi_ent* et);
  int (*is_discrete_ent)(struct gmi_model* m, struct gmi_ent* e);
  void (*destroy)(struct gmi_model* m);
};

struct gmi_model {
  struct gmi_model_ops const* ops;
  int n[4];
};

typedef struct gmi_model* (*gmi_creator)(const char* filename);

#ifdef __cplusplus
extern "C" {
#endif

struct gmi_set* gmi_make_set(int n);
void gmi_free_set(struct gmi_set* s);

struct gmi_iter* gmi_begin(struct gmi_model* m, int dim);
struct gmi_ent* gmi_next(struct gmi_model* m, struct gmi_iter* i);
void gmi_end(struct gmi_model* m, struct gmi_iter* i);
int gmi_dim(struct gmi_model* m, struct gmi_ent* e);
int gmi_tag(struct gmi_model* m, struct gmi_ent* e);
struct gmi_ent* gmi_find(struct gmi_model* m, int dim, int tag);
struct gmi_set* gmi_adjacent(struct gmi_model* m, struct gmi_ent* e, int dim);
int gmi_can_eval(struct gmi_model* m);
int gmi_can_get_closest_point(struct gmi_model* m);
int gmi_has_normal(struct gmi_model* m);
void gmi_eval(struct gmi_model* m, struct gmi_ent* e,
    double const p[2], double x[3]);
void gmi_reparam(struct gmi_model* m, struct gmi_ent* from,
    double const from_p[2], struct gmi_ent* to, double to_p[2]);
int gmi_periodic(struct gmi_model* m, struct gmi_ent* e, int dim);
void gmi_range(struct gmi_model* m, struct gmi_ent* e, int dim,
    double r[2]);
void gmi_closest_point(struct gmi_model* m, struct gmi_ent* e,
    double const from[3], double to[3], double to_p[2]);
void gmi_normal(struct gmi_model* m, struct gmi_ent* e,
    double const p[2], double n[3]);
void gmi_first_derivative(struct gmi_model* m, struct gmi_ent* e,
    double const p[2], double t0[3], double t1[3]);
int gmi_is_point_in_region(struct gmi_model* m, struct gmi_ent* e,
    double point[3]);
int gmi_is_in_closure_of(struct gmi_model* m, struct gmi_ent* e,
    struct gmi_ent* et);
void gmi_bbox(struct gmi_model* m, struct gmi_ent* e,
    double bmin[3], double bmax[3]);
int gmi_is_discrete_ent(struct gmi_model* m, struct gmi_ent* e);
void gmi_destroy(struct gmi_model* m);

void gmi_register(gmi_creator f, const char* ext);
struct gmi_model* gmi_load(const char* filename);

void gmi_write_dmg(struct gmi_model* m, const char* filename);

void gmi_fail(const char* why) __attribute__((noreturn));

void gmi_fscanf(FILE* f, int n, const char* format, ...);

int gmi_getline(char** line, size_t* cap, FILE* f);
int gmi_has_ext(const char* filename, const char* ext);

#ifdef __cplusplus
}
#endif

#endif