warp-lang 1.0.1__py3-none-manylinux2014_x86_64.whl → 1.1.0__py3-none-manylinux2014_x86_64.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 +108 -97
- warp/__init__.pyi +1 -1
- warp/bin/warp-clang.so +0 -0
- warp/bin/warp.so +0 -0
- warp/build.py +115 -113
- warp/build_dll.py +383 -375
- warp/builtins.py +3425 -3354
- warp/codegen.py +2878 -2792
- warp/config.py +40 -36
- warp/constants.py +45 -45
- warp/context.py +5194 -5102
- warp/dlpack.py +442 -442
- warp/examples/__init__.py +16 -16
- warp/examples/assets/bear.usd +0 -0
- warp/examples/assets/bunny.usd +0 -0
- warp/examples/assets/cartpole.urdf +110 -110
- warp/examples/assets/crazyflie.usd +0 -0
- warp/examples/assets/cube.usd +0 -0
- warp/examples/assets/nv_ant.xml +92 -92
- warp/examples/assets/nv_humanoid.xml +183 -183
- warp/examples/assets/quadruped.urdf +267 -267
- warp/examples/assets/rocks.nvdb +0 -0
- warp/examples/assets/rocks.usd +0 -0
- warp/examples/assets/sphere.usd +0 -0
- warp/examples/benchmarks/benchmark_api.py +383 -383
- warp/examples/benchmarks/benchmark_cloth.py +278 -279
- warp/examples/benchmarks/benchmark_cloth_cupy.py +88 -88
- warp/examples/benchmarks/benchmark_cloth_jax.py +97 -100
- warp/examples/benchmarks/benchmark_cloth_numba.py +146 -142
- warp/examples/benchmarks/benchmark_cloth_numpy.py +77 -77
- warp/examples/benchmarks/benchmark_cloth_pytorch.py +86 -86
- warp/examples/benchmarks/benchmark_cloth_taichi.py +112 -112
- warp/examples/benchmarks/benchmark_cloth_warp.py +146 -146
- warp/examples/benchmarks/benchmark_launches.py +295 -295
- warp/examples/browse.py +29 -28
- warp/examples/core/example_dem.py +234 -221
- warp/examples/core/example_fluid.py +293 -267
- warp/examples/core/example_graph_capture.py +144 -129
- warp/examples/core/example_marching_cubes.py +188 -176
- warp/examples/core/example_mesh.py +174 -154
- warp/examples/core/example_mesh_intersect.py +205 -193
- warp/examples/core/example_nvdb.py +176 -169
- warp/examples/core/example_raycast.py +105 -89
- warp/examples/core/example_raymarch.py +199 -178
- warp/examples/core/example_render_opengl.py +185 -141
- warp/examples/core/example_sph.py +405 -389
- warp/examples/core/example_torch.py +222 -181
- warp/examples/core/example_wave.py +263 -249
- warp/examples/fem/bsr_utils.py +378 -380
- warp/examples/fem/example_apic_fluid.py +407 -391
- warp/examples/fem/example_convection_diffusion.py +182 -168
- warp/examples/fem/example_convection_diffusion_dg.py +219 -209
- warp/examples/fem/example_convection_diffusion_dg0.py +204 -194
- warp/examples/fem/example_deformed_geometry.py +177 -159
- warp/examples/fem/example_diffusion.py +201 -173
- warp/examples/fem/example_diffusion_3d.py +177 -152
- warp/examples/fem/example_diffusion_mgpu.py +221 -214
- warp/examples/fem/example_mixed_elasticity.py +244 -222
- warp/examples/fem/example_navier_stokes.py +259 -243
- warp/examples/fem/example_stokes.py +220 -192
- warp/examples/fem/example_stokes_transfer.py +265 -249
- warp/examples/fem/mesh_utils.py +133 -109
- warp/examples/fem/plot_utils.py +292 -287
- warp/examples/optim/example_bounce.py +260 -248
- warp/examples/optim/example_cloth_throw.py +222 -210
- warp/examples/optim/example_diffray.py +566 -535
- warp/examples/optim/example_drone.py +864 -835
- warp/examples/optim/example_inverse_kinematics.py +176 -169
- warp/examples/optim/example_inverse_kinematics_torch.py +185 -170
- warp/examples/optim/example_spring_cage.py +239 -234
- warp/examples/optim/example_trajectory.py +223 -201
- warp/examples/optim/example_walker.py +306 -292
- warp/examples/sim/example_cartpole.py +139 -128
- warp/examples/sim/example_cloth.py +196 -184
- warp/examples/sim/example_granular.py +124 -113
- warp/examples/sim/example_granular_collision_sdf.py +197 -185
- warp/examples/sim/example_jacobian_ik.py +236 -213
- warp/examples/sim/example_particle_chain.py +118 -106
- warp/examples/sim/example_quadruped.py +193 -179
- warp/examples/sim/example_rigid_chain.py +197 -189
- warp/examples/sim/example_rigid_contact.py +189 -176
- warp/examples/sim/example_rigid_force.py +127 -126
- warp/examples/sim/example_rigid_gyroscopic.py +109 -97
- warp/examples/sim/example_rigid_soft_contact.py +134 -124
- warp/examples/sim/example_soft_body.py +190 -178
- warp/fabric.py +337 -335
- warp/fem/__init__.py +60 -27
- warp/fem/cache.py +401 -388
- warp/fem/dirichlet.py +178 -179
- warp/fem/domain.py +262 -263
- warp/fem/field/__init__.py +100 -101
- warp/fem/field/field.py +148 -149
- warp/fem/field/nodal_field.py +298 -299
- warp/fem/field/restriction.py +22 -21
- warp/fem/field/test.py +180 -181
- warp/fem/field/trial.py +183 -183
- warp/fem/geometry/__init__.py +15 -19
- warp/fem/geometry/closest_point.py +69 -70
- warp/fem/geometry/deformed_geometry.py +270 -271
- warp/fem/geometry/element.py +744 -744
- warp/fem/geometry/geometry.py +184 -186
- warp/fem/geometry/grid_2d.py +380 -373
- warp/fem/geometry/grid_3d.py +441 -435
- warp/fem/geometry/hexmesh.py +953 -953
- warp/fem/geometry/partition.py +374 -376
- warp/fem/geometry/quadmesh_2d.py +532 -532
- warp/fem/geometry/tetmesh.py +840 -840
- warp/fem/geometry/trimesh_2d.py +577 -577
- warp/fem/integrate.py +1630 -1615
- warp/fem/operator.py +190 -191
- warp/fem/polynomial.py +214 -213
- warp/fem/quadrature/__init__.py +2 -2
- warp/fem/quadrature/pic_quadrature.py +243 -245
- warp/fem/quadrature/quadrature.py +295 -294
- warp/fem/space/__init__.py +294 -292
- warp/fem/space/basis_space.py +488 -489
- warp/fem/space/collocated_function_space.py +100 -105
- warp/fem/space/dof_mapper.py +236 -236
- warp/fem/space/function_space.py +148 -145
- warp/fem/space/grid_2d_function_space.py +267 -267
- warp/fem/space/grid_3d_function_space.py +305 -306
- warp/fem/space/hexmesh_function_space.py +350 -352
- warp/fem/space/partition.py +350 -350
- warp/fem/space/quadmesh_2d_function_space.py +368 -369
- warp/fem/space/restriction.py +158 -160
- warp/fem/space/shape/__init__.py +13 -15
- warp/fem/space/shape/cube_shape_function.py +738 -738
- warp/fem/space/shape/shape_function.py +102 -103
- warp/fem/space/shape/square_shape_function.py +611 -611
- warp/fem/space/shape/tet_shape_function.py +565 -567
- warp/fem/space/shape/triangle_shape_function.py +429 -429
- warp/fem/space/tetmesh_function_space.py +294 -292
- warp/fem/space/topology.py +297 -295
- warp/fem/space/trimesh_2d_function_space.py +223 -221
- warp/fem/types.py +77 -77
- warp/fem/utils.py +495 -495
- warp/jax.py +166 -141
- warp/jax_experimental.py +341 -339
- warp/native/array.h +1072 -1025
- warp/native/builtin.h +1560 -1560
- warp/native/bvh.cpp +398 -398
- warp/native/bvh.cu +525 -525
- warp/native/bvh.h +429 -429
- warp/native/clang/clang.cpp +495 -464
- warp/native/crt.cpp +31 -31
- warp/native/crt.h +334 -334
- warp/native/cuda_crt.h +1049 -1049
- warp/native/cuda_util.cpp +549 -540
- warp/native/cuda_util.h +288 -203
- warp/native/cutlass_gemm.cpp +34 -34
- warp/native/cutlass_gemm.cu +372 -372
- warp/native/error.cpp +66 -66
- warp/native/error.h +27 -27
- warp/native/fabric.h +228 -228
- warp/native/hashgrid.cpp +301 -278
- warp/native/hashgrid.cu +78 -77
- warp/native/hashgrid.h +227 -227
- warp/native/initializer_array.h +32 -32
- warp/native/intersect.h +1204 -1204
- warp/native/intersect_adj.h +365 -365
- warp/native/intersect_tri.h +322 -322
- warp/native/marching.cpp +2 -2
- warp/native/marching.cu +497 -497
- warp/native/marching.h +2 -2
- warp/native/mat.h +1498 -1498
- warp/native/matnn.h +333 -333
- warp/native/mesh.cpp +203 -203
- warp/native/mesh.cu +293 -293
- warp/native/mesh.h +1887 -1887
- warp/native/nanovdb/NanoVDB.h +4782 -4782
- warp/native/nanovdb/PNanoVDB.h +2553 -2553
- warp/native/nanovdb/PNanoVDBWrite.h +294 -294
- warp/native/noise.h +850 -850
- warp/native/quat.h +1084 -1084
- warp/native/rand.h +299 -299
- warp/native/range.h +108 -108
- warp/native/reduce.cpp +156 -156
- warp/native/reduce.cu +348 -348
- warp/native/runlength_encode.cpp +61 -61
- warp/native/runlength_encode.cu +46 -46
- warp/native/scan.cpp +30 -30
- warp/native/scan.cu +36 -36
- warp/native/scan.h +7 -7
- warp/native/solid_angle.h +442 -442
- warp/native/sort.cpp +94 -94
- warp/native/sort.cu +97 -97
- warp/native/sort.h +14 -14
- warp/native/sparse.cpp +337 -337
- warp/native/sparse.cu +544 -544
- warp/native/spatial.h +630 -630
- warp/native/svd.h +562 -562
- warp/native/temp_buffer.h +30 -30
- warp/native/vec.h +1132 -1132
- warp/native/volume.cpp +297 -297
- warp/native/volume.cu +32 -32
- warp/native/volume.h +538 -538
- warp/native/volume_builder.cu +425 -425
- warp/native/volume_builder.h +19 -19
- warp/native/warp.cpp +1057 -1052
- warp/native/warp.cu +2943 -2828
- warp/native/warp.h +313 -305
- warp/optim/__init__.py +9 -9
- warp/optim/adam.py +120 -120
- warp/optim/linear.py +1104 -939
- warp/optim/sgd.py +104 -92
- warp/render/__init__.py +10 -10
- warp/render/render_opengl.py +3217 -3204
- warp/render/render_usd.py +768 -749
- warp/render/utils.py +152 -150
- warp/sim/__init__.py +52 -59
- warp/sim/articulation.py +685 -685
- warp/sim/collide.py +1594 -1590
- warp/sim/import_mjcf.py +489 -481
- warp/sim/import_snu.py +220 -221
- warp/sim/import_urdf.py +536 -516
- warp/sim/import_usd.py +887 -881
- warp/sim/inertia.py +316 -317
- warp/sim/integrator.py +234 -233
- warp/sim/integrator_euler.py +1956 -1956
- warp/sim/integrator_featherstone.py +1910 -1991
- warp/sim/integrator_xpbd.py +3294 -3312
- warp/sim/model.py +4473 -4314
- warp/sim/particles.py +113 -112
- warp/sim/render.py +417 -403
- warp/sim/utils.py +413 -410
- warp/sparse.py +1227 -1227
- warp/stubs.py +2109 -2469
- warp/tape.py +1162 -225
- warp/tests/__init__.py +1 -1
- warp/tests/__main__.py +4 -4
- warp/tests/assets/torus.usda +105 -105
- warp/tests/aux_test_class_kernel.py +26 -26
- warp/tests/aux_test_compile_consts_dummy.py +10 -10
- warp/tests/aux_test_conditional_unequal_types_kernels.py +21 -21
- warp/tests/aux_test_dependent.py +22 -22
- warp/tests/aux_test_grad_customs.py +23 -23
- warp/tests/aux_test_reference.py +11 -11
- warp/tests/aux_test_reference_reference.py +10 -10
- warp/tests/aux_test_square.py +17 -17
- warp/tests/aux_test_unresolved_func.py +14 -14
- warp/tests/aux_test_unresolved_symbol.py +14 -14
- warp/tests/disabled_kinematics.py +239 -239
- warp/tests/run_coverage_serial.py +31 -31
- warp/tests/test_adam.py +157 -157
- warp/tests/test_arithmetic.py +1124 -1124
- warp/tests/test_array.py +2417 -2326
- warp/tests/test_array_reduce.py +150 -150
- warp/tests/test_async.py +668 -656
- warp/tests/test_atomic.py +141 -141
- warp/tests/test_bool.py +204 -149
- warp/tests/test_builtins_resolution.py +1292 -1292
- warp/tests/test_bvh.py +164 -171
- warp/tests/test_closest_point_edge_edge.py +228 -228
- warp/tests/test_codegen.py +566 -553
- warp/tests/test_compile_consts.py +97 -101
- warp/tests/test_conditional.py +246 -246
- warp/tests/test_copy.py +232 -215
- warp/tests/test_ctypes.py +632 -632
- warp/tests/test_dense.py +67 -67
- warp/tests/test_devices.py +91 -98
- warp/tests/test_dlpack.py +530 -529
- warp/tests/test_examples.py +400 -378
- warp/tests/test_fabricarray.py +955 -955
- warp/tests/test_fast_math.py +62 -54
- warp/tests/test_fem.py +1277 -1278
- warp/tests/test_fp16.py +130 -130
- warp/tests/test_func.py +338 -337
- warp/tests/test_generics.py +571 -571
- warp/tests/test_grad.py +746 -640
- warp/tests/test_grad_customs.py +333 -336
- warp/tests/test_hash_grid.py +210 -164
- warp/tests/test_import.py +39 -39
- warp/tests/test_indexedarray.py +1134 -1134
- warp/tests/test_intersect.py +67 -67
- warp/tests/test_jax.py +307 -307
- warp/tests/test_large.py +167 -164
- warp/tests/test_launch.py +354 -354
- warp/tests/test_lerp.py +261 -261
- warp/tests/test_linear_solvers.py +191 -171
- warp/tests/test_lvalue.py +421 -493
- warp/tests/test_marching_cubes.py +65 -65
- warp/tests/test_mat.py +1801 -1827
- warp/tests/test_mat_lite.py +115 -115
- warp/tests/test_mat_scalar_ops.py +2907 -2889
- warp/tests/test_math.py +126 -193
- warp/tests/test_matmul.py +500 -499
- warp/tests/test_matmul_lite.py +410 -410
- warp/tests/test_mempool.py +188 -190
- warp/tests/test_mesh.py +284 -324
- warp/tests/test_mesh_query_aabb.py +228 -241
- warp/tests/test_mesh_query_point.py +692 -702
- warp/tests/test_mesh_query_ray.py +292 -303
- warp/tests/test_mlp.py +276 -276
- warp/tests/test_model.py +110 -110
- warp/tests/test_modules_lite.py +39 -39
- warp/tests/test_multigpu.py +163 -163
- warp/tests/test_noise.py +248 -248
- warp/tests/test_operators.py +250 -250
- warp/tests/test_options.py +123 -125
- warp/tests/test_peer.py +133 -137
- warp/tests/test_pinned.py +78 -78
- warp/tests/test_print.py +54 -54
- warp/tests/test_quat.py +2086 -2086
- warp/tests/test_rand.py +288 -288
- warp/tests/test_reload.py +217 -217
- warp/tests/test_rounding.py +179 -179
- warp/tests/test_runlength_encode.py +190 -190
- warp/tests/test_sim_grad.py +243 -0
- warp/tests/test_sim_kinematics.py +91 -97
- warp/tests/test_smoothstep.py +168 -168
- warp/tests/test_snippet.py +305 -266
- warp/tests/test_sparse.py +468 -460
- warp/tests/test_spatial.py +2148 -2148
- warp/tests/test_streams.py +486 -473
- warp/tests/test_struct.py +710 -675
- warp/tests/test_tape.py +173 -148
- warp/tests/test_torch.py +743 -743
- warp/tests/test_transient_module.py +87 -87
- warp/tests/test_types.py +556 -659
- warp/tests/test_utils.py +490 -499
- warp/tests/test_vec.py +1264 -1268
- warp/tests/test_vec_lite.py +73 -73
- warp/tests/test_vec_scalar_ops.py +2099 -2099
- warp/tests/test_verify_fp.py +94 -94
- warp/tests/test_volume.py +737 -736
- warp/tests/test_volume_write.py +255 -265
- warp/tests/unittest_serial.py +37 -37
- warp/tests/unittest_suites.py +363 -359
- warp/tests/unittest_utils.py +603 -578
- warp/tests/unused_test_misc.py +71 -71
- warp/tests/walkthrough_debug.py +85 -85
- warp/thirdparty/appdirs.py +598 -598
- warp/thirdparty/dlpack.py +143 -143
- warp/thirdparty/unittest_parallel.py +566 -561
- warp/torch.py +321 -295
- warp/types.py +4504 -4450
- warp/utils.py +1008 -821
- {warp_lang-1.0.1.dist-info → warp_lang-1.1.0.dist-info}/LICENSE.md +126 -126
- {warp_lang-1.0.1.dist-info → warp_lang-1.1.0.dist-info}/METADATA +338 -400
- warp_lang-1.1.0.dist-info/RECORD +352 -0
- warp/examples/assets/cube.usda +0 -42
- warp/examples/assets/sphere.usda +0 -56
- warp/examples/assets/torus.usda +0 -105
- warp_lang-1.0.1.dist-info/RECORD +0 -352
- {warp_lang-1.0.1.dist-info → warp_lang-1.1.0.dist-info}/WHEEL +0 -0
- {warp_lang-1.0.1.dist-info → warp_lang-1.1.0.dist-info}/top_level.txt +0 -0
warp/fem/cache.py
CHANGED
|
@@ -1,388 +1,401 @@
|
|
|
1
|
-
|
|
2
|
-
|
|
3
|
-
import
|
|
4
|
-
import
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
13
|
-
|
|
14
|
-
|
|
15
|
-
|
|
16
|
-
|
|
17
|
-
|
|
18
|
-
|
|
19
|
-
|
|
20
|
-
|
|
21
|
-
|
|
22
|
-
|
|
23
|
-
key
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
|
|
27
|
-
|
|
28
|
-
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
33
|
-
|
|
34
|
-
|
|
35
|
-
|
|
36
|
-
|
|
37
|
-
|
|
38
|
-
|
|
39
|
-
|
|
40
|
-
|
|
41
|
-
|
|
42
|
-
|
|
43
|
-
|
|
44
|
-
|
|
45
|
-
|
|
46
|
-
|
|
47
|
-
|
|
48
|
-
|
|
49
|
-
kernel_options
|
|
50
|
-
|
|
51
|
-
|
|
52
|
-
|
|
53
|
-
|
|
54
|
-
|
|
55
|
-
#
|
|
56
|
-
|
|
57
|
-
module_name =
|
|
58
|
-
|
|
59
|
-
module
|
|
60
|
-
module.options.
|
|
61
|
-
|
|
62
|
-
|
|
63
|
-
|
|
64
|
-
|
|
65
|
-
|
|
66
|
-
|
|
67
|
-
|
|
68
|
-
|
|
69
|
-
|
|
70
|
-
|
|
71
|
-
|
|
72
|
-
|
|
73
|
-
|
|
74
|
-
|
|
75
|
-
|
|
76
|
-
|
|
77
|
-
|
|
78
|
-
|
|
79
|
-
|
|
80
|
-
|
|
81
|
-
|
|
82
|
-
|
|
83
|
-
|
|
84
|
-
|
|
85
|
-
|
|
86
|
-
|
|
87
|
-
|
|
88
|
-
|
|
89
|
-
|
|
90
|
-
|
|
91
|
-
|
|
92
|
-
|
|
93
|
-
|
|
94
|
-
|
|
95
|
-
|
|
96
|
-
|
|
97
|
-
|
|
98
|
-
|
|
99
|
-
|
|
100
|
-
|
|
101
|
-
|
|
102
|
-
|
|
103
|
-
|
|
104
|
-
|
|
105
|
-
|
|
106
|
-
|
|
107
|
-
|
|
108
|
-
|
|
109
|
-
|
|
110
|
-
|
|
111
|
-
|
|
112
|
-
|
|
113
|
-
|
|
114
|
-
|
|
115
|
-
|
|
116
|
-
|
|
117
|
-
|
|
118
|
-
|
|
119
|
-
|
|
120
|
-
|
|
121
|
-
|
|
122
|
-
|
|
123
|
-
|
|
124
|
-
|
|
125
|
-
|
|
126
|
-
|
|
127
|
-
|
|
128
|
-
|
|
129
|
-
|
|
130
|
-
|
|
131
|
-
|
|
132
|
-
|
|
133
|
-
|
|
134
|
-
|
|
135
|
-
|
|
136
|
-
|
|
137
|
-
|
|
138
|
-
|
|
139
|
-
|
|
140
|
-
|
|
141
|
-
|
|
142
|
-
|
|
143
|
-
|
|
144
|
-
|
|
145
|
-
|
|
146
|
-
|
|
147
|
-
|
|
148
|
-
|
|
149
|
-
|
|
150
|
-
|
|
151
|
-
|
|
152
|
-
|
|
153
|
-
|
|
154
|
-
|
|
155
|
-
|
|
156
|
-
|
|
157
|
-
|
|
158
|
-
|
|
159
|
-
|
|
160
|
-
|
|
161
|
-
|
|
162
|
-
|
|
163
|
-
|
|
164
|
-
|
|
165
|
-
|
|
166
|
-
|
|
167
|
-
|
|
168
|
-
|
|
169
|
-
|
|
170
|
-
|
|
171
|
-
|
|
172
|
-
|
|
173
|
-
|
|
174
|
-
|
|
175
|
-
|
|
176
|
-
|
|
177
|
-
|
|
178
|
-
|
|
179
|
-
|
|
180
|
-
|
|
181
|
-
|
|
182
|
-
|
|
183
|
-
|
|
184
|
-
|
|
185
|
-
|
|
186
|
-
|
|
187
|
-
|
|
188
|
-
|
|
189
|
-
|
|
190
|
-
|
|
191
|
-
|
|
192
|
-
|
|
193
|
-
|
|
194
|
-
|
|
195
|
-
|
|
196
|
-
|
|
197
|
-
|
|
198
|
-
|
|
199
|
-
|
|
200
|
-
|
|
201
|
-
|
|
202
|
-
|
|
203
|
-
|
|
204
|
-
|
|
205
|
-
|
|
206
|
-
if
|
|
207
|
-
self.
|
|
208
|
-
|
|
209
|
-
|
|
210
|
-
|
|
211
|
-
|
|
212
|
-
|
|
213
|
-
|
|
214
|
-
|
|
215
|
-
|
|
216
|
-
|
|
217
|
-
|
|
218
|
-
|
|
219
|
-
|
|
220
|
-
|
|
221
|
-
|
|
222
|
-
|
|
223
|
-
|
|
224
|
-
|
|
225
|
-
|
|
226
|
-
|
|
227
|
-
|
|
228
|
-
|
|
229
|
-
|
|
230
|
-
|
|
231
|
-
|
|
232
|
-
|
|
233
|
-
|
|
234
|
-
|
|
235
|
-
|
|
236
|
-
|
|
237
|
-
|
|
238
|
-
|
|
239
|
-
|
|
240
|
-
|
|
241
|
-
|
|
242
|
-
|
|
243
|
-
|
|
244
|
-
|
|
245
|
-
|
|
246
|
-
|
|
247
|
-
|
|
248
|
-
|
|
249
|
-
|
|
250
|
-
|
|
251
|
-
|
|
252
|
-
|
|
253
|
-
|
|
254
|
-
|
|
255
|
-
|
|
256
|
-
|
|
257
|
-
|
|
258
|
-
|
|
259
|
-
|
|
260
|
-
|
|
261
|
-
|
|
262
|
-
|
|
263
|
-
|
|
264
|
-
|
|
265
|
-
|
|
266
|
-
|
|
267
|
-
|
|
268
|
-
|
|
269
|
-
|
|
270
|
-
|
|
271
|
-
|
|
272
|
-
|
|
273
|
-
|
|
274
|
-
|
|
275
|
-
|
|
276
|
-
|
|
277
|
-
|
|
278
|
-
|
|
279
|
-
|
|
280
|
-
size
|
|
281
|
-
|
|
282
|
-
|
|
283
|
-
|
|
284
|
-
|
|
285
|
-
|
|
286
|
-
|
|
287
|
-
|
|
288
|
-
|
|
289
|
-
|
|
290
|
-
|
|
291
|
-
|
|
292
|
-
|
|
293
|
-
|
|
294
|
-
|
|
295
|
-
|
|
296
|
-
|
|
297
|
-
|
|
298
|
-
|
|
299
|
-
|
|
300
|
-
|
|
301
|
-
self
|
|
302
|
-
|
|
303
|
-
|
|
304
|
-
|
|
305
|
-
|
|
306
|
-
|
|
307
|
-
|
|
308
|
-
|
|
309
|
-
|
|
310
|
-
|
|
311
|
-
|
|
312
|
-
|
|
313
|
-
|
|
314
|
-
|
|
315
|
-
|
|
316
|
-
|
|
317
|
-
|
|
318
|
-
|
|
319
|
-
|
|
320
|
-
|
|
321
|
-
|
|
322
|
-
|
|
323
|
-
|
|
324
|
-
|
|
325
|
-
|
|
326
|
-
|
|
327
|
-
|
|
328
|
-
|
|
329
|
-
|
|
330
|
-
|
|
331
|
-
|
|
332
|
-
|
|
333
|
-
|
|
334
|
-
|
|
335
|
-
|
|
336
|
-
|
|
337
|
-
|
|
338
|
-
|
|
339
|
-
|
|
340
|
-
|
|
341
|
-
|
|
342
|
-
|
|
343
|
-
"""
|
|
344
|
-
|
|
345
|
-
|
|
346
|
-
|
|
347
|
-
|
|
348
|
-
|
|
349
|
-
|
|
350
|
-
|
|
351
|
-
|
|
352
|
-
|
|
353
|
-
|
|
354
|
-
|
|
355
|
-
|
|
356
|
-
|
|
357
|
-
|
|
358
|
-
|
|
359
|
-
|
|
360
|
-
|
|
361
|
-
|
|
362
|
-
|
|
363
|
-
|
|
364
|
-
|
|
365
|
-
|
|
366
|
-
|
|
367
|
-
|
|
368
|
-
|
|
369
|
-
|
|
370
|
-
temporary_store:
|
|
371
|
-
|
|
372
|
-
|
|
373
|
-
|
|
374
|
-
|
|
375
|
-
|
|
376
|
-
|
|
377
|
-
|
|
378
|
-
|
|
379
|
-
|
|
380
|
-
|
|
381
|
-
|
|
382
|
-
|
|
383
|
-
|
|
384
|
-
|
|
385
|
-
|
|
386
|
-
|
|
387
|
-
|
|
388
|
-
|
|
1
|
+
import bisect
|
|
2
|
+
import re
|
|
3
|
+
from copy import copy
|
|
4
|
+
from typing import Any, Callable, Dict, Optional, Tuple, Union
|
|
5
|
+
|
|
6
|
+
import warp as wp
|
|
7
|
+
|
|
8
|
+
_kernel_cache = {}
|
|
9
|
+
_struct_cache = {}
|
|
10
|
+
_func_cache = {}
|
|
11
|
+
|
|
12
|
+
_key_re = re.compile("[^0-9a-zA-Z_]+")
|
|
13
|
+
|
|
14
|
+
|
|
15
|
+
def _make_key(obj, suffix: str, use_qualified_name):
|
|
16
|
+
base_name = f"{obj.__module__}.{obj.__qualname__}" if use_qualified_name else obj.__name__
|
|
17
|
+
return _key_re.sub("", f"{base_name}_{suffix}")
|
|
18
|
+
|
|
19
|
+
|
|
20
|
+
def get_func(func, suffix: str, use_qualified_name: bool = False):
|
|
21
|
+
key = _make_key(func, suffix, use_qualified_name)
|
|
22
|
+
|
|
23
|
+
if key not in _func_cache:
|
|
24
|
+
_func_cache[key] = wp.Function(
|
|
25
|
+
func=func,
|
|
26
|
+
key=key,
|
|
27
|
+
namespace="",
|
|
28
|
+
module=wp.get_module(
|
|
29
|
+
func.__module__,
|
|
30
|
+
),
|
|
31
|
+
)
|
|
32
|
+
|
|
33
|
+
return _func_cache[key]
|
|
34
|
+
|
|
35
|
+
|
|
36
|
+
def dynamic_func(suffix: str, use_qualified_name=False):
|
|
37
|
+
def wrap_func(func: Callable):
|
|
38
|
+
return get_func(func, suffix=suffix, use_qualified_name=use_qualified_name)
|
|
39
|
+
|
|
40
|
+
return wrap_func
|
|
41
|
+
|
|
42
|
+
|
|
43
|
+
def get_kernel(
|
|
44
|
+
func,
|
|
45
|
+
suffix: str,
|
|
46
|
+
use_qualified_name: bool = False,
|
|
47
|
+
kernel_options: Dict[str, Any] = None,
|
|
48
|
+
):
|
|
49
|
+
if kernel_options is None:
|
|
50
|
+
kernel_options = {}
|
|
51
|
+
|
|
52
|
+
key = _make_key(func, suffix, use_qualified_name)
|
|
53
|
+
|
|
54
|
+
if key not in _kernel_cache:
|
|
55
|
+
# Avoid creating too long file names -- can lead to issues on Windows
|
|
56
|
+
# We could hash the key, but prefer to keep it human-readable
|
|
57
|
+
module_name = f"{func.__module__}.dyn.{key}"
|
|
58
|
+
module_name = module_name[:128] if len(module_name) > 128 else module_name
|
|
59
|
+
module = wp.get_module(module_name)
|
|
60
|
+
module.options = copy(wp.get_module(func.__module__).options)
|
|
61
|
+
module.options.update(kernel_options)
|
|
62
|
+
_kernel_cache[key] = wp.Kernel(func=func, key=key, module=module)
|
|
63
|
+
return _kernel_cache[key]
|
|
64
|
+
|
|
65
|
+
|
|
66
|
+
def dynamic_kernel(suffix: str, use_qualified_name=False, kernel_options: Dict[str, Any] = None):
|
|
67
|
+
if kernel_options is None:
|
|
68
|
+
kernel_options = {}
|
|
69
|
+
|
|
70
|
+
def wrap_kernel(func: Callable):
|
|
71
|
+
return get_kernel(func, suffix=suffix, use_qualified_name=use_qualified_name, kernel_options=kernel_options)
|
|
72
|
+
|
|
73
|
+
return wrap_kernel
|
|
74
|
+
|
|
75
|
+
|
|
76
|
+
def get_struct(struct: type, suffix: str, use_qualified_name: bool = False):
|
|
77
|
+
key = _make_key(struct, suffix, use_qualified_name)
|
|
78
|
+
# used in codegen
|
|
79
|
+
struct.__qualname__ = key
|
|
80
|
+
|
|
81
|
+
if key not in _struct_cache:
|
|
82
|
+
module = wp.get_module(struct.__module__)
|
|
83
|
+
_struct_cache[key] = wp.codegen.Struct(
|
|
84
|
+
cls=struct,
|
|
85
|
+
key=key,
|
|
86
|
+
module=module,
|
|
87
|
+
)
|
|
88
|
+
|
|
89
|
+
return _struct_cache[key]
|
|
90
|
+
|
|
91
|
+
|
|
92
|
+
def dynamic_struct(suffix: str, use_qualified_name=False):
|
|
93
|
+
def wrap_struct(struct: type):
|
|
94
|
+
return get_struct(struct, suffix=suffix, use_qualified_name=use_qualified_name)
|
|
95
|
+
|
|
96
|
+
return wrap_struct
|
|
97
|
+
|
|
98
|
+
|
|
99
|
+
def get_integrand_function(
|
|
100
|
+
integrand: "warp.fem.operator.Integrand", # noqa: F821
|
|
101
|
+
suffix: str,
|
|
102
|
+
func=None,
|
|
103
|
+
annotations=None,
|
|
104
|
+
code_transformers=None,
|
|
105
|
+
):
|
|
106
|
+
if code_transformers is None:
|
|
107
|
+
code_transformers = []
|
|
108
|
+
|
|
109
|
+
key = _make_key(integrand.func, suffix, use_qualified_name=True)
|
|
110
|
+
|
|
111
|
+
if key not in _func_cache:
|
|
112
|
+
_func_cache[key] = wp.Function(
|
|
113
|
+
func=integrand.func if func is None else func,
|
|
114
|
+
key=key,
|
|
115
|
+
namespace="",
|
|
116
|
+
module=integrand.module,
|
|
117
|
+
overloaded_annotations=annotations,
|
|
118
|
+
code_transformers=code_transformers,
|
|
119
|
+
)
|
|
120
|
+
|
|
121
|
+
return _func_cache[key]
|
|
122
|
+
|
|
123
|
+
|
|
124
|
+
def get_integrand_kernel(
|
|
125
|
+
integrand: "warp.fem.operator.Integrand", # noqa: F821
|
|
126
|
+
suffix: str,
|
|
127
|
+
kernel_fn: Optional[Callable] = None,
|
|
128
|
+
kernel_options: Dict[str, Any] = None,
|
|
129
|
+
code_transformers=None,
|
|
130
|
+
):
|
|
131
|
+
if kernel_options is None:
|
|
132
|
+
kernel_options = {}
|
|
133
|
+
|
|
134
|
+
if code_transformers is None:
|
|
135
|
+
code_transformers = []
|
|
136
|
+
|
|
137
|
+
key = _make_key(integrand.func, suffix, use_qualified_name=True)
|
|
138
|
+
|
|
139
|
+
if key not in _kernel_cache:
|
|
140
|
+
if kernel_fn is None:
|
|
141
|
+
return None
|
|
142
|
+
|
|
143
|
+
module = wp.get_module(f"{integrand.module.name}.{integrand.name}")
|
|
144
|
+
module.options = copy(integrand.module.options)
|
|
145
|
+
module.options.update(kernel_options)
|
|
146
|
+
|
|
147
|
+
_kernel_cache[key] = wp.Kernel(func=kernel_fn, key=key, module=module, code_transformers=code_transformers)
|
|
148
|
+
return _kernel_cache[key]
|
|
149
|
+
|
|
150
|
+
|
|
151
|
+
def cached_arg_value(func: Callable):
|
|
152
|
+
"""Decorator to be applied to member methods assembling Arg structs, so that the result gets
|
|
153
|
+
automatically cached for the lifetime of the parent object
|
|
154
|
+
"""
|
|
155
|
+
|
|
156
|
+
cache_attr = f"_{func.__name__}_cache"
|
|
157
|
+
|
|
158
|
+
def get_arg(obj, device):
|
|
159
|
+
if not hasattr(obj, cache_attr):
|
|
160
|
+
setattr(obj, cache_attr, {})
|
|
161
|
+
|
|
162
|
+
cache = getattr(obj, cache_attr, {})
|
|
163
|
+
|
|
164
|
+
device = wp.get_device(device)
|
|
165
|
+
if device.ordinal not in cache:
|
|
166
|
+
cache[device.ordinal] = func(obj, device)
|
|
167
|
+
|
|
168
|
+
return cache[device.ordinal]
|
|
169
|
+
|
|
170
|
+
return get_arg
|
|
171
|
+
|
|
172
|
+
|
|
173
|
+
_cached_vec_types = {}
|
|
174
|
+
_cached_mat_types = {}
|
|
175
|
+
|
|
176
|
+
|
|
177
|
+
def cached_vec_type(length, dtype):
|
|
178
|
+
key = (length, dtype)
|
|
179
|
+
if key not in _cached_vec_types:
|
|
180
|
+
_cached_vec_types[key] = wp.vec(length=length, dtype=dtype)
|
|
181
|
+
|
|
182
|
+
return _cached_vec_types[key]
|
|
183
|
+
|
|
184
|
+
|
|
185
|
+
def cached_mat_type(shape, dtype):
|
|
186
|
+
key = (*shape, dtype)
|
|
187
|
+
if key not in _cached_mat_types:
|
|
188
|
+
_cached_mat_types[key] = wp.mat(shape=shape, dtype=dtype)
|
|
189
|
+
|
|
190
|
+
return _cached_mat_types[key]
|
|
191
|
+
|
|
192
|
+
|
|
193
|
+
class Temporary:
|
|
194
|
+
"""Handle over a temporary array from a :class:`TemporaryStore`.
|
|
195
|
+
|
|
196
|
+
The array will be automatically returned to the temporary pool for reuse upon destruction of this object, unless
|
|
197
|
+
the temporary is explicitly detached from the pool using :meth:`detach`.
|
|
198
|
+
The temporary may also be explicitly returned to the pool before destruction using :meth:`release`.
|
|
199
|
+
"""
|
|
200
|
+
|
|
201
|
+
def __init__(self, array: wp.array, pool: Optional["TemporaryStore.Pool"] = None, shape=None, dtype=None):
|
|
202
|
+
self._raw_array = array
|
|
203
|
+
self._array_view = array
|
|
204
|
+
self._pool = pool
|
|
205
|
+
|
|
206
|
+
if shape is not None or dtype is not None:
|
|
207
|
+
self._view_as(shape=shape, dtype=dtype)
|
|
208
|
+
|
|
209
|
+
def detach(self) -> wp.array:
|
|
210
|
+
"""Detaches the temporary so it is never returned to the pool"""
|
|
211
|
+
if self._pool is not None:
|
|
212
|
+
self._pool.detach(self._raw_array)
|
|
213
|
+
|
|
214
|
+
self._pool = None
|
|
215
|
+
return self._array_view
|
|
216
|
+
|
|
217
|
+
def release(self):
|
|
218
|
+
"""Returns the temporary array to the pool"""
|
|
219
|
+
if self._pool is not None:
|
|
220
|
+
self._pool.redeem(self._raw_array)
|
|
221
|
+
|
|
222
|
+
self._pool = None
|
|
223
|
+
|
|
224
|
+
@property
|
|
225
|
+
def array(self) -> wp.array:
|
|
226
|
+
"""View of the array with desired shape and data type."""
|
|
227
|
+
return self._array_view
|
|
228
|
+
|
|
229
|
+
def _view_as(self, shape, dtype) -> "Temporary":
|
|
230
|
+
def _view_reshaped_truncated(array):
|
|
231
|
+
return wp.types.array(
|
|
232
|
+
ptr=array.ptr,
|
|
233
|
+
dtype=dtype,
|
|
234
|
+
shape=shape,
|
|
235
|
+
device=array.device,
|
|
236
|
+
pinned=array.pinned,
|
|
237
|
+
capacity=array.capacity,
|
|
238
|
+
copy=False,
|
|
239
|
+
grad=None if array.grad is None else _view_reshaped_truncated(array.grad),
|
|
240
|
+
)
|
|
241
|
+
|
|
242
|
+
self._array_view = _view_reshaped_truncated(self._raw_array)
|
|
243
|
+
return self
|
|
244
|
+
|
|
245
|
+
def __del__(self):
|
|
246
|
+
self.release()
|
|
247
|
+
|
|
248
|
+
|
|
249
|
+
class TemporaryStore:
|
|
250
|
+
"""
|
|
251
|
+
Shared pool of temporary arrays that will be persisted and reused across invocations of ``warp.fem`` functions.
|
|
252
|
+
|
|
253
|
+
A :class:`TemporaryStore` instance may either be passed explicitly to ``warp.fem`` functions that accept such an argument, for instance :func:`.integrate.integrate`,
|
|
254
|
+
or can be set globally as the default store using :func:`set_default_temporary_store`.
|
|
255
|
+
|
|
256
|
+
By default, there is no default temporary store, so that temporary allocations are not persisted.
|
|
257
|
+
"""
|
|
258
|
+
|
|
259
|
+
_default_store: "TemporaryStore" = None
|
|
260
|
+
|
|
261
|
+
class Pool:
|
|
262
|
+
def __init__(self, dtype, device, pinned: bool):
|
|
263
|
+
self.dtype = dtype
|
|
264
|
+
self.device = device
|
|
265
|
+
self.pinned = pinned
|
|
266
|
+
|
|
267
|
+
self._pool = [] # Currently available arrays for borrowing, ordered by size
|
|
268
|
+
self._pool_sizes = [] # Sizes of available arrays for borrowing, ascending
|
|
269
|
+
self._allocs = {} # All allocated arrays, including borrowed ones
|
|
270
|
+
|
|
271
|
+
def borrow(self, shape, dtype, requires_grad: bool):
|
|
272
|
+
size = 1
|
|
273
|
+
if isinstance(shape, int):
|
|
274
|
+
shape = (shape,)
|
|
275
|
+
for d in shape:
|
|
276
|
+
size *= d
|
|
277
|
+
|
|
278
|
+
index = bisect.bisect_left(
|
|
279
|
+
a=self._pool_sizes,
|
|
280
|
+
x=size,
|
|
281
|
+
)
|
|
282
|
+
if index < len(self._pool):
|
|
283
|
+
# Big enough array found, remove from pool
|
|
284
|
+
array = self._pool.pop(index)
|
|
285
|
+
self._pool_sizes.pop(index)
|
|
286
|
+
if requires_grad and array.grad is None:
|
|
287
|
+
array.requires_grad = True
|
|
288
|
+
return Temporary(pool=self, array=array, shape=shape, dtype=dtype)
|
|
289
|
+
|
|
290
|
+
# No big enough array found, allocate new one
|
|
291
|
+
if len(self._pool) > 0:
|
|
292
|
+
grow_factor = 1.5
|
|
293
|
+
size = max(int(self._pool_sizes[-1] * grow_factor), size)
|
|
294
|
+
|
|
295
|
+
array = wp.empty(
|
|
296
|
+
shape=(size,), dtype=self.dtype, pinned=self.pinned, device=self.device, requires_grad=requires_grad
|
|
297
|
+
)
|
|
298
|
+
self._allocs[array.ptr] = array
|
|
299
|
+
return Temporary(pool=self, array=array, shape=shape, dtype=dtype)
|
|
300
|
+
|
|
301
|
+
def redeem(self, array):
|
|
302
|
+
# Insert back array into available pool
|
|
303
|
+
index = bisect.bisect_left(
|
|
304
|
+
a=self._pool_sizes,
|
|
305
|
+
x=array.size,
|
|
306
|
+
)
|
|
307
|
+
self._pool.insert(index, array)
|
|
308
|
+
self._pool_sizes.insert(index, array.size)
|
|
309
|
+
|
|
310
|
+
def detach(self, array):
|
|
311
|
+
del self._allocs[array.ptr]
|
|
312
|
+
|
|
313
|
+
def __init__(self):
|
|
314
|
+
self.clear()
|
|
315
|
+
|
|
316
|
+
def clear(self):
|
|
317
|
+
self._temporaries = {}
|
|
318
|
+
|
|
319
|
+
def borrow(self, shape, dtype, pinned: bool = False, device=None, requires_grad: bool = False) -> Temporary:
|
|
320
|
+
dtype = wp.types.type_to_warp(dtype)
|
|
321
|
+
device = wp.get_device(device)
|
|
322
|
+
|
|
323
|
+
type_length = wp.types.type_length(dtype)
|
|
324
|
+
key = (dtype._type_, type_length, pinned, device.ordinal)
|
|
325
|
+
|
|
326
|
+
pool = self._temporaries.get(key, None)
|
|
327
|
+
if pool is None:
|
|
328
|
+
value_type = (
|
|
329
|
+
cached_vec_type(length=type_length, dtype=wp.types.type_scalar_type(dtype))
|
|
330
|
+
if type_length > 1
|
|
331
|
+
else dtype
|
|
332
|
+
)
|
|
333
|
+
pool = TemporaryStore.Pool(value_type, device, pinned=pinned)
|
|
334
|
+
self._temporaries[key] = pool
|
|
335
|
+
|
|
336
|
+
return pool.borrow(dtype=dtype, shape=shape, requires_grad=requires_grad)
|
|
337
|
+
|
|
338
|
+
|
|
339
|
+
def set_default_temporary_store(temporary_store: Optional[TemporaryStore]):
|
|
340
|
+
"""Globally sets the default :class:`TemporaryStore` instance to use for temporary allocations in ``warp.fem`` functions.
|
|
341
|
+
|
|
342
|
+
If the default temporary store is set to ``None``, temporary allocations are not persisted unless a :class:`TemporaryStore` is provided at a per-function granularity.
|
|
343
|
+
"""
|
|
344
|
+
|
|
345
|
+
TemporaryStore._default_store = temporary_store
|
|
346
|
+
|
|
347
|
+
|
|
348
|
+
def borrow_temporary(
|
|
349
|
+
temporary_store: Optional[TemporaryStore],
|
|
350
|
+
shape: Union[int, Tuple[int]],
|
|
351
|
+
dtype: type,
|
|
352
|
+
pinned: bool = False,
|
|
353
|
+
requires_grad: bool = False,
|
|
354
|
+
device=None,
|
|
355
|
+
) -> Temporary:
|
|
356
|
+
"""
|
|
357
|
+
Borrows and returns a temporary array with specified attributes from a shared pool.
|
|
358
|
+
|
|
359
|
+
If an array with sufficient capacity and matching desired attributes is already available in the pool, it will be returned.
|
|
360
|
+
Otherwise, a new allocation will be performed.
|
|
361
|
+
|
|
362
|
+
Args:
|
|
363
|
+
temporary_store: the shared pool to borrow the temporary from. If `temporary_store` is ``None``, the global default temporary store, if set, will be used.
|
|
364
|
+
shape: desired dimensions for the temporary array
|
|
365
|
+
dtype: desired data type for the temporary array
|
|
366
|
+
pinned: whether a pinned allocation is desired
|
|
367
|
+
device: device on which the memory should be allocated; if ``None``, the current device will be used.
|
|
368
|
+
"""
|
|
369
|
+
|
|
370
|
+
if temporary_store is None:
|
|
371
|
+
temporary_store = TemporaryStore._default_store
|
|
372
|
+
|
|
373
|
+
if temporary_store is None:
|
|
374
|
+
return Temporary(
|
|
375
|
+
array=wp.empty(shape=shape, dtype=dtype, pinned=pinned, device=device, requires_grad=requires_grad)
|
|
376
|
+
)
|
|
377
|
+
|
|
378
|
+
return temporary_store.borrow(shape=shape, dtype=dtype, device=device, pinned=pinned, requires_grad=requires_grad)
|
|
379
|
+
|
|
380
|
+
|
|
381
|
+
def borrow_temporary_like(
|
|
382
|
+
array: Union[wp.array, Temporary],
|
|
383
|
+
temporary_store: Optional[TemporaryStore],
|
|
384
|
+
) -> Temporary:
|
|
385
|
+
"""
|
|
386
|
+
Borrows and returns a temporary array with the same attributes as another array or temporary.
|
|
387
|
+
|
|
388
|
+
Args:
|
|
389
|
+
array: Warp or temporary array to read the desired attributes from
|
|
390
|
+
temporary_store: the shared pool to borrow the temporary from. If `temporary_store` is ``None``, the global default temporary store, if set, will be used.
|
|
391
|
+
"""
|
|
392
|
+
if isinstance(array, Temporary):
|
|
393
|
+
array = array.array
|
|
394
|
+
return borrow_temporary(
|
|
395
|
+
temporary_store=temporary_store,
|
|
396
|
+
shape=array.shape,
|
|
397
|
+
dtype=array.dtype,
|
|
398
|
+
pinned=array.pinned,
|
|
399
|
+
device=array.device,
|
|
400
|
+
requires_grad=array.requires_grad,
|
|
401
|
+
)
|