MatematicaVect.h

#ifndef MATEMATICAVECT_H
#define MATEMATICAVECT_H
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

//-----------------Vettore-----------------
class Vettore{
 public:
  Vettore(int N);
  Vettore();
  ~Vettore();
  Vettore(double *Pos,int N);
  Vettore(double x,double y,double z);
  Vettore(double x,double y);
  int NDim;
  double *x;
  //  Vettore *CrossProduct(Vettore *u,Vettore *v);
  double Abs();
  double Norm();
  double Normalize();
  //Vettore Normal(const Vettore u,const Vettore v);
  void Normal(const Vettore *u,const Vettore *v);
  void NormalSurf(const Vettore *u,const Vettore *v,const Vettore *w);
  double ProjOnSurf(Vettore *S1,Vettore *S2,Vettore *S3,Vettore *P);
  void Mult(double Fact);
  void Subs(const Vettore *u,const Vettore *v);
  double ScalS(const Vettore *u,const Vettore *v);
  double CosAngle(Vettore *u);
  double CosAngle(Vettore *u,Vettore *v);
  double SinAngle(Vettore *u,Vettore *v);
  double SinAngle(Vettore *u);
  double Angle(Vettore *u,Vettore *v);
  double Angle(Vettore *u);
  double Col(int N);
  double Val(int N){return Col(N);};
  void Set(double Val,int Col);
  void Axis(Vettore *u,Vettore *v);
  void ScalV(const Vettore *u,const Vettore *v);
  double VetV(const Vettore *u,const Vettore *v);
  double VetV3(const Vettore *u,const Vettore *v);
  double VetV(const Vettore *u);
  double ProjOnAxis(Vettore *a);
  double ProjOnAxis(Vettore *Pos,Vettore *Axis);
  void ApplyOn(Vettore *o);
  void Copy(Vettore *o);
  void Export(double *x);
  void PerpTo(Vettore *o);
  double PerpTo(Vettore *Pos,Vettore *Axis);
  double PerpTo3(Vettore *Pos,Vettore *Axis);
  void Print();
  void Rescale(double NewLength);
  Vettore operator+(const Vettore &Vet);
  Vettore operator+(const Vettore &Vet) const;
  Vettore& operator+=(const Vettore &Vet);
  Vettore operator-(const Vettore &Vet);
  Vettore operator-(const Vettore &Vet) const;
  Vettore& operator-=(const Vettore &Vet);
  Vettore operator*(const Vettore &Vet); 
  Vettore& operator*=( const Vettore& vec );
  Vettore operator*(const double Fact);
  Vettore& operator*=( const double f ); 
  double operator%(const Vettore &Vet); 
  double& operator%=(const Vettore &Vet); 
  Vettore operator^(const Vettore& vec);
  Vettore& operator^=(const Vettore& vec);
  Vettore operator/(const Vettore& vec);
  Vettore& operator/=(const Vettore& vec);
  Vettore operator=(const Vettore &Vet);
  double * getPtr() {
    return (double*)&x;
  }
  const double * getPtr() const {
    return (const double *)&x;
  }
  double& operator[]( int col ){
    //return getPtr()[col];
    return x[col];
  }
  double operator[]( int col ) const {
    //return getPtr()[col];
    return x[col];
  }
};
  Vettore operator*(double Fact, const Vettore& vec);
#endif //MATEMATICA_H