miss.cu

#include <optix_world.h>
#include "datadef.h"

rtDeclareVariable(optix::Ray, ray, rtCurrentRay, );
rtDeclareVariable(intersection_point, payload, rtPayload, );
rtDeclareVariable(uint, launch_index_in, rtLaunchIndex, );
rtDeclareVariable(unsigned,  trace_type, , );
rtDeclareVariable(unsigned,  outer_cell, , );
rtBuffer<unsigned,1>          rxn_buffer;
rtBuffer<spatial_data,1>      positions_buffer;
rtBuffer<unsigned,1>      	  matnum_buffer;
rtBuffer<unsigned,1>      	  cellnum_buffer;
rtBuffer<unsigned,1>      	  remap_buffer;

RT_PROGRAM void miss()
{
	unsigned launch_index;
	
	// remap if transport
	if(trace_type==2){
		launch_index=remap_buffer[launch_index_in];
		//rtPrintf("!!!MISS IN TRANSPORT!!! rxn = %u outer_cell = %d launch_index = %d remapped index %u trace %d (x,y,z,xhat,yhat,zhat-source)= (% 10.9E % 10.9E % 10.9E)(% 10.9E % 10.9E % 10.9E)\n", rxn_buffer[launch_index_in], outer_cell, launch_index_in, launch_index, trace_type,positions_buffer[launch_index].x,positions_buffer[launch_index].y,positions_buffer[launch_index].z,positions_buffer[launch_index].xhat,positions_buffer[launch_index].yhat,positions_buffer[launch_index].zhat);
		rtPrintf("!!!MISS(x,y,z,xhat,yhat,zhat-source) = % 10.9E % 10.9E % 10.9E % 10.9E % 10.9E % 10.9E\n", positions_buffer[launch_index].x,positions_buffer[launch_index].y,positions_buffer[launch_index].z,positions_buffer[launch_index].xhat,positions_buffer[launch_index].yhat,positions_buffer[launch_index].zhat);
	}
	else{
		launch_index = launch_index_in;  // misses in fissile query and geometry plotting are expected for out-of-bounds regions in non-rectangular geometries
	}
	
	payload.sense 				= -9;
	rxn_buffer[launch_index_in]	=  997;     //miss code, same as leak basically
	payload.surf_dist 			= -1.0;
	payload.cell 				=  3000;
	payload.mat  				=  3000;
	payload.fiss 				=  0;


}