3.2.0

examples/grids/Wake_Flow.ipynb

@@ -15,13 +15,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "%pip install phiflow\n",
-    "from phi.jax.flow import *\n",
-    "# from phi.flow import *  # If JAX is not installed. You can use phi.torch or phi.tf as well."
+    "from phi.torch.flow import *\n",
+    "# from phi.flow import *  # If JAX is not installed. You can use phi.torch or phi.tf as well.\n",
+    "from tqdm.notebook import trange"
    ]
   },
   {
@@ -75,7 +76,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 18,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -86,7 +87,8 @@
     "\n",
     "boundary = {'x-': vec(x=2, y=0, z=0), 'x+': ZERO_GRADIENT, 'y': PERIODIC, 'z': PERIODIC}\n",
     "v0 = StaggeredGrid((8., 0, 0), boundary, x=128, y=64, z=8, bounds=Box(x=200, y=100, z=5))\n",
-    "v_trj, p_trj = iterate(step, batch(time=200), v0, None)"
+    "v0, p0 = fluid.make_incompressible(v0, cylinder, Solve('scipy-direct'))\n",
+    "v_trj, p_trj = iterate(step, batch(time=200), v0, p0, range=trange)"
    ]
   },
   {

phi/VERSION

@@ -1 +1 @@
-3.1.0
+3.2.0

phi/field/_angular_velocity.py

@@ -3,7 +3,7 @@
 
 from phi import math
 from ._field import FieldInitializer, get_sample_points
-from ..geom import Geometry
+from ..geom import Geometry, cross
 from ..math import Shape, spatial, instance, Tensor, wrap, Extrapolation
 
 
@@ -36,6 +36,6 @@ def _sample(self, geometry: Geometry, at: str, boundaries: Extrapolation, **kwar
         points = get_sample_points(geometry, at, boundaries)
         distances = points - self.location
         strength = self.strength if self.falloff is None else self.strength * self.falloff(distances)
-        velocity = math.cross_product(strength, distances)
+        velocity = cross(strength, distances)
         velocity = math.sum(velocity, self.location.shape.batch.without(points.shape))
         return velocity

phi/field/_embed.py

@@ -53,7 +53,7 @@ def is_face_valid(self, key) -> bool:
         return False
 
     def pad_values(self, value: Tensor, width: int, dim: str, upper_edge: bool, bounds: Box = None, already_padded: dict = None, **kwargs) -> Tensor:
-        assert bounds is not None, f"{type(self)}.pad() requires 'bounds' argument"
+        assert bounds is not None or (already_padded and not value.shape.is_non_uniform), f"{type(self)}.pad() requires 'bounds' argument"
         if value.shape.is_non_uniform:
             unstacked = unstack(value, value.shape.non_uniform_shape)
             indices = value.shape.non_uniform_shape.meshgrid(names=True)