From 237b64d0a429b3c724902ae366df27f8a32d126c Mon Sep 17 00:00:00 2001
From: "Hutschenreuter, Sebastian (hutsch)" <hutsch@mpa-garching.mpg.de>
Date: Wed, 3 Oct 2018 11:29:16 +0200
Subject: [PATCH] fixed power linearization
---
nifty5/linearization.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/nifty5/linearization.py b/nifty5/linearization.py
index 17675aade..566c6801b 100644
--- a/nifty5/linearization.py
+++ b/nifty5/linearization.py
@@ -1,3 +1,4 @@
+
from __future__ import absolute_import, division, print_function
import numpy as np
@@ -104,7 +105,7 @@ class Linearization(object):
def __pow__(self, power):
if not np.isscalar(power):
return NotImplemented
- return self.new(self._val**power, (self._jac**(power-1)).scale(power))
+ return self.new(self._val**power, makeOp(self._val**(power-1)).scale(power)(self._jac))
def inverse(self):
return self.new(1./self._val, makeOp(-1./(self._val**2))(self._jac))
--
GitLab