setup.hpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include <iostream>
	#include <string.h>
	//#include <cuda_runtime.h>
	#include <math.h>
	#include <sys/time.h>
	#include <fstream>
	#ifndef __xlC__
	#include <immintrin.h>
	#endif

	//#include "simd_math.h"

	#ifdef __WITH_LENSTOOL
	#include "structure.h"
	#include "structure_hpc.hpp"

	extern struct pot lens[NLMAX];
	extern struct g_grille G;
	void
	//setup_jauzac_LT(struct pot** lens, int* nlenses, double* x, double* y, double* sol_grad_x, double* sol_grad_y)
	convert_to_LT(Potential_SOA* lenses, int nlenses)
	{

	G.xmax = 50.0;
	G.xmin = -50.0;
	G.ymax = 50.0;
	G.ymin = -50.0;
	G.dx = 100./1000.;
	G.dy = 100./1000.;
	G.nlens = 6;
	//
	//lens = (struct pot) malloc(sizeof(struct pot)(nlenses));
	//
	for (int i = 0; i < nlenses; ++i)
	{
	//
	lens[i].C.x = lenses->position_x[i];
	lens[i].C.y = lenses->position_y[i];
	//
	lens[i].sigma = 1.; // disp
	if (lenses->type[i] == 5) {
	lens[i].type = 1;
	}
	else{
	lens[i].type = lenses->type[i];
	}
	lens[i].emass = 0.11;
	lens[i].epot = lenses->ellipticity_potential[i];
	//printf("%f %f\n", lenses->ellipticity_potential[i], lenses->ellipticity[i]);
	//lens[i].epot = lenses->ellipticity[i];
	lens[i].theta = lenses->ellipticity_angle[i];
	lens[i].rcut = lenses->rcut[i];
	lens[i].rc = lenses->rcore[i];
	lens[i].b0 = lenses->b0[i];
	lens[i].z = lenses->z[i];
	lens[i].masse = 0; // weight
	//lens[i].rc = 0; // rscale
	// (&lens)[i].exponent = 0;
	lens[i].alpha = 0.;
	// (&lens)[i].einasto_kappacritic = 0;
	lens[i].z = 0.4;
	}
	}
	#endif