Commit 6cc6aab0 authored by Martin Reinecke's avatar Martin Reinecke

fixes

parent fc68b05e
...@@ -5,4 +5,4 @@ def ApplyData(data, var, model_data): ...@@ -5,4 +5,4 @@ def ApplyData(data, var, model_data):
from ..sugar import sqrt from ..sugar import sqrt
sqrt_n = DiagonalOperator(sqrt(var)) sqrt_n = DiagonalOperator(sqrt(var))
data = Constant(model_data.position, data) data = Constant(model_data.position, data)
return sqrt_n.inverse_times(model_data - data) return sqrt_n.inverse(model_data - data)
...@@ -28,7 +28,7 @@ class PointSources(Model): ...@@ -28,7 +28,7 @@ class PointSources(Model):
@property @property
@memo @memo
def gradient(self): def jacobian(self):
u = self.position['points'].local_data u = self.position['points'].local_data
inner = norm.pdf(u) inner = norm.pdf(u)
outer_inv = invgamma.pdf(invgamma.ppf(norm.cdf(u), outer_inv = invgamma.pdf(invgamma.ppf(norm.cdf(u),
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment