warp-lang 1.0.2__py3-none-manylinux2014_x86_64.whl → 1.2.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 (356) 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 +88 -113
  6. warp/build_dll.py +383 -375
  7. warp/builtins.py +3693 -3354
  8. warp/codegen.py +2925 -2792
  9. warp/config.py +40 -36
  10. warp/constants.py +49 -45
  11. warp/context.py +5409 -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 +381 -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 +145 -146
  34. warp/examples/benchmarks/benchmark_launches.py +293 -295
  35. warp/examples/browse.py +29 -29
  36. warp/examples/core/example_dem.py +232 -219
  37. warp/examples/core/example_fluid.py +291 -267
  38. warp/examples/core/example_graph_capture.py +142 -126
  39. warp/examples/core/example_marching_cubes.py +186 -174
  40. warp/examples/core/example_mesh.py +172 -155
  41. warp/examples/core/example_mesh_intersect.py +203 -193
  42. warp/examples/core/example_nvdb.py +174 -170
  43. warp/examples/core/example_raycast.py +103 -90
  44. warp/examples/core/example_raymarch.py +197 -178
  45. warp/examples/core/example_render_opengl.py +183 -141
  46. warp/examples/core/example_sph.py +403 -387
  47. warp/examples/core/example_torch.py +219 -181
  48. warp/examples/core/example_wave.py +261 -248
  49. warp/examples/fem/bsr_utils.py +378 -380
  50. warp/examples/fem/example_apic_fluid.py +432 -389
  51. warp/examples/fem/example_burgers.py +262 -0
  52. warp/examples/fem/example_convection_diffusion.py +180 -168
  53. warp/examples/fem/example_convection_diffusion_dg.py +217 -209
  54. warp/examples/fem/example_deformed_geometry.py +175 -159
  55. warp/examples/fem/example_diffusion.py +199 -173
  56. warp/examples/fem/example_diffusion_3d.py +178 -152
  57. warp/examples/fem/example_diffusion_mgpu.py +219 -214
  58. warp/examples/fem/example_mixed_elasticity.py +242 -222
  59. warp/examples/fem/example_navier_stokes.py +257 -243
  60. warp/examples/fem/example_stokes.py +218 -192
  61. warp/examples/fem/example_stokes_transfer.py +263 -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 +258 -246
  65. warp/examples/optim/example_cloth_throw.py +220 -209
  66. warp/examples/optim/example_diffray.py +564 -536
  67. warp/examples/optim/example_drone.py +862 -835
  68. warp/examples/optim/example_inverse_kinematics.py +174 -168
  69. warp/examples/optim/example_inverse_kinematics_torch.py +183 -169
  70. warp/examples/optim/example_spring_cage.py +237 -231
  71. warp/examples/optim/example_trajectory.py +221 -199
  72. warp/examples/optim/example_walker.py +304 -293
  73. warp/examples/sim/example_cartpole.py +137 -129
  74. warp/examples/sim/example_cloth.py +194 -186
  75. warp/examples/sim/example_granular.py +122 -111
  76. warp/examples/sim/example_granular_collision_sdf.py +195 -186
  77. warp/examples/sim/example_jacobian_ik.py +234 -214
  78. warp/examples/sim/example_particle_chain.py +116 -105
  79. warp/examples/sim/example_quadruped.py +191 -180
  80. warp/examples/sim/example_rigid_chain.py +195 -187
  81. warp/examples/sim/example_rigid_contact.py +187 -177
  82. warp/examples/sim/example_rigid_force.py +125 -125
  83. warp/examples/sim/example_rigid_gyroscopic.py +107 -95
  84. warp/examples/sim/example_rigid_soft_contact.py +132 -122
  85. warp/examples/sim/example_soft_body.py +188 -177
  86. warp/fabric.py +337 -335
  87. warp/fem/__init__.py +61 -27
  88. warp/fem/cache.py +403 -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 +16 -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 +748 -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 +437 -435
  104. warp/fem/geometry/hexmesh.py +953 -953
  105. warp/fem/geometry/nanogrid.py +455 -0
  106. warp/fem/geometry/partition.py +374 -376
  107. warp/fem/geometry/quadmesh_2d.py +532 -532
  108. warp/fem/geometry/tetmesh.py +840 -840
  109. warp/fem/geometry/trimesh_2d.py +577 -577
  110. warp/fem/integrate.py +1684 -1615
  111. warp/fem/operator.py +190 -191
  112. warp/fem/polynomial.py +214 -213
  113. warp/fem/quadrature/__init__.py +2 -2
  114. warp/fem/quadrature/pic_quadrature.py +243 -245
  115. warp/fem/quadrature/quadrature.py +295 -294
  116. warp/fem/space/__init__.py +179 -292
  117. warp/fem/space/basis_space.py +522 -489
  118. warp/fem/space/collocated_function_space.py +100 -105
  119. warp/fem/space/dof_mapper.py +236 -236
  120. warp/fem/space/function_space.py +148 -145
  121. warp/fem/space/grid_2d_function_space.py +148 -267
  122. warp/fem/space/grid_3d_function_space.py +167 -306
  123. warp/fem/space/hexmesh_function_space.py +253 -352
  124. warp/fem/space/nanogrid_function_space.py +202 -0
  125. warp/fem/space/partition.py +350 -350
  126. warp/fem/space/quadmesh_2d_function_space.py +261 -369
  127. warp/fem/space/restriction.py +161 -160
  128. warp/fem/space/shape/__init__.py +90 -15
  129. warp/fem/space/shape/cube_shape_function.py +728 -738
  130. warp/fem/space/shape/shape_function.py +102 -103
  131. warp/fem/space/shape/square_shape_function.py +611 -611
  132. warp/fem/space/shape/tet_shape_function.py +565 -567
  133. warp/fem/space/shape/triangle_shape_function.py +429 -429
  134. warp/fem/space/tetmesh_function_space.py +224 -292
  135. warp/fem/space/topology.py +297 -295
  136. warp/fem/space/trimesh_2d_function_space.py +153 -221
  137. warp/fem/types.py +77 -77
  138. warp/fem/utils.py +495 -495
  139. warp/jax.py +166 -141
  140. warp/jax_experimental.py +341 -339
  141. warp/native/array.h +1081 -1025
  142. warp/native/builtin.h +1603 -1560
  143. warp/native/bvh.cpp +402 -398
  144. warp/native/bvh.cu +533 -525
  145. warp/native/bvh.h +430 -429
  146. warp/native/clang/clang.cpp +496 -464
  147. warp/native/crt.cpp +42 -32
  148. warp/native/crt.h +352 -335
  149. warp/native/cuda_crt.h +1049 -1049
  150. warp/native/cuda_util.cpp +549 -540
  151. warp/native/cuda_util.h +288 -203
  152. warp/native/cutlass_gemm.cpp +34 -34
  153. warp/native/cutlass_gemm.cu +372 -372
  154. warp/native/error.cpp +66 -66
  155. warp/native/error.h +27 -27
  156. warp/native/exports.h +187 -0
  157. warp/native/fabric.h +228 -228
  158. warp/native/hashgrid.cpp +301 -278
  159. warp/native/hashgrid.cu +78 -77
  160. warp/native/hashgrid.h +227 -227
  161. warp/native/initializer_array.h +32 -32
  162. warp/native/intersect.h +1204 -1204
  163. warp/native/intersect_adj.h +365 -365
  164. warp/native/intersect_tri.h +322 -322
  165. warp/native/marching.cpp +2 -2
  166. warp/native/marching.cu +497 -497
  167. warp/native/marching.h +2 -2
  168. warp/native/mat.h +1545 -1498
  169. warp/native/matnn.h +333 -333
  170. warp/native/mesh.cpp +203 -203
  171. warp/native/mesh.cu +292 -293
  172. warp/native/mesh.h +1887 -1887
  173. warp/native/nanovdb/GridHandle.h +366 -0
  174. warp/native/nanovdb/HostBuffer.h +590 -0
  175. warp/native/nanovdb/NanoVDB.h +6624 -4782
  176. warp/native/nanovdb/PNanoVDB.h +3390 -2553
  177. warp/native/noise.h +850 -850
  178. warp/native/quat.h +1112 -1085
  179. warp/native/rand.h +303 -299
  180. warp/native/range.h +108 -108
  181. warp/native/reduce.cpp +156 -156
  182. warp/native/reduce.cu +348 -348
  183. warp/native/runlength_encode.cpp +61 -61
  184. warp/native/runlength_encode.cu +46 -46
  185. warp/native/scan.cpp +30 -30
  186. warp/native/scan.cu +36 -36
  187. warp/native/scan.h +7 -7
  188. warp/native/solid_angle.h +442 -442
  189. warp/native/sort.cpp +94 -94
  190. warp/native/sort.cu +97 -97
  191. warp/native/sort.h +14 -14
  192. warp/native/sparse.cpp +337 -337
  193. warp/native/sparse.cu +544 -544
  194. warp/native/spatial.h +630 -630
  195. warp/native/svd.h +562 -562
  196. warp/native/temp_buffer.h +30 -30
  197. warp/native/vec.h +1177 -1133
  198. warp/native/volume.cpp +529 -297
  199. warp/native/volume.cu +58 -32
  200. warp/native/volume.h +960 -538
  201. warp/native/volume_builder.cu +446 -425
  202. warp/native/volume_builder.h +34 -19
  203. warp/native/volume_impl.h +61 -0
  204. warp/native/warp.cpp +1057 -1052
  205. warp/native/warp.cu +2949 -2828
  206. warp/native/warp.h +321 -305
  207. warp/optim/__init__.py +9 -9
  208. warp/optim/adam.py +120 -120
  209. warp/optim/linear.py +1104 -939
  210. warp/optim/sgd.py +104 -92
  211. warp/render/__init__.py +10 -10
  212. warp/render/render_opengl.py +3356 -3204
  213. warp/render/render_usd.py +768 -749
  214. warp/render/utils.py +152 -150
  215. warp/sim/__init__.py +52 -59
  216. warp/sim/articulation.py +685 -685
  217. warp/sim/collide.py +1594 -1590
  218. warp/sim/import_mjcf.py +489 -481
  219. warp/sim/import_snu.py +220 -221
  220. warp/sim/import_urdf.py +536 -516
  221. warp/sim/import_usd.py +887 -881
  222. warp/sim/inertia.py +316 -317
  223. warp/sim/integrator.py +234 -233
  224. warp/sim/integrator_euler.py +1956 -1956
  225. warp/sim/integrator_featherstone.py +1917 -1991
  226. warp/sim/integrator_xpbd.py +3288 -3312
  227. warp/sim/model.py +4473 -4314
  228. warp/sim/particles.py +113 -112
  229. warp/sim/render.py +417 -403
  230. warp/sim/utils.py +413 -410
  231. warp/sparse.py +1289 -1227
  232. warp/stubs.py +2192 -2469
  233. warp/tape.py +1162 -225
  234. warp/tests/__init__.py +1 -1
  235. warp/tests/__main__.py +4 -4
  236. warp/tests/assets/test_index_grid.nvdb +0 -0
  237. warp/tests/assets/torus.usda +105 -105
  238. warp/tests/aux_test_class_kernel.py +26 -26
  239. warp/tests/aux_test_compile_consts_dummy.py +10 -10
  240. warp/tests/aux_test_conditional_unequal_types_kernels.py +21 -21
  241. warp/tests/aux_test_dependent.py +20 -22
  242. warp/tests/aux_test_grad_customs.py +21 -23
  243. warp/tests/aux_test_reference.py +9 -11
  244. warp/tests/aux_test_reference_reference.py +8 -10
  245. warp/tests/aux_test_square.py +15 -17
  246. warp/tests/aux_test_unresolved_func.py +14 -14
  247. warp/tests/aux_test_unresolved_symbol.py +14 -14
  248. warp/tests/disabled_kinematics.py +237 -239
  249. warp/tests/run_coverage_serial.py +31 -31
  250. warp/tests/test_adam.py +155 -157
  251. warp/tests/test_arithmetic.py +1088 -1124
  252. warp/tests/test_array.py +2415 -2326
  253. warp/tests/test_array_reduce.py +148 -150
  254. warp/tests/test_async.py +666 -656
  255. warp/tests/test_atomic.py +139 -141
  256. warp/tests/test_bool.py +212 -149
  257. warp/tests/test_builtins_resolution.py +1290 -1292
  258. warp/tests/test_bvh.py +162 -171
  259. warp/tests/test_closest_point_edge_edge.py +227 -228
  260. warp/tests/test_codegen.py +562 -553
  261. warp/tests/test_compile_consts.py +217 -101
  262. warp/tests/test_conditional.py +244 -246
  263. warp/tests/test_copy.py +230 -215
  264. warp/tests/test_ctypes.py +630 -632
  265. warp/tests/test_dense.py +65 -67
  266. warp/tests/test_devices.py +89 -98
  267. warp/tests/test_dlpack.py +528 -529
  268. warp/tests/test_examples.py +403 -378
  269. warp/tests/test_fabricarray.py +952 -955
  270. warp/tests/test_fast_math.py +60 -54
  271. warp/tests/test_fem.py +1298 -1278
  272. warp/tests/test_fp16.py +128 -130
  273. warp/tests/test_func.py +336 -337
  274. warp/tests/test_generics.py +596 -571
  275. warp/tests/test_grad.py +885 -640
  276. warp/tests/test_grad_customs.py +331 -336
  277. warp/tests/test_hash_grid.py +208 -164
  278. warp/tests/test_import.py +37 -39
  279. warp/tests/test_indexedarray.py +1132 -1134
  280. warp/tests/test_intersect.py +65 -67
  281. warp/tests/test_jax.py +305 -307
  282. warp/tests/test_large.py +169 -164
  283. warp/tests/test_launch.py +352 -354
  284. warp/tests/test_lerp.py +217 -261
  285. warp/tests/test_linear_solvers.py +189 -171
  286. warp/tests/test_lvalue.py +419 -493
  287. warp/tests/test_marching_cubes.py +63 -65
  288. warp/tests/test_mat.py +1799 -1827
  289. warp/tests/test_mat_lite.py +113 -115
  290. warp/tests/test_mat_scalar_ops.py +2905 -2889
  291. warp/tests/test_math.py +124 -193
  292. warp/tests/test_matmul.py +498 -499
  293. warp/tests/test_matmul_lite.py +408 -410
  294. warp/tests/test_mempool.py +186 -190
  295. warp/tests/test_mesh.py +281 -324
  296. warp/tests/test_mesh_query_aabb.py +226 -241
  297. warp/tests/test_mesh_query_point.py +690 -702
  298. warp/tests/test_mesh_query_ray.py +290 -303
  299. warp/tests/test_mlp.py +274 -276
  300. warp/tests/test_model.py +108 -110
  301. warp/tests/test_module_hashing.py +111 -0
  302. warp/tests/test_modules_lite.py +36 -39
  303. warp/tests/test_multigpu.py +161 -163
  304. warp/tests/test_noise.py +244 -248
  305. warp/tests/test_operators.py +248 -250
  306. warp/tests/test_options.py +121 -125
  307. warp/tests/test_peer.py +131 -137
  308. warp/tests/test_pinned.py +76 -78
  309. warp/tests/test_print.py +52 -54
  310. warp/tests/test_quat.py +2084 -2086
  311. warp/tests/test_rand.py +324 -288
  312. warp/tests/test_reload.py +207 -217
  313. warp/tests/test_rounding.py +177 -179
  314. warp/tests/test_runlength_encode.py +188 -190
  315. warp/tests/test_sim_grad.py +241 -0
  316. warp/tests/test_sim_kinematics.py +89 -97
  317. warp/tests/test_smoothstep.py +166 -168
  318. warp/tests/test_snippet.py +303 -266
  319. warp/tests/test_sparse.py +466 -460
  320. warp/tests/test_spatial.py +2146 -2148
  321. warp/tests/test_special_values.py +362 -0
  322. warp/tests/test_streams.py +484 -473
  323. warp/tests/test_struct.py +708 -675
  324. warp/tests/test_tape.py +171 -148
  325. warp/tests/test_torch.py +741 -743
  326. warp/tests/test_transient_module.py +85 -87
  327. warp/tests/test_types.py +554 -659
  328. warp/tests/test_utils.py +488 -499
  329. warp/tests/test_vec.py +1262 -1268
  330. warp/tests/test_vec_lite.py +71 -73
  331. warp/tests/test_vec_scalar_ops.py +2097 -2099
  332. warp/tests/test_verify_fp.py +92 -94
  333. warp/tests/test_volume.py +961 -736
  334. warp/tests/test_volume_write.py +338 -265
  335. warp/tests/unittest_serial.py +38 -37
  336. warp/tests/unittest_suites.py +367 -359
  337. warp/tests/unittest_utils.py +434 -578
  338. warp/tests/unused_test_misc.py +69 -71
  339. warp/tests/walkthrough_debug.py +85 -85
  340. warp/thirdparty/appdirs.py +598 -598
  341. warp/thirdparty/dlpack.py +143 -143
  342. warp/thirdparty/unittest_parallel.py +563 -561
  343. warp/torch.py +321 -295
  344. warp/types.py +4941 -4450
  345. warp/utils.py +1008 -821
  346. {warp_lang-1.0.2.dist-info → warp_lang-1.2.0.dist-info}/LICENSE.md +126 -126
  347. {warp_lang-1.0.2.dist-info → warp_lang-1.2.0.dist-info}/METADATA +365 -400
  348. warp_lang-1.2.0.dist-info/RECORD +359 -0
  349. warp/examples/assets/cube.usda +0 -42
  350. warp/examples/assets/sphere.usda +0 -56
  351. warp/examples/assets/torus.usda +0 -105
  352. warp/examples/fem/example_convection_diffusion_dg0.py +0 -194
  353. warp/native/nanovdb/PNanoVDBWrite.h +0 -295
  354. warp_lang-1.0.2.dist-info/RECORD +0 -352
  355. {warp_lang-1.0.2.dist-info → warp_lang-1.2.0.dist-info}/WHEEL +0 -0
  356. {warp_lang-1.0.2.dist-info → warp_lang-1.2.0.dist-info}/top_level.txt +0 -0
warp/fem/cache.py CHANGED
@@ -1,388 +1,403 @@
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
+ view = 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
+ view._ref = array
242
+ return view
243
+
244
+ self._array_view = _view_reshaped_truncated(self._raw_array)
245
+ return self
246
+
247
+ def __del__(self):
248
+ self.release()
249
+
250
+
251
+ class TemporaryStore:
252
+ """
253
+ Shared pool of temporary arrays that will be persisted and reused across invocations of ``warp.fem`` functions.
254
+
255
+ A :class:`TemporaryStore` instance may either be passed explicitly to ``warp.fem`` functions that accept such an argument, for instance :func:`.integrate.integrate`,
256
+ or can be set globally as the default store using :func:`set_default_temporary_store`.
257
+
258
+ By default, there is no default temporary store, so that temporary allocations are not persisted.
259
+ """
260
+
261
+ _default_store: "TemporaryStore" = None
262
+
263
+ class Pool:
264
+ def __init__(self, dtype, device, pinned: bool):
265
+ self.dtype = dtype
266
+ self.device = device
267
+ self.pinned = pinned
268
+
269
+ self._pool = [] # Currently available arrays for borrowing, ordered by size
270
+ self._pool_sizes = [] # Sizes of available arrays for borrowing, ascending
271
+ self._allocs = {} # All allocated arrays, including borrowed ones
272
+
273
+ def borrow(self, shape, dtype, requires_grad: bool):
274
+ size = 1
275
+ if isinstance(shape, int):
276
+ shape = (shape,)
277
+ for d in shape:
278
+ size *= d
279
+
280
+ index = bisect.bisect_left(
281
+ a=self._pool_sizes,
282
+ x=size,
283
+ )
284
+ if index < len(self._pool):
285
+ # Big enough array found, remove from pool
286
+ array = self._pool.pop(index)
287
+ self._pool_sizes.pop(index)
288
+ if requires_grad and array.grad is None:
289
+ array.requires_grad = True
290
+ return Temporary(pool=self, array=array, shape=shape, dtype=dtype)
291
+
292
+ # No big enough array found, allocate new one
293
+ if len(self._pool) > 0:
294
+ grow_factor = 1.5
295
+ size = max(int(self._pool_sizes[-1] * grow_factor), size)
296
+
297
+ array = wp.empty(
298
+ shape=(size,), dtype=self.dtype, pinned=self.pinned, device=self.device, requires_grad=requires_grad
299
+ )
300
+ self._allocs[array.ptr] = array
301
+ return Temporary(pool=self, array=array, shape=shape, dtype=dtype)
302
+
303
+ def redeem(self, array):
304
+ # Insert back array into available pool
305
+ index = bisect.bisect_left(
306
+ a=self._pool_sizes,
307
+ x=array.size,
308
+ )
309
+ self._pool.insert(index, array)
310
+ self._pool_sizes.insert(index, array.size)
311
+
312
+ def detach(self, array):
313
+ del self._allocs[array.ptr]
314
+
315
+ def __init__(self):
316
+ self.clear()
317
+
318
+ def clear(self):
319
+ self._temporaries = {}
320
+
321
+ def borrow(self, shape, dtype, pinned: bool = False, device=None, requires_grad: bool = False) -> Temporary:
322
+ dtype = wp.types.type_to_warp(dtype)
323
+ device = wp.get_device(device)
324
+
325
+ type_length = wp.types.type_length(dtype)
326
+ key = (dtype._type_, type_length, pinned, device.ordinal)
327
+
328
+ pool = self._temporaries.get(key, None)
329
+ if pool is None:
330
+ value_type = (
331
+ cached_vec_type(length=type_length, dtype=wp.types.type_scalar_type(dtype))
332
+ if type_length > 1
333
+ else dtype
334
+ )
335
+ pool = TemporaryStore.Pool(value_type, device, pinned=pinned)
336
+ self._temporaries[key] = pool
337
+
338
+ return pool.borrow(dtype=dtype, shape=shape, requires_grad=requires_grad)
339
+
340
+
341
+ def set_default_temporary_store(temporary_store: Optional[TemporaryStore]):
342
+ """Globally sets the default :class:`TemporaryStore` instance to use for temporary allocations in ``warp.fem`` functions.
343
+
344
+ 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.
345
+ """
346
+
347
+ TemporaryStore._default_store = temporary_store
348
+
349
+
350
+ def borrow_temporary(
351
+ temporary_store: Optional[TemporaryStore],
352
+ shape: Union[int, Tuple[int]],
353
+ dtype: type,
354
+ pinned: bool = False,
355
+ requires_grad: bool = False,
356
+ device=None,
357
+ ) -> Temporary:
358
+ """
359
+ Borrows and returns a temporary array with specified attributes from a shared pool.
360
+
361
+ If an array with sufficient capacity and matching desired attributes is already available in the pool, it will be returned.
362
+ Otherwise, a new allocation will be performed.
363
+
364
+ Args:
365
+ 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.
366
+ shape: desired dimensions for the temporary array
367
+ dtype: desired data type for the temporary array
368
+ pinned: whether a pinned allocation is desired
369
+ device: device on which the memory should be allocated; if ``None``, the current device will be used.
370
+ """
371
+
372
+ if temporary_store is None:
373
+ temporary_store = TemporaryStore._default_store
374
+
375
+ if temporary_store is None:
376
+ return Temporary(
377
+ array=wp.empty(shape=shape, dtype=dtype, pinned=pinned, device=device, requires_grad=requires_grad)
378
+ )
379
+
380
+ return temporary_store.borrow(shape=shape, dtype=dtype, device=device, pinned=pinned, requires_grad=requires_grad)
381
+
382
+
383
+ def borrow_temporary_like(
384
+ array: Union[wp.array, Temporary],
385
+ temporary_store: Optional[TemporaryStore],
386
+ ) -> Temporary:
387
+ """
388
+ Borrows and returns a temporary array with the same attributes as another array or temporary.
389
+
390
+ Args:
391
+ array: Warp or temporary array to read the desired attributes from
392
+ 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.
393
+ """
394
+ if isinstance(array, Temporary):
395
+ array = array.array
396
+ return borrow_temporary(
397
+ temporary_store=temporary_store,
398
+ shape=array.shape,
399
+ dtype=array.dtype,
400
+ pinned=array.pinned,
401
+ device=array.device,
402
+ requires_grad=array.requires_grad,
403
+ )