summaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
2011-08-15woops. dot product should be clamped to [-1.0,1.0].Anthony LaTorre
2011-08-15fix nan bug by clamping dot() of surface normal and -photon direction to ↵Anthony LaTorre
[0.0,1.0]
2011-08-14Rewrite the color_solid function in gpu.GPU (and associated CUDA code)Stan Seibert
to make it 100x faster. Instead of having each CUDA thread loop over the full triangle list, we give each thread a single triangle and ask it to look up the hit status and color for that triangle. The hit array and color array are small enough (approx 30,000 entries) to fit into the cache, so this goes much faster. Now the event viewer is quite snappy!
2011-08-14Add a helper function to the ROOT Channel class that will populate 3Stan Seibert
arrays with hit information: hit or not, time and charge. Handy for fetching all the hits in an event into numpy arrays.
2011-08-14Oops, view.py is still in use.Stan Seibert
2011-08-13Remove unused render and view files. They have both been superceded by ↵Stan Seibert
camera.py.
2011-08-13A faulty optical process can make the position or direction of theStan Seibert
photon into NaN on the GPU. Now we abort these photons rather than let them lock up the intersect_mesh() method. There is a new history bit (#31) that indicates when a NAN_ABORT has occurred, and this bit is checked for by GPU.propagate(). If set for any of the photons, a warning message is printed. While not as good as preventing the NaN problem in the first place, this at least ensures we are aware of the problem.
2011-08-13For consistency, we should keep the water index of refraction down toStan Seibert
200 nm, like all the other optical properties. It no longer is a speed benefit, anyway.
2011-08-13Finally get GEANT4 to stop printing debug information about all the physics ↵Stan Seibert
processes
2011-08-12mergeAnthony LaTorre
2011-08-12do not check child nodes of a node at which the distance to the bounding box ↵Anthony LaTorre
is further than a triangle the ray/photon has already intersected
2011-08-12MergeStan Seibert
2011-08-12Use an input and output photon queue in order to consolidate all theStan Seibert
photons that didn't die during propagation into the beginning of the list. This speeds up propagation by reducing the number of partially filled CUDA warps on the next propagation step. 2.2 million photons/sec on LBNE!
2011-08-12mergeAnthony LaTorre
2011-08-12speedup Mesh.remove_duplicate_vertices()Anthony LaTorre
2011-08-12Widen the field of view in camera to make it easier to look at Cherenkov ↵Stan Seibert
rings from inside the detector.
2011-08-12Allow multiple GEANT4 generator threads to keep up with the speed of photon ↵Stan Seibert
propagation. Defaults to 4.
2011-08-12Segmented cylinder FOR MOAR SPEED! 1.8 Mphotons/sec!Stan Seibert
2011-08-12G4Generator and sim.py can now generate boosted pi0 decays.Stan Seibert
The ROOT data structure has been extended to allow storage of the individual gamma rays as subtracks.
2011-08-12Refactor ROOT file writing into fileio.root.RootWriter classStan Seibert
2011-08-12Storage for subtracks in ROOT data structure.Stan Seibert
2011-08-11Make the generator process a "daemon" process, so it will be cleanedStan Seibert
up when the parent exits, and also start the generator process earlier.
2011-08-11mergeStan Seibert
2011-08-11Tell CUDA to enlarge L1 cache.Stan Seibert
2011-08-11merge headsAnthony LaTorre
2011-08-11add consume() to itertoolset. checksum the number of bits shifted in each ↵Anthony LaTorre
iteration of the bounding volume hierarchy construction. only cache the bounding volume hierarchy information.
2011-08-11No need for __noinline__ now that kernel caching works.Stan Seibert
2011-08-11Show number of registers used in CUDA kernelsStan Seibert
2011-08-11Switch from texture to float3 array for upper and lower bounds. 10% speed ↵Stan Seibert
boost!
2011-08-10Set the GEANT4 and CUDA RNG seeds using current time and process ID if not ↵Stan Seibert
set on command line.
2011-08-10Use WCSim properties for water with the lbne detectorStan Seibert
2011-08-10Fix misuse of Material.set(), so now wavelengths and values of materialStan Seibert
properties are not transposed.
2011-08-10Using gzip compression level 1, the BVH at 10 bit is 7x smaller. AddsStan Seibert
10 seconds to cache load time if the BVH file is in the memory cache, otherwise this is faster if the BVH file must be read fresh from disk.
2011-08-10Rename chroma.io to chroma.fileio to avoid collision with Python package ↵Stan Seibert
named io
2011-08-10merge headsAnthony LaTorre
2011-08-10add the ability to cache a geometry along with its bounding volume ↵Anthony LaTorre
hierarchy. cached files are stored in $HOME/.chroma. fixed the timeit() decorator so that the decorated function is still able to pass back a return value.
2011-08-09Automatically detect kernprof and profile main() in sim.py. Also ↵Stan Seibert
synchronize CUDA calls.
2011-08-09The --use-fast-math gives the simulation another 25% speed boost.Stan Seibert
2011-08-09Autodetect presence of Space Navigator rather than require an option.Stan Seibert
2011-08-09Raise exception if unable to spnav_x11_open()Stan Seibert
2011-08-09Bug fix that caused sim.py to only record hits on the first 3 PMTsStan Seibert
2011-08-09fix variable name in GPU.color_solids()Anthony LaTorre
2011-08-09merge headsAnthony LaTorre
2011-08-09switch to indexing child nodes by start and stop indices instead of start ↵Anthony LaTorre
and length; this reduces a bit of arithmetic when traversing the bounding volume hierarchy and makes the code in Geometry.build() more concise. add an untested cuda kernel to interleave the bits in three uint64 arrays.
2011-08-09Put number of detected photons into charge value for channel.Stan Seibert
2011-08-09Store a photon history for each hit channel. If multiple photons hit theStan Seibert
same channel, their history bits are OR'ed together.
2011-08-09Update camera to use SpaceNav over X11 protocol. Now events pass through ↵Stan Seibert
pygame loop.
2011-08-09Update spnav ctypes wrapper to support X11-based communication.Stan Seibert
2011-08-08propagate() takes an array of photon offsets and a range ofStan Seibert
offsets to load. Now events with more photons than RNG states can be propagated through multiple kernel calls. Also lays the groundwork for consolidating photons between steps to reduce the amount of propagation work required.
2011-08-08Record generated particle information in data structureStan Seibert