From 7d4e30a5afc77a8c50a20c70294fb14fdce4ad87 Mon Sep 17 00:00:00 2001
From: Martin Reinecke <martin@mpa-garching.mpg.de>
Date: Wed, 25 Jul 2018 16:17:51 +0200
Subject: [PATCH] reintroduce some clip() calls
---
nifty5/library/point_sources.py | 4 ++++
1 file changed, 4 insertions(+)
diff --git a/nifty5/library/point_sources.py b/nifty5/library/point_sources.py
index 7b45e9571..fa121f675 100644
--- a/nifty5/library/point_sources.py
+++ b/nifty5/library/point_sources.py
@@ -43,6 +43,8 @@ class PointSources(Model):
@memo
def value(self):
points = self.position['points'].local_data
+ # MR FIXME?!
+ points = np.clip(points, None, 8.2)
points = Field.from_local_data(self.position['points'].domain, points)
return self.IG(points, self._alpha, self._q)
@@ -55,6 +57,8 @@ class PointSources(Model):
self._alpha,
scale=self._q),
self._alpha, scale=self._q)
+ # FIXME
+ outer_inv = np.clip(outer_inv, 1e-20, None)
outer = 1/outer_inv
grad = Field.from_local_data(self.position['points'].domain,
inner*outer)
--
GitLab