VarDatFile.h

#ifndef VAR_DAT_FILE_H
#define VAR_DAT_FILE_H
#include "../include/Matematica.h"
#define DIS_IF_TYPE(c,t) ( ((c)&(t))==(t) )
#define DIS_ADD_TYPE(c,t) ( (c)|=(t) )
#define DIS_REM_TYPE(c,t) ( (c)&=(~t) )
enum VisLog{
  DIS_NOLOG = 0x001,
  DIS_LOGX  = 0x002,
  DIS_LOGY  = 0x004,
  DIS_LOGLOG= 0x006,
};

enum RefType{
  REF_SEQ  = -1,
  REF_ASC  = -2,
};
class VarDatFile:public Matematica{
 private:
  char XFormula[60];
  char YFormula[60];
  char Header[STRSIZE];
  void Init();
  void Allocate();
  int NMax;
  int NMaxPunti;
  int NVar;
  double **st;
  double *sMin;
  double *sMax;
  char **sc;
  int *SetNMax;
  int *sType;
  int *sColor;
  int *RefAbsc;
  double *Punti;
  double *PuntiErr;
  double *Punti1;
  double *dInter;
  double *dError;
  double *dInter1;
  double GlobMax;
  double GlobMin;
  double xGlobMax;
  double xGlobMin;
  int IfPuntiAlloc;
  int IfPuntiFree;
  //int Incr;
 public:
  VarDatFile(char *file,int NBin);
  VarDatFile(int ExtNMax,int ExtNVar,int ExtNBin);
  VarDatFile(double **ExtSt,int ExtNMax,int ExtNVar,int ExtNBin);
  VarDatFile(char **FileList,int *Pos,int NFile,int Values);
  ~VarDatFile();
  int Aggiungi(char *file);
  int FileInfo(FILE *InFIle,int *NCol);
  int ReadLine(char *cLine,double *Value);
  int  ReadLine(FILE *InFile, const char *Format, ...);
  int ElabSegnale(int NElMin,int NElMax);
  MOMENTI DistrSegnale(int NElMin,int NElMax,int IfNorm);
  MOMENTI DistrLogSegnale(int NElMin,int NElMax,int IfNorm);
  MOMENTI DistrExpSegnale(int NElMin,int NElMax,int IfNorm);
  MOMENTI DistrSegnale(int NElMin,int NElMax,double *Border,int IfNorm);
  MOMENTI DistrSignErr(int NElMin,int NElMax,double *Border,int IfNorm);
  MOMENTI DistrGaussSegnale(int ElMin,int ElMax,int IfNorm);
  MOMENTI WeightAverageSet(int CoordY,int ElMin, int ElMax);
  void DistrSignSample(int ElMin, int ElMax,double **Distr,int NSample,int IfNorm,double *xBound);
  int SpettroSegnale(int NElMin,int NElMax);
  void SpeLine(int ElMin,int ElMax,int NBin,double *Spe);
  int NormalizzaSegnale(int NElMin,int NElMax);
  int NormalizzaInter();
  bool RadiceSegnale();
  int AutocorSegnale(int NElMin,int NElMax);
  double SumSegnale(int CoordY,int NElMin,int NElMax);
  double IntSegnale();
  void SommaSegnali();
  void AverageOrdinate(int ElMin, int ElMax,double *Distr,double *xBound);
  void Punta(int n);
  void Punta(double **ExtSt,int n);
  void Punta(double *sp,int n);
  void PuntaInt(double *sp);
  void DerivataSegnale(){Elab = &Matematica::DerO4;};
  void VarieSegnale(){Elab = &Matematica::SquareGradient;};
  void ModuloSegnale(){Elab = &Matematica::Modulo;};
  double VarieSegnale(int NElMin,int NElMax);
  RETTA InterRettSegnale(int CoordY,int NElMin,int NElMax,int LogLog);
  RETTA InterExpSegnale(int CoordY,int NElMin,int NElMax,int LogLog);
  MOMENTI InterGaussSegnale(int CoordY,int NElMin,int NElMax,int LogLog);
  PARABOLA ParabolaSegnale(int CoordY,int NElMin,int NElMax,int LogLog);
  int MediaMobSegnale(int);
  int WeightHistoSign(int NHisto);
  int CorrelaADuePunti(int dist);
  void AutosimilaritaSegnale(int);
  void CambiaNBin(int);
  void CambiaPunti();
  void Print();
  void Sort();
  void ScriviPunti(char *file);
  void ScriviFile(char *file,int CoordY,int LogLog,int NVisMin,int NVisMax);
  void ScriviTutto(char *file,int LogLog,int NVisMin,int NVisMax);
  void RescaleToBulk(char *FName);
  void ExportTxvl(char *FName,int NElMin,int NElMax);
  void TecPlot(char *FName);
  void PuntiAlloc();
  void PuntiFree();
  int NBin;
  int pNVar(){ return NVar;};
  int pNMax(){return NMax;};
  int pNRow(int CoordY);
  void pGlobBorder(double *xMin,double *xMax,double *yMin,double *yMax){
    *xMin = xGlobMin;
    *xMax = xGlobMax;
    *yMin = GlobMin;
    *yMax = GlobMax;
  };
  double Val(int CoordY,int n);
  double Abscissa(int CoordY,int n);
  double pPunti(int n);
  double pPuntiErr(int n);
  void pMinMaxGlob(int NVisMin,int NVisMax);
  double pMaxGlob(int NVisMin,int NVisMax);//{return GlobMax;};
  double pMinGlob(int NVisMin,int NVisMax);//{return GlobMin;};
  double pxMinGlob(int NVisMin,int NVisMax);//{return GlobMin;};
  double pxMaxGlob(int NVisMin,int NVisMax);//{return GlobMax;};
  double pMax(int CoordY,int NVisMin,int NVisMax);//{return sMax[n];};
  double pMin(int CoordY,int NVisMin,int NVisMax);//{return sMin[n];};
  double pMaxGlobLog(int NVisMin,int NVisMax);//{return GlobMax;};
  double pMinGlobLog(int NVisMin,int NVisMax);//{return GlobMin;};
  double pMaxLog(int CoordY,int NVisMin,int NVisMax);//{return sMax[n];};
  double pMinLog(int CoordY,int NVisMin,int NVisMax);//{return sMin[n];};
  double pInter(int n){return dInter[n];};
  double pError(int n){return dError[n];};
  double pInter1(int n){return dInter1[n];};
  double PuntiMin();
  double PuntiMax();
  int IsAbscissa(int Col){
    if(RefAbsc[Col] == Col) return 1; return 0;};
  int pRefAbsc(int Col){return RefAbsc[Col];};
  int pSetNMax(int Col){return SetNMax[Col];};
  int IsSequence(int Col){return RefAbsc[Col]==REF_SEQ?1:0;};
  void ImpSequence(int Col);
  void ImpSequence();
  void ImpCoordX(int vAbs);
  void ImpCoordX(int vSet,int vAbs);
  void ImpCoordY(int Ext);
  void setXFormula(char *str);
  void setYFormula(char *str);
  void Reverse();
  int Smooth(double Fact,int CoordY,int NVisMin,int NVisMax);
  void SmoothGauss(double Fact,int CoordY,int NVisMin,int NVisMax);
  void DoubleDistFluct();
  void WriteFormula(char *Exit);
  char *PrintHeader();
  double *sp;
};
#endif //VARELEMENTIGRAFICI_H