summaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
2011-08-03Merge with TonyStan Seibert
2011-07-30reduce number of lookup calls.Anthony LaTorre
2011-07-30you can rotate just the camera by holding the control key and take movies by ↵Anthony LaTorre
pressing the m key.
2011-07-30when throwing photons from the light source out onto the scene, photons are ↵Anthony LaTorre
now thrown randomly across each triangle instead of only at the center of each triangle. all of the rendering kernels have been rewritten so that they operate additively; for example, you may now throw photons from the light source onto the scene, render from the camera to the scene, then throw more photons and render again.
2011-07-29Add SNO acrylic and LAB+PPO scintillator to detector model.Stan Seibert
2011-07-28Alternate "real_sno" detector configuration using the CAD model of theStan Seibert
acrylic vessel. There are serious issues with this model which probably make it unsuitable for photon propagation at the moment. (Example: Why do you need to traverse 8 surfaces to cross the acrylic vessel? There should only be 4.)
2011-07-28Oops, forgot to include the reduced 8" PMT profile.Stan Seibert
2011-07-28Increase alpha depth for complex models. Remove intermediate roundingStan Seibert
to integers in alpha blending until the last step.
2011-07-28Support for reading bz2 compressed STL files.Stan Seibert
2011-07-27Merge changesStan Seibert
2011-07-27Use the more realistic SNO database and add PMT concentrators.Stan Seibert
2011-07-27camera automatically increases render quality when it is not being movedAnthony LaTorre
2011-07-26Preliminary SNO detector.Stan Seibert
PMTs in the right place and acrylic vessel. No new material properties yet.
2011-07-26added a camera class which is able to render using the simple ray tracer or ↵Anthony LaTorre
the hybrid monte carlo ray tracer in a separate thread. the camera object is initialized by passing a pycuda.driver.Context object and a threading.Lock object; you can then run kernels and copy arrays to and from the device in the same context as that used by the camera by acquiring the lock, calling context.push(), executing the kernel and/or copying arrays, calling context.pop(), and then releasing the lock. fixed mistakes in both build_rgb_lookup() and render() where I had accidently switched the green and blue photons. updated the lbne geometry with the latest specifications. added profiles for the sno pmt and concentrator.
2011-07-25moved triangle colors to a separate global device array so that the ray ↵Anthony LaTorre
tracer and photon simulation can be run in the same context. added the ability to add alpha channel to triangle color so that triangles can be made transparent. added __noinline__ modifier to certain device functions to speed up kernel compilation.
2011-07-21pull random number generator states out of the photon struct; this allows ↵Anthony LaTorre
you to copy photon information within a thread but still keep a single random number generator throughout the thread.
2011-07-21hybrid monte carlo render now distinguishes between the two different sides ↵Anthony LaTorre
of each triangle. reduced the number of runs to average when propagating photons from each pixel in render.py from 5 to 1; the speed improvement outweighs any small improvement in the quality of the rendered image.
2011-07-20in the previous commit i wrongly edited the code to reflect/transmit across ↵Anthony LaTorre
a boundary
2011-07-20pulled a lot of the photon propagation code out of src/kernel.cu into ↵Anthony LaTorre
src/photon.h so that photon propagation by propagate() in kernel.cu and the hybrid monte carlo ray tracing use the same code. instead of a single state, photons now carry the history of the processes they've undergone. this history is stored as a bitmask; see src/photon.h. start_node and first_node of the mesh are now stored as global variables in mesh.h instead of being passed to kernel functions.
2011-07-19removed STL pmt models; pmt models are now built by calling rotate_extrude() ↵Anthony LaTorre
on a profile of the PMT model (see build_pmt() in solids/pmts.py). triangle intersection now allows one of the two coefficients multiplying the vectors which span the triangle to float slightly negative (up to -EPSILON; EPSILON is defined in src/intersect.h) in order to eliminate rays passing through the line between two triangles. cleaned up a lot of unused code. pulled duplicate code in view() and render() into functions in view.py. in order to allow view.py and render.py to search pre-defined geometries, solids, meshes, etc. without requiring them to be pre-built, pre-defined geometries, solids, meshes, etc. should be returned by a function tagged by the decorator @buildable(identifier) defined in view.py, where identifier is a string used to identify the object as an argument to either view.py or render.py. optical materials and surfaces are now defined in optics.py. added an image directory to save cool screenshots.
2011-07-10added a hybrid monte carlo ray tracing rendering algorithmAnthony LaTorre
2011-06-24argument '-j' to threadtest.py now specifies a list of device ids to be ↵Anthony LaTorre
used. GPUThread objects now shallow copy the geometry so that threads are not writing to the same memory when the geometry is loaded onto the gpu. the model number for the 12" Hamamatsu PMT is R11708, not r7081 (which is the model for the 10" PMT); all references to the 12" PMT were changed accordingly. only allocate space for 20 materials and 20 surfaces on the gpu instead of 100 to save some space. started to modify track.py to build its own photons and module since the GPUThread object only copies photon hit times back from the gpu (not track information), but I am waiting to find out if pycuda GPUArrays can be used with vector types.
2011-06-24speedup mesh intersection by skipping directly to the child node of nodes ↵Anthony LaTorre
with only one child.
2011-06-23move photon initialization to the gpu. it's unclear if this is a speed ↵Anthony LaTorre
improvement.
2011-06-23Replace calculation of likelihood variance from uncertainties becauseStan Seibert
it is too slow. Now we do it directly.
2011-06-22Rescale the number of PMTs and dimensions of LBNE to match the latestStan Seibert
Water Cherenkov case study.
2011-06-22Start jittering the PMT times in the DAQ simulation.Stan Seibert
2011-06-22Switch from float4 to float3 for vertices.Stan Seibert
2011-06-22Fix accidental reuse of absorption_length when computing scattering length.Stan Seibert
2011-06-21speedup geometry.build() by using more iterators and keeping a list of ↵Anthony LaTorre
unique materials and surfaces in each solid.
2011-06-21Merge ancient headStan Seibert
2011-06-21Merge speedup of build()Stan Seibert
2011-06-21added rotate_extrude() and a number of functions to produce meshes of useful ↵Anthony LaTorre
solid shapes. tried to speed up geometry.build() by replacing some for loops with np.fromiter() and imap().
2011-06-21Set a detector object variable near the top of the main function inStan Seibert
threadtest.py to make it easier to flip between different models.
2011-06-21Switch vertex list on GPU from a texture to a standard device array.Stan Seibert
This avoids the texture size limitation.
2011-06-21Create a r7081_cut mesh that uses the cut PMT model.Stan Seibert
Modify LBNE class to take a cut_pmt option to switch to cut PMT model. Add corresponding lbne_cut, minilbne_cut, microlbne_cut detectors.
2011-06-21Fix PMT models to work with OpenSCAD, add reduced models cut justStan Seibert
behind the PMT equator to save memory. Note: To repair a PMT produced by Sketchup, run Netfabb, start repair mode, and run the "Remove degenerate faces" action, and apply repair. Now you can export the part and OpenSCAD will be happy.
2011-06-21A GPU-side "DAQ" implementation that identifies the first photon on each ↵Stan Seibert
channel, which is presumed to trigger that channel. Major speed up in conversion of detection times to time PDFs.
2011-06-21Add a microlbne detector for debugging purposes.Stan Seibert
2011-06-20pack material and surface indices into the fourth byte of the triangle array ↵Anthony LaTorre
on the GPU. you can now take a screenshot of an image rendered with view.py() by pressing the f12 key.
2011-06-19reverse the face orientation on all triangles in the stl files exported from ↵Anthony LaTorre
sketchup.
2011-06-18stupid mistake. event_times should not be converted to an array when it is ↵Anthony LaTorre
returned from generate_event()
2011-06-18change curandstate array length to 100,000. gputhread.run() now correctly ↵Anthony LaTorre
initializes ALL of the states.
2011-06-18forgot to initialize gpu random number states in gputhread.run()Anthony LaTorre
2011-06-18moved class definitions for Solid, Mesh, Material, and Surface into ↵Anthony LaTorre
geometry.py and moved instances of these classes into separate folders. the Solid object no longer contains a rotation, displacement, or id variable; instead, they are passed to a geometry object when calling add_solid().
2011-06-17visually tested optics code. added models of the inner and outer meshes for ↵Anthony LaTorre
the 12" hamamatsu and sno pmts. ratdb.py is able to parse ratdb files. chromaticity.py provides a function to map wavelength -> rgb color. lbne detector model now includes an outer black cylinder and pmts with a glass layer and photocathode/reflective surfaces.
2011-06-12added some fun models; added some untested code to implement absorption, ↵Anthony LaTorre
scattering, reflection, and refraction
2011-06-07Switch triangle texture to device array, use int32 and float32Stan Seibert
datatypes everywhere, and build final mesh without concatenation of lists. This allows for very large detectors, like full size LBNE.
2011-06-06Thin lens with radius of curvature of 5 units, and thickness of 2 units.Stan Seibert
2011-06-03load material/surface index lookup arrays to the gpu and bind them to ↵Anthony LaTorre
textures. also, forgot to include the gpu code for material/surface structures