From c6791711cd1e20a53362cc45649c1e5252efa33e Mon Sep 17 00:00:00 2001
From: Philipp Arras <parras@mpa-garching.mpg.de>
Date: Fri, 6 Mar 2020 11:17:35 +0100
Subject: [PATCH] Fixup

---
 nifty6/operators/energy_operators.py | 2 +-
 test/test_energy_gradients.py        | 1 +
 2 files changed, 2 insertions(+), 1 deletion(-)

diff --git a/nifty6/operators/energy_operators.py b/nifty6/operators/energy_operators.py
index c715fc90e..597d66bab 100644
--- a/nifty6/operators/energy_operators.py
+++ b/nifty6/operators/energy_operators.py
@@ -138,7 +138,7 @@ class VariableCovarianceGaussianEnergy(EnergyOperator):
             return res
         mf = {self._r: x.val[self._icov], self._icov: .5*x.val[self._icov]**(-2)}
         metric = makeOp(MultiField.from_dict(mf))
-        return res.add_metric(SandwichOperator(x.jac, metric))
+        return res.add_metric(SandwichOperator.make(x.jac, metric))
 
 
 class GaussianEnergy(EnergyOperator):
diff --git a/test/test_energy_gradients.py b/test/test_energy_gradients.py
index adea87054..1b0972c22 100644
--- a/test/test_energy_gradients.py
+++ b/test/test_energy_gradients.py
@@ -48,6 +48,7 @@ def test_variablecovariancegaussian(field):
     mf = ift.MultiField.from_dict(dc)
     energy = ift.VariableCovarianceGaussianEnergy(field.domain, 'a', 'b')
     ift.extra.check_jacobian_consistency(energy, mf, tol=1e-6)
+    energy(ift.Linearization.make_var(mf, want_metric=True)).metric.draw_sample()
 
 
 def test_gaussian(field):
-- 
GitLab