warp-lang 1.6.1__py3-none-macosx_10_13_universal2.whl → 1.7.0__py3-none-macosx_10_13_universal2.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Potentially problematic release.
This version of warp-lang might be problematic. Click here for more details.
- warp/__init__.py +21 -7
- warp/autograd.py +14 -6
- warp/bin/libwarp-clang.dylib +0 -0
- warp/bin/libwarp.dylib +0 -0
- warp/build.py +424 -6
- warp/build_dll.py +20 -20
- warp/builtins.py +467 -368
- warp/codegen.py +193 -125
- warp/config.py +56 -12
- warp/constants.py +14 -6
- warp/context.py +524 -277
- warp/dlpack.py +22 -12
- warp/examples/__init__.py +14 -6
- warp/examples/assets/nonuniform.usd +0 -0
- warp/examples/assets/nvidia_logo.png +0 -0
- warp/examples/benchmarks/benchmark_api.py +14 -6
- warp/examples/benchmarks/benchmark_cloth.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_cupy.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_jax.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_numba.py +15 -0
- warp/examples/benchmarks/benchmark_cloth_numpy.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_paddle.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_pytorch.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_taichi.py +14 -6
- warp/examples/benchmarks/benchmark_cloth_warp.py +14 -6
- warp/examples/benchmarks/benchmark_gemm.py +82 -48
- warp/examples/benchmarks/benchmark_interop_paddle.py +14 -6
- warp/examples/benchmarks/benchmark_interop_torch.py +14 -6
- warp/examples/benchmarks/benchmark_launches.py +14 -6
- warp/examples/benchmarks/benchmark_tile_load_store.py +103 -0
- warp/examples/browse.py +14 -6
- warp/examples/core/example_cupy.py +14 -6
- warp/examples/core/example_dem.py +14 -6
- warp/examples/core/example_fluid.py +14 -6
- warp/examples/core/example_graph_capture.py +14 -6
- warp/examples/core/example_marching_cubes.py +14 -6
- warp/examples/core/example_mesh.py +14 -6
- warp/examples/core/example_mesh_intersect.py +14 -6
- warp/examples/core/example_nvdb.py +14 -6
- warp/examples/core/example_raycast.py +14 -6
- warp/examples/core/example_raymarch.py +14 -6
- warp/examples/core/example_render_opengl.py +14 -6
- warp/examples/core/example_sample_mesh.py +300 -0
- warp/examples/core/example_sph.py +14 -6
- warp/examples/core/example_torch.py +14 -6
- warp/examples/core/example_wave.py +14 -6
- warp/examples/fem/example_adaptive_grid.py +14 -6
- warp/examples/fem/example_apic_fluid.py +15 -7
- warp/examples/fem/example_burgers.py +16 -8
- warp/examples/fem/example_convection_diffusion.py +14 -6
- warp/examples/fem/example_convection_diffusion_dg.py +14 -6
- warp/examples/fem/example_deformed_geometry.py +15 -7
- warp/examples/fem/example_diffusion.py +14 -6
- warp/examples/fem/example_diffusion_3d.py +14 -6
- warp/examples/fem/example_diffusion_mgpu.py +14 -6
- warp/examples/fem/example_distortion_energy.py +15 -7
- warp/examples/fem/example_magnetostatics.py +20 -12
- warp/examples/fem/example_mixed_elasticity.py +14 -6
- warp/examples/fem/example_navier_stokes.py +14 -6
- warp/examples/fem/example_nonconforming_contact.py +14 -6
- warp/examples/fem/example_stokes.py +14 -6
- warp/examples/fem/example_stokes_transfer.py +14 -6
- warp/examples/fem/example_streamlines.py +14 -6
- warp/examples/fem/utils.py +24 -3
- warp/examples/interop/example_jax_callable.py +116 -0
- warp/examples/interop/example_jax_ffi_callback.py +132 -0
- warp/examples/interop/example_jax_kernel.py +205 -0
- warp/examples/optim/example_bounce.py +14 -6
- warp/examples/optim/example_cloth_throw.py +14 -6
- warp/examples/optim/example_diffray.py +14 -6
- warp/examples/optim/example_drone.py +14 -6
- warp/examples/optim/example_fluid_checkpoint.py +497 -0
- warp/examples/optim/example_inverse_kinematics.py +14 -6
- warp/examples/optim/example_inverse_kinematics_torch.py +14 -6
- warp/examples/optim/example_softbody_properties.py +14 -6
- warp/examples/optim/example_spring_cage.py +14 -6
- warp/examples/optim/example_trajectory.py +14 -6
- warp/examples/sim/example_cartpole.py +14 -6
- warp/examples/sim/example_cloth.py +14 -6
- warp/examples/sim/example_cloth_self_contact.py +14 -6
- warp/examples/sim/example_granular.py +14 -6
- warp/examples/sim/example_granular_collision_sdf.py +14 -6
- warp/examples/sim/example_jacobian_ik.py +14 -6
- warp/examples/sim/example_particle_chain.py +14 -6
- warp/examples/sim/example_quadruped.py +14 -6
- warp/examples/sim/example_rigid_chain.py +14 -6
- warp/examples/sim/example_rigid_contact.py +14 -6
- warp/examples/sim/example_rigid_force.py +14 -6
- warp/examples/sim/example_rigid_gyroscopic.py +14 -6
- warp/examples/sim/example_rigid_soft_contact.py +14 -6
- warp/examples/sim/example_soft_body.py +14 -6
- warp/examples/tile/example_tile_cholesky.py +14 -6
- warp/examples/tile/example_tile_convolution.py +14 -6
- warp/examples/tile/example_tile_fft.py +14 -6
- warp/examples/tile/example_tile_filtering.py +14 -6
- warp/examples/tile/example_tile_matmul.py +16 -10
- warp/examples/tile/example_tile_mlp.py +14 -6
- warp/examples/tile/example_tile_nbody.py +14 -6
- warp/examples/tile/example_tile_walker.py +14 -6
- warp/fabric.py +15 -0
- warp/fem/__init__.py +26 -1
- warp/fem/adaptivity.py +19 -4
- warp/fem/cache.py +15 -0
- warp/fem/dirichlet.py +15 -0
- warp/fem/domain.py +15 -0
- warp/fem/field/__init__.py +15 -0
- warp/fem/field/field.py +15 -0
- warp/fem/field/nodal_field.py +37 -68
- warp/fem/field/restriction.py +15 -0
- warp/fem/field/virtual.py +77 -23
- warp/fem/geometry/__init__.py +15 -0
- warp/fem/geometry/adaptive_nanogrid.py +24 -10
- warp/fem/geometry/closest_point.py +16 -1
- warp/fem/geometry/deformed_geometry.py +20 -2
- warp/fem/geometry/element.py +15 -0
- warp/fem/geometry/geometry.py +20 -0
- warp/fem/geometry/grid_2d.py +27 -12
- warp/fem/geometry/grid_3d.py +27 -15
- warp/fem/geometry/hexmesh.py +20 -7
- warp/fem/geometry/nanogrid.py +24 -11
- warp/fem/geometry/partition.py +15 -0
- warp/fem/geometry/quadmesh.py +28 -13
- warp/fem/geometry/tetmesh.py +18 -4
- warp/fem/geometry/trimesh.py +18 -8
- warp/fem/integrate.py +277 -93
- warp/fem/linalg.py +20 -5
- warp/fem/operator.py +15 -0
- warp/fem/polynomial.py +15 -0
- warp/fem/quadrature/__init__.py +15 -0
- warp/fem/quadrature/pic_quadrature.py +52 -22
- warp/fem/quadrature/quadrature.py +209 -25
- warp/fem/space/__init__.py +16 -1
- warp/fem/space/basis_function_space.py +19 -2
- warp/fem/space/basis_space.py +40 -18
- warp/fem/space/dof_mapper.py +15 -0
- warp/fem/space/function_space.py +15 -0
- warp/fem/space/grid_2d_function_space.py +15 -0
- warp/fem/space/grid_3d_function_space.py +15 -0
- warp/fem/space/hexmesh_function_space.py +17 -2
- warp/fem/space/nanogrid_function_space.py +15 -0
- warp/fem/space/partition.py +21 -2
- warp/fem/space/quadmesh_function_space.py +23 -8
- warp/fem/space/restriction.py +15 -0
- warp/fem/space/shape/__init__.py +15 -0
- warp/fem/space/shape/cube_shape_function.py +38 -23
- warp/fem/space/shape/shape_function.py +15 -0
- warp/fem/space/shape/square_shape_function.py +27 -12
- warp/fem/space/shape/tet_shape_function.py +15 -0
- warp/fem/space/shape/triangle_shape_function.py +16 -1
- warp/fem/space/tetmesh_function_space.py +18 -3
- warp/fem/space/topology.py +15 -0
- warp/fem/space/trimesh_function_space.py +17 -2
- warp/fem/types.py +15 -0
- warp/fem/utils.py +27 -6
- warp/jax.py +28 -7
- warp/jax_experimental/__init__.py +16 -0
- warp/{jax_experimental.py → jax_experimental/custom_call.py} +28 -33
- warp/jax_experimental/ffi.py +698 -0
- warp/jax_experimental/xla_ffi.py +602 -0
- warp/math.py +103 -6
- warp/native/array.h +28 -6
- warp/native/builtin.h +44 -9
- warp/native/bvh.cpp +18 -7
- warp/native/bvh.cu +57 -20
- warp/native/bvh.h +17 -7
- warp/native/clang/clang.cpp +45 -9
- warp/native/coloring.cpp +15 -6
- warp/native/crt.cpp +15 -6
- warp/native/crt.h +15 -6
- warp/native/cuda_crt.h +15 -6
- warp/native/cuda_util.cpp +29 -6
- warp/native/cuda_util.h +17 -6
- warp/native/error.cpp +15 -6
- warp/native/error.h +15 -6
- warp/native/exports.h +85 -63
- warp/native/fabric.h +15 -6
- warp/native/hashgrid.cpp +15 -6
- warp/native/hashgrid.cu +15 -6
- warp/native/hashgrid.h +15 -6
- warp/native/initializer_array.h +15 -6
- warp/native/intersect.h +41 -32
- warp/native/intersect_adj.h +48 -39
- warp/native/intersect_tri.h +17 -0
- warp/native/marching.cpp +16 -0
- warp/native/marching.cu +16 -7
- warp/native/marching.h +17 -0
- warp/native/mat.h +528 -15
- warp/native/mathdx.cpp +15 -6
- warp/native/matnn.h +15 -6
- warp/native/mesh.cpp +15 -6
- warp/native/mesh.cu +15 -6
- warp/native/mesh.h +25 -16
- warp/native/noise.h +15 -6
- warp/native/quat.h +114 -17
- warp/native/rand.h +21 -6
- warp/native/range.h +15 -6
- warp/native/reduce.cpp +15 -6
- warp/native/reduce.cu +15 -6
- warp/native/runlength_encode.cpp +15 -6
- warp/native/runlength_encode.cu +15 -6
- warp/native/scan.cpp +15 -6
- warp/native/scan.cu +15 -6
- warp/native/scan.h +15 -6
- warp/native/solid_angle.h +17 -0
- warp/native/sort.cpp +137 -65
- warp/native/sort.cu +167 -21
- warp/native/sort.h +23 -7
- warp/native/sparse.cpp +58 -28
- warp/native/sparse.cu +67 -23
- warp/native/spatial.h +15 -6
- warp/native/svd.h +131 -6
- warp/native/temp_buffer.h +15 -6
- warp/native/tile.h +316 -111
- warp/native/tile_reduce.h +61 -9
- warp/native/vec.h +83 -13
- warp/native/volume.cpp +100 -119
- warp/native/volume.cu +15 -6
- warp/native/volume.h +15 -6
- warp/native/volume_builder.cu +40 -16
- warp/native/volume_builder.h +21 -6
- warp/native/volume_impl.h +15 -6
- warp/native/warp.cpp +20 -12
- warp/native/warp.cu +114 -16
- warp/native/warp.h +34 -16
- warp/optim/__init__.py +14 -6
- warp/optim/adam.py +14 -6
- warp/optim/linear.py +25 -10
- warp/optim/sgd.py +14 -6
- warp/paddle.py +14 -6
- warp/render/__init__.py +14 -6
- warp/render/render_opengl.py +14 -6
- warp/render/render_usd.py +14 -6
- warp/render/utils.py +14 -6
- warp/sim/__init__.py +14 -7
- warp/sim/articulation.py +18 -10
- warp/sim/collide.py +35 -16
- warp/sim/graph_coloring.py +14 -6
- warp/sim/import_mjcf.py +463 -162
- warp/sim/import_snu.py +14 -7
- warp/sim/import_urdf.py +46 -18
- warp/sim/import_usd.py +14 -7
- warp/sim/inertia.py +14 -6
- warp/sim/integrator.py +14 -6
- warp/sim/integrator_euler.py +19 -11
- warp/sim/integrator_featherstone.py +17 -16
- warp/sim/integrator_vbd.py +222 -8
- warp/sim/integrator_xpbd.py +19 -11
- warp/sim/model.py +56 -19
- warp/sim/particles.py +14 -6
- warp/sim/render.py +14 -6
- warp/sim/utils.py +17 -2
- warp/sparse.py +657 -555
- warp/stubs.py +231 -19
- warp/tape.py +14 -6
- warp/tests/aux_test_class_kernel.py +14 -6
- warp/tests/aux_test_compile_consts_dummy.py +14 -6
- warp/tests/aux_test_conditional_unequal_types_kernels.py +14 -6
- warp/tests/aux_test_dependent.py +14 -6
- warp/tests/aux_test_grad_customs.py +14 -6
- warp/tests/aux_test_instancing_gc.py +14 -6
- warp/tests/aux_test_module_unload.py +14 -6
- warp/tests/aux_test_name_clash1.py +14 -6
- warp/tests/aux_test_name_clash2.py +14 -6
- warp/tests/aux_test_unresolved_func.py +14 -6
- warp/tests/aux_test_unresolved_symbol.py +14 -6
- warp/tests/cuda/__init__.py +0 -0
- warp/tests/{test_async.py → cuda/test_async.py} +14 -6
- warp/tests/{test_ipc.py → cuda/test_ipc.py} +14 -6
- warp/tests/{test_mempool.py → cuda/test_mempool.py} +53 -6
- warp/tests/{test_multigpu.py → cuda/test_multigpu.py} +14 -6
- warp/tests/{test_peer.py → cuda/test_peer.py} +14 -6
- warp/tests/{test_pinned.py → cuda/test_pinned.py} +14 -6
- warp/tests/{test_streams.py → cuda/test_streams.py} +85 -6
- warp/tests/geometry/__init__.py +0 -0
- warp/tests/{test_bvh.py → geometry/test_bvh.py} +14 -6
- warp/tests/{test_hash_grid.py → geometry/test_hash_grid.py} +14 -6
- warp/tests/{test_marching_cubes.py → geometry/test_marching_cubes.py} +14 -6
- warp/tests/{test_mesh.py → geometry/test_mesh.py} +14 -6
- warp/tests/{test_mesh_query_aabb.py → geometry/test_mesh_query_aabb.py} +14 -6
- warp/tests/{test_mesh_query_point.py → geometry/test_mesh_query_point.py} +80 -69
- warp/tests/{test_mesh_query_ray.py → geometry/test_mesh_query_ray.py} +15 -7
- warp/tests/{test_volume.py → geometry/test_volume.py} +55 -12
- warp/tests/{test_volume_write.py → geometry/test_volume_write.py} +14 -6
- warp/tests/interop/__init__.py +0 -0
- warp/tests/{test_dlpack.py → interop/test_dlpack.py} +42 -11
- warp/tests/{test_jax.py → interop/test_jax.py} +14 -6
- warp/tests/{test_paddle.py → interop/test_paddle.py} +14 -6
- warp/tests/{test_torch.py → interop/test_torch.py} +14 -6
- warp/tests/run_coverage_serial.py +14 -6
- warp/tests/sim/__init__.py +0 -0
- warp/tests/{disabled_kinematics.py → sim/disabled_kinematics.py} +23 -16
- warp/tests/{flaky_test_sim_grad.py → sim/flaky_test_sim_grad.py} +14 -6
- warp/tests/{test_collision.py → sim/test_collision.py} +16 -8
- warp/tests/{test_coloring.py → sim/test_coloring.py} +14 -7
- warp/tests/{test_model.py → sim/test_model.py} +55 -7
- warp/tests/{test_sim_grad_bounce_linear.py → sim/test_sim_grad_bounce_linear.py} +14 -6
- warp/tests/{test_sim_kinematics.py → sim/test_sim_kinematics.py} +16 -7
- warp/tests/sim/test_vbd.py +597 -0
- warp/tests/test_adam.py +14 -6
- warp/tests/test_arithmetic.py +14 -6
- warp/tests/test_array.py +14 -6
- warp/tests/test_array_reduce.py +14 -6
- warp/tests/test_assert.py +14 -6
- warp/tests/test_atomic.py +14 -6
- warp/tests/test_bool.py +15 -7
- warp/tests/test_builtins_resolution.py +14 -6
- warp/tests/test_closest_point_edge_edge.py +14 -6
- warp/tests/test_codegen.py +14 -6
- warp/tests/test_codegen_instancing.py +14 -6
- warp/tests/test_compile_consts.py +14 -6
- warp/tests/test_conditional.py +14 -6
- warp/tests/test_context.py +14 -6
- warp/tests/test_copy.py +14 -6
- warp/tests/test_ctypes.py +14 -6
- warp/tests/test_dense.py +14 -6
- warp/tests/test_devices.py +14 -6
- warp/tests/test_examples.py +42 -42
- warp/tests/test_fabricarray.py +14 -6
- warp/tests/test_fast_math.py +14 -6
- warp/tests/test_fem.py +37 -10
- warp/tests/test_fp16.py +14 -6
- warp/tests/test_func.py +14 -6
- warp/tests/test_future_annotations.py +14 -6
- warp/tests/test_generics.py +14 -6
- warp/tests/test_grad.py +14 -6
- warp/tests/test_grad_customs.py +14 -6
- warp/tests/test_grad_debug.py +14 -6
- warp/tests/test_implicit_init.py +14 -6
- warp/tests/test_import.py +14 -6
- warp/tests/test_indexedarray.py +14 -6
- warp/tests/test_intersect.py +14 -6
- warp/tests/test_iter.py +14 -6
- warp/tests/test_large.py +14 -6
- warp/tests/test_launch.py +14 -6
- warp/tests/test_lerp.py +14 -6
- warp/tests/test_linear_solvers.py +15 -11
- warp/tests/test_lvalue.py +14 -6
- warp/tests/test_mat.py +247 -85
- warp/tests/test_mat_lite.py +14 -6
- warp/tests/test_mat_scalar_ops.py +18 -10
- warp/tests/test_math.py +14 -6
- warp/tests/test_mlp.py +14 -6
- warp/tests/test_module_hashing.py +14 -6
- warp/tests/test_modules_lite.py +14 -6
- warp/tests/test_noise.py +14 -6
- warp/tests/test_operators.py +14 -6
- warp/tests/test_options.py +14 -6
- warp/tests/test_overwrite.py +15 -60
- warp/tests/test_print.py +14 -6
- warp/tests/test_quat.py +81 -52
- warp/tests/test_rand.py +58 -43
- warp/tests/test_reload.py +14 -6
- warp/tests/test_rounding.py +14 -6
- warp/tests/test_runlength_encode.py +14 -6
- warp/tests/test_scalar_ops.py +14 -6
- warp/tests/test_smoothstep.py +14 -6
- warp/tests/test_snippet.py +15 -0
- warp/tests/test_sparse.py +61 -12
- warp/tests/test_spatial.py +89 -6
- warp/tests/test_special_values.py +14 -6
- warp/tests/test_static.py +15 -7
- warp/tests/test_struct.py +14 -6
- warp/tests/test_tape.py +14 -6
- warp/tests/test_transient_module.py +14 -6
- warp/tests/test_triangle_closest_point.py +14 -6
- warp/tests/test_types.py +14 -6
- warp/tests/test_utils.py +98 -10
- warp/tests/test_vec.py +60 -40
- warp/tests/test_vec_lite.py +14 -6
- warp/tests/test_vec_scalar_ops.py +14 -6
- warp/tests/test_verify_fp.py +14 -6
- warp/tests/tile/__init__.py +0 -0
- warp/tests/{test_tile.py → tile/test_tile.py} +150 -57
- warp/tests/{test_tile_load.py → tile/test_tile_load.py} +15 -7
- warp/tests/{test_tile_mathdx.py → tile/test_tile_mathdx.py} +23 -12
- warp/tests/{test_tile_mlp.py → tile/test_tile_mlp.py} +39 -20
- warp/tests/{test_tile_reduce.py → tile/test_tile_reduce.py} +74 -7
- warp/tests/{test_tile_shared_memory.py → tile/test_tile_shared_memory.py} +14 -6
- warp/tests/{test_tile_view.py → tile/test_tile_view.py} +15 -7
- warp/tests/unittest_serial.py +15 -6
- warp/tests/unittest_suites.py +59 -65
- warp/tests/unittest_utils.py +16 -7
- warp/tests/walkthrough_debug.py +14 -6
- warp/thirdparty/unittest_parallel.py +15 -8
- warp/torch.py +14 -6
- warp/types.py +124 -664
- warp/utils.py +151 -78
- {warp_lang-1.6.1.dist-info → warp_lang-1.7.0.dist-info}/METADATA +39 -12
- warp_lang-1.7.0.dist-info/RECORD +429 -0
- {warp_lang-1.6.1.dist-info → warp_lang-1.7.0.dist-info}/WHEEL +1 -1
- warp_lang-1.7.0.dist-info/licenses/LICENSE.md +202 -0
- warp/examples/optim/example_walker.py +0 -309
- warp/native/cutlass_gemm.cpp +0 -34
- warp/native/cutlass_gemm.cu +0 -373
- warp/tests/test_matmul.py +0 -503
- warp/tests/test_matmul_lite.py +0 -403
- warp/tests/test_vbd.py +0 -378
- warp/tests/unused_test_misc.py +0 -69
- warp_lang-1.6.1.dist-info/LICENSE.md +0 -126
- warp_lang-1.6.1.dist-info/RECORD +0 -419
- {warp_lang-1.6.1.dist-info → warp_lang-1.7.0.dist-info}/top_level.txt +0 -0
|
@@ -0,0 +1,300 @@
|
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2025 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
15
|
+
|
|
16
|
+
###########################################################################
|
|
17
|
+
# Example Sample Mesh
|
|
18
|
+
#
|
|
19
|
+
# Shows how to sample points on a mesh's surface using
|
|
20
|
+
# a Cumulative Distribution Function (CDF).
|
|
21
|
+
#
|
|
22
|
+
# The CDF enables uniform sampling of points across the mesh's surface,
|
|
23
|
+
# even when the density of triangles varies. It represents the cumulative
|
|
24
|
+
# probability of selecting a triangle from the mesh, with each triangle
|
|
25
|
+
# weighted by its area relative to the total surface area of the mesh.
|
|
26
|
+
#
|
|
27
|
+
###########################################################################
|
|
28
|
+
|
|
29
|
+
import numpy as np
|
|
30
|
+
|
|
31
|
+
import warp as wp
|
|
32
|
+
import warp.render
|
|
33
|
+
|
|
34
|
+
# fmt: off
|
|
35
|
+
POINTS = np.array(
|
|
36
|
+
(
|
|
37
|
+
(-0.986598, -0.400638, -0.175759), (-0.81036 , -0.482105, -0.541125),
|
|
38
|
+
(-1.079616, 0.022652, -0.023381), (-0.894468, -0.080795, -0.618379),
|
|
39
|
+
(-0.607365, -0.702012, -0.556551), (-0.366107, -0.800096, -0.620734),
|
|
40
|
+
(-0.801777, -0.690991, -0.239593), (-0.553576, -0.871746, -0.335518),
|
|
41
|
+
(-0.309133, -0.370805, -0.965784), (-0.288299, -0.956987, -0.402091),
|
|
42
|
+
(-0.051878, -0.894342, -0.597583), (-0.386774, -1.003107, -0.145116),
|
|
43
|
+
(-0.19062 , -1.061165, 0.012418), (-0.176053, -1.044838, -0.217194),
|
|
44
|
+
( 0.001479, -1.020045, -0.356905), (-0.105375, -0.655117, -0.861365),
|
|
45
|
+
(-0.542102, -0.517255, -0.795259), (-0.476599, -0.105709, -0.981171),
|
|
46
|
+
(-1.047915, -0.121584, 0.322098), (-0.527852, 0.137252, 0.501813),
|
|
47
|
+
(-0.721762, -0.803275, 0.117162), (-0.904992, -0.573281, 0.168408),
|
|
48
|
+
(-0.796762, -0.473428, 0.569649), (-0.606446, -0.753374, 0.492938),
|
|
49
|
+
(-0.466481, -0.576566, 0.802562), (-0.50476 , -0.908596, 0.300064),
|
|
50
|
+
(-0.337425, -1.008902, 0.170911), (-0.048676, -1.055594, 0.246732),
|
|
51
|
+
(-0.212871, -0.760442, 0.738447), (-0.281356, -0.9322 , 0.474965),
|
|
52
|
+
(-0.560476, 0.062512, -0.561019), (-0.003252, 0.083237, -1.049784),
|
|
53
|
+
(-0.009392, 0.593703, -0.522479), (-0.530465, 0.577231, 0.007172),
|
|
54
|
+
(-0.02106 , 0.064189, 1.066722), (-0.003512, 0.59714 , 0.516904),
|
|
55
|
+
( 0.000194, 1.093899, 0.001113), ( 0.256861, -0.955856, -0.445325),
|
|
56
|
+
( 0.251205, -1.038759, -0.174212), ( 0.170201, -0.800019, -0.712158),
|
|
57
|
+
( 0.364385, -0.560298, -0.866843), ( 0.092809, -0.269437, -1.058467),
|
|
58
|
+
( 0.628127, -0.12359 , -0.9012 ), ( 0.507433, -0.930658, -0.215908),
|
|
59
|
+
( 0.496448, -0.800205, -0.545904), ( 0.757415, -0.527449, -0.565395),
|
|
60
|
+
( 0.908704, -0.596257, 0.028995), ( 0.754069, -0.731365, -0.256687),
|
|
61
|
+
( 0.921362, -0.09028 , -0.546421), ( 1.017846, -0.335787, -0.263017),
|
|
62
|
+
( 0.016768, -1.080014, -0.058473), ( 0.204245, -1.056388, 0.078346),
|
|
63
|
+
( 0.260892, -1.001704, 0.322104), ( 0.16608 , -0.739172, 0.788097),
|
|
64
|
+
( 0.021091, -0.931327, 0.557789), (-0.046158, -0.408417, 1.011046),
|
|
65
|
+
( 0.429623, -0.987237, 0.088537), ( 0.704993, -0.739396, 0.386838),
|
|
66
|
+
( 0.37277 , -0.825639, 0.591102), ( 0.493947, -0.896091, 0.339163),
|
|
67
|
+
( 0.321112, -0.540547, 0.890161), ( 0.654753, -0.520495, 0.690104),
|
|
68
|
+
( 0.922472, -0.124429, 0.530498), ( 0.662544, -0.85601 , 0.054375),
|
|
69
|
+
( 0.950976, -0.422783, 0.327726), ( 0.536849, 0.109943, -0.52279 ),
|
|
70
|
+
( 0.517242, 0.120634, 0.535708), ( 0.532707, 0.598943, -0.000767),
|
|
71
|
+
( 1.086691, 0.048722, 0.032517), ( 0.528734, -0.109809, 0.96863 ),
|
|
72
|
+
(-0.581832, -0.916941, -0.027829), (-0.625071, -0.14445 , 0.906538),
|
|
73
|
+
),
|
|
74
|
+
dtype=np.float32,
|
|
75
|
+
)
|
|
76
|
+
|
|
77
|
+
FACE_VERTEX_INDICES = np.array(
|
|
78
|
+
(
|
|
79
|
+
6, 0, 1, 6, 21, 0, 2, 0, 18, 0, 3, 1, 2, 3, 0, 5,
|
|
80
|
+
7, 4, 70, 7, 11, 4, 6, 1, 16, 1, 3, 7, 6, 4, 4, 1,
|
|
81
|
+
16, 9, 7, 5, 3, 17, 16, 16, 17, 8, 41, 8, 17, 30, 17, 3,
|
|
82
|
+
10, 14, 9, 5, 10, 9, 10, 37, 14, 15, 10, 5, 7, 9, 11, 11,
|
|
83
|
+
9, 13, 11, 13, 12, 50, 12, 13, 9, 14, 13, 15, 16, 8, 15, 8,
|
|
84
|
+
41, 16, 5, 4, 16, 15, 5, 17, 31, 41, 21, 22, 18, 20, 21, 6,
|
|
85
|
+
18, 0, 21, 20, 25, 23, 20, 70, 25, 70, 11, 26, 26, 25, 70, 25,
|
|
86
|
+
29, 23, 21, 20, 23, 21, 23, 22, 23, 24, 22, 24, 71, 22, 26, 29,
|
|
87
|
+
25, 26, 11, 12, 12, 27, 26, 26, 27, 29, 27, 54, 29, 27, 12, 50,
|
|
88
|
+
28, 29, 54, 54, 53, 28, 23, 28, 24, 29, 28, 23, 28, 55, 24, 28,
|
|
89
|
+
53, 55, 53, 60, 55, 24, 55, 71, 55, 34, 71, 30, 3, 2, 2, 33,
|
|
90
|
+
30, 17, 30, 31, 32, 31, 30, 33, 36, 32, 19, 33, 2, 19, 35, 33,
|
|
91
|
+
19, 71, 34, 35, 19, 34, 34, 66, 35, 35, 36, 33, 35, 67, 36, 15,
|
|
92
|
+
39, 10, 10, 39, 37, 44, 37, 39, 14, 50, 13, 14, 38, 50, 14, 37,
|
|
93
|
+
38, 37, 43, 38, 40, 15, 41, 40, 39, 15, 41, 42, 40, 44, 39, 40,
|
|
94
|
+
31, 42, 41, 38, 43, 56, 44, 43, 37, 44, 47, 43, 47, 63, 43, 44,
|
|
95
|
+
40, 45, 42, 45, 40, 46, 63, 47, 45, 47, 44, 65, 48, 42, 46, 47,
|
|
96
|
+
49, 49, 47, 45, 48, 45, 42, 45, 48, 49, 68, 49, 48, 27, 52, 54,
|
|
97
|
+
50, 51, 27, 27, 51, 52, 50, 38, 51, 38, 56, 51, 51, 56, 52, 54,
|
|
98
|
+
52, 58, 52, 59, 58, 53, 54, 58, 60, 69, 55, 55, 69, 34, 43, 63,
|
|
99
|
+
56, 59, 52, 56, 63, 59, 56, 63, 57, 59, 58, 60, 53, 57, 58, 59,
|
|
100
|
+
58, 57, 61, 60, 58, 61, 57, 64, 61, 62, 61, 64, 60, 61, 69, 62,
|
|
101
|
+
69, 61, 46, 57, 63, 64, 57, 46, 46, 49, 64, 68, 64, 49, 62, 64,
|
|
102
|
+
68, 32, 65, 31, 65, 32, 67, 32, 36, 67, 65, 42, 31, 67, 68, 65,
|
|
103
|
+
48, 65, 68, 34, 69, 66, 67, 35, 66, 68, 66, 62, 66, 69, 62, 67,
|
|
104
|
+
66, 68, 33, 32, 30, 19, 2, 18, 20, 6, 70, 7, 70, 6, 18, 71,
|
|
105
|
+
19, 22, 71, 18,
|
|
106
|
+
),
|
|
107
|
+
dtype=np.int32,
|
|
108
|
+
)
|
|
109
|
+
# fmt: on
|
|
110
|
+
|
|
111
|
+
|
|
112
|
+
@wp.kernel(enable_backward=False)
|
|
113
|
+
def compute_tri_areas(
|
|
114
|
+
points: wp.array(dtype=wp.vec3),
|
|
115
|
+
face_vertex_indices: wp.array(dtype=wp.int32),
|
|
116
|
+
out_tri_areas: wp.array(dtype=wp.float32),
|
|
117
|
+
out_total_area: wp.array(dtype=wp.float32),
|
|
118
|
+
):
|
|
119
|
+
tri = wp.tid()
|
|
120
|
+
|
|
121
|
+
# Retrieve the indices of the three vertices that form the current triangle.
|
|
122
|
+
vtx_0 = face_vertex_indices[tri * 3]
|
|
123
|
+
vtx_1 = face_vertex_indices[tri * 3 + 1]
|
|
124
|
+
vtx_2 = face_vertex_indices[tri * 3 + 2]
|
|
125
|
+
|
|
126
|
+
# Retrieve their 3D position.
|
|
127
|
+
pt_0 = points[vtx_0]
|
|
128
|
+
pt_1 = points[vtx_1]
|
|
129
|
+
pt_2 = points[vtx_2]
|
|
130
|
+
|
|
131
|
+
# Calculate the cross product of two edges of the triangle,
|
|
132
|
+
# which gives a vector whose magnitude is twice the area of the triangle.
|
|
133
|
+
cross = wp.cross((pt_1 - pt_0), (pt_2 - pt_0))
|
|
134
|
+
area = wp.length(cross) * 0.5
|
|
135
|
+
|
|
136
|
+
# Store the result.
|
|
137
|
+
out_tri_areas[tri] = area
|
|
138
|
+
wp.atomic_add(out_total_area, 0, area)
|
|
139
|
+
|
|
140
|
+
|
|
141
|
+
@wp.kernel(enable_backward=False)
|
|
142
|
+
def compute_probability_distribution(
|
|
143
|
+
tri_areas: wp.array(dtype=wp.float32),
|
|
144
|
+
total_area: wp.array(dtype=wp.float32),
|
|
145
|
+
out_probabilities: wp.array(dtype=wp.float32),
|
|
146
|
+
):
|
|
147
|
+
tri = wp.tid()
|
|
148
|
+
|
|
149
|
+
# Calculate the probability of selecting this triangle,
|
|
150
|
+
# which is proportional to the triangle's area relative to total mesh area.
|
|
151
|
+
out_probabilities[tri] = tri_areas[tri] / total_area[0]
|
|
152
|
+
|
|
153
|
+
|
|
154
|
+
@wp.kernel(enable_backward=False)
|
|
155
|
+
def accumulate_cdf(
|
|
156
|
+
tri_count: wp.int32,
|
|
157
|
+
out_cdf: wp.array(dtype=wp.float32),
|
|
158
|
+
):
|
|
159
|
+
# Transform probability values into a Cumulative Distribution Function (CDF).
|
|
160
|
+
for tri in range(1, tri_count):
|
|
161
|
+
out_cdf[tri] += out_cdf[tri - 1]
|
|
162
|
+
|
|
163
|
+
|
|
164
|
+
@wp.kernel(enable_backward=False)
|
|
165
|
+
def sample_mesh(
|
|
166
|
+
mesh: wp.uint64,
|
|
167
|
+
cdf: wp.array(dtype=wp.float32),
|
|
168
|
+
seed: wp.int32,
|
|
169
|
+
out_points: wp.array(dtype=wp.vec3),
|
|
170
|
+
):
|
|
171
|
+
tid = wp.tid()
|
|
172
|
+
|
|
173
|
+
rng = wp.rand_init(seed, tid)
|
|
174
|
+
|
|
175
|
+
# Sample the triangle index using the CDF.
|
|
176
|
+
sample = wp.randf(rng)
|
|
177
|
+
tri = wp.lower_bound(cdf, sample)
|
|
178
|
+
|
|
179
|
+
# Sample the location in that triangle using random barycentric cordinates.
|
|
180
|
+
ru = wp.randf(rng)
|
|
181
|
+
rv = wp.randf(rng)
|
|
182
|
+
tri_u = 1.0 - wp.sqrt(ru)
|
|
183
|
+
tri_v = wp.sqrt(ru) * (1.0 - rv)
|
|
184
|
+
pos = wp.mesh_eval_position(mesh, tri, tri_u, tri_v)
|
|
185
|
+
|
|
186
|
+
# Store the result.
|
|
187
|
+
out_points[tid] = pos
|
|
188
|
+
|
|
189
|
+
|
|
190
|
+
class Example:
|
|
191
|
+
def __init__(self, stage_path="example_sample_mesh.usd"):
|
|
192
|
+
self.mesh = wp.Mesh(
|
|
193
|
+
points=wp.array(POINTS, dtype=wp.vec3),
|
|
194
|
+
indices=wp.array(FACE_VERTEX_INDICES, dtype=wp.int32),
|
|
195
|
+
)
|
|
196
|
+
self.tri_count = len(FACE_VERTEX_INDICES) // 3
|
|
197
|
+
|
|
198
|
+
# Compute the area of each triangle and the total area of the mesh.
|
|
199
|
+
tri_areas = wp.empty(shape=(self.tri_count,), dtype=wp.float32)
|
|
200
|
+
total_area = wp.zeros(shape=(1,), dtype=wp.float32)
|
|
201
|
+
wp.launch(
|
|
202
|
+
compute_tri_areas,
|
|
203
|
+
dim=tri_areas.shape,
|
|
204
|
+
inputs=(
|
|
205
|
+
self.mesh.points,
|
|
206
|
+
self.mesh.indices,
|
|
207
|
+
),
|
|
208
|
+
outputs=(
|
|
209
|
+
tri_areas,
|
|
210
|
+
total_area,
|
|
211
|
+
),
|
|
212
|
+
)
|
|
213
|
+
|
|
214
|
+
# Build a Cumulative Distribution Function (CDF) where the probability
|
|
215
|
+
# of sampling a given triangle is proportional to its area.
|
|
216
|
+
self.cdf = wp.empty(shape=(self.tri_count,), dtype=wp.float32)
|
|
217
|
+
wp.launch(
|
|
218
|
+
compute_probability_distribution,
|
|
219
|
+
dim=self.cdf.shape,
|
|
220
|
+
inputs=(
|
|
221
|
+
tri_areas,
|
|
222
|
+
total_area,
|
|
223
|
+
),
|
|
224
|
+
outputs=(self.cdf,),
|
|
225
|
+
)
|
|
226
|
+
wp.launch(
|
|
227
|
+
accumulate_cdf,
|
|
228
|
+
dim=(1,),
|
|
229
|
+
inputs=(self.tri_count,),
|
|
230
|
+
outputs=(self.cdf,),
|
|
231
|
+
)
|
|
232
|
+
|
|
233
|
+
# Array to store the sampled points.
|
|
234
|
+
self.points = wp.empty(shape=(100,), dtype=wp.vec3)
|
|
235
|
+
|
|
236
|
+
self.fps = 4
|
|
237
|
+
self.frame = 0
|
|
238
|
+
|
|
239
|
+
if stage_path:
|
|
240
|
+
self.renderer = wp.render.UsdRenderer(stage_path, fps=self.fps)
|
|
241
|
+
else:
|
|
242
|
+
self.renderer = None
|
|
243
|
+
|
|
244
|
+
def step(self):
|
|
245
|
+
with wp.ScopedTimer("step"):
|
|
246
|
+
# Sample new points on the mesh using the CDF and the current frame
|
|
247
|
+
# number as seed to ensure different samples each frame.
|
|
248
|
+
wp.launch(
|
|
249
|
+
sample_mesh,
|
|
250
|
+
dim=self.points.shape,
|
|
251
|
+
inputs=(
|
|
252
|
+
self.mesh.id,
|
|
253
|
+
self.cdf,
|
|
254
|
+
self.frame,
|
|
255
|
+
),
|
|
256
|
+
outputs=(self.points,),
|
|
257
|
+
)
|
|
258
|
+
|
|
259
|
+
self.frame += 1
|
|
260
|
+
|
|
261
|
+
def render(self):
|
|
262
|
+
if self.renderer is None:
|
|
263
|
+
return
|
|
264
|
+
|
|
265
|
+
with wp.ScopedTimer("render"):
|
|
266
|
+
self.renderer.begin_frame(self.frame / self.fps)
|
|
267
|
+
self.renderer.render_mesh(
|
|
268
|
+
name="mesh",
|
|
269
|
+
points=self.mesh.points.numpy(),
|
|
270
|
+
indices=self.mesh.indices.numpy(),
|
|
271
|
+
colors=(0.35, 0.55, 0.9),
|
|
272
|
+
)
|
|
273
|
+
self.renderer.render_points(name="points", points=self.points.numpy(), radius=0.05, colors=(0.8, 0.3, 0.2))
|
|
274
|
+
self.renderer.end_frame()
|
|
275
|
+
|
|
276
|
+
|
|
277
|
+
if __name__ == "__main__":
|
|
278
|
+
import argparse
|
|
279
|
+
|
|
280
|
+
parser = argparse.ArgumentParser(formatter_class=argparse.ArgumentDefaultsHelpFormatter)
|
|
281
|
+
parser.add_argument("--device", type=str, default=None, help="Override the default Warp device.")
|
|
282
|
+
parser.add_argument(
|
|
283
|
+
"--stage_path",
|
|
284
|
+
type=lambda x: None if x == "None" else str(x),
|
|
285
|
+
default="example_sample_mesh.usd",
|
|
286
|
+
help="Path to the output USD file.",
|
|
287
|
+
)
|
|
288
|
+
parser.add_argument("--num_frames", type=int, default=16, help="Total number of frames.")
|
|
289
|
+
|
|
290
|
+
args = parser.parse_known_args()[0]
|
|
291
|
+
|
|
292
|
+
with wp.ScopedDevice(args.device):
|
|
293
|
+
example = Example(stage_path=args.stage_path)
|
|
294
|
+
|
|
295
|
+
for _ in range(args.num_frames):
|
|
296
|
+
example.step()
|
|
297
|
+
example.render()
|
|
298
|
+
|
|
299
|
+
if example.renderer:
|
|
300
|
+
example.renderer.save()
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Smoothed Particle Hydrodynamics
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Torch
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Wave
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2024 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2024 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Adaptive Grid
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example APIC Fluid Simulation
|
|
@@ -109,7 +117,7 @@ def divergence_form(s: Sample, domain: Domain, u: Field, psi: Field):
|
|
|
109
117
|
def invert_volume_kernel(values: wp.array(dtype=float)):
|
|
110
118
|
i = wp.tid()
|
|
111
119
|
m = values[i]
|
|
112
|
-
values[i] = wp.
|
|
120
|
+
values[i] = wp.where(m == 0.0, 0.0, 1.0 / m)
|
|
113
121
|
|
|
114
122
|
|
|
115
123
|
@wp.kernel
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Burgers
|
|
@@ -67,7 +75,7 @@ def cell_transport_form(s: fem.Sample, domain: fem.Domain, u: fem.Field, v: fem.
|
|
|
67
75
|
def initial_condition(s: fem.Sample, domain: fem.Domain):
|
|
68
76
|
x = domain(s)[0] * 2.0
|
|
69
77
|
wave = wp.sin(x * wp.pi)
|
|
70
|
-
return wp.vec2(wp.
|
|
78
|
+
return wp.vec2(wp.where(x <= 1.0, wave, 0.0), 0.0)
|
|
71
79
|
|
|
72
80
|
|
|
73
81
|
@fem.integrand
|
|
@@ -79,7 +87,7 @@ def velocity_norm(s: fem.Sample, u: fem.Field):
|
|
|
79
87
|
def minmod(a: float, b: float):
|
|
80
88
|
sa = wp.sign(a)
|
|
81
89
|
sb = wp.sign(b)
|
|
82
|
-
return wp.
|
|
90
|
+
return wp.where(sa == sb, sa * wp.min(wp.abs(a), wp.abs(b)), 0.0)
|
|
83
91
|
|
|
84
92
|
|
|
85
93
|
@fem.integrand
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Convection Diffusion
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Convection Diffusion DG
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Deformed Geometry
|
|
@@ -49,7 +57,7 @@ def boundary_projector_form(
|
|
|
49
57
|
Bilinear boundary condition projector form, non-zero on radial edges
|
|
50
58
|
"""
|
|
51
59
|
nor = fem.normal(domain, s)
|
|
52
|
-
active = wp.
|
|
60
|
+
active = wp.where(nor[0] < -0.9999 or nor[1] < -0.9999, 1.0, 0.0)
|
|
53
61
|
return active * u(s) * v(s)
|
|
54
62
|
|
|
55
63
|
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Diffusion
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Diffusion 3D
|
|
@@ -1,9 +1,17 @@
|
|
|
1
|
-
# Copyright (c) 2022 NVIDIA CORPORATION.
|
|
2
|
-
#
|
|
3
|
-
#
|
|
4
|
-
#
|
|
5
|
-
#
|
|
6
|
-
#
|
|
1
|
+
# SPDX-FileCopyrightText: Copyright (c) 2022 NVIDIA CORPORATION & AFFILIATES. All rights reserved.
|
|
2
|
+
# SPDX-License-Identifier: Apache-2.0
|
|
3
|
+
#
|
|
4
|
+
# Licensed under the Apache License, Version 2.0 (the "License");
|
|
5
|
+
# you may not use this file except in compliance with the License.
|
|
6
|
+
# You may obtain a copy of the License at
|
|
7
|
+
#
|
|
8
|
+
# http://www.apache.org/licenses/LICENSE-2.0
|
|
9
|
+
#
|
|
10
|
+
# Unless required by applicable law or agreed to in writing, software
|
|
11
|
+
# distributed under the License is distributed on an "AS IS" BASIS,
|
|
12
|
+
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
|
13
|
+
# See the License for the specific language governing permissions and
|
|
14
|
+
# limitations under the License.
|
|
7
15
|
|
|
8
16
|
###########################################################################
|
|
9
17
|
# Example Diffusion MGPU
|