From 1b35443a948dac91edd0b900fccc19b237c6d20e Mon Sep 17 00:00:00 2001 From: Philipp Holl Date: Wed, 13 Mar 2024 15:37:27 +0100 Subject: [PATCH] [field] PointCloud no longer simplifies sparse values on construction --- phi/field/_point_cloud.py | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/phi/field/_point_cloud.py b/phi/field/_point_cloud.py index e15befaaf..72430eb31 100644 --- a/phi/field/_point_cloud.py +++ b/phi/field/_point_cloud.py @@ -9,7 +9,7 @@ from ..math.extrapolation import Extrapolation, ConstantExtrapolation, PERIODIC -def PointCloud(elements: Union[Tensor, Geometry], values: Any = 1., extrapolation: Union[Extrapolation, float] = 0., bounds: Box = None): +def PointCloud(elements: Union[Tensor, Geometry], values: Any = 1., extrapolation: Union[Extrapolation, float] = 0., bounds: Box = None) -> Field: """ A `PointCloud` comprises: @@ -41,12 +41,12 @@ def PointCloud(elements: Union[Tensor, Geometry], values: Any = 1., extrapolatio """ if bounds is not None: warnings.warn("bounds argument is deprecated since 2.5 and will be ignored.") - if dual(values): - assert dual(values).rank == 1, f"PointCloud cannot convert values with more than 1 dual dimension." - non_dual_name = dual(values).name[1:] - indices = math.stored_indices(values)[non_dual_name] - values = math.stored_values(values) - elements = elements[{non_dual_name: indices}] + # if dual(values): + # assert dual(values).rank == 1, f"PointCloud cannot convert values with more than 1 dual dimension." + # non_dual_name = dual(values).name[1:] + # indices = math.stored_indices(values)[non_dual_name] + # values = math.stored_values(values) + # elements = elements[{non_dual_name: indices}] if isinstance(elements, Tensor): elements = geom.Point(elements) result = Field(elements, values, extrapolation)