ElPoly.h

#ifndef ELPOLY_H
#define ELPOLY_H
#include "../include/VarData.h"
#include "../include/Cubo.h"
//#include "../include/Draw.h"
#ifdef USE_GL
#include <GL/glut.h>
#endif//USE_GL
#include "../include/Draw.h"
#include "../include/SingProc.h"
enum ElVisualize{
  EL_PART  =  0x0000,
  EL_QUAD    ,
  EL_CHAIN   ,
  EL_VECTORS   ,
  EL_DENS    ,
  EL_QUAD1   ,
  EL_POTENTIAL   ,
  EL_SURF    ,
  EL_CROSS   ,
  EL_ISOIPSE ,
  EL_SAMPLE  ,
  EL_SPECTRUM,
  EL_COLOR,
  EL_VORONOI,
  EL_SKIN,
  EL_STALK,
  EL_TRIA,
  EL_MESH,
  EL_ISOLEVEL,
  EL_POLYGON,
  EL_SQUAREMESH,
  EL_OPENGL,
  EL_POVRAY
};
class ElPoly : public VarData{
 private:
  char **cFile;
 public:
  ElPoly(int argc,char **argv,int *FilePos);
  ~ElPoly();
  //-----------------ElPolyProfile.cpp------------------
  int DensProf(int NBin,int NSample,int Coord);
  void DensProfNormalSlab(int NBin,int NSample,int Coord);
  int Diff2Files(int NBin,int How);
  void RestPress(int NBin);
  void SlabProf(int NBin,int nNano,int Coord1);
  void PrintDens(FILE *FileToWrite,double **Plot,double *LatDim,int NBin);
  void RadDens2Thick(int NBin);
  void RadDens2Thick2d(int NBin);
  int Temperature(int NBin,int Coord);
  void CartDens(int NBin,int nNano);
  void RadDistrF(int NBin,int How,int nNano);
  void BondDistr(int NSample);
  void E2EDistr(int NSample);
  void SplayDistr(int NSample);
  int RadialShell(int NBin);
  int NanoParticle(int NBin);
  int WormF(int Partition,int NBin);
  void StalkF(int NSample);
  void ThickFromDens(int NBin);
  void AreaDistrF(int NBin);
  void RadNormPos(int NBin,int NGrid);
  int SurfTens(int NBin);
  void SumTens();
  int PressTrace();
  int PressRadial();
  int Tens2dCartRad();
  void StalkLineProfF(int NBin);
  void Prova();
  //-----------------ElPolyMeasure.cpp------------------
  void PairCorr(int NBin,int NDim);
  int PairCorrelationF(int NBin,int How);
  int Diffusivity();
  void DiffSlab(int NSlab);
  int ScatteringF(int NBin,int How);
  int NChainPSquareF();
  void AreaCompr(int NSample);
  int SpectrumF(int NBin);
  void Midplane(int NBin);
  void SpectrumMidplane(int NBin);
  void HeaderAverage(int nNano);
  void WidomOut(char *InFile,int NBin);
  void WidomOut();
  void WidomIn(char *InFile,int NBin);
  void WidomIn();
  void End2EndDistr(char *OutFile);
  void BondDistr(char *OutFile,int NBin);
  void Decoupling(int What);
  void ElasticCoupling(int NSample);
  void ElasticCouplingNVT();
  void BilayerDistance(char *OutFile,int NBin);
  void EndToEndDist();
   //-----------------ElPolyRepr.cpp---------------------
  int ProjectionF(int NBin,int Coord);
  int CoreF(int NBin,int How);
  int Surface(int NBin,int Coord);
  int From3To2d(int Coord,double Param);
  int From2To1d(int Coord);
  int From3To1d(int Coord);
  void IsoSurf(int NSample,double *IsoValue,int NIso);
  void IsoLine(int NSample,double *IsoValue,int NIso,int How);
  void IsoLine(FILE *F2Write,double *Plot,int NSample,double *IsoLevel,int NIso);
  void FetchStalk();
  void FetchPore();
  void StalkArea();
  void AvSnap();
  void SlabAngleProfs(int NBin,int NAngle,int Coord);
  void Sample(int NSample);
  //-----------------ElPolyOutput.cpp-------------------
  void Conv2Tecplot(int NBin,int How);
  void Conv2Vmd();  
  void Conv2Povray();  
  void Conv2rzd(int NSample);
  void Conv2xyzd(int NSample);
  void DrField(int NGrid,double IsoLevel,int nNano,FILE *FWrite);
  void DrBondPovRay(double *Pos1,double *Pos2,float *Color);
  void HeaderPovRay();
  void DrNanoPovRay(int n);
  void DrPartPovRay(int p);
  void ConvLattice(int NSample,char *FName);
  //-----------------ElPolyEl.cpp-----------------------
  void RemoveChains();
  int Angle(int NBin);
  double ContactAngle(double x);
  int CenterOfMass(int Coord);
  int ChangeFile();
  int SpecifyCoord();
  int OpenFile(int f);
  int OpenFile(char *FileName);
  int PropertiesF();
  char *ChooseDraw(int ExtWhat2Draw);
  void ChooseDraw(char *String);
  void DivideLayers(int How);
  void SetBoundFile(int InitFile,int EndFile);
  void Processing(int f);
  void Shift2Center();
  void SetBackFold(int Bf){
    NBackFold = Bf;
  }
  void SetNVisSkip(int NSkip){
    if(NSkip > 0){
      NVisSkip = NSkip;
    }
  }
#ifdef __glut_h__
  //---------------ElPolyDraw.cpp----------------------
  void RenderPart(void);
  void DrRunTime();
  void DrColor();
  void DrPartList();
  void DrVector(Vettore v,Vettore Origin);
  void DrVectors();
  int DrIntorno(int p,double Blue);
  void DrCrossLinks();
  int Graphics(int argc,char **argv);
  void ESlide();
  void ESlide1();
  void DrProtein(const char *FileName,int nBlock);
  void DrNano();
  void DrCrossLinks(char *FileName);
  typedef void(ElPoly::*DRAW_PART)(int p);
  DRAW_PART Draw_Part;
  void DrawPart(int p){(*this.*Draw_Part)(p);}
  void DrPartOpenGl(int p);
  DRAW_PART Draw_Nano;
  void DrawNano(int n){(*this.*Draw_Nano)(n);}
  void DrNanoOpenGl(int n);
  void DrawFuncHeader();
  void DrawFuncFooter();
  typedef void(ElPoly::*DRAW_BOND)(double *Pos1,double *Pos2,float *Color);
  DRAW_BOND Draw_Bond;
  void DrBond(int p);
  void DrawBond(double *Pos1,double *Pos2,float *Color){(*this.*Draw_Bond)(Pos1,Pos2,Color);}
  void DrBondNo(double *Pos1,double *Pos2,float *Color){};
  void DrBondOpenGl(double *Pos1,double *Pos2,float *Color);
  //---------------ElPolyDrawControl.cpp----------------------
  void keyboard(unsigned char key,int x, int y);
  /* /// Definition of the key bindings */
  /* void keyboardDraw(unsigned char key); */
  void ElDrawMouse(int button,int state,int x,int y);
  void CompileList();
  void DrPosCol();
  void Menu();
  void ElMenuChoise(int option);
  void ElMenuVisual(int option);
  void CreateMenu();
  //---------------ElPolyDrawSurf.cpp----------------------
  void DrIsoipse(int NBin,int NIsoipse,int CoordN);
  void DrSurface();
  void DrSmooth(double *Plot,int NSample,double Min,double Max);
  void DrSample(int NSample);
  void DrSpectrum();
  void DrDerivative();
  void DrChains();
  void DrPolygon();
  void DrQuad();
  void DrDensity();
  void DrQuad1();
  void DrPotential();
  void DrShell();
  void DrVoronoi();
  void DrInterpSurface();
  void DrIsolevel(int NSample,double IsoLevel);
  void DrField(int NGrid,double IsoLevel,int nNano);
  void DrSquareMesh();
  void DrNormalPoint(int p,int NEdge);
  //-----------------------------------------------
  GLuint Quad,Point,*Cylinder,MetalCylinder,Hexagon,Cube,Arrow,GlWall;
  double Saturation;
  char Block2Draw[20];
  void DrTria(Vettore *v00,Vettore *v01,Vettore *v11,Vettore *vN);
  void DrTriaContour(Vettore *v00,Vettore *v01,Vettore *v11);
  void DrDoubleTria(Vettore *v00,Vettore *v01,Vettore *v11,Vettore *v10,Vettore *vN);
  void DrCreateStalk();
  void DrStalk();
  void Tile();
  //----------------------------ElPolyDrawCGAL----------------
  void DrTriangulate();
  void DrMesh();
  void DrGenMesh();
  void DrCells();
  void DefineSkin(int NSample);
  void DefineSurf();
#endif //__glut_h__
  int NFile[2];
  int NFileTot;
  int NPro ;
  int quando;
  //  int NBin;
  int LineSize;
  int IfIntorno;
  //int IfColour;
  //int IfChains;
  int IfLine;
  int IfChType;
  int NBackFold;
  int  What2Draw;
  int NVisSkip;
  //int *ChainTypes;
  int DrawOutput;
  //double Diameter;
  //double StepDiameter;
  //double ExtraDiam;
  double Vicinanze;
  double InvScaleUn;
  double ExtParam;
  double ScaleFact;
  FILE *DrawOutFile;
#ifdef USE_MPI
  SingProc *Proc;
#endif
};
static int IfImage;

extern void Slide();
#ifdef __glut_h__
#endif//__glut_h__
#endif //ELPOLY_H