Homec4science
Diffusion Lenstool-HPC (lenstool-mixed) Browse

R1448/Benchmarks/GridGradientBenchmarkMixed/structure_mixed.hpp
3f0b763e9149lenstool-mixed

Last Change
Raw File

structure_mixed.hpp
View Options

#pragma once
template<typename fp_type = double>
struct Potential_Mixed
{
int* type; // 1=PIEMD ; 2=NFW; 3=SIES, 4=point mass
char type_name[10]; // PIEMD, NFW, SIES, point
fp_type* name; // name of the clump (e.g. name of the galaxy) : not compulsory
//struct point position; // position of the center of the halo
fp_type* position_x; // position of the center of the halo
fp_type* position_y; // position of the center of the halo
fp_type* weight; // weight of the clump (the projected mass sigma0 for PIEMD, the density rhoscale for NFW)
fp_type* b0; // Impact parameter
fp_type* vdisp; //Dispersion velocity
fp_type* ellipticity_angle; // orientation of the clump
fp_type* ellipticity; // ellipticity of the mass distribition
fp_type* ellipticity_potential; //ellipticity of the potential
fp_type* rcore; // core radius
fp_type* rcut; // cut radius
fp_type* rscale; // scale radius for NFW, Einasto
fp_type* exponent; // exponent for Einasto
fp_type* alpha; // exponent for general NFW
fp_type* einasto_kappacritic; // critical kappa for Einasto profile
fp_type* z; // redshift of the clump
fp_type* mag; //magnitude
fp_type* lum; //luminosity
fp_type* theta; //theta
fp_type* anglecos; //theta precomputation of cosinus and sinus values
fp_type* anglesin; //theta
fp_type* sigma; // sigma
};
//
//
//
template<typename fp_type = double>
void
alloc(struct Potential_Mixed<fp_type> &lens_mixed, const int Nlens)
{
lens_mixed.type = (int*) malloc(Nlens*sizeof(int));
lens_mixed.name = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.position_x = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.position_y = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.weight = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.b0 = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.vdisp = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.ellipticity_angle = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.ellipticity = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.ellipticity_potential = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.rcore = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.rcut = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.rscale = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.exponent = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.alpha = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.einasto_kappacritic = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.z = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.mag = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.lum = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.theta = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.anglecos = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.anglesin = (fp_type*) malloc(Nlens*sizeof(fp_type));
lens_mixed.sigma = (fp_type*) malloc(Nlens*sizeof(fp_type));
};
//
//
//
template<typename fp_type = double>
void
copy(struct Potential_Mixed<fp_type> &lens_mixed, const struct Potential_SOA *lens, int Nlenses)
{
#if 1
for (int i = 0; i < Nlenses; ++i)
{
lens_mixed.type[i] = (int) lens->type[i];
lens_mixed.name[i] = (fp_type) lens->name[i];
lens_mixed.position_x[i] = (fp_type) lens->position_x[i];
lens_mixed.position_y[i] = (fp_type) lens->position_y[i];
lens_mixed.b0[i] = (fp_type) lens->b0[i];
lens_mixed.ellipticity[i] = (fp_type) lens->ellipticity[i];
lens_mixed.ellipticity_angle[i] = (fp_type) lens->ellipticity_angle[i];
lens_mixed.ellipticity_potential[i] = (fp_type) lens->ellipticity_potential[i];
lens_mixed.rcore[i] = (fp_type) lens->rcore[i];
lens_mixed.rcut[i] = (fp_type) lens->rcut[i];
lens_mixed.z[i] = (fp_type) lens->z[i];
lens_mixed.anglecos[i] = (fp_type) lens->anglecos[i];
lens_mixed.anglesin[i] = (fp_type) lens->anglesin[i];
}
#else
memcpy(&(lens_mixed.type)[0] , &(lens->type)[0] , Nlenses*sizeof(int));
memcpy(&(lens_mixed.position_x)[0] , &(lens->position_x)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.position_y)[0] , &(lens->position_y)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.b0)[0] , &(lens->b0)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.ellipticity_angle)[0] , &(lens->ellipticity_angle)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.ellipticity)[0] , &(lens->ellipticity)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.ellipticity_potential)[0] , &(lens->ellipticity_potential)[0], Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.rcore)[0] , &(lens->rcore)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.rcut)[0] , &(lens->rcut)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.z)[0] , &(lens->z)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.anglecos)[0] , &(lens->anglecos)[0] , Nlenses*sizeof(fp_type));
memcpy(&(lens_mixed.anglesin)[0] , &(lens->anglesin)[0] , Nlenses*sizeof(fp_type));
#endif
}
c4science · Help