#include <iostream>
using namespace std;

#include <stdio.h>

#include <JANA/JCalibrationFile.h>
using namespace jana;

#include <HDGEOMETRY/DMagneticFieldMapParameterized.h>

// Declarations of wrapper routines (prefixed with gpu_)
void gpu_VectAdd(float *A, float *B, float *C);

int main(int narg, char *argv[])
{

	// Create JCalibration object
	const char *url = getenv("JANA_CALIB_URL");
	if(!url){
		cerr<<"JANA_CALIB_URL not set!!!"<<endl;
		exit(-1);
	}
	JCalibration *jcalib = new JCalibrationFile(string(url), 1);
	
	// Create DMagenticField object
	DMagneticFieldMapParameterized *bfield = new DMagneticFieldMapParameterized(jcalib);

	// Add two 3-D vectors on GPU
	float A[3], B[3], C[3];
	A[0] = 1;	B[0] = 27;
	A[1] = 7;	B[1] = 31;
	A[2] = 13;	B[2] = 47;
	
	gpu_VectAdd(A, B, C);
	
	for(int i=0; i<3; i++){
		cout<<"A["<<i<<"] + B["<<i<<"] = "<<A[i]<<" + "<<B[i]<<" = "<<C[i]<<endl;
	}

	return 0;
}