From 91cb0f37c829bcc0e708283b423859c0dfbf029f Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Thu, 17 Feb 2022 12:16:53 -0500 Subject: [PATCH] favor numpy intrinsic over list-comprehension --- python/damask/_rotation.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/python/damask/_rotation.py b/python/damask/_rotation.py index 1626312d1..e99d2c1b5 100644 --- a/python/damask/_rotation.py +++ b/python/damask/_rotation.py @@ -1567,7 +1567,7 @@ class Rotation: +0.000059719705868660826, -0.00001980756723965647, +0.000003953714684212874, -0.00000036555001439719544]) hmag_squared = np.sum(ho**2.,axis=-1,keepdims=True) - s = np.sum(np.array([t*hmag_squared**i for i,t in enumerate(tfit)]),0) + s = np.sum(tfit*hmag_squared**np.arange(len(tfit)),axis=-1,keepdims=True) with np.errstate(invalid='ignore'): ax = np.where(np.broadcast_to(np.abs(hmag_squared)<1.e-8,ho.shape[:-1]+(4,)), [ 0.0, 0.0, 1.0, 0.0 ],