forceTree.h

This file has been depricated and still exists only to recycle code.
 
/*
 * forceTree.h
 *
 *  Created on: Oct 14, 2011
 *      Author: bmetcalf
 */
 
#ifndef FORCE_TREE_H_
#define FORCE_TREE_H_
 
#include "simpleTree.h"
#include "lens_halos.h"
 
//enum PartProf {gaussian,powerlaw};
 
class TreeForce : public TreeSimple {
public:
 TreeForce(PosType **xp,IndexType Npoints,float *masses,float *rsph,bool Multimass,bool Multisize
   ,PosType my_kappa_background = 0,int bucket = 5,int dimensions = 2,bool median = false,PosType theta = 0.1
   );
 
 TreeForce(PosType **xp,IndexType Npoints,LensHalo *my_halos
   ,bool Multisize = true,PosType my_kappa_bk=0.0,int bucket = 5,int dimensions = 2
   ,bool median = false,PosType theta = 0.1
   );
 
 ~TreeForce();
 
 float * CalculateSPHsmoothing(int N);
 void force2D(PosType const *ray,PosType *alpha,KappaType *kappa,KappaType *gamma
                 ,KappaType *phi);
 // void force2D(PosType *ray,PosType *alpha,KappaType *kappa,KappaType *gamma);
 //void ChangeParticleProfile(PartProf partprof);
 
 
protected:
 
 bool init;
 
 bool MultiMass;
 bool MultiRadius;
 float *masses;
 float *rsph;
 PosType kappa_background;
 
 PosType force_theta;
 
 //PosType (*alpha_particle)(PosType r,float xmax);
 //PosType (*kappa_particle)(PosType r,float xmax);
 //PosType (*gamma_particle)(PosType r,float xmax);
 
 //PosType (*alpha_halo)(PosType r,HaloInternal &par);
 //PosType (*kappa_halo)(PosType r,HaloInternal &par);
 //PosType (*gamma_halo)(PosType r,HaloInternal &par);
 
 void CalcMoments();
 void rotate_coordinates(PosType **coord);
 //void spread_particles();
 
 inline virtual PosType alpha_h(PosType r2s2,PosType sigma){
   return (sigma > 0.0 ) ? ( exp(-0.5*r2s2) - 1.0 ) : -1.0;
 }
 inline virtual PosType kappa_h(PosType r2s2,PosType sigma){
   return 0.5*r2s2*exp(-0.5*r2s2);
 }
 inline virtual PosType gamma_h(PosType r2s2,PosType sigma){
   return (sigma > 0.0 ) ? (-2.0 + (2.0 + r2s2)*exp(-0.5*r2s2) ) : -2.0;
 }
 inline virtual PosType phi_o(PosType r2,PosType sigma){
  ERROR_MESSAGE();  // not yet written
  std::cout << "time delay has not been fixed fot this profile yet." << std::endl;
  exit(1);
  return 0;
 }
 
 bool haloON;
 LensHalo *halos;
};
 
typedef TreeForce *TreeForceHndl;
 
#endif /* FORCE_TREE_H_ */