cpolytope.h

Go to the documentation of this file.

#ifndef CPOLYTOPE
#define CPOLYTOPE


#include "boolean.h"
#include "box.h"
#include "parameters.h"
#include "scpolytope.h" /* To get some common constants */
#include "wcs.h"

#define POLYTOPE_R_ENLARGEMENT 1.01

typedef struct {
  unsigned int k;   
  double r;         
  double v;         
  boolean emptyPolytope;     
  unsigned int nFaces;       
  unsigned int maxFaces;     
  double **face;             
  unsigned int *ID;          
  unsigned int maxVertices;  
  unsigned int nVertices;    
  double **vertex;           
  boolean *expandible;       
  unsigned int nExpandible;  
  unsigned int freeVertex;   
  unsigned int *nIndices;    
  unsigned int *maxIndices;  
  unsigned int **indices;    
  Tbox bb;                   
  boolean open;              
} Tcpolytope;


void InitEmptyPolytope(unsigned int k,double r,Tcpolytope *mp);

void DefinePolytope(Tcpolytope *mp);

void SPolytope2Polytope(Tparameters *pr,Tscpolytope *sp,Tcpolytope *p);


void Polytope2SPolytope(Tparameters *pr,Tcpolytope *p,Tscpolytope *sp);


void CopyPolytope(Tcpolytope *mp_dst,Tcpolytope *mp_src);

boolean CutPolytope(Tparameters *pr,double *t,double r,unsigned int id,void *wcs,void *c,
                    unsigned int m,unsigned int *tp,Tbox *ambient,Tcpolytope *mp);

boolean CutPolytopeWithFace(Tparameters *pr,double *t,double offset,unsigned int id,void *wcs,
                            void *c,unsigned int m,unsigned int *tp,Tbox *ambient,Tcpolytope *mp);

void PolytopeCenter(double *t,Tcpolytope *mp);

boolean InsidePolytope(double *t,Tcpolytope *mp);


double GetPolytopeRadius(Tcpolytope *mp);

unsigned int GetPolytopeDim(Tcpolytope *mp);

Tbox *GetPolytopeBB(Tcpolytope *mp);

void SetPolytopeBB(Tbox *bb,Tcpolytope *mp);


boolean ExpandiblePolytope(Tcpolytope *mp);

void WrongPolytopeCorner(unsigned int nv,Tcpolytope *mp);

boolean GetPolytopeInteriorPoint(double rSample,double *t,Tcpolytope *mp);

boolean PolytopeBoundaryPointFromExternalCorner(double rSample,boolean rand,
                                                unsigned int *nv,double *t,Tcpolytope *mp);


void PolytopeBoundaryPointsFromExternalCorners(double rSample,unsigned int *n,unsigned int **nv,
                                               double ***t,Tcpolytope *mp);

boolean PolytopeRandomPointOnBoundary(double rSample,double *t,Tcpolytope *mp);

boolean RandomPointInPolytope(double *t,Tcpolytope *mp);

double PolytopeMaxVolume(Tcpolytope *mp);

double PolytopeVolume(Tcpolytope *mp);

unsigned int PolytopeNumNeighbours(Tcpolytope *mp);

unsigned int PolytopeNeighbourID(unsigned int n,Tcpolytope *mp);

void GetPolytopeVertices(unsigned int *nv,double ***v,Tcpolytope *mp);

void GetPolytopeNeighboursFromVertices(unsigned int *nv,unsigned int **cID1,unsigned int **cID2,Tcpolytope *mp);

void GetPolytopeEdges(unsigned int *ne,unsigned int **vID1,unsigned int **vID2,
                      Tcpolytope *mp);


unsigned int PolytopeMemSize(Tcpolytope *mp);

void SavePolytope(FILE *f,Tcpolytope *mp);

void LoadPolytope(FILE *f,Tcpolytope *mp);

void PrintPolytopeInfo(Tcpolytope *mp);

void DeletePolytope(Tcpolytope *mp);

#endif