grid.h

Go to the documentation of this file.

#ifndef __SEAM_GRID
#define __SEAM_GRID

/* tolerance parameter - in C++ would get
   from numeric limits. use: two grid locations
   x and y are identified if abs(x-y) < TOL*d.
*/
#define TOL 100.0*REAL_EPS

// define boundary index between external, internal extended axis indices
#define EXTINT 100

#include "utils.h"
#include "except.hh"

typedef struct s_axis {
  size_t n; 
  ireal d;
  ireal o;
  int id;
} axis;

void init_default_axis(axis * a);
  
void init_axis(axis * a, size_t n, ireal d, ireal o);
  
void copy_axis(axis * tgt, const axis * src);

void print_axis(axis a);
  
void fprint_axis(FILE * fp, axis a);
void fprint_num_axis(FILE * fp, int i, axis a);
int compare_axis(axis a1, axis a2);
typedef struct {
  int dim;
  int gdim;
  axis axes[RARR_MAX_NDIM];
} grid;
  
void init_default_grid(grid * g);

int init_grid(grid * g, int dim, int gdim);

void copy_grid(grid * tgt, const grid * src);

void print_grid(grid a);
  
void fprint_grid(FILE * fp, grid a);
  
int compare_grid(const grid g1, const grid g2);

int compatible_grid(const grid g1, const grid g2);

int get_dimension_grid(grid g);

int get_datasize_grid(grid g);

int get_extended_datasize_grid(grid g);

size_t get_global_datasize_grid(grid g);

ireal get_cellvol_grid(grid g);

ireal get_extended_cellvol_grid(grid g);

ireal get_global_cellvol_grid(grid g);

int get_panelnum_grid(grid g);

void get_n(IPNT n, grid g);
void get_d(_RPNT d, grid g);
void get_o(_RPNT o, grid g);

void get_gs(IPNT gs, grid g);

void get_ge(IPNT ge, grid g);

void get_id(IPNT id, grid g);

void get_ord(IPNT a, grid g);

bool grid_union(grid * g, axis const * ax);

bool init_step(grid g, IPNT step, bool fwd);

bool next_step(grid g, IPNT step);

#endif /* __SEAM_GRID */

---