warp-lang 1.0.2__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.

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