Robert Maynard Kitware, Inc. Kenneth Moreland Sandia National Laboratories

Transcription

1 Robert Maynard Kitware, Inc. Kenneth Moreland Sandia National Laboratories

2 Data Analysis at Extreme Data Analysis and Visualization library Multi- and Many-core processor support Allows users to create visualization algorithms

3 So What Can It Do? Threshold Tetrahedralize Cut Marching Cubes Slice And many more

4 How Can It Do All That? 1. Device Adapter 2. Worklet 3. Scheduler 4. Array Handle

5 Device Adapter

6 Device Adapter Allows Dax Worklets to be executed on any backend: CUDA TBB OpenMP Serial CPU Users can add new Device Adapters.

7 Device Adapter Transfer Control Environment Transfer Execu0on Environment Schedule worklet Schedule worklet worklet worklet worklet worklet worklet worklet worklet Scan Sort Compute Compute Other Support Algorithms Stream Compact, Copy, Parallel Find, Unique

8 Worklet

9 struct Sine: public dax::exec::workletmapfield { typedef void ControlSignature(Field(In), Field(Out)); typedef _2 ExecutionSignature(_1); DAX_EXEC_EXPORT dax::scalar operator()(dax::scalar v) const { return dax::math::sin(v); } };

10 Scheduler and Array Handle

11 struct Sine: public dax::exec::workletmapfield { typedef void ControlSignature(Field(In), Field(Out)); typedef _2 ExecutionSignature(_1); DAX_EXEC_EXPORT dax::scalar operator()(dax::scalar v) const { return dax::math::sin(v); } }; Control Environment Execution Environment dax::cont::arrayhandle<dax::scalar> inputhandle = dax::cont::make_arrayhandle(input); dax::cont::arrayhandle<dax::scalar> result; dax::cont::scheduler<> scheduler; scheduler.invoke(sine(), inputhandle, result);

12 Sign Me Up! Code available from these git repositories

13 Prerequisites Git CMake (2.8.8 or newer) Boost ( headers only) Optionally: CUDA Toolkit 4+ (for CUDA Device Adapter) Thrust ( for CUDA and OpenMP Device Adapter) Intel TBB ( for TBB Device Adapter)

14 Let s Get Dax git clone mkdir daxbuild cd daxbuild ccmake../daxtoolkit //build tests make ctest

15 Let s Make Something

16 Determine the task size Worklets operate on fields, cells, or points Think about decomposing the problem into multiple worklets

17 Our task Surface Normals of a Triangle Mesh

18 struct Normals: public dax::exec::workletmapcell { typedef void ControlSignature(Topology, Field(In,Point), Field(Out)); typedef _3 ExecutionSignature(_2); DAX_EXEC_EXPORT dax::vector3 operator()(const dax::cellfield< dax::celltagtriangle, dax::vector3 >& coordinates) const { const dax::vector3 a = coordinates[1] coordinates[0]; const dax::vector3 b = coordinates[2] coordinates[0]; const dax::vector3 c = dax::math::cross(a,b); return dax::math::normalize(c); } };

19 struct Normals: public dax::exec::workletmapcell { typedef void ControlSignature(Topology, Field(In,Point), Field(Out)); typedef _3 ExecutionSignature(_2); DAX_EXEC_EXPORT dax::vector3 operator()(const dax::cellfield< dax::celltagtriangle, dax::vector3 >& coordinates) const { const dax::vector3 a = coordinates[1] coordinates[0]; const dax::vector3 b = coordinates[2] coordinates[0]; const dax::vector3 c = dax::math::cross(a,b); return dax::math::normalize(c); } };

20 struct Normals: public dax::exec::workletmapcell { typedef void ControlSignature(Topology, Field(In,Point), Field(Out)); typedef _3 ExecutionSignature(_2); DAX_EXEC_EXPORT dax::vector3 operator()(const dax::cellfield< dax::celltagtriangle, dax::vector3 >& coordinates) const { const dax::vector3 a = coordinates[1] coordinates[0]; const dax::vector3 b = coordinates[2] coordinates[0]; const dax::vector3 c = dax::math::cross(a,b); return dax::math::normalize(c); } };

21 struct Normals: public dax::exec::workletmapcell { typedef void ControlSignature(Topology, Field(In,Point), Field(Out)); typedef _3 ExecutionSignature(_2); DAX_EXEC_EXPORT dax::vector3 operator()(const dax::cellfield< dax::celltagtriangle, dax::vector3 >& coordinates) const { const dax::vector3 a = coordinates[1] coordinates[0]; const dax::vector3 b = coordinates[2] coordinates[0]; const dax::vector3 c = dax::math::cross(a,b); return dax::math::normalize(c); } };

22 Let s Run The Worklet

23 struct Normals: public dax::exec::workletmapcell { typedef void ControlSignature(Topology,Field(In,Point),Field(Out)); typedef _3 ExecutionSignature(_2); DAX_EXEC_EXPORT dax::vector3 operator()(const dax::cellfield<dax::celltagtriangle,dax::vector3>& coordinates) const { const dax::vector3 a = coordinates[1] coordinates[0]; Execution Environment Control Environment dax::cont::unstructuredgrid<dax::celltagtriangle> grid; dax::cont::arrayhandle<dax::vector3> normals_result; dax::cont::scheduler<> scheduler; scheduler.invoke(normals(), grid, normals_result);

24 What About Performance?

25 Example: Threshold Implemented as two worklets. First worklet determines which cells will pass Second worklet generates new topology and resolves duplicate point IDs.

26

27 Now Without Point Resolution

28

29

30 Acknowledgements This work was supported by the DOE Office of Science, Advanced Scientific Computing Research, under award number , program manager Lucy Nowell. Additional support by the Director, Office of Advanced Scientific Computing Research, Office of Science, of the U.S. Department of Energy under Contract No , through the Scientific Discovery through Advanced Computing (SciDAC) Institute of Scalable Data Management, Analysis and Visualization.