elliptic.h

//
//  elliptic.h
//  GLAMER
//
//  Created by bmetcalf on 8/5/14.
//
//
 
#ifndef __GLAMER__elliptic__
#define __GLAMER__elliptic__
 
#include <iostream>
#include "standard.h"
//#include "slsimlib.h"
#include "lens.h"
 
class Elliptic{
public:
  Elliptic(
           LensHalo *isohalo    
           ,PosType fratio      
           ,PosType inclination 
           )
  : inclination(inclination),isohalo(isohalo){
    a = isohalo->getRsize()*sqrt(fratio);
    b = isohalo->getRsize()/sqrt(fratio);
  };
  
  void alpha(PosType x[],PosType alpha[]);
private:
 
  PosType a;
  PosType b;
  PosType inclination;
  LensHalo *isohalo;
 
  struct DALPHAXDM{
    DALPHAXDM(PosType lambda,PosType a2,PosType b2,PosType X[],LensHalo* isohalo)
    :lambda(lambda),a2(a2),b2(b2),x(X),isohalo(isohalo){};
    
    PosType operator()(PosType m);
  private:
    PosType lambda;
    PosType a2;
    PosType b2;
    PosType *x;
    LensHalo* isohalo;
  };
  struct DALPHAYDM{
    DALPHAYDM(PosType lambda,PosType a2,PosType b2,PosType X[],LensHalo* isohalo)
    :lambda(lambda),a2(a2),b2(b2),x(X),isohalo(isohalo){};
    
    PosType operator()(PosType m);
  private:
    PosType lambda;
    PosType a2;
    PosType b2;
    PosType *x;
    LensHalo* isohalo;
  };
  
};
 
 
#endif /* defined(__GLAMER__elliptic__) */